CR23683: Geom_BSplineSurface incorrectly determines continuity for periodic cases
[occt.git] / src / QABugs / QABugs_19.cxx
old mode 100644 (file)
new mode 100755 (executable)
index cdb55f4..36c214d
@@ -18,7 +18,6 @@
 // and conditions governing the rights and limitations under the License.
 
 
-
 #include <QABugs.hxx>
 
 #include <Draw_Interpretor.hxx>
 #include <TopoDS_Shape.hxx>
 
 #include <gp_Pnt2d.hxx>
+#include <gp_Ax1.hxx>
 #include <GCE2d_MakeSegment.hxx>
 #include <Geom2d_TrimmedCurve.hxx>
 #include <DrawTrSurf.hxx>
 
+#include <Precision.hxx>
+
 #include <PCollection_HAsciiString.hxx>
 
-//static Standard_Integer OCC230 (Draw_Interpretor& /*di*/, Standard_Integer /*argc*/, const char ** /*argv*/)
+#include <cstdio>
+#include <cmath>
+#include <iostream>
+#include <OSD_PerfMeter.hxx>
+#include <OSD_Timer.hxx>
+#include <BRepPrimAPI_MakeBox.hxx>
+#include <BRepPrimAPI_MakeSphere.hxx>
+#include <BRepAlgo_Cut.hxx>
+
 static Standard_Integer OCC230 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
 {
   if ( argc != 4) {
@@ -70,11 +80,346 @@ static Standard_Integer OCC142 (Draw_Interpretor& di, Standard_Integer /*argc*/,
   return 0;
 }
 
+static Standard_Integer OCC23361 (Draw_Interpretor& di, Standard_Integer /*argc*/, const char ** /*argv*/)
+{
+  gp_Pnt p(0, 0, 2);
+  
+  gp_Trsf t1, t2;
+  t1.SetRotation(gp_Ax1(p, gp_Dir(0, 1, 0)), -0.49328285294022267);
+  t2.SetRotation(gp_Ax1(p, gp_Dir(0, 0, 1)), 0.87538474718473880);
+
+  gp_Trsf tComp = t2 * t1;
+
+  gp_Pnt p1(10, 3, 4);
+  gp_Pnt p2 = p1.Transformed(tComp);
+  gp_Pnt p3 = p1.Transformed(t1);
+  p3.Transform(t2);
+
+  // points must be equal
+  if ( ! p2.IsEqual(p3, Precision::Confusion()) )
+    di << "ERROR OCC23361: equivalent transformations does not produce equal points" << "\n";
+  else 
+    di << "OCC23361: OK" << "\n";
+
+  return 0;
+}
+
+static Standard_Integer OCC23237 (Draw_Interpretor& di, Standard_Integer /*argc*/, const char** /*argv*/)
+{
+  OSD_PerfMeter aPM("TestMeter",0);
+  OSD_Timer aTM;
+  
+  // run some operation in cycle for about 2 seconds to have good values of times to compare
+  int count = 0;
+  printf("OSD_PerfMeter test.\nRunning Boolean operation on solids in loop.\n");
+  for (; aTM.ElapsedTime() < 2.; count++)
+  {
+    aPM.Start();
+    aTM.Start();
+
+    // do some operation that will take considerable time compared with time or starting / stopping timers
+    BRepPrimAPI_MakeBox aBox (10., 10., 10.);
+    BRepPrimAPI_MakeSphere aSphere (10.);
+    BRepAlgo_Cut (aBox.Shape(), aSphere.Shape());
+
+    aTM.Stop();
+    aPM.Stop();
+  }
+  int aNbEnters = 0;
+  Standard_Real aPerfMeter_CPUtime = 0., aTimer_ElapsedTime = aTM.ElapsedTime();
+
+  perf_get_meter("TestMeter", &aNbEnters, &aPerfMeter_CPUtime);
+
+  Standard_Real aTimeDiff = (fabs(aTimer_ElapsedTime - aPerfMeter_CPUtime) / aTimer_ElapsedTime);
+
+  printf("\nMeasurement results (%d cycles):\n", count);
+  printf("\nOSD_PerfMeter CPU time: %lf\nOSD_Timer elapsed time: %lf\n", aPerfMeter_CPUtime, aTimer_ElapsedTime);
+  printf("Time delta is: %.3lf %%\n", aTimeDiff * 100);
+
+  if (aTimeDiff > 0.2)
+    di << "OCC23237: Error: too much difference between CPU and elapsed times";
+  else if (aNbEnters != count)
+    di << "OCC23237: Error: counter reported by PerfMeter (" << aNbEnters << ") does not correspond to actual number of cycles";
+  else
+    di << "OCC23237: OK";
+
+  return 0;
+}
+
+#ifdef HAVE_TBB
+
+#include <Standard_Atomic.hxx>
+#include <tbb/blocked_range.h>
+#include <tbb/parallel_for.h>
+
+class IncrementerDecrementer
+{
+public:
+    IncrementerDecrementer (Standard_Integer* theVal, Standard_Boolean thePositive) : myVal (theVal), myPositive (thePositive)
+    {}
+    void operator() (const tbb::blocked_range<size_t>& r) const
+    {
+        if (myPositive)
+            for (size_t i = r.begin(); i != r.end(); ++i)
+                Standard_Atomic_Increment (myVal);
+        else
+            for (size_t i = r.begin(); i != r.end(); ++i)
+                Standard_Atomic_Decrement (myVal);
+    }
+private:
+    Standard_Integer*   myVal;
+    Standard_Boolean   myPositive;
+};
+#endif
+
+#define QCOMPARE(val1, val2) \
+  di << "Checking " #val1 " == " #val2 << \
+        ((val1) == (val2) ? ": OK\n" : ": Error\n")
+
+#ifdef HAVE_TBB
+static Standard_Integer OCC22980 (Draw_Interpretor& di, Standard_Integer /*argc*/, const char ** /*argv*/)
+{
+  int aSum = 0;
+
+  //check returned value
+  QCOMPARE (Standard_Atomic_Decrement (&aSum), -1);
+  QCOMPARE (Standard_Atomic_Increment (&aSum), 0);
+  QCOMPARE (Standard_Atomic_Increment (&aSum), 1);
+  QCOMPARE (Standard_Atomic_Increment (&aSum), 2);
+//  QCOMPARE (Standard_Atomic_DecrementTest (&aSum), 0);
+//  QCOMPARE (Standard_Atomic_DecrementTest (&aSum), 1);
+
+  //check atomicity 
+  aSum = 0;
+  const int N = 1 << 24; //big enough to ensure concurrency
+
+  //increment
+  tbb::parallel_for (tbb::blocked_range<size_t> (0, N), IncrementerDecrementer (&aSum, true));
+  QCOMPARE (aSum, N);
+
+  //decrement
+  tbb::parallel_for (tbb::blocked_range<size_t> (0, N), IncrementerDecrementer (&aSum, false));
+  QCOMPARE (aSum, 0);
+
+  return 0;
+}
+
+#else /* HAVE_TBB */
+
+static Standard_Integer OCC22980 (Draw_Interpretor& di, Standard_Integer /*argc*/, const char **argv)
+{
+  di << "Test skipped: command " << argv[0] << " requires TBB library\n";
+  return 0;
+}
+
+#endif /* HAVE_TBB */
+
+#include <TDocStd_Application.hxx>
+#include <XCAFApp_Application.hxx>
+#include <TDocStd_Document.hxx>
+#include <XCAFDoc_ShapeTool.hxx>
+#include <XCAFDoc_DocumentTool.hxx>
+#include <TDF_Label.hxx>
+#include <TDataStd_Name.hxx>
+
+static Standard_Integer OCC23595 (Draw_Interpretor& di, Standard_Integer /*argc*/, const char **argv)
+{
+  const Handle(TDocStd_Application)& anApp = XCAFApp_Application::GetApplication();
+  Handle(TDocStd_Document) aDoc;
+  anApp->NewDocument ("XmlXCAF", aDoc);
+  QCOMPARE (!aDoc.IsNull(), Standard_True);
+
+  Handle(XCAFDoc_ShapeTool) aShTool = XCAFDoc_DocumentTool::ShapeTool (aDoc->Main());
+
+  //check default value
+  Standard_Boolean aValue = XCAFDoc_ShapeTool::AutoNaming();
+  QCOMPARE (aValue, Standard_True);
+
+  //true
+  XCAFDoc_ShapeTool::SetAutoNaming (Standard_True);
+  TopoDS_Shape aShape = BRepPrimAPI_MakeBox (100., 200., 300.).Shape();
+  TDF_Label aLabel = aShTool->AddShape (aShape);
+  Handle(TDataStd_Name) anAttr;
+  QCOMPARE (aLabel.FindAttribute (TDataStd_Name::GetID(), anAttr), Standard_True);
+
+  //false
+  XCAFDoc_ShapeTool::SetAutoNaming (Standard_False);
+  aShape = BRepPrimAPI_MakeBox (300., 200., 100.).Shape();
+  aLabel = aShTool->AddShape (aShape);
+  QCOMPARE (!aLabel.FindAttribute (TDataStd_Name::GetID(), anAttr), Standard_True);
+
+  //restore
+  XCAFDoc_ShapeTool::SetAutoNaming (aValue);
+
+  return 0;
+}
+
+#include <ExprIntrp_GenExp.hxx>
+Standard_Integer OCC22611 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
+{
+
+  if (argc != 3) {
+    di << "Usage : " << argv[0] << " string nb\n";
+    return 1;
+  }
+
+  TCollection_AsciiString aToken = argv[1];
+  Standard_Integer aNb = atoi(argv[2]);
+
+  Handle(ExprIntrp_GenExp) aGen = ExprIntrp_GenExp::Create();
+  for (Standard_Integer i=0; i < aNb; i++)
+  {
+    aGen->Process(aToken);
+    Handle(Expr_GeneralExpression) aExpr = aGen->Expression();
+  }
+
+  return 0;
+}
+
+Standard_Integer OCC22595 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
+{
+  gp_Mat M0;
+  di << "M0 = "
+  << " {" << M0(1,1) << "} {" << M0(1,2) << "} {" << M0(1,3) <<"}"
+  << " {" << M0(2,1) << "} {" << M0(2,2) << "} {" << M0(2,3) <<"}"
+  << " {" << M0(1,1) << "} {" << M0(1,2) << "} {" << M0(1,3) <<"}";
+  return 0;
+}
+
+#include <TopoDS_Face.hxx>
+#include <TopoDS_Face.hxx>
+#include <TopoDS.hxx>
+#include <BRepBuilderAPI_Transform.hxx>
+#include <BRepExtrema_DistShapeShape.hxx>
+#include <BRepTools.hxx>
+
+Standard_Boolean static OCC23774Test(const TopoDS_Face& grossPlateFace, const TopoDS_Shape& originalWire, Draw_Interpretor& di)
+{
+  BRepExtrema_DistShapeShape distShapeShape(grossPlateFace,originalWire,Extrema_ExtFlag_MIN);
+  if(!distShapeShape.IsDone()) {
+    di <<"Distance ShapeShape is Not Done\n";
+    return Standard_False;
+  }
+
+  if(distShapeShape.Value() > 0.01) {
+    di << "Wrong Dist = " <<distShapeShape.Value() << "\n";
+    return Standard_False;
+  } else
+    di << "Dist0 = " <<distShapeShape.Value() <<"\n";
+
+  //////////////////////////////////////////////////////////////////////////
+  /// First Flip Y
+  const gp_Pnt2d axis1P1(1474.8199035519228,1249.9995745636970);
+  const gp_Pnt2d axis1P2(1474.8199035519228,1250.9995745636970);
+
+  gp_Vec2d mirrorVector1(axis1P1,axis1P2);
+
+  gp_Trsf2d mirror1;
+  mirror1.SetMirror(gp_Ax2d(axis1P1,mirrorVector1));
+
+  BRepBuilderAPI_Transform transformer1(mirror1);
+  transformer1.Perform(originalWire);
+  if(!transformer1.IsDone()) {
+    di << "Not Done1 " << "\n";
+    return Standard_False;
+  }
+  const TopoDS_Shape& step1ModifiedShape = transformer1.ModifiedShape(originalWire);
+  
+  BRepExtrema_DistShapeShape distShapeShape1(grossPlateFace,step1ModifiedShape,Extrema_ExtFlag_MIN);
+  if(!distShapeShape1.IsDone())
+    return Standard_False;
+  if(distShapeShape1.Value() > 0.01) {
+    di << "Dist = " <<distShapeShape1.Value() <<"\n";
+    return Standard_False;
+  } else
+    di << "Dist1 = " <<distShapeShape1.Value() <<"\n";
+
+  //////////////////////////////////////////////////////////////////////////
+  /// Second flip Y
+  transformer1.Perform(step1ModifiedShape);
+  if(!transformer1.IsDone()) {
+    di << "Not Done1 \n";
+    return Standard_False;
+  }
+  const TopoDS_Shape& step2ModifiedShape = transformer1.ModifiedShape(step1ModifiedShape);
+  BRepTools::Write(step2ModifiedShape, Standard_CString ("SecondMirrorWire.brep"));
+
+  //This is identity matrix for values but for type is gp_Rotation ?!
+  gp_Trsf2d mirror11 = mirror1;
+  mirror11.PreMultiply(mirror1);
+
+  BRepExtrema_DistShapeShape distShapeShape2(grossPlateFace,step2ModifiedShape);//,Extrema_ExtFlag_MIN);
+  if(!distShapeShape2.IsDone())
+    return Standard_False;
+
+  //This last test case give error (the value is 1008.8822038689706)
+  if(distShapeShape2.Value() > 0.01) {
+    di  << "Wrong Dist2 = " <<distShapeShape2.Value() <<"\n";
+    Standard_Integer N = distShapeShape2.NbSolution();
+    di << "Nb = " <<N <<"\n";
+    for (Standard_Integer i=1;i <= N;i++)
+        di <<"Sol(" <<i<<") = " <<distShapeShape2.PointOnShape1(i).Distance(distShapeShape2.PointOnShape2(i)) <<"\n";
+    return Standard_False;
+  }
+  di << "Distance2 = " <<distShapeShape2.Value() <<"\n";
+  return Standard_True;
+}
+static Standard_Integer OCC23774(Draw_Interpretor& di, Standard_Integer n, const char** a)
+{ 
+
+  if (n != 3) {
+       di <<"OCC23774: invalid number of input parameters\n";
+       return 1;
+  }
+
+  const char *ns1 = (a[1]), *ns2 = (a[2]);
+  TopoDS_Shape S1(DBRep::Get(ns1)), S2(DBRep::Get(ns2));
+  if (S1.IsNull() || S2.IsNull()) {
+       di <<"OCC23774: Null input shapes\n";
+       return 1;
+  }
+  const TopoDS_Face& aFace  = TopoDS::Face(S1);
+  if(!OCC23774Test(aFace, S2, di))
+       di << "Something is wrong\n";
+
+ return 0;
+}
+
+static Standard_Integer OCC23683 (Draw_Interpretor& di, Standard_Integer argc,const char ** argv)
+{
+  if (argc < 2) {
+    di<<"Usage: " << argv[0] << " invalid number of arguments"<<"\n";
+    return 1;
+  }
+
+  Standard_Integer ucontinuity = 1;
+  Standard_Integer vcontinuity = 1;
+  Standard_Boolean iscnu = false;
+  Standard_Boolean iscnv = false;
+  
+  Handle(Geom_Surface) aSurf = DrawTrSurf::GetSurface(argv[1]);
+
+  QCOMPARE (aSurf->IsCNu (ucontinuity), iscnu);
+  QCOMPARE (aSurf->IsCNv (vcontinuity), iscnv);
+
+  return 0;
+}
+
 void QABugs::Commands_19(Draw_Interpretor& theCommands) {
-  char *group = "QABugs";
+  const char *group = "QABugs";
 
   theCommands.Add ("OCC230", "OCC230 TrimmedCurve Pnt2d Pnt2d", __FILE__, OCC230, group);
   theCommands.Add ("OCC142", "OCC142", __FILE__, OCC142, group);
+  theCommands.Add ("OCC23361", "OCC23361", __FILE__, OCC23361, group);
+  theCommands.Add ("OCC23237", "OCC23237", __FILE__, OCC23237, group); 
+  theCommands.Add ("OCC22980", "OCC22980", __FILE__, OCC22980, group);
+  theCommands.Add ("OCC23595", "OCC23595", __FILE__, OCC23595, group);
+  theCommands.Add ("OCC22611", "OCC22611 string nb", __FILE__, OCC22611, group);
+  theCommands.Add ("OCC22595", "OCC22595", __FILE__, OCC22595, group);
+  theCommands.Add ("OCC23774", "OCC23774 shape1 shape2", __FILE__, OCC23774, group);
+  theCommands.Add ("OCC23683", "OCC23683 shape", __FILE__, OCC23683, group);
 
   return;
 }