0028840: Data Exchange - rewrite the STL Reader/Writer
[occt.git] / src / QABugs / QABugs_2.cxx
1 // Created on: 2002-07-18
2 // Created by: QA Admin
3 // Copyright (c) 2002-2014 OPEN CASCADE SAS
4 //
5 // This file is part of Open CASCADE Technology software library.
6 //
7 // This library is free software; you can redistribute it and/or modify it under
8 // the terms of the GNU Lesser General Public License version 2.1 as published
9 // by the Free Software Foundation, with special exception defined in the file
10 // OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
11 // distribution for complete text of the license and disclaimer of any warranty.
12 //
13 // Alternatively, this file may be used under the terms of Open CASCADE
14 // commercial license or contractual agreement.
15
16 #include <QABugs.hxx>
17
18 #include <Draw.hxx>
19 #include <Draw_Interpretor.hxx>
20 #include <DBRep.hxx>
21 #include <DrawTrSurf.hxx>
22 #include <AIS_InteractiveContext.hxx>
23 #include <ViewerTest.hxx>
24 #include <AIS_Shape.hxx>
25 #include <TopoDS_Shape.hxx>
26 #include <TopoDS_Vertex.hxx>
27
28 #include <TopExp_Explorer.hxx>
29 #include <TopoDS_Face.hxx>
30 #include <TopoDS.hxx>
31 #include <BRepBndLib.hxx>
32 #include <gp_Pln.hxx>
33 #include <BRep_Tool.hxx>
34 #include <BRepAlgoAPI_Section.hxx>
35 #include <BRepMesh_IncrementalMesh.hxx>
36 #include <Precision.hxx>
37 #include <Standard_ErrorHandler.hxx>
38
39 #include <stdio.h>
40
41 //=======================================================================
42 //function : OCC527 
43 //purpose  : 
44 //=======================================================================
45 static Standard_Integer OCC527(Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
46 {
47   try
48   {
49     OCC_CATCH_SIGNALS
50     // 1. Verify amount of arguments of the command
51     if (argc != 2) {
52       di << "Usage : " << argv[0] << "OCC527 shape\n";
53       return 1;
54     }
55
56     // 2. Get selected shape
57     TopoDS_Shape aShape = DBRep::Get(argv[1]);
58     if(aShape.IsNull()) { di << "OCC527 FAULTY. Entry shape is NULL"; return 0;} 
59
60     // 3. Explode entry shape on faces and build sections from Zmin to Zmax with step aStep
61     const Standard_Real Zmin = -40.228173882121, Zmax = 96.408126285268, aStep = 1.0;
62     char str[100]; str[0] = 0; Sprintf(str,"Test range: [%f, %f] with step %f\n",Zmin,Zmax,aStep); di << str;
63     int nbf = 0;
64     TopExp_Explorer aExp1;
65     for (aExp1.Init(aShape,TopAbs_FACE); aExp1.More(); aExp1.Next())
66     {
67       // Process one face
68       str[0] = 0; Sprintf(str,"Face #%d: \t",nbf++); di << str;
69       TopoDS_Face aFace = TopoDS::Face(aExp1.Current());
70       
71       // Build BndBox in order to avoid try of building section 
72       // if plane of the one does not intersect BndBox of the face
73       Bnd_Box aFaceBox;
74       BRepBndLib::Add(aFace,aFaceBox);
75       Standard_Real X1,X2,Y1,Y2,Z1,Z2;
76       aFaceBox.Get(X1,Y1,Z1,X2,Y2,Z2);
77
78       // Build sections from Zmin to Zmax with step aStep
79       double gmaxdist = 0.0, gzmax = Zmax;
80       for (double zcur = Zmax; zcur > Zmin; zcur -= aStep)
81       {
82         // If plane of the section does not intersect BndBox of the face do nothing
83         if(zcur < Z1 || zcur > Z2 ) continue;
84                 
85         // Build current section
86         gp_Pln pl(0,0,1,-zcur);
87
88         //
89         di << "BRepAlgoAPI_Section aSection(aFace,pl,Standard_False)\n";
90         BRepAlgoAPI_Section aSection(aFace, pl, Standard_False);
91         aSection.Approximation(Standard_True);
92         aSection.Build();
93         Standard_Boolean IsDone = aSection.IsDone();
94
95         if (IsDone)
96         {
97           const TopoDS_Shape& aResult = aSection.Shape();
98           if (!aResult.IsNull())
99           {
100             double lmaxdist = 0.0;
101             TopExp_Explorer aExp2;
102             for (aExp2.Init(aResult, TopAbs_VERTEX); aExp2.More(); aExp2.Next())
103             {
104               TopoDS_Vertex aV = TopoDS::Vertex(aExp2.Current());
105               Standard_Real  toler = BRep_Tool::Tolerance(aV);
106               double dist = pl.Distance(BRep_Tool::Pnt(aV));
107               if (dist > lmaxdist) lmaxdist = dist;
108               // If section was built check distance beetwen vertexes and plane of the one
109               str[0] = 0;
110               if (dist > toler)
111                 Sprintf(str, "Dist=%f, Toler=%f, Param=%f FAULTY\n", dist, toler, gzmax);
112               else
113                 Sprintf(str, "Dist=%f, Toler=%f, Param=%f\n", dist, toler, gzmax);
114               di << str;
115             }
116             if (lmaxdist > gmaxdist)
117             {
118               gmaxdist = lmaxdist;
119               gzmax = zcur;
120             }
121           }
122         }
123       }
124     }
125   }
126   catch (Standard_Failure) {di << "OCC527 Exception \n" ;return 0;}
127   
128   return 0;
129 }
130
131 void QABugs::Commands_2(Draw_Interpretor& theCommands) {
132   const char *group = "QABugs";
133
134   theCommands.Add("OCC527", "OCC527 shape", __FILE__, OCC527, group);
135   return;
136 }