Small correction
[occt.git] / src / QABugs / QABugs_19.cxx
1 // Created on: 2002-05-21
2 // Created by: QA Admin
3 // Copyright (c) 2002-2012 OPEN CASCADE SAS
4 //
5 // The content of this file is subject to the Open CASCADE Technology Public
6 // License Version 6.5 (the "License"). You may not use the content of this file
7 // except in compliance with the License. Please obtain a copy of the License
8 // at http://www.opencascade.org and read it completely before using this file.
9 //
10 // The Initial Developer of the Original Code is Open CASCADE S.A.S., having its
11 // main offices at: 1, place des Freres Montgolfier, 78280 Guyancourt, France.
12 //
13 // The Original Code and all software distributed under the License is
14 // distributed on an "AS IS" basis, without warranty of any kind, and the
15 // Initial Developer hereby disclaims all such warranties, including without
16 // limitation, any warranties of merchantability, fitness for a particular
17 // purpose or non-infringement. Please see the License for the specific terms
18 // and conditions governing the rights and limitations under the License.
19
20
21 #include <QABugs.hxx>
22
23 #include <Draw_Interpretor.hxx>
24 #include <DBRep.hxx>
25 #include <DrawTrSurf.hxx>
26 #include <AIS_InteractiveContext.hxx>
27 #include <ViewerTest.hxx>
28 #include <AIS_Shape.hxx>
29 #include <TopoDS_Shape.hxx>
30
31 #include <gp_Pnt2d.hxx>
32 #include <gp_Ax1.hxx>
33 #include <GCE2d_MakeSegment.hxx>
34 #include <Geom2d_TrimmedCurve.hxx>
35 #include <DrawTrSurf.hxx>
36
37 #include <Precision.hxx>
38
39 #include <PCollection_HAsciiString.hxx>
40
41 #include <cstdio>
42 #include <cmath>
43 #include <iostream>
44 #include <OSD_PerfMeter.hxx>
45 #include <OSD_Timer.hxx>
46 #include <BRepPrimAPI_MakeBox.hxx>
47 #include <BRepPrimAPI_MakeSphere.hxx>
48 #include <BRepAlgo_Cut.hxx>
49
50
51 static Standard_Integer OCC230 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
52 {
53   if ( argc != 4) {
54     di << "ERROR OCC230: Usage : " << argv[0] << " TrimmedCurve Pnt2d Pnt2d" << "\n";
55     return 1;
56   }
57
58   gp_Pnt2d P1, P2;
59   if ( !DrawTrSurf::GetPoint2d(argv[2],P1)) {
60     di << "ERROR OCC230: " << argv[2] << " is not Pnt2d" << "\n";
61     return 1;
62   }
63   if ( !DrawTrSurf::GetPoint2d(argv[3],P2)) {
64     di << "ERROR OCC230: " << argv[3] << " is not Pnt2d" << "\n";
65     return 1;
66   }
67
68   GCE2d_MakeSegment MakeSegment(P1,P2);
69   Handle(Geom2d_TrimmedCurve) TrimmedCurve = MakeSegment.Value();
70   DrawTrSurf::Set(argv[1], TrimmedCurve);
71   return 0;
72 }
73
74 static Standard_Integer OCC142 (Draw_Interpretor& di, Standard_Integer /*argc*/, const char ** /*argv*/)
75 {
76   for(Standard_Integer i= 0;i <= 20;i++){
77     Handle(PCollection_HAsciiString) pstr = new PCollection_HAsciiString("TEST");
78     pstr->Clear();
79   }
80   di << "OCC142: OK" << "\n";
81   return 0;
82 }
83
84 static Standard_Integer OCC23361 (Draw_Interpretor& di, Standard_Integer /*argc*/, const char ** /*argv*/)
85 {
86   gp_Pnt p(0, 0, 2);
87   
88   gp_Trsf t1, t2;
89   t1.SetRotation(gp_Ax1(p, gp_Dir(0, 1, 0)), -0.49328285294022267);
90   t2.SetRotation(gp_Ax1(p, gp_Dir(0, 0, 1)), 0.87538474718473880);
91
92   gp_Trsf tComp = t2 * t1;
93
94   gp_Pnt p1(10, 3, 4);
95   gp_Pnt p2 = p1.Transformed(tComp);
96   gp_Pnt p3 = p1.Transformed(t1);
97   p3.Transform(t2);
98
99   // points must be equal
100   if ( ! p2.IsEqual(p3, Precision::Confusion()) )
101     di << "ERROR OCC23361: equivalent transformations does not produce equal points" << "\n";
102   else 
103     di << "OCC23361: OK" << "\n";
104
105   return 0;
106 }
107
108 static Standard_Integer OCC23237 (Draw_Interpretor& di, Standard_Integer /*argc*/, const char** /*argv*/)
109 {
110   OSD_PerfMeter aPM("TestMeter",0);
111   OSD_Timer aTM;
112   
113   // run some operation in cycle for about 2 seconds to have good values of times to compare
114   int count = 0;
115   printf("OSD_PerfMeter test.\nRunning Boolean operation on solids in loop.\n");
116   for (; aTM.ElapsedTime() < 2.; count++)
117   {
118     aPM.Start();
119     aTM.Start();
120
121     // do some operation that will take considerable time compared with time or starting / stopping timers
122     BRepPrimAPI_MakeBox aBox (10., 10., 10.);
123     BRepPrimAPI_MakeSphere aSphere (10.);
124     BRepAlgo_Cut (aBox.Shape(), aSphere.Shape());
125
126     aTM.Stop();
127     aPM.Stop();
128   }
129  
130   int aNbEnters = 0;
131   Standard_Real aPerfMeter_CPUtime = 0., aTimer_ElapsedTime = aTM.ElapsedTime();
132
133   perf_get_meter("TestMeter", &aNbEnters, &aPerfMeter_CPUtime);
134
135   Standard_Real aTimeDiff = (fabs(aTimer_ElapsedTime - aPerfMeter_CPUtime) / aTimer_ElapsedTime);
136
137   printf("\nMeasurement results (%d cycles):\n", count);
138   printf("\nOSD_PerfMeter CPU time: %lf\nOSD_Timer elapsed time: %lf\n", aPerfMeter_CPUtime, aTimer_ElapsedTime);
139   printf("Time delta is: %.3lf %%\n", aTimeDiff * 100);
140
141   if (aTimeDiff > 0.2)
142     di << "OCC23237: Error: too much difference between CPU and elapsed times";
143   else if (aNbEnters != count)
144     di << "OCC23237: Error: counter reported by PerfMeter (" << aNbEnters << ") does not correspond to actual number of cycles";
145   else
146     di << "OCC23237: OK";
147
148   return 0;
149 }
150
151 #ifdef HAVE_TBB
152
153 #include <Standard_Atomic.hxx>
154 #include <tbb/blocked_range.h>
155 #include <tbb/parallel_for.h>
156
157 class IncrementerDecrementer
158 {
159 public:
160     IncrementerDecrementer (Standard_Integer* theVal, Standard_Boolean thePositive) : myVal (theVal), myPositive (thePositive)
161     {}
162     void operator() (const tbb::blocked_range<size_t>& r) const
163     {
164         if (myPositive)
165             for (size_t i = r.begin(); i != r.end(); ++i)
166                 Standard_Atomic_Increment (myVal);
167         else
168             for (size_t i = r.begin(); i != r.end(); ++i)
169                 Standard_Atomic_Decrement (myVal);
170     }
171 private:
172     Standard_Integer*   myVal;
173     Standard_Boolean   myPositive;
174 };
175 #endif
176
177 #define QCOMPARE(val1, val2) \
178   di << "Checking " #val1 " == " #val2 << \
179         ((val1) == (val2) ? ": OK\n" : ": Error\n")
180
181 #ifdef HAVE_TBB
182 static Standard_Integer OCC22980 (Draw_Interpretor& di, Standard_Integer /*argc*/, const char ** /*argv*/)
183 {
184   int aSum = 0;
185
186   //check returned value
187   QCOMPARE (Standard_Atomic_Decrement (&aSum), -1);
188   QCOMPARE (Standard_Atomic_Increment (&aSum), 0);
189   QCOMPARE (Standard_Atomic_Increment (&aSum), 1);
190   QCOMPARE (Standard_Atomic_Increment (&aSum), 2);
191 //  QCOMPARE (Standard_Atomic_DecrementTest (&aSum), 0);
192 //  QCOMPARE (Standard_Atomic_DecrementTest (&aSum), 1);
193
194   //check atomicity 
195   aSum = 0;
196   const int N = 1 << 24; //big enough to ensure concurrency
197
198   //increment
199   tbb::parallel_for (tbb::blocked_range<size_t> (0, N), IncrementerDecrementer (&aSum, true));
200   QCOMPARE (aSum, N);
201
202   //decrement
203   tbb::parallel_for (tbb::blocked_range<size_t> (0, N), IncrementerDecrementer (&aSum, false));
204   QCOMPARE (aSum, 0);
205
206   return 0;
207 }
208
209 #else /* HAVE_TBB */
210
211 static Standard_Integer OCC22980 (Draw_Interpretor& di, Standard_Integer /*argc*/, const char **argv)
212 {
213   di << "Test skipped: command " << argv[0] << " requires TBB library\n";
214   return 0;
215 }
216
217 #endif /* HAVE_TBB */
218
219 #include <TDocStd_Application.hxx>
220 #include <XCAFApp_Application.hxx>
221 #include <TDocStd_Document.hxx>
222 #include <XCAFDoc_ShapeTool.hxx>
223 #include <XCAFDoc_DocumentTool.hxx>
224 #include <TDF_Label.hxx>
225 #include <TDataStd_Name.hxx>
226
227 static Standard_Integer OCC23595 (Draw_Interpretor& di, Standard_Integer /*argc*/, const char **argv)
228 {
229   const Handle(TDocStd_Application)& anApp = XCAFApp_Application::GetApplication();
230   Handle(TDocStd_Document) aDoc;
231   anApp->NewDocument ("XmlXCAF", aDoc);
232   QCOMPARE (!aDoc.IsNull(), Standard_True);
233
234   Handle(XCAFDoc_ShapeTool) aShTool = XCAFDoc_DocumentTool::ShapeTool (aDoc->Main());
235
236   //check default value
237   Standard_Boolean aValue = XCAFDoc_ShapeTool::AutoNaming();
238   QCOMPARE (aValue, Standard_True);
239
240   //true
241   XCAFDoc_ShapeTool::SetAutoNaming (Standard_True);
242   TopoDS_Shape aShape = BRepPrimAPI_MakeBox (100., 200., 300.).Shape();
243   TDF_Label aLabel = aShTool->AddShape (aShape);
244   Handle(TDataStd_Name) anAttr;
245   QCOMPARE (aLabel.FindAttribute (TDataStd_Name::GetID(), anAttr), Standard_True);
246
247   //false
248   XCAFDoc_ShapeTool::SetAutoNaming (Standard_False);
249   aShape = BRepPrimAPI_MakeBox (300., 200., 100.).Shape();
250   aLabel = aShTool->AddShape (aShape);
251   QCOMPARE (!aLabel.FindAttribute (TDataStd_Name::GetID(), anAttr), Standard_True);
252
253   //restore
254   XCAFDoc_ShapeTool::SetAutoNaming (aValue);
255
256   return 0;
257 }
258
259 #include <ExprIntrp_GenExp.hxx>
260 Standard_Integer OCC22611 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
261 {
262
263   if (argc != 3) {
264     di << "Usage : " << argv[0] << " string nb\n";
265     return 1;
266   }
267
268   TCollection_AsciiString aToken = argv[1];
269   Standard_Integer aNb = atoi(argv[2]);
270
271   Handle(ExprIntrp_GenExp) aGen = ExprIntrp_GenExp::Create();
272   for (Standard_Integer i=0; i < aNb; i++)
273   {
274     aGen->Process(aToken);
275     Handle(Expr_GeneralExpression) aExpr = aGen->Expression();
276   }
277
278   return 0;
279 }
280
281 Standard_Integer OCC22595 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
282 {
283   gp_Mat M0;
284   di << "M0 = "
285   << " {" << M0(1,1) << "} {" << M0(1,2) << "} {" << M0(1,3) <<"}"
286   << " {" << M0(2,1) << "} {" << M0(2,2) << "} {" << M0(2,3) <<"}"
287   << " {" << M0(1,1) << "} {" << M0(1,2) << "} {" << M0(1,3) <<"}";
288   return 0;
289 }
290
291 #include <TopoDS_Face.hxx>
292 #include <TopoDS_Face.hxx>
293 #include <TopoDS.hxx>
294 #include <BRepBuilderAPI_Transform.hxx>
295 #include <BRepExtrema_DistShapeShape.hxx>
296 #include <BRepTools.hxx>
297
298 Standard_Boolean static OCC23774Test(const TopoDS_Face& grossPlateFace, const TopoDS_Shape& originalWire, Draw_Interpretor& di)
299 {
300   BRepExtrema_DistShapeShape distShapeShape(grossPlateFace,originalWire,Extrema_ExtFlag_MIN);
301   if(!distShapeShape.IsDone()) {
302     di <<"Distance ShapeShape is Not Done\n";
303     return Standard_False;
304   }
305
306   if(distShapeShape.Value() > 0.01) {
307     di << "Wrong Dist = " <<distShapeShape.Value() << "\n";
308     return Standard_False;
309   } else
310     di << "Dist0 = " <<distShapeShape.Value() <<"\n";
311
312   //////////////////////////////////////////////////////////////////////////
313   /// First Flip Y
314   const gp_Pnt2d axis1P1(1474.8199035519228,1249.9995745636970);
315   const gp_Pnt2d axis1P2(1474.8199035519228,1250.9995745636970);
316
317   gp_Vec2d mirrorVector1(axis1P1,axis1P2);
318
319   gp_Trsf2d mirror1;
320   mirror1.SetMirror(gp_Ax2d(axis1P1,mirrorVector1));
321
322   BRepBuilderAPI_Transform transformer1(mirror1);
323   transformer1.Perform(originalWire);
324   if(!transformer1.IsDone()) {
325     di << "Not Done1 " << "\n";
326     return Standard_False;
327   }
328   const TopoDS_Shape& step1ModifiedShape = transformer1.ModifiedShape(originalWire);
329   
330   BRepExtrema_DistShapeShape distShapeShape1(grossPlateFace,step1ModifiedShape,Extrema_ExtFlag_MIN);
331   if(!distShapeShape1.IsDone())
332     return Standard_False;
333   if(distShapeShape1.Value() > 0.01) {
334     di << "Dist = " <<distShapeShape1.Value() <<"\n";
335     return Standard_False;
336   } else
337     di << "Dist1 = " <<distShapeShape1.Value() <<"\n";
338
339   //////////////////////////////////////////////////////////////////////////
340   /// Second flip Y
341   transformer1.Perform(step1ModifiedShape);
342   if(!transformer1.IsDone()) {
343     di << "Not Done1 \n";
344     return Standard_False;
345   }
346   const TopoDS_Shape& step2ModifiedShape = transformer1.ModifiedShape(step1ModifiedShape);
347   BRepTools::Write(step2ModifiedShape, Standard_CString ("SecondMirrorWire.brep"));
348
349   //This is identity matrix for values but for type is gp_Rotation ?!
350   gp_Trsf2d mirror11 = mirror1;
351   mirror11.PreMultiply(mirror1);
352
353   BRepExtrema_DistShapeShape distShapeShape2(grossPlateFace,step2ModifiedShape);//,Extrema_ExtFlag_MIN);
354   if(!distShapeShape2.IsDone())
355     return Standard_False;
356
357   //This last test case give error (the value is 1008.8822038689706)
358   if(distShapeShape2.Value() > 0.01) {
359     di  << "Wrong Dist2 = " <<distShapeShape2.Value() <<"\n";
360     Standard_Integer N = distShapeShape2.NbSolution();
361     di << "Nb = " <<N <<"\n";
362     for (Standard_Integer i=1;i <= N;i++)
363         di <<"Sol(" <<i<<") = " <<distShapeShape2.PointOnShape1(i).Distance(distShapeShape2.PointOnShape2(i)) <<"\n";
364     return Standard_False;
365   }
366   di << "Distance2 = " <<distShapeShape2.Value() <<"\n";
367  
368   return Standard_True;
369 }
370 static Standard_Integer OCC23774(Draw_Interpretor& di, Standard_Integer n, const char** a)
371
372
373   if (n != 3) {
374         di <<"OCC23774: invalid number of input parameters\n";
375         return 1;
376   }
377
378   const char *ns1 = (a[1]), *ns2 = (a[2]);
379   TopoDS_Shape S1(DBRep::Get(ns1)), S2(DBRep::Get(ns2));
380   if (S1.IsNull() || S2.IsNull()) {
381         di <<"OCC23774: Null input shapes\n";
382         return 1;
383   }
384   const TopoDS_Face& aFace  = TopoDS::Face(S1);
385   if(!OCC23774Test(aFace, S2, di))
386         di << "Something is wrong\n";
387
388  return 0;
389 }
390
391 #include <Geom_SurfaceOfRevolution.hxx> 
392
393 static Standard_Integer OCC23683 (Draw_Interpretor& di, Standard_Integer argc,const char ** argv)
394 {
395   if (argc < 2) {
396     di<<"Usage: " << argv[0] << " invalid number of arguments"<<"\n";
397     return 1;
398   }
399
400   Standard_Integer ucontinuity = 1;
401   Standard_Integer vcontinuity = 1;
402   Standard_Boolean iscnu = false;
403   Standard_Boolean iscnv = false;
404   
405   Handle(Geom_Surface) aSurf = DrawTrSurf::GetSurface(argv[1]);
406
407   QCOMPARE (aSurf->IsCNu (ucontinuity), iscnu);
408   QCOMPARE (aSurf->IsCNv (vcontinuity), iscnv);
409
410   return 0;
411 }
412
413 void QABugs::Commands_19(Draw_Interpretor& theCommands) {
414   const char *group = "QABugs";
415
416   theCommands.Add ("OCC230", "OCC230 TrimmedCurve Pnt2d Pnt2d", __FILE__, OCC230, group);
417   theCommands.Add ("OCC142", "OCC142", __FILE__, OCC142, group);
418   theCommands.Add ("OCC23361", "OCC23361", __FILE__, OCC23361, group);
419   theCommands.Add ("OCC23237", "OCC23237", __FILE__, OCC23237, group); 
420   theCommands.Add ("OCC22980", "OCC22980", __FILE__, OCC22980, group);
421   theCommands.Add ("OCC23595", "OCC23595", __FILE__, OCC23595, group);
422   theCommands.Add ("OCC22611", "OCC22611 string nb", __FILE__, OCC22611, group);
423   theCommands.Add ("OCC22595", "OCC22595", __FILE__, OCC22595, group);
424   theCommands.Add ("OCC23774", "OCC23774 shape1 shape2", __FILE__, OCC23774, group);
425   theCommands.Add ("OCC23683", "OCC23683 shape", __FILE__, OCC23683, group);
426
427   return;
428 }