#include <TColStd_Array2OfReal.hxx>
#include <TColStd_Array1OfReal.hxx>
#include <TColStd_Array1OfInteger.hxx>
+#include <Geom_BezierCurve.hxx>
#include <Geom_BSplineSurface.hxx>
+#include <GeomConvert.hxx>
#include <Geom2d_Curve.hxx>
#include <Geom2d_Line.hxx>
+#include <GeomFill_BSplineCurves.hxx>
#include <Draw.hxx>
#include <DrawTrSurf.hxx>
#include <ShapeConstruct_ProjectCurveOnSurface.hxx>
#include <XCAFDoc_ShapeTool.hxx>
#include <HLRAppli_ReflectLines.hxx>
+#include <HLRBRep_PolyHLRToShape.hxx>
+#include <HLRBRep_PolyAlgo.hxx>
+
+#include <Standard_Failure.hxx>
+
+#include <Bnd_OBB.hxx>
+#include <BRepBndLib.hxx>
+#include <OSD_Timer.hxx>
+
+#include <limits>
//=======================================================================
//function : SurfaceGenOCC26675_1
// The test is based of occupying of all available virtual memory.
// Obviously it has no sense on 64-bit platforms.
- enum Status
+ enum AllocTestStatus
{
NotApplicable = 0x1,
OUMCatchOK = 0x2,
for (;;)
aList.Append(Standard::Allocate(aBlockSizes[i]));
}
- catch (Standard_Failure)
+ catch (Standard_Failure const&)
{
aStatus |= OUMCatchOK;
}
{
if (theNArg != 2)
{
- cout << "Use: " << theArgVal[0] << " shape" << endl;
+ std::cout << "Use: " << theArgVal[0] << " shape" << std::endl;
return 1;
}
std::pair<gp_Pnt, gp_Pnt> vertices;
clock_t t = clock();
- theDI << "\nRetrieving " << iterations << " vertices using approach A)...";
+ theDI << "\nRetrieving " << iterations << " vertices using approach A)...\n";
for (int i = 0; i < iterations; ++i)
{
vertices = getVerticesA(edge);
}
- theDI << "done in " << (clock() - t) / (double)CLOCKS_PER_SEC << " seconds\n";
+ theDI << "COUNTER RetrievingVertA" << ": " << (clock() - t) / (double)CLOCKS_PER_SEC << "\n";
t = clock();
- theDI << "\nRetrieving " << iterations << " vertices using approach B)...";
+ theDI << "\nRetrieving " << iterations << " vertices using approach B)...\n";
for (int i = 0; i < iterations; ++i)
{
vertices = getVerticesB(edge);
}
- theDI << "done in " << (clock() - t) / (double)CLOCKS_PER_SEC << " seconds\n";
+ theDI << "COUNTER RetrievingVertB" << ": " << (clock() - t) / (double)CLOCKS_PER_SEC << "\n";
return 0;
}
{
if (theNArg != 5)
{
- cout << "Use: " << theArgVal[0] <<" surface curve start end" << endl;
+ std::cout << "Use: " << theArgVal[0] <<" surface curve start end" << std::endl;
return 1;
}
return 0;
}
-//=======================================================================
-//function : OCC27341
-//purpose : check exact HLR algorighm's work
-//=======================================================================
-static Standard_Integer OCC27341 (Draw_Interpretor& , Standard_Integer n, const char** a)
-{
- if (n != 4)
- {
- cout << "Use: OCC27341 res shape axo/top/bottom/front/back/left/right" << endl;
- return 1;
- }
-
- TopoDS_Shape aShape = DBRep::Get(a[2]);
- if (aShape.IsNull())
- return 1;
-
- gp_Pnt anOrigin(0.,0.,0.);
- gp_Dir aNormal(0.57735026918962573, -0.57735026918962573, 0.57735026918962573);
- gp_Ax2 anAxes(anOrigin, aNormal);
- gp_Dir aDX = anAxes.XDirection();
-
- HLRAppli_ReflectLines Reflector(aShape);
-
- if (strcmp(a[3],"axo") == 0)
- {
- aNormal.SetCoord(0.57735026918962573, -0.57735026918962573, 0.57735026918962573);
- aDX.SetCoord(-0.40824829046386307, 0.40824829046386307, 0.81649658092772615);
- }
- else if (strcmp(a[3],"top") == 0)
- {
- aNormal.SetCoord(0,0,1);
- aDX.SetCoord(0,1,0);
- }
- else if (strcmp(a[3],"bottom") == 0)
- {
- aNormal.SetCoord(0,0,-1);
- aDX.SetCoord(0,-1,0);
- }
- else if (strcmp(a[3],"front") == 0)
- {
- aNormal.SetCoord(0,-1,0);
- aDX.SetCoord(0,0,1);
- }
- else if (strcmp(a[3],"back") == 0)
- {
- aNormal.SetCoord(0,1,0);
- aDX.SetCoord(0,0,1);
- }
- else if (strcmp(a[3],"left") == 0)
- {
- aNormal.SetCoord(-1,0,0);
- aDX.SetCoord(0,0,1);
- }
- else if (strcmp(a[3],"right") == 0)
- {
- aNormal.SetCoord(1,0,0);
- aDX.SetCoord(0,0,1);
- }
-
- Reflector.SetAxes(aNormal.X(), aNormal.Y(), aNormal.Z(),
- anOrigin.X(), anOrigin.Y(), anOrigin.Z(),
- aDX.X(), aDX.Y(), aDX.Z());
-
- Reflector.Perform();
-
- TopoDS_Compound Result;
- BRep_Builder BB;
- BB.MakeCompound(Result);
-
- TopoDS_Shape SharpEdges = Reflector.GetCompoundOf3dEdges(HLRBRep_Sharp, Standard_True, Standard_False);
- if (!SharpEdges.IsNull())
- BB.Add(Result, SharpEdges);
- TopoDS_Shape OutLines = Reflector.GetCompoundOf3dEdges(HLRBRep_OutLine, Standard_True, Standard_False);
- if (!OutLines.IsNull())
- BB.Add(Result, OutLines);
- TopoDS_Shape SmoothEdges = Reflector.GetCompoundOf3dEdges(HLRBRep_Rg1Line, Standard_True, Standard_False);
- if (!SmoothEdges.IsNull())
- BB.Add(Result, SmoothEdges);
-
- DBRep::Set(a[1], Result);
-
- return 0;
-}
-
//=======================================================================
//function : OCC27466
//purpose :
{
if (theNArg != 4)
{
- cout << "Use: " << theArgVal[0] << " face point start_pnt2d" << endl;
+ std::cout << "Use: " << theArgVal[0] << " face point start_pnt2d" << std::endl;
return 1;
}
Geom2dAdaptor_Curve(normalLine),1e-9,g1,g1,g3);
aDuumyList.Append(aCircleBuilder.NbSolutions());
}
- catch(Standard_Failure)
+ catch(Standard_Failure const&)
{
theDI << "Exception was caught\n";
}
DBRep::Set(theArgVal[2], bsp_face1);
}
}
- catch (Standard_Failure)
+ catch (Standard_Failure const&)
{
theDI << "ERROR: Exception in GeomFill_NSections\n";
}
}
+#include <TDF_Tool.hxx>
+#include <XCAFDoc_View.hxx>
+#include <XCAFDoc_ViewTool.hxx>
+#include <XCAFView_Object.hxx>
+#include <XCAFView_ProjectionType.hxx>
+static Standard_Integer OCC28389(Draw_Interpretor& di, Standard_Integer argc, const char** argv)
+{
+ if (argc < 20) {
+ di << "Use: OCC28389 Doc label nb_shapes nb_GDT nb_planes name type pp_x pp_y pp_z vd_x vd_y vd_z ud_x ud_y ud_z zoom width height";
+ return 1;
+ }
+ Handle(TDocStd_Document) aDoc;
+ DDocStd::GetDocument(argv[1], aDoc);
+ if (aDoc.IsNull()) {
+ di << "Error: Wrong document";
+ return 1;
+ }
+ Handle(XCAFDoc_ViewTool) aViewTool = XCAFDoc_DocumentTool::ViewTool(aDoc->Main());
+
+ TDF_Label aLabel;
+ TDF_Tool::Label(aDoc->GetData(), argv[2], aLabel);
+ if (aLabel.IsNull()) {
+ di << "Error: Wrong label";
+ return 1;
+ }
+ Handle(XCAFDoc_View) aView;
+ if (!aLabel.FindAttribute(XCAFDoc_View::GetID(), aView)) {
+ di << "Error: Wrong label";
+ return 1;
+ }
+ Handle(XCAFView_Object) anObj = aView->GetObject();
+ if (anObj.IsNull()) {
+ di << "Error: Wrong label";
+ return 1;
+ }
+
+ Standard_Boolean isOK = Standard_True;
+ // check links
+ Standard_Integer nbShapes = Draw::Atoi(argv[3]);
+ Standard_Integer nbGDTs = Draw::Atoi(argv[4]);
+ Standard_Integer nbPlanes = Draw::Atoi(argv[5]);
+ TDF_LabelSequence aSequence;
+ aViewTool->GetRefShapeLabel(aLabel, aSequence);
+ if (aSequence.Length() != nbShapes)
+ isOK = Standard_False;
+ aSequence.Clear();
+ aViewTool->GetRefGDTLabel(aLabel, aSequence);
+ if (aSequence.Length() != nbGDTs)
+ isOK = Standard_False;
+ aSequence.Clear();
+ aViewTool->GetRefClippingPlaneLabel(aLabel, aSequence);
+ if (aSequence.Length() != nbPlanes)
+ isOK = Standard_False;
+ if (!isOK) {
+ di << "Error: Wrong references";
+ return 1;
+ }
+
+ if (anObj->Name()->IsDifferent(new TCollection_HAsciiString(argv[6]))) {
+ di << "Error: Wrong name";
+ return 1;
+ }
+
+ XCAFView_ProjectionType aType = XCAFView_ProjectionType_NoCamera;
+ if (argv[7][0] == 'p')
+ aType = XCAFView_ProjectionType_Parallel;
+ else if (argv[7][0] == 'c')
+ aType = XCAFView_ProjectionType_Central;
+
+ if (anObj->Type()!= aType) {
+ di << "Error: Wrong type";
+ return 1;
+ }
+
+ gp_Pnt aPP(Draw::Atof(argv[8]), Draw::Atof(argv[9]), Draw::Atof(argv[10]));
+ if (aPP.Distance(anObj->ProjectionPoint()) > Precision::Confusion()) {
+ di << "Error: Wrong projection point";
+ return 1;
+ }
+
+ gp_Dir aVD(Draw::Atof(argv[11]), Draw::Atof(argv[12]), Draw::Atof(argv[13]));
+ if (!aVD.IsEqual(anObj->ViewDirection(), Precision::Angular())) {
+ di << "Error: Wrong view direction";
+ return 1;
+ }
+
+ gp_Dir aUD(Draw::Atof(argv[14]), Draw::Atof(argv[15]), Draw::Atof(argv[16]));
+ if (!aUD.IsEqual(anObj->UpDirection(), Precision::Angular())) {
+ di << "Error: Wrong up direction";
+ return 1;
+ }
+
+ if (fabs(anObj->ZoomFactor() - Draw::Atof(argv[17])) > Precision::Confusion()) {
+ di << "Error: Wrong zoom factor";
+ return 1;
+ }
+
+ if (fabs(anObj->WindowHorizontalSize() - Draw::Atof(argv[18])) > Precision::Confusion())
+ isOK = Standard_False;
+ if (fabs(anObj->WindowVerticalSize() - Draw::Atof(argv[19])) > Precision::Confusion())
+ isOK = Standard_False;
+ if (!isOK) {
+ di << "Error: Wrong Window size";
+ return 1;
+ }
+
+ di << argv[2] << " OK";
+ return 0;
+}
+
+#include <TColgp_HArray1OfPnt2d.hxx>
+#include <TColgp_Array1OfVec2d.hxx>
+#include <TColStd_HArray1OfBoolean.hxx>
+#include <Geom2d_BSplineCurve.hxx>
+#include <Geom2dAPI_Interpolate.hxx>
+#include <GeomAPI.hxx>
+#include <BRepBuilderAPI_MakeEdge2d.hxx>
+
+static Standard_Integer OCC28594(Draw_Interpretor& di, Standard_Integer argc, const char** argv)
+{
+ if (argc != 3)
+ {
+ di << "Usage :" << argv[0] << " curve_with_scale curve_without_scale\n";
+ return 0;
+ }
+ Handle(TColgp_HArray1OfPnt2d) points_2d = new TColgp_HArray1OfPnt2d(1, 6);
+ (*points_2d)(1) = gp_Pnt2d(-30.4, 8);
+ (*points_2d)(2) = gp_Pnt2d(-16.689912, 17.498217);
+ (*points_2d)(3) = gp_Pnt2d(-23.803064, 24.748543);
+ (*points_2d)(4) = gp_Pnt2d(-16.907466, 32.919615);
+ (*points_2d)(5) = gp_Pnt2d(-8.543829, 26.549421);
+ (*points_2d)(6) = gp_Pnt2d(0, 39.200000);
+
+ TColgp_Array1OfVec2d tangent_2d(1, 6);
+ (tangent_2d)(1) = gp_Vec2d(0.3, 0.4);
+ (tangent_2d)(2) = gp_Vec2d(0, 0);
+ (tangent_2d)(3) = gp_Vec2d(0, 0);
+ (tangent_2d)(4) = gp_Vec2d(0, 0);
+ (tangent_2d)(5) = gp_Vec2d(0, 0);
+ (tangent_2d)(6) = gp_Vec2d(1, 0);
+
+ Handle(TColStd_HArray1OfBoolean) tangent_flags = new TColStd_HArray1OfBoolean(1, 6);
+ (*tangent_flags)(1) = true;
+ (*tangent_flags)(2) = false;
+ (*tangent_flags)(3) = false;
+ (*tangent_flags)(4) = false;
+ (*tangent_flags)(5) = false;
+ (*tangent_flags)(6) = true;
+
+ Geom2dAPI_Interpolate interp_2d_with_scale(points_2d, Standard_False, Precision::Confusion());
+ interp_2d_with_scale.Load(tangent_2d, tangent_flags);
+ interp_2d_with_scale.Perform();
+ Handle(Geom2d_BSplineCurve) curve_2d_with_scale = interp_2d_with_scale.Curve();
+
+ Geom2dAPI_Interpolate interp_2d_without_scale(points_2d, Standard_False, Precision::Confusion());
+ interp_2d_without_scale.Load(tangent_2d, tangent_flags, Standard_False);
+ interp_2d_without_scale.Perform();
+ Handle(Geom2d_BSplineCurve) curve_2d_without_scale = interp_2d_without_scale.Curve();
+
+ DrawTrSurf::Set(argv[1], curve_2d_with_scale);
+ DrawTrSurf::Set(argv[2], curve_2d_without_scale);
+ return 0;
+}
+
+static Standard_Integer OCC28784(Draw_Interpretor&, Standard_Integer argc, const char** argv)
+{
+ if (argc < 3)
+ return 1;
+
+ TopoDS_Shape aShape = DBRep::Get(argv[2]);
+ if (aShape.IsNull())
+ return 1;
+
+ gp_Ax2 aPlane (gp::Origin(), gp::DX(), -gp::DZ());
+ HLRAlgo_Projector aProjector(aPlane);
+
+ Handle(HLRBRep_PolyAlgo) aHLR = new HLRBRep_PolyAlgo(aShape);
+ aHLR->Projector(aProjector);
+ aHLR->Update();
+
+ HLRBRep_PolyHLRToShape aHLRtoShape;
+ aHLRtoShape.Update(aHLR);
+
+ TopoDS_Shape aHidden = aHLRtoShape.HCompound();
+
+ DBRep::Set(argv[1], aHidden);
+
+ return 0;
+}
+
+static Standard_Integer OCC28829 (Draw_Interpretor&, Standard_Integer, const char**)
+{
+ // do something that causes FPE exception
+ std::cout << "sqrt(-1) = " << sqrt (-1.) << std::endl;
+ return 0;
+}
+
+#include <NCollection_Buffer.hxx>
+#include <DDocStd_DrawDocument.hxx>
+#include <OSD_OpenFile.hxx>
+#include <Standard_ArrayStreamBuffer.hxx>
+#include <TDataStd_Name.hxx>
+#include <TDocStd_Application.hxx>
+
+#ifdef max
+ #undef max
+#endif
+
+static Standard_Integer OCC28887 (Draw_Interpretor&, Standard_Integer theNbArgs, const char** theArgVec)
+{
+ if (theNbArgs < 3)
+ {
+ std::cout << "Syntax error: wrong number of arguments!\n";
+ return 1;
+ }
+
+ const TCollection_AsciiString aFilePath (theArgVec[1]);
+ const TCollection_AsciiString aName (theArgVec[2]);
+ Handle(NCollection_Buffer) aBuffer;
+ {
+ std::ifstream aFile;
+ OSD_OpenStream (aFile, aFilePath.ToCString(), std::ios::binary | std::ios::in);
+ if (!aFile.is_open())
+ {
+ std::cout << "Error: input file '" << aFilePath << "' cannot be read\n";
+ return 1;
+ }
+ aFile.seekg (0, std::ios_base::end);
+ const int64_t aFileLength = int64_t (aFile.tellg());
+ if (aFileLength > int64_t (std::numeric_limits<ptrdiff_t>::max())
+ || aFileLength < 1)
+ {
+ std::cout << "Error: input file '" << aFilePath << "' is too large\n";
+ return 1;
+ }
+ aFile.seekg (0, std::ios_base::beg);
+
+ aBuffer = new NCollection_Buffer (NCollection_BaseAllocator::CommonBaseAllocator());
+ if (!aBuffer->Allocate (size_t(aFileLength)))
+ {
+ std::cout << "Error: memory allocation (" << aFileLength << ") has failed\n";
+ return 1;
+ }
+
+ aFile.read ((char* )aBuffer->ChangeData(), aBuffer->Size());
+ if (!aFile.good())
+ {
+ std::cout << "Error: input file '" << aFilePath << "' reading failure\n";
+ return 1;
+ }
+ }
+
+ Standard_ArrayStreamBuffer aStreamBuffer ((const char* )aBuffer->ChangeData(), aBuffer->Size());
+ std::istream aStream (&aStreamBuffer);
+ // just play with seeking
+ aStream.seekg (0, std::ios_base::end);
+ aStream.seekg (0, std::ios_base::beg);
+ if (aFilePath.EndsWith (".brep")
+ || aFilePath.EndsWith (".rle"))
+ {
+ TopoDS_Shape aShape;
+ BRep_Builder aBuilder;
+ BRepTools::Read (aShape, aStream, aBuilder);
+ DBRep::Set (aName.ToCString(), aShape);
+ }
+ else
+ {
+ Handle(TDocStd_Document) aDoc;
+ Handle(TDocStd_Application) anApp = DDocStd::GetApplication();
+ Standard_CString aNameVar = aName.ToCString();
+ if (DDocStd::GetDocument (aNameVar, aDoc, Standard_False))
+ {
+ std::cout << "Error: document with name " << aName << " already exists\n";
+ return 1;
+ }
+
+ if (anApp->Open (aStream, aDoc) != PCDM_RS_OK)
+ {
+ std::cout << "Error: cannot open XDE document\n";
+ return 1;
+ }
+
+ Handle(DDocStd_DrawDocument) aDrawDoc = new DDocStd_DrawDocument (aDoc);
+ TDataStd_Name::Set (aDoc->GetData()->Root(), aName.ToCString());
+ Draw::Set (aName.ToCString(), aDrawDoc);
+ }
+
+ return 0;
+}
+
+static Standard_Integer OCC28131 (Draw_Interpretor&, Standard_Integer theNbArgs, const char** theArgVec)
+{
+ if (theNbArgs != 2)
+ {
+ std::cerr << "Error: wrong number of arguments" << std::endl;
+ return 1;
+ }
+
+ double height = 8.5;
+ gp_Pnt JiZhunXian2_v0 = gp_Pnt(-17.6, 0.0, 0.0);
+ gp_Pnt JiZhunXian2_v1 = gp_Pnt(0, 32.8, 0.0);
+
+ // Outline
+ TColgp_Array1OfPnt outer_e_bzr_geom_v(1, 4);
+ {
+ outer_e_bzr_geom_v(1) = JiZhunXian2_v0;
+ outer_e_bzr_geom_v(4) = JiZhunXian2_v1;
+
+ Standard_Real ratio1 = 5.4 / 13.2;
+ outer_e_bzr_geom_v(2) = gp_Pnt(outer_e_bzr_geom_v(1).X(), ratio1*outer_e_bzr_geom_v(4).Y(), 0);
+ Standard_Real ratio2 = 6.0 / 6.8;
+ outer_e_bzr_geom_v(3) = gp_Pnt(ratio2*outer_e_bzr_geom_v(1).X(), outer_e_bzr_geom_v(4).Y(), 0);
+ }
+
+ Handle(Geom_BezierCurve) outer_e_bzr_geom = new Geom_BezierCurve(outer_e_bzr_geom_v);
+ Handle(Geom_BSplineCurve) outer_e_bsp_geom = GeomConvert::CurveToBSplineCurve(outer_e_bzr_geom);
+ TopoDS_Edge outer_e = BRepBuilderAPI_MakeEdge(outer_e_bsp_geom);
+
+ Handle(Geom_BSplineCurve) curve1;
+ {
+ Handle(TColgp_HArray1OfPnt2d) harray = new TColgp_HArray1OfPnt2d(1, 2); // sizing harray
+ harray->SetValue(1, gp_Pnt2d(-JiZhunXian2_v1.Y(), 0));
+ harray->SetValue(2, gp_Pnt2d(0, height + height / 2));
+
+ Geom2dAPI_Interpolate anInterpolation(harray, Standard_False, 1e-6);
+
+ gp_Vec2d vtangent1(0, 1);
+ gp_Vec2d vtangent2(1, 0);
+ anInterpolation.Load(vtangent1, vtangent2);
+ anInterpolation.Perform();
+
+ Handle(Geom2d_BSplineCurve) c = anInterpolation.Curve();
+
+ gp_Pln pln(gp_Ax3(gp_Pnt(), gp_Dir(1, 0, 0), gp_Dir(0, -1, 0)));
+
+ Handle(Geom_BSplineCurve) c3d = Handle(Geom_BSplineCurve)::DownCast(GeomAPI::To3d(c, pln));
+ curve1 = c3d;
+ }
+
+ Handle(Geom_BSplineCurve) curve2;
+ {
+ Handle(TColgp_HArray1OfPnt2d) harray = new TColgp_HArray1OfPnt2d(1, 3); // sizing harray
+ harray->SetValue(1, gp_Pnt2d(-JiZhunXian2_v0.X(), 0));
+ harray->SetValue(2, gp_Pnt2d(-JiZhunXian2_v0.X() - 2.6, height));
+ harray->SetValue(3, gp_Pnt2d(0, height + height / 2));
+
+ Geom2dAPI_Interpolate anInterpolation(harray, Standard_False, 1e-6);
+ anInterpolation.Perform();
+
+ Handle(Geom2d_BSplineCurve) c = anInterpolation.Curve();
+ gp_Pln pln(gp_Ax3(gp_Pnt(), gp_Dir(0, -1, 0), gp_Dir(-1, 0, 0)));
+ Handle(Geom_BSplineCurve) c3d = Handle(Geom_BSplineCurve)::DownCast(GeomAPI::To3d(c, pln));
+ curve2 = c3d;
+ }
+
+ //////////////////////////////////////
+ GeomFill_BSplineCurves fill2;
+ fill2.Init(outer_e_bsp_geom, curve1, curve2, GeomFill_CoonsStyle);
+
+ const Handle(Geom_BSplineSurface)& surf_geom = fill2.Surface();
+
+ TopoDS_Shape filled_face = BRepBuilderAPI_MakeFace(surf_geom, 0);
+
+ DBRep::Set (theArgVec[1], filled_face);
+
+/*
+ ///////////////////////////////////////////////////////////////////////
+ TopoDS_Solid first_solid;
+ {
+ BRepOffset_MakeOffset myOffsetShape(filled_face, -offset_thick, 1e-4,
+ BRepOffset_Skin, //Mode
+ Standard_False, //Intersection
+ Standard_False, //SelfInter
+ GeomAbs_Intersection, //Join
+ Standard_True, //Thickening
+ Standard_False //RemoveIntEdges
+ ); //RemoveInvalidFaces
+ first_solid = TopoDS::Solid(myOffsetShape.Shape());
+ }
+*/
+ return 0;
+}
+#include <math_NewtonFunctionRoot.hxx>
+#include <math_TrigonometricFunctionRoots.hxx>
+#include <math_TrigonometricEquationFunction.hxx>
+#include <gp_Elips2d.hxx>
+#include <Geom2d_Ellipse.hxx>
+#include <Geom2dAPI_InterCurveCurve.hxx>
+static Standard_Integer OCC29289(Draw_Interpretor&, Standard_Integer , const char** )
+{
+ gp_Elips2d e1(gp_Ax2d(gp_Pnt2d(0., 0.), gp_Dir2d(1., 0)), 2., 1.);
+ Handle(Geom2d_Ellipse) Ge1 = new Geom2d_Ellipse(e1);
+ gp_Elips2d e2(gp_Ax2d(gp_Pnt2d(0.5, 0.5), gp_Dir2d(1., 1.)), 2., 1.);
+ Handle(Geom2d_Ellipse) Ge2 = new Geom2d_Ellipse(e2);
+
+ Standard_Integer err = 0;
+ Geom2dAPI_InterCurveCurve Intersector;
+ Intersector.Init(Ge1, Ge2, 1.e-7);
+ if (Intersector.NbPoints() == 0)
+ {
+ std::cout << "Error: intersector is not done \n";
+ err = 1;
+ }
+
+
+ Standard_Real A, B, C, D, E;
+ A = 1.875;
+ B = -.75;
+ C = -.5;
+ D = -.25;
+ E = -.25;
+ math_TrigonometricEquationFunction MyF(A, B, C, D, E);
+ Standard_Real X, Tol1, Eps, Teta, TetaNewton;
+ Tol1 = 1.e-15;
+ Eps = 1.5e-12;
+ Standard_Integer Nit[] = { 5, 6, 7, 6 };
+
+ Standard_Real TetaPrev = 0.;
+ Standard_Integer i;
+ for (i = 1; i <= Intersector.NbPoints(); i++) {
+ Teta = Intersector.Intersector().Point(i).ParamOnFirst();
+ X = Teta - 0.1 * (Teta - TetaPrev);
+ TetaPrev = Teta;
+ math_NewtonFunctionRoot Resol(MyF, X, Tol1, Eps, Nit[i-1]);
+ if (Resol.IsDone()) {
+ TetaNewton = Resol.Root();
+ if (Abs(Teta - TetaNewton) > 1.e-7)
+ {
+ std::cout << "Error: Newton root is wrong for " << Teta << " \n";
+ err = 1;
+ }
+ }
+ else
+ {
+ std::cout << "Error: Newton is not done for " << Teta << " \n";
+ err = 1;
+ }
+ }
+
+ return err;
+}
+
+//===============================================================================================
+Standard_Boolean IsSameGuid (const Standard_GUID& aGuidNull, const Standard_GUID& aGuid2)
+{
+ Standard_Boolean isSame (Standard_False);
+ if(Standard_GUID::IsEqual(aGuidNull, aGuid2)) {
+ aGuid2.ShallowDump(std::cout);
+ isSame = Standard_True;
+ } else {
+ aGuid2.ShallowDump(std::cout);
+ std::cout <<std::endl;
+ }
+ return isSame;
+}
+
+#include <TDataStd_AsciiString.hxx>
+#include <TDataStd_BooleanArray.hxx>
+#include <TDataStd_BooleanList.hxx>
+#include <TDataStd_ByteArray.hxx>
+#include <TDataStd_ExtStringArray.hxx>
+#include <TDataStd_ExtStringList.hxx>
+#include <TDataStd_Integer.hxx>
+#include <TDataStd_IntegerArray.hxx>
+#include <TDataStd_IntegerList.hxx>
+#include <TDataStd_Name.hxx>
+#include <TDataStd_Real.hxx>
+#include <TDataStd_RealArray.hxx>
+#include <TDataStd_RealList.hxx>
+#include <TDataStd_ReferenceArray.hxx>
+#include <TDataStd_ReferenceList.hxx>
+
+#define QCOMPARE(val1, val2) \
+ di << "Checking " #val1 " == " #val2 << \
+ ((val1) == (val2) ? ": OK\n" : ": Error\n")
+
+static Standard_Integer OCC29371 (Draw_Interpretor& di, Standard_Integer n, const char** a)
+{
+ if (n != 1)
+ {
+ std::cout << "Usage : " << a[0] << "\n";
+ return 1;
+ }
+
+ Handle(TDocStd_Application) anApp = DDocStd::GetApplication();
+ Handle(TDocStd_Document) aDoc;
+ anApp->NewDocument ("BinOcaf", aDoc);
+ TDF_Label aLab = aDoc->Main();
+ Standard_GUID aNullGuid("00000000-0000-0000-0000-000000000000");
+ Standard_Boolean IsNullGuid(Standard_False);
+
+ try {
+ //1. Set TDataStd_AsciiString
+ Handle(TDataStd_AsciiString) aStrAtt = new TDataStd_AsciiString();
+ aLab.AddAttribute(aStrAtt);
+ if(!aStrAtt.IsNull()) {
+ Standard_GUID aGuid = aStrAtt->ID();
+ IsNullGuid = IsSameGuid(aNullGuid, aGuid);
+ }
+
+ //2. Set TDataStd_BooleanArray
+ Handle(TDataStd_BooleanArray) aBArAtt = new TDataStd_BooleanArray();
+ aLab.AddAttribute(aBArAtt);
+ if(!aBArAtt.IsNull()) {
+ Standard_GUID aGuid = aBArAtt->ID();
+ IsNullGuid = IsSameGuid(aNullGuid, aGuid);
+ }
+
+ //3. Set TDataStd_BooleanList
+ Handle(TDataStd_BooleanList) aBListAtt = new TDataStd_BooleanList();
+ aLab.AddAttribute(aBListAtt);
+ if(!aBListAtt.IsNull()) {
+ Standard_GUID aGuid = aBListAtt->ID();
+ IsNullGuid = IsSameGuid(aNullGuid, aGuid);
+ }
+
+ //4. Set TDataStd_ByteArray
+ Handle(TDataStd_ByteArray) aByteArAtt = new TDataStd_ByteArray();
+ aLab.AddAttribute(aByteArAtt);
+ if(!aByteArAtt.IsNull()) {
+ Standard_GUID aGuid = aByteArAtt->ID();
+ IsNullGuid = IsSameGuid(aNullGuid, aGuid);
+ }
+
+ //5. Set TDataStd_ExtStringArray
+ Handle(TDataStd_ExtStringArray) anExtStrArAtt = new TDataStd_ExtStringArray();
+ aLab.AddAttribute(anExtStrArAtt);
+ if(!anExtStrArAtt.IsNull()) {
+ Standard_GUID aGuid = anExtStrArAtt->ID();
+ IsNullGuid = IsSameGuid(aNullGuid, aGuid);
+ }
+
+ //6. Set TDataStd_ExtStringList
+ Handle(TDataStd_ExtStringList) anExtStrListAtt = new TDataStd_ExtStringList();
+ aLab.AddAttribute(anExtStrListAtt);
+ if(!anExtStrListAtt.IsNull()) {
+ Standard_GUID aGuid = anExtStrListAtt->ID();
+ IsNullGuid = IsSameGuid(aNullGuid, aGuid);
+ }
+
+ //7. Set TDataStd_Integer
+ Handle(TDataStd_Integer) anIntAtt = new TDataStd_Integer();
+ aLab.AddAttribute(anIntAtt);
+ if(!anIntAtt.IsNull()) {
+ Standard_GUID aGuid = anIntAtt->ID();
+ IsNullGuid = IsSameGuid(aNullGuid, aGuid);
+ }
+
+ //8. Set TDataStd_IntegerArray
+ Handle(TDataStd_IntegerArray) anIntArrAtt = new TDataStd_IntegerArray();
+ aLab.AddAttribute(anIntArrAtt);
+ if(!anIntArrAtt.IsNull()) {
+ Standard_GUID aGuid = anIntArrAtt->ID();
+ IsNullGuid = IsSameGuid(aNullGuid, aGuid);
+ }
+
+ //9. Set TDataStd_IntegerList
+ Handle(TDataStd_IntegerList) anIntListAtt = new TDataStd_IntegerList();
+ aLab.AddAttribute(anIntListAtt);
+ if(!anIntListAtt.IsNull()) {
+ Standard_GUID aGuid = anIntListAtt->ID();
+ IsNullGuid = IsSameGuid(aNullGuid, aGuid);
+ }
+
+ //10. Set TDataStd_Name
+ Handle(TDataStd_Name) aNameAtt = new TDataStd_Name();
+ aLab.AddAttribute(aNameAtt);
+ if(!aNameAtt.IsNull()) {
+ Standard_GUID aGuid = aNameAtt->ID();
+ IsNullGuid = IsSameGuid(aNullGuid, aGuid);
+ }
+
+ //11. Set TDataStd_Real
+ Handle(TDataStd_Real) aRealAtt = new TDataStd_Real();
+ aLab.AddAttribute(aRealAtt);
+ if(!aRealAtt.IsNull()) {
+ Standard_GUID aGuid = aRealAtt->ID();
+ IsNullGuid = IsSameGuid(aNullGuid, aGuid);
+ }
+
+ //12. Set TDataStd_RealArray
+ Handle(TDataStd_RealArray) aRealArrAtt = new TDataStd_RealArray();
+ aLab.AddAttribute(aRealArrAtt);
+ if(!aRealArrAtt.IsNull()) {
+ Standard_GUID aGuid = aRealArrAtt->ID();
+ IsNullGuid = IsSameGuid(aNullGuid, aGuid);
+ }
+
+ //13. Set TDataStd_RealList
+ Handle(TDataStd_RealList) aRealListAtt = new TDataStd_RealList();
+ aLab.AddAttribute(aRealListAtt);
+ if(!aRealListAtt.IsNull()) {
+ Standard_GUID aGuid = aRealListAtt->ID();
+ IsNullGuid = IsSameGuid(aNullGuid, aGuid);
+ }
+
+ //14. Set TDataStd_ReferenceArray
+ Handle(TDataStd_ReferenceArray) aRefArrAtt = new TDataStd_ReferenceArray();
+ aLab.AddAttribute(aRefArrAtt);
+ if(!aRefArrAtt.IsNull()) {
+ Standard_GUID aGuid = aRefArrAtt->ID();
+ IsNullGuid = IsSameGuid(aNullGuid, aGuid);
+ }
+
+ //15. Set TDataStd_ReferenceList
+ Handle(TDataStd_ReferenceList) aRefListAtt = new TDataStd_ReferenceList();
+ aLab.AddAttribute(aRefListAtt);
+ if(!aRefListAtt.IsNull()) {
+ Standard_GUID aGuid = aRefListAtt->ID();
+ IsNullGuid = IsSameGuid(aNullGuid, aGuid);
+ }
+ } catch (...)
+ {
+ IsNullGuid = Standard_True;
+ }
+ QCOMPARE (IsNullGuid, Standard_False);
+ anApp->Close(aDoc);
+ return 0;
+}
+
+#include <NCollection_DoubleMap.hxx>
+#include <NCollection_IndexedMap.hxx>
+#include <NCollection_DataMap.hxx>
+#include <NCollection_IndexedDataMap.hxx>
+#include <OSD_MemInfo.hxx>
+
+// check that copying of empty maps does not allocate extra memory
+template<typename T> void AllocDummyArr (Draw_Interpretor& theDI, int theN1, int theN2)
+{
+ NCollection_Array1<T> aMapArr1(0, theN1), aMapArr2(0, theN2);
+
+ OSD_MemInfo aMemTool;
+ Standard_Size aMem0 = aMemTool.Value (OSD_MemInfo::MemHeapUsage);
+
+ for (int i = 1; i < theN1; i++)
+ aMapArr1(i) = aMapArr1(i-1);
+ for (int i = 1; i < theN2; i++)
+ aMapArr2(i) = aMapArr2(0);
+
+ aMemTool.Update();
+ Standard_Size aMem1 = aMemTool.Value (OSD_MemInfo::MemHeapUsage);
+
+ theDI << "Heap usage before copy = " << (int)aMem0 << ", after = " << (int)aMem1 << "\n";
+
+ if (aMem1 > aMem0)
+ theDI << "Error: memory increased by " << (int)(aMem1 - aMem0) << " bytes\n";
+};
+
+static Standard_Integer OCC29064 (Draw_Interpretor& theDI, Standard_Integer theArgc, const char** theArgv)
+{
+ if (theArgc < 2)
+ {
+ std::cout << "Error: give argument indicating type of map (map, doublemap, datamap, indexedmap, indexeddatamap)" << std::endl;
+ return 1;
+ }
+
+ const int nbm1 = 10000, nbm2 = 10000;
+ if (strcasecmp (theArgv[1], "map") == 0)
+ AllocDummyArr<NCollection_Map<int> > (theDI, nbm1, nbm2);
+ else if (strcasecmp (theArgv[1], "doublemap") == 0)
+ AllocDummyArr<NCollection_DoubleMap<int, int> > (theDI, nbm1, nbm2);
+ else if (strcasecmp (theArgv[1], "datamap") == 0)
+ AllocDummyArr<NCollection_DataMap<int, int> > (theDI, nbm1, nbm2);
+ else if (strcasecmp (theArgv[1], "indexedmap") == 0)
+ AllocDummyArr<NCollection_IndexedMap<int> > (theDI, nbm1, nbm2);
+ else if (strcasecmp (theArgv[1], "indexeddatamap") == 0)
+ AllocDummyArr<NCollection_IndexedDataMap<int, int> > (theDI, nbm1, nbm2);
+ else
+ {
+ std::cout << "Error: unrecognized argument " << theArgv[1] << std::endl;
+ return 1;
+ }
+ return 0;
+}
+
+#include <BRepOffsetAPI_MakePipeShell.hxx>
+#include <GC_MakeArcOfCircle.hxx>
+#include <BRepAdaptor_CompCurve.hxx>
+#include <gp_Circ.hxx>
+//=======================================================================
+//function : OCC29430
+//purpose :
+//=======================================================================
+static Standard_Integer OCC29430(Draw_Interpretor& theDI,
+ Standard_Integer /*theNArg*/,
+ const char** theArgVal)
+{
+ const Standard_Real r45 = M_PI / 4.0, r225 = 3.0*M_PI / 4.0;
+
+ GC_MakeArcOfCircle arcMaker(gp_Circ(gp_Ax2(gp_Pnt(0.0, 0.0, 0.0), gp_Dir(0.0, 0.0, 1.0), gp_Dir(1.0, 0.0, 0.0)), 1.0), r45, r225, Standard_True);
+ BRepBuilderAPI_MakeEdge edgeMaker(arcMaker.Value());
+ BRepBuilderAPI_MakeWire wireMaker(edgeMaker.Edge());
+ const TopoDS_Wire circle = wireMaker.Wire();
+
+ DBRep::Set(theArgVal[1], circle);
+
+ BRepAdaptor_CompCurve curve(circle);
+ theDI << "Curve.FirstParameter() = " << curve.FirstParameter() << "\n";
+ theDI << "Curve.LastParameter() = " << curve.LastParameter() << "\n";
+ theDI << "Curve.Period() = " << (curve.IsPeriodic()? curve.Period() : 0.0) << "\n";
+ const gp_Pnt aStartPt = curve.Value(curve.FirstParameter());
+ const gp_Pnt anEndPt = curve.Value(curve.LastParameter());
+
+ DrawTrSurf::Set(theArgVal[2], aStartPt);
+ DrawTrSurf::Set(theArgVal[3], anEndPt);
+
+ return 0;
+}
+
+#include <STEPCAFControl_Reader.hxx>
+
+//=======================================================================
+//function : OCC29531
+//purpose :
+//=======================================================================
+static Standard_Integer OCC29531(Draw_Interpretor&, Standard_Integer, const char** theArgV)
+{
+ Handle(TDocStd_Application) anApp = DDocStd::GetApplication();
+ Handle(TDocStd_Document) aDoc;
+ anApp->NewDocument("BinOcaf", aDoc);
+ aDoc->SetUndoLimit(1);
+
+ STEPCAFControl_Reader Reader;
+ Reader.ReadFile(theArgV[1]);
+ Reader.Transfer(aDoc);
+ TDF_Label aShL, aDL;
+ TDF_Tool::Label(aDoc->GetData(), "0:1:1:2:672", aShL);
+ TDF_Tool::Label(aDoc->GetData(), "0:1:4:10", aDL);
+
+ aDoc->OpenCommand();
+
+ Handle(XCAFDoc_DimTolTool) aDimTolTool = XCAFDoc_DocumentTool::DimTolTool(aDoc->Main());
+ aDimTolTool->SetDimension(aShL, aDL);
+
+ aDoc->CommitCommand();
+
+ aDoc->Undo();
+ aDoc->Redo();
+ return 0;
+}
+
+//=======================================================================
+//function : OCC29807
+//purpose :
+//=======================================================================
+#include <GeomAdaptor_HSurface.hxx>
+#include <IntPatch_PointLine.hxx>
+#include <IntSurf_PntOn2S.hxx>
+static Standard_Integer OCC29807(Draw_Interpretor& theDI, Standard_Integer theNArg, const char** theArgV)
+{
+ if (theNArg != 7)
+ {
+ theDI << "Use: " << theArgV[0] << "surface1 surface2 u1 v1 u2 v2\n";
+ return 1;
+ }
+
+ const Handle(Geom_Surface) aS1 = DrawTrSurf::GetSurface(theArgV[1]);
+ const Handle(Geom_Surface) aS2 = DrawTrSurf::GetSurface(theArgV[2]);
+
+ if (aS1.IsNull() || aS2.IsNull())
+ {
+ theDI << "Error. Null surface is not supported.\n";
+ return 1;
+ }
+
+ const Standard_Real aU1 = Draw::Atof(theArgV[3]);
+ const Standard_Real aV1 = Draw::Atof(theArgV[4]);
+ const Standard_Real aU2 = Draw::Atof(theArgV[5]);
+ const Standard_Real aV2 = Draw::Atof(theArgV[6]);
+
+ const Handle(GeomAdaptor_HSurface) anAS1 = new GeomAdaptor_HSurface(aS1);
+ const Handle(GeomAdaptor_HSurface) anAS2 = new GeomAdaptor_HSurface(aS2);
+
+ const gp_Pnt aP1 = anAS1->Value(aU1, aV1);
+ const gp_Pnt aP2 = anAS2->Value(aU2, aV2);
+
+ if (aP1.SquareDistance(aP2) > Precision::SquareConfusion())
+ {
+ theDI << "Error. True intersection point must be specified. "
+ "Please check parameters: u1 v1 u2 v2.\n";
+ return 1;
+ }
+
+ IntSurf_PntOn2S aPOn2S;
+ aPOn2S.SetValue(0.5*(aP1.XYZ() + aP2.XYZ()), aU1, aV1, aU2, aV2);
+
+ const Standard_Real aCurvatureRadius = IntPatch_PointLine::CurvatureRadiusOfIntersLine(anAS1, anAS2, aPOn2S);
+ theDI << "Radius of curvature is " << aCurvatureRadius << "\n";
+ return 0;
+}
+
+//=======================================================================
+//function : OCC29925
+//purpose : check safety of functions like IsSpace(), LowerCase(), etc. for all chars
+//=======================================================================
+static Standard_Integer OCC29925 (Draw_Interpretor& theDI, Standard_Integer, const char**)
+{
+ // iterate by all valid ASCII chars (including extended)
+ for (int i = 0; i < 256; i++)
+ {
+ Standard_Character c = (char)(unsigned char)i;
+// if (c != i) theDI << c << " != " << i << "\n";
+ const char* anOp = "";
+ try {
+ anOp = "IsAlphabetic";
+ IsAlphabetic (c);
+ anOp = "IsDigit";
+ IsDigit (c);
+ anOp = "IsXDigit";
+ IsXDigit (c);
+ anOp = "IsAlphanumeric";
+ IsAlphanumeric (c);
+ anOp = "IsControl";
+ IsControl (c);
+ anOp = "IsGraphic";
+ IsGraphic (c);
+ anOp = "IsLowerCase";
+ IsLowerCase (c);
+ anOp = "IsPrintable";
+ IsPrintable (c);
+ anOp = "IsPunctuation";
+ IsPunctuation (c);
+ anOp = "IsSpace";
+ IsSpace (c);
+ anOp = "IsUpperCase";
+ IsUpperCase (c);
+ anOp = "LowerCase";
+ LowerCase (c);
+ anOp = "UpperCase";
+ UpperCase (c);
+ }
+ catch (const Handle(Standard_Failure)& e)
+ {
+ theDI << anOp << "() fails for " << c << " (" << e->DynamicType()->Name() << ")\n";
+ }
+ }
+
+ return 0;
+}
+
+//=======================================================================
+//function : OCC29311
+//purpose : check performance of OBB calculations
+//=======================================================================
+static Standard_Integer OCC29311 (Draw_Interpretor& theDI, Standard_Integer theArgc, const char** theArgv)
+{
+ if (theArgc < 4)
+ {
+ std::cerr << "Use: " << theArgv[0] << " shape counter_name nb_iterations" << std::endl;
+ return 1;
+ }
+
+ TopoDS_Shape aShape = DBRep::Get (theArgv[1]);
+ Standard_Integer aNbIter = Draw::Atoi (theArgv[3]);
+
+ Bnd_OBB anOBB;
+ OSD_Timer aTimer;
+ aTimer.Start ();
+ for (Standard_Integer aN = aNbIter; aN > 0; --aN)
+ {
+ anOBB.SetVoid ();
+ BRepBndLib::AddOBB (aShape, anOBB, Standard_False, Standard_False, Standard_False);
+ }
+ aTimer.Stop ();
+
+ theDI << "COUNTER " << theArgv[2] << ": " << aTimer.ElapsedTime() << "\n";
+
+ return 0;
+}
+
+//=======================================================================
+//function : OCC30391
+//purpose :
+//=======================================================================
+#include <BRepOffset_Tool.hxx>
+static Standard_Integer OCC30391(Draw_Interpretor& theDI,
+ Standard_Integer theNArg,
+ const char** theArgV)
+{
+ if (theNArg < 7)
+ {
+ theDI << "Use: " << theArgV[0] << "result face LenBeforeUfirst LenAfterUlast LenBeforeVfirst LenAfterVlast\n";
+ return 1;
+ }
+
+ TopoDS_Shape aShape = DBRep::Get(theArgV[2], TopAbs_FACE);
+ if (aShape.IsNull())
+ return 1;
+
+ const TopoDS_Face& aFace = TopoDS::Face(aShape);
+
+ Standard_Real aLenBeforeUfirst = atof(theArgV[3]);
+ Standard_Real aLenAfterUlast = atof(theArgV[4]);
+ Standard_Real aLenBeforeVfirst = atof(theArgV[5]);
+ Standard_Real aLenAfterVlast = atof(theArgV[6]);
+
+ TopoDS_Face Result;
+ BRepOffset_Tool::EnLargeFace(aFace, Result,
+ Standard_True,Standard_True,Standard_True,Standard_True,Standard_True,1,
+ aLenBeforeUfirst, aLenAfterUlast,
+ aLenBeforeVfirst, aLenAfterVlast);
+
+ DBRep::Set(theArgV[1], Result);
+ return 0;
+}
+
+//=======================================================================
+//function : QAStartsWith string startstring
+//=======================================================================
+static Standard_Integer QAStartsWith(Draw_Interpretor& di, Standard_Integer n, const char** a)
+{
+ if (n == 3)
+ {
+ TCollection_ExtendedString str = a[1];
+ TCollection_ExtendedString startstr = a[2];
+ if (str.StartsWith(startstr))
+ di << "Yes";
+ else
+ di << "No";
+ return 0;
+ }
+ std::cerr << "Syntax error\n";
+ return 1;
+}
+
+//=======================================================================
+//function : QAEndsWith string endstring
+//=======================================================================
+static Standard_Integer QAEndsWith(Draw_Interpretor& di, Standard_Integer n, const char** a)
+{
+ if (n == 3)
+ {
+ TCollection_ExtendedString str = a[1];
+ TCollection_ExtendedString endstr = a[2];
+ if (str.EndsWith(endstr))
+ di << "Yes";
+ else
+ di << "No";
+ return 0;
+ }
+ std::cerr << "Syntax error\n";
+ return 1;
+}
+
+//Class is used in OCC30435
+#include <AppCont_Function.hxx>
+#include <Adaptor3d_HCurve.hxx>
+class CurveEvaluator : public AppCont_Function
+
+{
+
+public:
+ Handle(Adaptor3d_HCurve) myCurve;
+
+ CurveEvaluator(const Handle(Adaptor3d_HCurve)& C)
+ : myCurve(C)
+ {
+ myNbPnt = 1;
+ myNbPnt2d = 0;
+ }
+
+ Standard_Real FirstParameter() const
+ {
+ return myCurve->FirstParameter();
+ }
+
+ Standard_Real LastParameter() const
+ {
+ return myCurve->LastParameter();
+ }
+
+ Standard_Boolean Value(const Standard_Real theT,
+ NCollection_Array1<gp_Pnt2d>& /*thePnt2d*/,
+ NCollection_Array1<gp_Pnt>& thePnt) const
+ {
+ thePnt(1) = myCurve->Value(theT);
+ return Standard_True;
+ }
+
+ Standard_Boolean D1(const Standard_Real theT,
+ NCollection_Array1<gp_Vec2d>& /*theVec2d*/,
+ NCollection_Array1<gp_Vec>& theVec) const
+ {
+ gp_Pnt aDummyPnt;
+ myCurve->D1(theT, aDummyPnt, theVec(1));
+ return Standard_True;
+ }
+};
+
+#include <GeomAdaptor_HCurve.hxx>
+#include <Approx_FitAndDivide.hxx>
+#include <Convert_CompBezierCurvesToBSplineCurve.hxx>
+static Standard_Integer OCC30435(Draw_Interpretor& di, Standard_Integer, const char** a)
+{
+
+ Handle(Geom_Curve) GC;
+ GC = DrawTrSurf::GetCurve(a[2]);
+ if (GC.IsNull())
+ return 1;
+
+ Standard_Integer Dmin = 3;
+ Standard_Integer Dmax = 12;
+ Standard_Real Tol3d = 1.e-7;
+ Standard_Boolean inverse = Standard_True;
+
+
+ Standard_Integer inv = atoi(a[3]);
+ if (inv > 0)
+ {
+ inverse = Standard_True;
+ }
+ else
+ {
+ inverse = Standard_False;
+ }
+
+ Standard_Integer maxit = atoi(a[4]);
+
+ Handle(GeomAdaptor_HCurve) aGAC = new GeomAdaptor_HCurve(GC);
+
+ CurveEvaluator aCE(aGAC);
+
+ Approx_FitAndDivide anAppro(Dmin, Dmax, Tol3d, 0., Standard_True);
+ anAppro.SetInvOrder(inverse);
+ Standard_Integer i;
+ for (i = 1; i <= maxit; ++i)
+ anAppro.Perform(aCE);
+
+ if (!anAppro.IsAllApproximated())
+ {
+ di << "Approximation failed \n";
+ return 1;
+ }
+ Standard_Integer NbCurves = anAppro.NbMultiCurves();
+
+ Convert_CompBezierCurvesToBSplineCurve Conv;
+
+ Standard_Real tol3d, tol2d, tolreached = 0.;
+ for (i = 1; i <= NbCurves; i++) {
+ anAppro.Error(i, tol3d, tol2d);
+ tolreached = Max(tolreached, tol3d);
+ AppParCurves_MultiCurve MC = anAppro.Value(i);
+ TColgp_Array1OfPnt Poles(1, MC.Degree() + 1);
+ MC.Curve(1, Poles);
+ Conv.AddCurve(Poles);
+ }
+ Conv.Perform();
+ Standard_Integer NbPoles = Conv.NbPoles();
+ Standard_Integer NbKnots = Conv.NbKnots();
+
+ TColgp_Array1OfPnt NewPoles(1, NbPoles);
+ TColStd_Array1OfReal NewKnots(1, NbKnots);
+ TColStd_Array1OfInteger NewMults(1, NbKnots);
+
+ Conv.KnotsAndMults(NewKnots, NewMults);
+ Conv.Poles(NewPoles);
+
+ BSplCLib::Reparametrize(GC->FirstParameter(),
+ GC->LastParameter(),
+ NewKnots);
+ Handle(Geom_BSplineCurve) TheCurve = new Geom_BSplineCurve(NewPoles, NewKnots, NewMults, Conv.Degree());
+
+ DrawTrSurf::Set(a[1], TheCurve);
+ di << a[1] << ": tolreached = " << tolreached << "\n";
+
+ return 0;
+
+}
+
+//=======================================================================
+//function : OCC30708_1
+//purpose : Tests initialization of the TopoDS_Iterator with null shape
+//=======================================================================
+static Standard_Integer OCC30708_1 (Draw_Interpretor& di, Standard_Integer, const char**)
+{
+ TopoDS_Iterator it;
+ try
+ {
+ OCC_CATCH_SIGNALS
+
+ TopoDS_Shape empty;
+ it.Initialize (empty);
+
+ }
+ catch (Standard_Failure)
+ {
+ di << "Cannot initialize TopoDS_Iterator with null shape\n";
+ return 0;
+ }
+
+ if (it.More())
+ di << "Incorrect Iterator initialization: method More() returns true on null shape\n";
+
+ return 0;
+}
+
+//=======================================================================
+//function : OCC30708_2
+//purpose : Tests initialization of the BRepLib_MakeWire with null wire
+//=======================================================================
+static Standard_Integer OCC30708_2 (Draw_Interpretor& di, Standard_Integer, const char**)
+{
+ try
+ {
+ OCC_CATCH_SIGNALS
+
+ TopoDS_Wire empty;
+ BRepLib_MakeWire aWBuilder (empty);
+ }
+ catch (Standard_Failure)
+ {
+ di << "Cannot initialize BRepLib_MakeWire with null wire\n";
+ }
+
+ return 0;
+}
+
+//=======================================================================
+//function : OCC30869
+//purpose :
+//=======================================================================
+static Standard_Integer OCC30869 (Draw_Interpretor& theDI, Standard_Integer theArgc, const char** theArgv)
+{
+ if (theArgc != 2)
+ {
+ theDI.PrintHelp (theArgv[0]);
+ return 1;
+ }
+
+ TopoDS_Shape aWire = DBRep::Get (theArgv[1]);
+ if (aWire.IsNull() || aWire.ShapeType() != TopAbs_WIRE)
+ {
+ theDI << theArgv[1] << " is not a wire.\n";
+ return 1;
+ }
+
+ BRepAdaptor_CompCurve aBACC (TopoDS::Wire (aWire));
+
+ Standard_Real aFirst = aBACC.FirstParameter();
+ Standard_Real aLast = aBACC.LastParameter();
+
+ gp_Pnt aPFirst, aPLast;
+ gp_Vec aVFirst, aVLast;
+
+ aBACC.D1 (aFirst, aPFirst, aVFirst);
+ aBACC.D1 (aLast, aPLast, aVLast);
+
+ if (aVFirst.SquareMagnitude() > gp::Resolution())
+ aVFirst.Normalize();
+ if (aVLast.SquareMagnitude() > gp::Resolution())
+ aVLast.Normalize();
+
+ theDI << aFirst << ": point " << aPFirst.X() << " "
+ << aPFirst.Y() << " "
+ << aPFirst.Z()
+ << ", tangent " << aVFirst.X() << " "
+ << aVFirst.Y() << " "
+ << aVFirst.Z() << "\n";
+
+ theDI << aLast << ": point " << aPLast.X() << " "
+ << aPLast.Y() << " "
+ << aPLast.Z()
+ << ", tangent " << aVLast.X() << " "
+ << aVLast.Y() << " "
+ << aVLast.Z() << "\n";
+
+ return 0;
+}
+
+#include <BRepPrimAPI_MakeBox.hxx>
+static Standard_Integer OCC30704(Draw_Interpretor& di, Standard_Integer, const char**)
+{
+ // Make a shape somewhere far from (0, 0, 0).
+ BRepPrimAPI_MakeBox mkBox(gp_Pnt(100, 100, 100), 100, 100, 100);
+ const TopoDS_Shape& box = mkBox.Shape();
+
+ // Add a bounding box of a shape to a void bounding box.
+ Bnd_OBB aVoidBox, aBox;
+ BRepBndLib::AddOBB(box, aBox, Standard_False, Standard_False, Standard_False);
+ aVoidBox.Add(aBox);
+
+ // Print the center point of the bounding box.
+ const gp_XYZ& center = aVoidBox.Center();
+ di << center.X() << " " << center.Y() << " " << center.Z();
+ return 0;
+}
+static Standard_Integer OCC30704_1(Draw_Interpretor& di, Standard_Integer, const char**)
+{
+ // A point.
+ gp_Pnt aP(100, 200, 300);
+
+ // Add the point to a void bounding box.
+ Bnd_OBB aVoidBox;
+ aVoidBox.Add(aP);
+
+ // Print the center point of the bounding box.
+ const gp_XYZ& center = aVoidBox.Center();
+ di << center.X() << " " << center.Y() << " " << center.Z();
+ return 0;
+}
+
void QABugs::Commands_20(Draw_Interpretor& theCommands) {
const char *group = "QABugs";
theCommands.Add("OCC27235", "OCC27235", __FILE__, OCC27235, group);
theCommands.Add("OCC26930", "OCC26930", __FILE__, OCC26930, group);
theCommands.Add("OCC27466", "OCC27466", __FILE__, OCC27466, group);
- theCommands.Add("OCC27341",
- "OCC27341 res shape axo/top/bottom/front/back/left/right",
- __FILE__, OCC27341, group);
theCommands.Add ("OCC26747_1", "OCC26747_1 result", __FILE__, OCC26747_1, group);
theCommands.Add ("OCC26747_2", "OCC26747_2 result", __FILE__, OCC26747_2, group);
theCommands.Add ("OCC26747_3", "OCC26747_3 result", __FILE__, OCC26747_3, group);
theCommands.Add("OCC26270", "OCC26270 shape result", __FILE__, OCC26270, group);
theCommands.Add ("OCC27552", "OCC27552", __FILE__, OCC27552, group);
theCommands.Add("OCC27875", "OCC27875 curve", __FILE__, OCC27875, group);
+ theCommands.Add("OCC28389", "OCC28389", __FILE__, OCC28389, group);
+ theCommands.Add("OCC28594", "OCC28594", __FILE__, OCC28594, group);
+ theCommands.Add("OCC28784", "OCC28784 result shape", __FILE__, OCC28784, group);
+ theCommands.Add("OCC28829", "OCC28829: perform invalid FPE operation", __FILE__, OCC28829, group);
+ theCommands.Add("OCC28887",
+ "OCC28887 filePath result"
+ "\n\t\t: Check interface for reading BRep from memory.",
+ __FILE__, OCC28887, group);
+ theCommands.Add("OCC28131", "OCC28131 name: creates face problematic for offset", __FILE__, OCC28131, group);
+ theCommands.Add("OCC29289", "OCC29289 : searching trigonometric root by Newton iterations", __FILE__, OCC29289, group);
+ theCommands.Add ("OCC29371", "OCC29371", __FILE__, OCC29371, group);
+ theCommands.Add("OCC29430", "OCC29430 <result wire> "
+ "<result first point> <result last point>",
+ __FILE__, OCC29430, group);
+ theCommands.Add("OCC29531", "OCC29531 <step file name>", __FILE__, OCC29531, group);
+
+ theCommands.Add ("OCC29064", "OCC29064: test memory usage by copying empty maps", __FILE__, OCC29064, group);
+ theCommands.Add ("OCC29925", "OCC29925: check safety of character classification functions", __FILE__, OCC29925, group);
+ theCommands.Add("OCC29807", "OCC29807 surface1 surface2 u1 v1 u2 v2", __FILE__, OCC29807, group);
+ theCommands.Add("OCC29311", "OCC29311 shape counter nbiter: check performance of OBB calculation", __FILE__, OCC29311, group);
+ theCommands.Add("OCC30391", "OCC30391 result face LenBeforeUfirst LenAfterUlast LenBeforeVfirst LenAfterVlast", __FILE__, OCC30391, group);
+ theCommands.Add("OCC30435", "OCC30435 result curve inverse nbit", __FILE__, OCC30435, group);
+
+ theCommands.Add("QAStartsWith",
+ "QAStartsWith string startstring",
+ __FILE__, QAStartsWith, group);
+
+ theCommands.Add("QAEndsWith",
+ "QAEndsWith string endstring",
+ __FILE__, QAEndsWith, group);
+
+ theCommands.Add ("OCC30708_1", "Tests initialization of the TopoDS_Iterator with null shape",
+ __FILE__, OCC30708_1, group);
+
+ theCommands.Add ("OCC30708_2", "Tests initialization of the BRepLib_MakeWire with null shape",
+ __FILE__, OCC30708_2, group);
+
+ theCommands.Add ("OCC30869", "Prints bounding points of the given wire and tangent vectors at these points.\n"
+ "Usage: OCC30869 wire",
+ __FILE__, OCC30869, group);
+
+ theCommands.Add("OCC30704", "OCC30704", __FILE__, OCC30704, group);
+ theCommands.Add("OCC30704_1", "OCC30704_1", __FILE__, OCC30704_1, group);
return;
}