#include <cstdio>
#include <cmath>
#include <iostream>
-#include <OSD_PerfMeter.hxx>
#include <OSD_Timer.hxx>
+#include <OSD_Parallel.hxx>
+#include <OSD_PerfMeter.hxx>
#include <BRepPrimAPI_MakeBox.hxx>
#include <BRepPrimAPI_MakeSphere.hxx>
#include <BRepAlgo_Cut.hxx>
#include <TCollection_HAsciiString.hxx>
#include <GeomFill_Trihedron.hxx>
#include <BRepOffsetAPI_MakePipe.hxx>
+#include <Standard_Atomic.hxx>
#include <Standard_Version.hxx>
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
+ void operator() (const size_t) 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);
+ if ( myPositive )
+ Standard_Atomic_Increment(myVal);
+ else
+ Standard_Atomic_Decrement(myVal);
}
private:
Standard_Integer* myVal;
- Standard_Boolean myPositive;
+ Standard_Boolean myPositive;
};
-#endif
-#ifdef HAVE_TBB
static Standard_Integer OCC22980 (Draw_Interpretor& di, Standard_Integer /*argc*/, const char ** /*argv*/)
{
int 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));
+ OSD_Parallel::For(0, N, IncrementerDecrementer (&aSum, true));
QCOMPARE (aSum, N);
//decrement
- tbb::parallel_for (tbb::blocked_range<size_t> (0, N), IncrementerDecrementer (&aSum, false));
+ OSD_Parallel::For(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 <Geom_BSplineSurface.hxx>
#include <Draw.hxx>
#include <OSD_Thread.hxx>
-static void GeomConvertTest (Draw_Interpretor& di, Standard_Integer theTargetNbUPoles, Standard_CString theFileName)
-{
- Handle(Geom_Surface) aSurf = DrawTrSurf::GetSurface(theFileName);
- GeomConvert_ApproxSurface aGAS (aSurf, 1e-4, GeomAbs_C1, GeomAbs_C1, 9, 9, 100, 1);
- if (!aGAS.IsDone()) {
- di << "ApproxSurface is not done!" << "\n";
- return;
- }
- const Handle(Geom_BSplineSurface)& aBSurf = aGAS.Surface();
- if (aBSurf.IsNull()) {
- di << "BSplineSurface is not created!" << "\n";
- return;
- }
- di << "Number of UPoles:" << aBSurf->NbUPoles() << "\n";
- QCOMPARE (aBSurf->NbUPoles(), theTargetNbUPoles);
-}
-struct aData {
- Draw_Interpretor* di;
- Standard_Integer nbupoles;
- Standard_CString filename;
+struct GeomConvertTest_Data
+{
+ Standard_Integer nbupoles;
+ Handle(Geom_Surface) surf;
};
-Standard_EXPORT Standard_Address convert(Standard_Address data)
+static Standard_Address GeomConvertTest (Standard_Address data)
{
- aData* info = (aData*) data;
- GeomConvertTest(*(info->di),info->nbupoles,info->filename);
- return NULL;
+ GeomConvertTest_Data* info = (GeomConvertTest_Data*)data;
+
+ GeomConvert_ApproxSurface aGAS (info->surf, 1e-4, GeomAbs_C1, GeomAbs_C1, 9, 9, 100, 1);
+ if (!aGAS.IsDone()) {
+ cout << "Error: ApproxSurface is not done!" << endl;
+ return Standard_False;
+ }
+ const Handle(Geom_BSplineSurface)& aBSurf = aGAS.Surface();
+ if (aBSurf.IsNull()) {
+ cout << "Error: BSplineSurface is not created!" << endl;
+ return Standard_False;
+ }
+ cout << "Number of UPoles:" << aBSurf->NbUPoles();
+ if (aBSurf->NbUPoles() == info->nbupoles)
+ {
+ cout << ": OK" << endl;
+ return data; // any non-null pointer
+ }
+ else
+ {
+ cout << ": Error, must be " << info->nbupoles << endl;
+ return 0;
+ }
}
static Standard_Integer OCC23952sweep (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
{
- if (argc != 3) {
- di << "Usage: " << argv[0] << " invalid number of arguments" << "\n";
- return 1;
- }
- struct aData aStorage;
- aStorage.di = &di;
- aStorage.nbupoles = Draw::Atoi(argv[1]);
- aStorage.filename = argv[2];
+ if (argc != 3) {
+ cout << "Error: invalid number of arguments" << endl;
+ return 1;
+ }
- OSD_Thread aThread1(convert);
- aThread1.Run(&aStorage);
- GeomConvertTest(di,aStorage.nbupoles,aStorage.filename);
- cout << "result of thread: " << aThread1.Wait() << endl;
+ struct GeomConvertTest_Data aStorage;
+ aStorage.nbupoles = Draw::Atoi(argv[1]);
+ aStorage.surf = DrawTrSurf::GetSurface(argv[2]);
+ if (aStorage.surf.IsNull())
+ {
+ cout << "Error: " << argv[2] << " is not a DRAW surface!" << endl;
+ return 0;
+ }
- return 0;
-}
+ // start conversion in several threads
+ const int NBTHREADS = 100;
+ OSD_Thread aThread[NBTHREADS];
+ for (int i=0; i < NBTHREADS; i++)
+ {
+ aThread[i].SetFunction (GeomConvertTest);
+ if (!aThread[i].Run(&aStorage))
+ di << "Error: Cannot start thread << " << i << "\n";
+ }
-#include <GeomInt_IntSS.hxx>
-static void GeomIntSSTest (Draw_Interpretor& di, Standard_Integer theNbSol, Standard_CString theFileName1, Standard_CString theFileName2)
-{
- Handle(Geom_Surface) aSurf1 = DrawTrSurf::GetSurface(theFileName1);
- Handle(Geom_Surface) aSurf2 = DrawTrSurf::GetSurface(theFileName2);
- GeomInt_IntSS anInter;
- anInter.Perform(aSurf1, aSurf2, Precision::Confusion(), Standard_True);
- if (!anInter.IsDone()) {
- di << "An intersection is not done!" << "\n";
- return;
- }
+ // check results
+ for (int i=0; i < NBTHREADS; i++)
+ {
+ Standard_Address aResult = 0;
+ if (!aThread[i].Wait(aResult))
+ di << "Error: Failed waiting for thread << " << i << "\n";
+ if (!aResult)
+ di << "Error: wrong number of poles in thread " << i << "!\n";
+ }
- di << "Number of Lines:" << anInter.NbLines() << "\n";
- QCOMPARE (anInter.NbLines(), theNbSol);
+ return 0;
}
-struct aNewData {
- Draw_Interpretor* di;
- Standard_Integer nbsol;
- Standard_CString filename1;
- Standard_CString filename2;
+#include <GeomInt_IntSS.hxx>
+
+struct GeomIntSSTest_Data
+{
+ Standard_Integer nbsol;
+ Handle(Geom_Surface) surf1, surf2;
};
-Standard_EXPORT Standard_Address convert_inter(Standard_Address data)
+
+static Standard_Address GeomIntSSTest (Standard_Address data)
{
- aNewData* info = (aNewData*) data;
- GeomIntSSTest(*(info->di),info->nbsol,info->filename1,info->filename2);
- return NULL;
+ GeomIntSSTest_Data* info = (GeomIntSSTest_Data*)data;
+ GeomInt_IntSS anInter;
+ anInter.Perform (info->surf1, info->surf2, Precision::Confusion(), Standard_True);
+ if (!anInter.IsDone()) {
+ cout << "An intersection is not done!" << endl;
+ return 0;
+ }
+
+ cout << "Number of Lines:" << anInter.NbLines();
+ if (anInter.NbLines() == info->nbsol)
+ {
+ cout << ": OK" << endl;
+ return data; // any non-null pointer
+ }
+ else
+ {
+ cout << ": Error, must be " << info->nbsol << endl;
+ return 0;
+ }
}
static Standard_Integer OCC23952intersect (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
{
- if (argc != 4) {
- di << "Usage: " << argv[0] << " invalid number of arguments" << "\n";
- return 1;
- }
- struct aNewData aStorage;
- aStorage.di = &di;
- aStorage.nbsol = Draw::Atoi(argv[1]);
- aStorage.filename1 = argv[2];
- aStorage.filename2 = argv[3];
+ if (argc != 4) {
+ cout << "Error: invalid number of arguments" << endl;
+ return 1;
+ }
- OSD_Thread aThread1(convert_inter);
- aThread1.Run(&aStorage);
- GeomIntSSTest(di,aStorage.nbsol,aStorage.filename1,aStorage.filename2);
- cout << "result of thread: " << aThread1.Wait() << endl;
+ struct GeomIntSSTest_Data aStorage;
+ aStorage.nbsol = Draw::Atoi(argv[1]);
+ aStorage.surf1 = DrawTrSurf::GetSurface(argv[2]);
+ aStorage.surf2 = DrawTrSurf::GetSurface(argv[3]);
+ if (aStorage.surf1.IsNull() || aStorage.surf2.IsNull())
+ {
+ cout << "Error: Either " << argv[2] << " or " << argv[3] << " is not a DRAW surface!" << endl;
+ return 0;
+ }
- return 0;
+ // start conversion in several threads
+ const int NBTHREADS = 100;
+ OSD_Thread aThread[NBTHREADS];
+ for (int i=0; i < NBTHREADS; i++)
+ {
+ aThread[i].SetFunction (GeomIntSSTest);
+ if (!aThread[i].Run(&aStorage))
+ di << "Error: Cannot start thread << " << i << "\n";
+ }
+
+ // check results
+ for (int i=0; i < NBTHREADS; i++)
+ {
+ Standard_Address aResult = 0;
+ if (!aThread[i].Wait(aResult))
+ di << "Error: Failed waiting for thread << " << i << "\n";
+ if (!aResult)
+ di << "Error: wrong number of intersections in thread " << i << "!\n";
+ }
+
+ return 0;
}
#include <Geom_SurfaceOfRevolution.hxx>
#include <STEPCAFControl_Writer.hxx>
static Standard_Integer OCC23951 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
{
- if (argc != 1) {
+ if (argc != 2) {
di << "Usage: " << argv[0] << " invalid number of arguments" << "\n";
return 1;
}
return 1;
}
- writer.Write("test_box.step");
+ writer.Write(argv[1]);
return 0;
}
#include <TNaming_Builder.hxx>
static Standard_Integer OCC24931 (Draw_Interpretor& di, Standard_Integer argc, const char** argv)
{
- if (argc != 1) {
+ if (argc != 2) {
di << "Usage: " << argv[0] << " invalid number of arguments"<<"\n";
return 1;
}
- TCollection_ExtendedString aFileName ("testdocument.xml");
+ TCollection_ExtendedString aFileName (argv[1]);
PCDM_StoreStatus aSStatus = PCDM_SS_Failure;
Handle(TDocStd_Application) anApp = new AppStd_Application;
return 0;
}
+#include <ShapeFix_Wireframe.hxx>
+//=======================================================================
+//function : OCC7570
+//purpose :
+//=======================================================================
+static Standard_Integer OCC7570 (Draw_Interpretor& di, Standard_Integer n, const char** a)
+{
+ if (n != 2) {
+ di<<"Usage: "<<a[0]<<" invalid number of arguments"<<"\n";
+ return 1;
+ }
+ TopoDS_Shape in_shape (DBRep::Get (a[1]));
+ ShapeFix_Wireframe fix_tool (in_shape);
+ fix_tool.ModeDropSmallEdges () = Standard_True;
+ fix_tool.SetPrecision (1.e+6);
+ fix_tool.SetLimitAngle (0.01);
+ fix_tool.FixSmallEdges ();
+ TopoDS_Shape new_shape = fix_tool.Shape ();
+ return 0;
+}
+
+#include <AIS_TypeFilter.hxx>
+//=======================================================================
+//function : OCC25340
+//purpose :
+//=======================================================================
+static Standard_Integer OCC25340 (Draw_Interpretor& /*theDI*/,
+ Standard_Integer /*theArgNb*/,
+ const char** /*theArgVec*/)
+{
+ Handle(AIS_InteractiveContext) aCtx = ViewerTest::GetAISContext();
+ if (aCtx.IsNull())
+ {
+ std::cerr << "Error: No opened viewer!\n";
+ return 1;
+ }
+ Handle(AIS_TypeFilter) aFilter = new AIS_TypeFilter (AIS_KOI_Shape);
+ aCtx->AddFilter (aFilter);
+ return 0;
+}
+
+//=======================================================================
+//function : OCC24826
+//purpose :
+//=======================================================================
+class ParallelTest_Saxpy
+{
+public:
+ typedef NCollection_Array1<Standard_Real> Vector;
+
+ //! Constructor
+ ParallelTest_Saxpy(const Vector& theX, Vector& theY, Standard_Real theScalar)
+ : myX(theX),
+ myY(theY),
+ myScalar(theScalar)
+ {
+ }
+
+ //! Dummy calculation
+ void operator() (const Standard_Integer theIndex) const
+ {
+ myY(theIndex) = myScalar * myX(theIndex) + myY(theIndex);
+ }
+
+private:
+ ParallelTest_Saxpy( const ParallelTest_Saxpy& );
+ ParallelTest_Saxpy& operator =( ParallelTest_Saxpy& );
+
+private:
+ const Vector& myX;
+ Vector& myY;
+ const Standard_Real myScalar;
+};
+
+//---------------------------------------------------------------------
+static Standard_Integer OCC24826(Draw_Interpretor& theDI,
+ Standard_Integer trheArgc,
+ const char** theArgv)
+{
+ if ( trheArgc != 2 )
+ {
+ theDI << "Usage: "
+ << theArgv[0]
+ << " vec_length\n";
+ return 1;
+ }
+
+ // Generate data;
+ Standard_Integer aLength = Draw::Atoi(theArgv[1]);
+
+ NCollection_Array1<Standard_Real> aX (0, aLength - 1);
+ NCollection_Array1<Standard_Real> anY(0, aLength - 1);
+
+ for ( Standard_Integer i = 0; i < aLength; ++i )
+ {
+ aX(i) = anY(i) = (Standard_Real) i;
+ }
+
+ OSD_Timer aTimer;
+
+ aTimer.Start();
+
+ //! Serial proccesing
+ for ( Standard_Integer i = 0; i < aLength; ++i )
+ {
+ anY(i) = 1e-6 * aX(i) + anY(i);
+ }
+
+ aTimer.Stop();
+ cout << "Processing time (sequential mode):\n";
+ aTimer.Show();
+
+ const ParallelTest_Saxpy aFunctor(aX, anY, 1e-6);
+
+ aTimer.Reset();
+ aTimer.Start();
+
+ // Parallel processing
+ OSD_Parallel::For(0, aLength, aFunctor);
+
+ aTimer.Stop();
+ cout << "Processing time (parallel mode):\n";
+ aTimer.Show();
+
+ return 0;
+}
+
/*****************************************************************************/
+#include <GeomAPI_IntSS.hxx>
+//=======================================================================
+//function : OCC25100
+//purpose :
+//=======================================================================
+static Standard_Integer OCC25100 (Draw_Interpretor& di, Standard_Integer argc, const char ** argv)
+{
+ if (argc < 2)
+ {
+ di << "the method requires a shape name\n";
+ return 1;
+ }
+
+ TopoDS_Shape S = DBRep::Get(argv[1]);
+ if ( S.IsNull() )
+ {
+ di << "Shape is empty" << "\n";
+ return 1;
+ }
+
+ TopExp_Explorer aFaceExp(S, TopAbs_FACE);
+ const Handle(Geom_Surface)& aSurf = BRep_Tool::Surface(TopoDS::Face(aFaceExp.Current()));
+
+ GeomAPI_IntSS anIntersector(aSurf, aSurf, Precision::Confusion());
+
+ if (!anIntersector.IsDone())
+ {
+ di << "Error. Intersection is not done\n";
+ return 1;
+ }
+
+ di << "Test complete\n";
+
+ return 0;
+}
+
+//=======================================================================
+//function : OCC25348
+//purpose :
+//=======================================================================
+static Standard_Integer OCC25348 (Draw_Interpretor& theDI,
+ Standard_Integer /*theArgNb*/,
+ const char** /*theArgVec*/)
+{
+ Handle(NCollection_IncAllocator) anAlloc1;
+ NCollection_List<int> aList1(anAlloc1);
+ for (int i=0; i < 10; i++)
+ {
+ Handle(NCollection_IncAllocator) anAlloc2;
+ NCollection_List<int> aList2(anAlloc2);
+ aList2.Append(i);
+ aList1.Assign(aList2);
+ }
+ theDI << "Test complete\n";
+ return 0;
+}
+
+#include <IntCurvesFace_ShapeIntersector.hxx>
+#include <BRepBndLib.hxx>
+//=======================================================================
+//function : OCC25413
+//purpose :
+//=======================================================================
+static Standard_Integer OCC25413 (Draw_Interpretor& di, Standard_Integer narg , const char** a)
+{
+ if (narg != 2) {
+ di << "Usage: " << a[0] << " invalid number of arguments" << "\n";
+ return 1;
+ }
+ TopoDS_Shape aShape = DBRep::Get (a[1]);
+
+ IntCurvesFace_ShapeIntersector Inter;
+ Inter.Load(aShape, Precision::Confusion());
+
+ Bnd_Box aBndBox;
+ BRepBndLib::Add(aShape, aBndBox);
+
+ gp_Dir aDir(0., 1., 0.);
+ const int N = 250;
+ Standard_Real xMin = aBndBox.CornerMin().X();
+ Standard_Real zMin = aBndBox.CornerMin().Z();
+ Standard_Real xMax = aBndBox.CornerMax().X();
+ Standard_Real zMax = aBndBox.CornerMax().Z();
+ Standard_Real xStep = (xMax - xMin) / N;
+ Standard_Real zStep = (zMax - zMin) / N;
+
+ for (Standard_Real x = xMin; x <= xMax; x += xStep)
+ for (Standard_Real z = zMin; z <= zMax; z += zStep)
+ {
+ gp_Pnt aPoint(x, 0.0, z);
+ gp_Lin aLine(aPoint, aDir);
+ Inter.PerformNearest(aLine, -100., 100.);
+ }
+ return 0;
+}
+
+
+#include <BOPAlgo_PaveFiller.hxx>
+//
+#include <BRepAlgoAPI_BooleanOperation.hxx>
+#include <BRepAlgoAPI_Common.hxx>
+#include <BRepAlgoAPI_Fuse.hxx>
+#include <BRepAlgoAPI_Cut.hxx>
+#include <BRepAlgoAPI_Section.hxx>
+//
+#include <BOPTools.hxx>
+//
+#include <BOPCol_MapOfShape.hxx>
+#include <BOPCol_ListOfShape.hxx>
+//=======================================================================
+//function : OCC25446
+//purpose :
+//=======================================================================
+static Standard_Integer OCC25446 (Draw_Interpretor& theDI,
+ Standard_Integer argc,
+ const char ** argv)
+{
+ if (argc != 5) {
+ theDI << "Usage: OCC25446 res b1 b2 op\n";
+ return 1;
+ }
+ //
+ TopoDS_Shape aS1 = DBRep::Get(argv[2]);
+ if (aS1.IsNull()) {
+ theDI << argv[2] << " shape is NULL\n";
+ return 1;
+ }
+ //
+ TopoDS_Shape aS2 = DBRep::Get(argv[3]);
+ if (aS2.IsNull()) {
+ theDI << argv[3] << " shape is NULL\n";
+ return 1;
+ }
+ //
+ Standard_Integer iOp;
+ BOPAlgo_Operation aOp;
+ //
+ iOp = Draw::Atoi(argv[4]);
+ if (iOp < 0 || iOp > 4) {
+ theDI << "Invalid operation type\n";
+ return 1;
+ }
+ aOp = (BOPAlgo_Operation)iOp;
+ //
+ Standard_Integer iErr;
+ BOPCol_ListOfShape aLS;
+ BOPAlgo_PaveFiller aPF;
+ //
+ aLS.Append(aS1);
+ aLS.Append(aS2);
+ aPF.SetArguments(aLS);
+ //
+ aPF.Perform();
+ iErr = aPF.ErrorStatus();
+ if (iErr) {
+ theDI << "Intersection failed with error status: " << iErr << "\n";
+ return 1;
+ }
+ //
+ BRepAlgoAPI_BooleanOperation* pBuilder = NULL;
+ //
+ switch (aOp) {
+ case BOPAlgo_COMMON:
+ pBuilder = new BRepAlgoAPI_Common(aS1, aS2, aPF);
+ break;
+ case BOPAlgo_FUSE:
+ pBuilder = new BRepAlgoAPI_Fuse(aS1, aS2, aPF);
+ break;
+ case BOPAlgo_CUT:
+ pBuilder = new BRepAlgoAPI_Cut (aS1, aS2, aPF);
+ break;
+ case BOPAlgo_CUT21:
+ pBuilder = new BRepAlgoAPI_Cut(aS1, aS2, aPF, Standard_False);
+ break;
+ case BOPAlgo_SECTION:
+ pBuilder = new BRepAlgoAPI_Section(aS1, aS2, aPF);
+ break;
+ default:
+ break;
+ }
+ //
+ iErr = pBuilder->ErrorStatus();
+ if (!pBuilder->IsDone()) {
+ theDI << "BOP failed with error status: " << iErr << "\n";
+ return 1;
+ }
+ //
+ const TopoDS_Shape& aRes = pBuilder->Shape();
+ DBRep::Set(argv[1], aRes);
+ //
+ BOPCol_MapOfShape aMapArgs, aMapShape;
+ BOPCol_MapIteratorOfMapOfShape aIt;
+ Standard_Boolean bIsDeletedHist, bIsDeletedMap;
+ TopAbs_ShapeEnum aType;
+ //
+ BOPTools::MapShapes(aS1, aMapArgs);
+ BOPTools::MapShapes(aS2, aMapArgs);
+ BOPTools::MapShapes(aRes, aMapShape);
+ //
+ aIt.Initialize(aMapArgs);
+ for (; aIt.More(); aIt.Next()) {
+ const TopoDS_Shape& aS = aIt.Value();
+ aType = aS.ShapeType();
+ if (!(aType==TopAbs_EDGE || aType==TopAbs_FACE ||
+ aType==TopAbs_VERTEX || aType==TopAbs_SOLID)) {
+ continue;
+ }
+ //
+ bIsDeletedHist = pBuilder->IsDeleted(aS);
+ bIsDeletedMap = !aMapShape.Contains(aS) &&
+ (pBuilder->Modified(aS).Extent() == 0);
+ //
+ if (bIsDeletedHist != bIsDeletedMap) {
+ theDI << "Error. Wrong value of IsDeleted flag.\n";
+ return 1;
+ }
+ }
+ //
+ theDI << "Test complete\n";
+ return 0;
+}
+
+//====================================================
+// Auxiliary functor class for the command OCC25545;
+// it gets access to a vertex with the given index and
+// checks that X coordinate of the point is equal to index;
+// if it is not so then a data race is reported.
+//====================================================
+struct OCC25545_Functor
+{
+ OCC25545_Functor(const std::vector<TopoDS_Shape>& theShapeVec)
+ : myShapeVec(&theShapeVec),
+ myIsRaceDetected(0)
+ {}
+
+ void operator()(size_t i) const
+ {
+ if (!myIsRaceDetected) {
+ const TopoDS_Vertex& aV = TopoDS::Vertex (myShapeVec->at(i));
+ gp_Pnt aP = BRep_Tool::Pnt (aV);
+ if (aP.X () != static_cast<double> (i)) {
+ Standard_Atomic_Increment(&myIsRaceDetected);
+ }
+ }
+ }
+
+ const std::vector<TopoDS_Shape>* myShapeVec;
+ mutable volatile int myIsRaceDetected;
+};
+
+//=======================================================================
+//function : OCC25545
+//purpose : Tests data race when concurrently accessing TopLoc_Location::Transformation()
+//=======================================================================
+#ifdef HAVE_TBB
+static Standard_Integer OCC25545 (Draw_Interpretor& di,
+ Standard_Integer,
+ const char **)
+{
+ // Place vertices in a vector, giving the i-th vertex the
+ // transformation that translates it on the vector (i,0,0) from the origin.
+ size_t n = 1000;
+ std::vector<TopoDS_Shape> aShapeVec (n);
+ std::vector<TopLoc_Location> aLocVec (n);
+ TopoDS_Shape aShape = BRepBuilderAPI_MakeVertex (gp::Origin ());
+ aShapeVec[0] = aShape;
+ for (size_t i = 1; i < n; ++i) {
+ gp_Trsf aT;
+ aT.SetTranslation (gp_Vec (1, 0, 0));
+ aLocVec[i] = aLocVec[i - 1] * aT;
+ aShapeVec[i] = aShape.Moved (aLocVec[i]);
+ }
+
+ // Evaluator function will access vertices geometry
+ // concurrently
+ OCC25545_Functor aFunc(aShapeVec);
+
+ //concurrently process
+ tbb::parallel_for (size_t (0), n, aFunc, tbb::simple_partitioner ());
+ QVERIFY (!aFunc.myIsRaceDetected);
+ return 0;
+}
+#else
+static Standard_Integer OCC25545 (Draw_Interpretor&,
+ Standard_Integer,
+ const char **argv)
+{
+ cout << "Test skipped: command " << argv[0] << " requires TBB library" << endl;
+ return 0;
+}
+#endif
+
+//=======================================================================
+//function : OCC25547
+//purpose :
+//=======================================================================
+#include <BRepMesh_GeomTool.hxx>
+#include <BRepAdaptor_Curve.hxx>
+#include <Geom_TrimmedCurve.hxx>
+#include <BRepBuilderAPI_MakeFace.hxx>
+#include <BRepAdaptor_HSurface.hxx>
+#include <BRepAdaptor_Surface.hxx>
+static Standard_Integer OCC25547(
+ Draw_Interpretor& theDI,
+ Standard_Integer /*argc*/,
+ const char ** /*argv*/)
+{
+ // The general aim of this test is to prevent linkage errors due to missed
+ // Standard_EXPORT attribute for static methods.
+
+ // However, start checking the main functionality at first.
+ const Standard_Real aFirstP = 0., aLastP = M_PI;
+ Handle(Geom_Circle) aCircle = new Geom_Circle(gp_Ax2(gp::Origin(), gp::DZ()), 10);
+ Handle(Geom_TrimmedCurve) aHalf = new Geom_TrimmedCurve(aCircle, aFirstP, aLastP);
+ TopoDS_Edge aEdge = BRepBuilderAPI_MakeEdge(aHalf);
+ BRepAdaptor_Curve aAdaptor(aEdge);
+ BRepMesh_GeomTool aGeomTool(aAdaptor, aFirstP, aLastP, 0.1, 0.5);
+
+ if (aGeomTool.NbPoints() == 0)
+ {
+ theDI << "Error. BRepMesh_GeomTool failed to discretize an arc.\n";
+ return 1;
+ }
+
+ // Test static methods.
+ TopoDS_Face aFace = BRepBuilderAPI_MakeFace(gp_Pln(gp::Origin(), gp::DZ()));
+ BRepAdaptor_Surface aSurf(aFace);
+ Handle(BRepAdaptor_HSurface) aHSurf = new BRepAdaptor_HSurface(aSurf);
+
+ gp_Pnt aPnt;
+ gp_Dir aNormal;
+ if (!BRepMesh_GeomTool::Normal(aHSurf, 10., 10., aPnt, aNormal))
+ {
+ theDI << "Error. BRepMesh_GeomTool failed to take a normal of surface.\n";
+ return 1;
+ }
+
+ gp_XY aRefPnts[4] = {
+ gp_XY(-10., -10.), gp_XY(10., 10.),
+ gp_XY(-10., 10.), gp_XY(10., -10.)
+ };
+
+ gp_Pnt2d aIntPnt;
+ Standard_Real aParams[2];
+ BRepMesh_GeomTool::IntFlag aIntFlag = BRepMesh_GeomTool::IntLinLin(
+ aRefPnts[0], aRefPnts[1], aRefPnts[2], aRefPnts[3],
+ aIntPnt.ChangeCoord(), aParams);
+
+ Standard_Real aDiff = aIntPnt.Distance(gp::Origin2d());
+ if (aIntFlag != BRepMesh_GeomTool::Cross || aDiff > Precision::PConfusion())
+ {
+ theDI << "Error. BRepMesh_GeomTool failed to intersect two lines.\n";
+ return 1;
+ }
+
+ aIntFlag = BRepMesh_GeomTool::IntSegSeg(
+ aRefPnts[0], aRefPnts[1], aRefPnts[2], aRefPnts[3],
+ Standard_False, Standard_False, aIntPnt);
+
+ aDiff = aIntPnt.Distance(gp::Origin2d());
+ if (aIntFlag != BRepMesh_GeomTool::Cross || aDiff > Precision::PConfusion())
+ {
+ theDI << "Error. BRepMesh_GeomTool failed to intersect two segments.\n";
+ return 1;
+ }
+
+
+ theDI << "Test complete\n";
+ return 0;
+}
+
+#include <TColStd_DataMapIteratorOfDataMapOfIntegerInteger.hxx>
+#include <TColStd_DataMapOfIntegerInteger.hxx>
+#include <OSD.hxx>
+#include <ShapeFix_Wire.hxx>
+#include <ShapeExtend_Status.hxx>
+#ifdef WNT
+#define EXCEPTION ...
+#else
+#define EXCEPTION Standard_Failure
+#endif
+
+static ShapeExtend_Status getStatusGap(const Handle(ShapeFix_Wire)& theFix,
+ const Standard_Boolean theIs3d)
+{
+ for (Standard_Integer i=ShapeExtend_OK; i<=ShapeExtend_FAIL; i++)
+ {
+ Standard_Boolean isFound;
+ if (theIs3d)
+ isFound = theFix->StatusGaps3d( (ShapeExtend_Status) i );
+ else
+ isFound = theFix->StatusGaps2d( (ShapeExtend_Status) i );
+ if (isFound) return ShapeExtend_Status(i);
+ }
+ return ShapeExtend_OK;
+}
+
+//===================
+//function : OCC24881
+//purpose :
+//===================
+static Standard_Integer OCC24881 (Draw_Interpretor& di, Standard_Integer narg , const char** a)
+{
+ if (narg < 2) {
+ di<<"Usage: "<<a[0]<<" invalid number of arguments"<<"\n";
+ return 1;
+ }
+// cout <<"FileName1: " << argv[1] <<endl;
+
+ TopoDS_Shape aShape = DBRep::Get (a[1]);
+
+ OSD::SetSignal();
+ Handle(ShapeFix_Wire) aWireFix = new ShapeFix_Wire;
+
+ // map FixStatus - NbSuchStatuses
+ TColStd_DataMapOfIntegerInteger aStatusNbDMap;
+ Standard_Integer nbFixed=0, nbOk=0;
+
+//Begin: STEP 7
+ ShapeExtend_Status aStatus=ShapeExtend_OK;
+ try {
+ TopExp_Explorer aFaceExplorer(aShape, TopAbs_FACE);
+ for (; aFaceExplorer.More(); aFaceExplorer.Next())
+ {
+ TopoDS_Shape aFace = aFaceExplorer.Current();
+ // loop on wires
+ TopoDS_Iterator aWireItr(aFace);
+ for (; aWireItr.More(); aWireItr.Next() )
+ {
+ Standard_Boolean wasOk = Standard_False;
+ TopoDS_Wire aSrcWire = TopoDS::Wire(aWireItr.Value());
+
+ aWireFix->Load (aSrcWire);
+ aWireFix->SetFace (TopoDS::Face(aFace));
+ aWireFix->FixReorder(); //correct order is a prerequisite
+ // fix 3d
+ if (!aWireFix->FixGaps3d())
+ {
+ // not fixed, why?
+ aStatus = getStatusGap(aWireFix, Standard_True);
+ if (aStatus == ShapeExtend_OK)
+ wasOk = Standard_True;
+ else
+ {
+ // keep 3d fail status
+ if (aStatusNbDMap.IsBound (aStatus))
+ aStatusNbDMap(aStatus)++;
+ else
+ aStatusNbDMap.Bind(aStatus,1);
+ continue;
+ }
+ }
+
+ // fix 2d
+ if (aWireFix->FixGaps2d())
+ nbFixed++;
+ else
+ {
+ aStatus = getStatusGap(aWireFix, Standard_False);
+ if (aStatus == ShapeExtend_OK)
+ {
+ if (wasOk)
+ {
+ nbOk++;
+ continue;
+ }
+ else
+ nbFixed++;
+ }
+ else
+ {
+ // keep 2d fail status
+ Standard_Integer aStatus2d = aStatus + ShapeExtend_FAIL;
+ if (aStatusNbDMap.IsBound (aStatus2d))
+ aStatusNbDMap(aStatus2d)++;
+ else
+ aStatusNbDMap.Bind(aStatus2d,1);
+ continue;
+ }
+ }
+ }
+ }
+//End: STEP 7
+ } catch (EXCEPTION) {
+ di << "Exception is raised = " <<aStatus << "\n";
+ return 1;
+
+ }
+// report what is done
+
+ if (nbFixed)
+ {
+ di <<"Fix_FillGaps_Fixed: nbFixed = "<<nbFixed <<"\n";
+
+ }
+ if (nbOk)
+ {
+ di << "Fix_FillGaps_NothingToDo" <<"\n";
+
+ }
+ TColStd_DataMapIteratorOfDataMapOfIntegerInteger aStatusItr(aStatusNbDMap);
+ for (; aStatusItr.More(); aStatusItr.Next())
+ {
+ switch ((ShapeExtend_Status) aStatusItr.Key())
+ {
+ // treat 3d status
+ case ShapeExtend_FAIL1:
+ di <<"Fix_FillGaps_3dNoCurveFail, nb failed = ";
+ break;
+ case ShapeExtend_FAIL2:
+ di <<"Fix_FillGaps_3dSomeGapsFail, nb failed = ";
+ break;
+ default:
+ // treat 2d status
+ switch ((ShapeExtend_Status) (aStatusItr.Key() - ShapeExtend_FAIL))
+ {
+ case ShapeExtend_FAIL1:
+ di <<"Fix_FillGaps_2dNoPCurveFail, nb failed = ";
+ break;
+ case ShapeExtend_FAIL2:
+ di <<"Fix_FillGaps_2dSomeGapsFail, nb failed = ";
+ break;
+ default:
+ break;
+ }
+ }
+ di <<aStatusItr.Value()<< "\n";
+ }
+ di << ("__________________________________") <<"\n";
+
+ return 0;
+}
+
+#include <IntTools_Context.hxx>
+#include <GeomAPI_ProjectPointOnSurf.hxx>
+
+//=======================================================================
+//function : xprojponf
+//purpose :
+//=======================================================================
+Standard_Integer xprojponf (Draw_Interpretor& di,
+ Standard_Integer n,
+ const char** a)
+{
+ if (n!=3) {
+ di<<" use xprojponf p f \n";
+ return 0;
+ }
+ //
+ gp_Pnt aP, aPS;
+ TopoDS_Shape aS;
+ TopoDS_Face aF;
+ Handle(IntTools_Context) aCtx;
+ //
+ DrawTrSurf::GetPoint(a[1], aP);
+ aS=DBRep::Get(a[2]);
+ //
+ if (aS.IsNull()) {
+ di<<" null shape is not allowed\n";
+ return 0;
+ }
+ //
+ if (aS.ShapeType()!=TopAbs_FACE) {
+ di << a[2] << " not a face\n";
+ return 0;
+ }
+ //
+ aCtx=new IntTools_Context;
+ //
+ aF=TopoDS::Face(aS);
+ GeomAPI_ProjectPointOnSurf& aPPS=aCtx->ProjPS(aF);
+ //
+ aPPS.Perform(aP);
+ if (!aPPS.IsDone()) {
+ di<<" projection failed\n";
+ return 0;
+ }
+ //
+ aPS=aPPS.NearestPoint();
+ di<< " point px " << aPS.X() << " " << aPS.Y() << " " << aPS.Z() << "\n";
+ //
+ return 0;
+}
+
void QABugs::Commands_19(Draw_Interpretor& theCommands) {
const char *group = "QABugs";
theCommands.Add ("OCC24755", "OCC24755", __FILE__, OCC24755, group);
theCommands.Add ("OCC24834", "OCC24834", __FILE__, OCC24834, group);
theCommands.Add ("OCC24889", "OCC24889", __FILE__, OCC24889, group);
- theCommands.Add ("OCC23951", "OCC23951", __FILE__, OCC23951, group);
- theCommands.Add ("OCC24931", "OCC24931", __FILE__, OCC24931, group);
+ theCommands.Add ("OCC23951", "OCC23951 path to saved step file", __FILE__, OCC23951, group);
+ theCommands.Add ("OCC24931", "OCC24931 path to saved xml file", __FILE__, OCC24931, group);
theCommands.Add ("OCC24945", "OCC24945", __FILE__, OCC24945, group);
theCommands.Add ("OCC23950", "OCC23950 step_file", __FILE__, OCC23950, group);
theCommands.Add ("OCC25004", "OCC25004", __FILE__, OCC25004, group);
__FILE__, OCC24925, group);
theCommands.Add ("OCC23010", "OCC23010 STEP_file", __FILE__, OCC23010, group);
theCommands.Add ("OCC25043", "OCC25043 shape", __FILE__, OCC25043, group);
+ theCommands.Add ("OCC24826,", "This test performs simple saxpy test.\n Usage: OCC24826 length", __FILE__, OCC24826, group);
theCommands.Add ("OCC24606", "OCC24606 : Tests ::FitAll for V3d view ('vfit' is for NIS view)", __FILE__, OCC24606, group);
theCommands.Add ("OCC25202", "OCC25202 res shape numF1 face1 numF2 face2", __FILE__, OCC25202, group);
+ theCommands.Add ("OCC7570", "OCC7570 shape", __FILE__, OCC7570, group);
+ theCommands.Add ("OCC25100", "OCC25100 shape", __FILE__, OCC25100, group);
+ theCommands.Add ("OCC25340", "OCC25340", __FILE__, OCC25340, group);
+ theCommands.Add ("OCC25348", "OCC25348", __FILE__, OCC25348, group);
+ theCommands.Add ("OCC25413", "OCC25413 shape", __FILE__, OCC25413, group);
+ theCommands.Add ("OCC25446", "OCC25446 res b1 b2 op", __FILE__, OCC25446, group);
+ theCommands.Add ("OCC25545",
+ "no args; tests data race when concurrently accessing \n"
+ "\t\tTopLoc_Location::Transformation()",
+ __FILE__, OCC25545, group);
+ theCommands.Add ("OCC25547", "OCC25547", __FILE__, OCC25547, group);
+ theCommands.Add ("OCC24881", "OCC24881 shape", __FILE__, OCC24881, group);
+ theCommands.Add ("xprojponf", "xprojponf p f", __FILE__, xprojponf, group);
return;
}
projects / occt.git / blobdiff