// File: ShapeAnalysis_Wire.cxx // Created: Thu Jan 20 14:46:55 2000 // Author: data exchange team // //:pdn 11.12.98: FixDegenerated improved //:pdn 05.01.99: renaming method CheckLittle to CheckSmall //:l0 abv 10.01.99: CATIA01 #1727: fix intersecting edges always if edge is lacking //:n2 abv 22.01.99: ma-test5.igs: IGES read (pref3d): remove degen edge with no pcurve //:o4 abv 17.02.99: r0301_db.stp #53082: adding parameter isClosed to CheckOrder // rln 03.03.99 S4135: using updated ShapeAnalysis_Surface for checking of singularities //:p9 abv 11.03.99: PRO7226 #489490: fix :i9 moved to allow fixing a set of degenerated edges //#77 rln 11.03.99: S4135: using singularity which has minimum gap between singular point and input 3D point //#84 rln 18.03.99: inserting degenerated edge between ends of pcurves //pdn 12.03.99 S4135 check degenerated applies minimal tolerance first. //pdn 16.03.99 S4135 adding check of non adjacent edjes. //#83 rln 19.03.99: processing segments in intersection as in BRepCheck //%15 pdn 15.03.99 checking of small area wire added //#2 smh 26.03.99 S4163 Zero divide //#4 szv S4163 optimizing //:r6 abv 08.04.99: protect FixIE against working out of curve range //:s1 abv 22.04.99: PRO7226 #489490: ensure fixing of degenerated edge //#9 smh 14.12.99 BUC60615 Using tolerance of verteces during checking degenerated edge. #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include //szvsh addition #include #include #include #include #include #include #include #include #include #include #include #include #include #include //======================================================================= //function : ShapeAnalysis_Wire //purpose : //======================================================================= ShapeAnalysis_Wire::ShapeAnalysis_Wire() { ClearStatuses(); myPrecision = ::Precision::Confusion(); } //======================================================================= //function : ShapeAnalysis_Wire //purpose : //======================================================================= ShapeAnalysis_Wire::ShapeAnalysis_Wire (const TopoDS_Wire& wire, const TopoDS_Face& face, const Standard_Real precision) { Init (wire, face, precision); } //======================================================================= //function : ShapeAnalysis_Wire //purpose : //======================================================================= ShapeAnalysis_Wire::ShapeAnalysis_Wire (const Handle(ShapeExtend_WireData)& sbwd, const TopoDS_Face& face, const Standard_Real precision) { Init (sbwd, face, precision); } //======================================================================= //function : Init //purpose : //======================================================================= void ShapeAnalysis_Wire::Init (const TopoDS_Wire& wire, const TopoDS_Face& face, const Standard_Real precision) { Init (new ShapeExtend_WireData (wire), face, precision); } //======================================================================= //function : Init //purpose : //======================================================================= void ShapeAnalysis_Wire::Init (const Handle(ShapeExtend_WireData)& sbwd, const TopoDS_Face& face, const Standard_Real precision) { Load (sbwd); SetFace (face); SetPrecision (precision); } //======================================================================= //function : Load //purpose : //======================================================================= void ShapeAnalysis_Wire::Load (const TopoDS_Wire& wire) { ClearStatuses(); myWire = new ShapeExtend_WireData (wire); } //======================================================================= //function : Load //purpose : //======================================================================= void ShapeAnalysis_Wire::Load (const Handle(ShapeExtend_WireData)& sbwd) { ClearStatuses(); myWire = sbwd; } //======================================================================= //function : SetFace //purpose : //======================================================================= void ShapeAnalysis_Wire::SetFace(const TopoDS_Face& face) { myFace = face; if(!face.IsNull()) mySurf = new ShapeAnalysis_Surface ( BRep_Tool::Surface ( myFace ) ); } //======================================================================= //function : SetSurface //purpose : //======================================================================= void ShapeAnalysis_Wire::SetSurface (const Handle(Geom_Surface)& surface) { SetSurface ( surface, TopLoc_Location() ); } //======================================================================= //function : SetSurface //purpose : //======================================================================= void ShapeAnalysis_Wire::SetSurface (const Handle(Geom_Surface)& surface, const TopLoc_Location& location) { BRep_Builder B; TopoDS_Face face; B.MakeFace ( face, surface, location, ::Precision::Confusion() ); SetFace ( face ); } //======================================================================= //function : SetPrecision //purpose : //======================================================================= void ShapeAnalysis_Wire::SetPrecision(const Standard_Real precision) { myPrecision = precision; } //======================================================================= //function : ClearStatuses //purpose : //======================================================================= void ShapeAnalysis_Wire::ClearStatuses() { myStatusOrder = myStatusConnected = myStatusEdgeCurves = myStatusDegenerated = myStatusClosed = myStatusLacking = myStatusSelfIntersection = myStatusSmall = myStatusGaps3d = myStatusGaps2d = myStatusCurveGaps = myStatusLoop = myStatus = 0; myMin3d = myMin2d = myMax3d = myMax2d = 0.; } //======================================================================= //function : Perform //purpose : //======================================================================= Standard_Boolean ShapeAnalysis_Wire::Perform() { Standard_Boolean result = Standard_False; result |= CheckOrder(); result |= CheckSmall(); result |= CheckConnected(); result |= CheckEdgeCurves(); result |= CheckDegenerated(); result |= CheckSelfIntersection(); result |= CheckLacking(); result |= CheckClosed(); return result; } //======================================================================= //function : CheckOrder //purpose : //======================================================================= Standard_Boolean ShapeAnalysis_Wire::CheckOrder (const Standard_Boolean isClosed, const Standard_Boolean mode3d) { ShapeAnalysis_WireOrder sawo; CheckOrder (sawo, isClosed, mode3d); myStatusOrder = myStatus; return StatusOrder (ShapeExtend_DONE); } //======================================================================= //function : CheckSmall //purpose : //======================================================================= Standard_Boolean ShapeAnalysis_Wire::CheckSmall(const Standard_Real precsmall) { for (Standard_Integer i = 1; i <= myWire->NbEdges(); i++) { CheckSmall (i, precsmall); myStatusSmall |= myStatus; } return StatusSmall (ShapeExtend_DONE); } //======================================================================= //function : CheckConnected //purpose : //======================================================================= Standard_Boolean ShapeAnalysis_Wire::CheckConnected(const Standard_Real prec) { for (Standard_Integer i = 1; i <= myWire->NbEdges(); i++) { CheckConnected ( i, prec ); myStatusConnected |= myStatus; } return StatusConnected (ShapeExtend_DONE); } //======================================================================= //function : CheckEdgeCurves //purpose : //======================================================================= Standard_Boolean ShapeAnalysis_Wire::CheckEdgeCurves() { myStatusEdgeCurves = ShapeExtend::EncodeStatus ( ShapeExtend_OK ); if ( ! IsReady() ) return Standard_False; Standard_Integer i, nb = myWire->NbEdges(); ShapeAnalysis_Edge SAE; for (i = 1; i <= nb; i++) { TopoDS_Edge E = myWire->Edge (i); SAE.CheckCurve3dWithPCurve (E, myFace); if (SAE.Status (ShapeExtend_DONE)) myStatusEdgeCurves |= ShapeExtend::EncodeStatus (ShapeExtend_DONE1); if (SAE.Status ( ShapeExtend_FAIL)) myStatusEdgeCurves |= ShapeExtend::EncodeStatus (ShapeExtend_FAIL1); SAE.CheckVerticesWithPCurve (E, myFace); if (SAE.Status (ShapeExtend_DONE)) myStatusEdgeCurves |= ShapeExtend::EncodeStatus (ShapeExtend_DONE2); if (SAE.Status ( ShapeExtend_FAIL)) myStatusEdgeCurves |= ShapeExtend::EncodeStatus (ShapeExtend_FAIL2); SAE.CheckVerticesWithCurve3d (E); if (SAE.Status (ShapeExtend_DONE)) myStatusEdgeCurves |= ShapeExtend::EncodeStatus (ShapeExtend_DONE3); if (SAE.Status ( ShapeExtend_FAIL)) myStatusEdgeCurves |= ShapeExtend::EncodeStatus (ShapeExtend_FAIL3); CheckSeam (i); if (LastCheckStatus (ShapeExtend_DONE)) myStatusEdgeCurves |= ShapeExtend::EncodeStatus (ShapeExtend_DONE4); if (LastCheckStatus (ShapeExtend_FAIL)) myStatusEdgeCurves |= ShapeExtend::EncodeStatus (ShapeExtend_FAIL4); CheckGap3d (i); if (LastCheckStatus (ShapeExtend_DONE)) myStatusEdgeCurves |= ShapeExtend::EncodeStatus (ShapeExtend_DONE5); if (LastCheckStatus (ShapeExtend_FAIL)) myStatusEdgeCurves |= ShapeExtend::EncodeStatus (ShapeExtend_FAIL5); CheckGap2d (i); if (LastCheckStatus (ShapeExtend_DONE)) myStatusEdgeCurves |= ShapeExtend::EncodeStatus (ShapeExtend_DONE6); if (LastCheckStatus (ShapeExtend_FAIL)) myStatusEdgeCurves |= ShapeExtend::EncodeStatus (ShapeExtend_FAIL6); Standard_Real maxdev = 0.0; SAE.CheckSameParameter (myWire->Edge (i), maxdev); if (SAE.Status (ShapeExtend_DONE)) myStatusEdgeCurves |= ShapeExtend::EncodeStatus (ShapeExtend_DONE7); if (SAE.Status ( ShapeExtend_FAIL)) myStatusEdgeCurves |= ShapeExtend::EncodeStatus (ShapeExtend_FAIL7); } return StatusEdgeCurves (ShapeExtend_DONE); } //======================================================================= //function : CheckDegenerated //purpose : //======================================================================= Standard_Boolean ShapeAnalysis_Wire::CheckDegenerated() { for (Standard_Integer i = 1; i <= myWire->NbEdges(); i++) { CheckDegenerated (i); myStatusDegenerated |= myStatus; } return StatusDegenerated (ShapeExtend_DONE); } //======================================================================= //function : CheckSelfIntersection //purpose : //======================================================================= Standard_Boolean ShapeAnalysis_Wire::CheckSelfIntersection() { myStatusSelfIntersection = ShapeExtend::EncodeStatus ( ShapeExtend_OK ); if (!IsReady()) return Standard_False; Standard_Integer i, nb = myWire->NbEdges(); for (i = 1; i <= nb; i++) { CheckSelfIntersectingEdge (i); if (LastCheckStatus (ShapeExtend_DONE)) myStatusSelfIntersection |= ShapeExtend::EncodeStatus (ShapeExtend_DONE1); if (LastCheckStatus (ShapeExtend_FAIL)) myStatusSelfIntersection |= ShapeExtend::EncodeStatus (ShapeExtend_FAIL1); CheckIntersectingEdges (i); if (LastCheckStatus (ShapeExtend_DONE)) myStatusSelfIntersection |= ShapeExtend::EncodeStatus (ShapeExtend_DONE2); if (LastCheckStatus (ShapeExtend_FAIL)) myStatusSelfIntersection |= ShapeExtend::EncodeStatus (ShapeExtend_FAIL2); } Bnd_Array1OfBox2d boxes(1,nb); TopLoc_Location L; const Handle(Geom_Surface)& S = BRep_Tool::Surface(Face(), L); Handle(Geom2d_Curve) c2d; Standard_Real cf,cl; ShapeAnalysis_Edge sae; Handle(ShapeExtend_WireData) sbwd = WireData(); for(i = 1; i <= nb; i++){ TopoDS_Edge E = sbwd->Edge (i); if(sae.PCurve (E,S,L,c2d,cf,cl,Standard_False)) { Bnd_Box2d box; Geom2dAdaptor_Curve gac(c2d,cf,cl); BndLib_Add2dCurve::Add(gac,::Precision::Confusion(),box); boxes(i) = box; } } Standard_Boolean isFail = Standard_False, isDone = Standard_False; for(Standard_Integer num1 = 1; num1 < nb-1; num1++) { Standard_Integer fin = (num1 == 1 ? nb-1 : nb); for(Standard_Integer num2 = num1+2; num2 <= fin; num2++) if(!boxes(num1).IsOut(boxes(num2))){ CheckIntersectingEdges(num1, num2); isFail |= LastCheckStatus ( ShapeExtend_FAIL1 ); isDone |= LastCheckStatus ( ShapeExtend_DONE1 ); } } if(isFail) myStatusSelfIntersection |= ShapeExtend::EncodeStatus ( ShapeExtend_FAIL3 ); if(isDone) myStatusSelfIntersection |= ShapeExtend::EncodeStatus ( ShapeExtend_DONE3 ); return StatusSelfIntersection (ShapeExtend_DONE); } //======================================================================= //function : CheckLacking //purpose : //======================================================================= Standard_Boolean ShapeAnalysis_Wire::CheckLacking() { if (!IsReady() || NbEdges() < 2) return Standard_False; for (Standard_Integer i = 1; i <= myWire->NbEdges(); i++) { CheckLacking (i); myStatusLacking |= myStatus; } return StatusLacking (ShapeExtend_DONE); } //======================================================================= //function : CheckClosed //purpose : //======================================================================= Standard_Boolean ShapeAnalysis_Wire::CheckClosed(const Standard_Real prec) { myStatusClosed = ShapeExtend::EncodeStatus ( ShapeExtend_OK ); if (!IsReady() || NbEdges() < 1) return Standard_False; CheckConnected (1, prec); if ( LastCheckStatus ( ShapeExtend_DONE ) ) myStatusClosed |= ShapeExtend::EncodeStatus ( ShapeExtend_DONE1 ); if ( LastCheckStatus ( ShapeExtend_FAIL ) ) myStatusClosed |= ShapeExtend::EncodeStatus ( ShapeExtend_FAIL1 ); CheckDegenerated ( 1 ); if ( LastCheckStatus ( ShapeExtend_DONE ) ) myStatusClosed |= ShapeExtend::EncodeStatus ( ShapeExtend_DONE2 ); if ( LastCheckStatus ( ShapeExtend_FAIL ) ) myStatusClosed |= ShapeExtend::EncodeStatus ( ShapeExtend_FAIL2 ); return StatusClosed ( ShapeExtend_DONE ); } //======================================================================= //function : CheckGaps3d //purpose : //======================================================================= Standard_Boolean ShapeAnalysis_Wire::CheckGaps3d () { myStatusGaps3d = ShapeExtend::EncodeStatus ( ShapeExtend_OK ); if (!IsLoaded() || NbEdges() < 1) return Standard_False; //gka IsLoaded Standard_Real dist, maxdist = 0.; for (Standard_Integer i = 1; i <= NbEdges(); i++) { CheckGap3d(i); myStatusGaps3d |= myStatus; if (!LastCheckStatus(ShapeExtend_FAIL1)) { dist = MinDistance3d(); if (maxdistNbEdges(); ShapeAnalysis_Edge EA; for (i = 1; i <= nb; i ++) { TopoDS_Edge E = myWire->Edge(i); if ( mode3d ) { TopoDS_Vertex V1 = EA.FirstVertex (E); TopoDS_Vertex V2 = EA.LastVertex (E); gp_Pnt p1 = BRep_Tool::Pnt (V1); gp_Pnt p2 = BRep_Tool::Pnt (V2); sawo.Add (p1.XYZ(),p2.XYZ()); } else { Standard_Real f,l; Handle(Geom2d_Curve) c2d; TopoDS_Shape tmpF = myFace.Oriented(TopAbs_FORWARD); if ( ! EA.PCurve(E,TopoDS::Face(tmpF),c2d,f,l) ) { myStatus = ShapeExtend::EncodeStatus (ShapeExtend_FAIL2); return Standard_False; } sawo.Add(c2d->Value(f).XY(),c2d->Value(l).XY()); } } sawo.Perform(isClosed); Standard_Integer stat = sawo.Status(); switch (stat) { case 0: myStatus = ShapeExtend::EncodeStatus (ShapeExtend_OK); break; case 1: myStatus = ShapeExtend::EncodeStatus (ShapeExtend_DONE1); break; case 2: myStatus = ShapeExtend::EncodeStatus (ShapeExtend_DONE2); break; case -1: myStatus = ShapeExtend::EncodeStatus (ShapeExtend_DONE3); break; case -2: myStatus = ShapeExtend::EncodeStatus (ShapeExtend_DONE4); break; case -10: myStatus = ShapeExtend::EncodeStatus (ShapeExtend_FAIL1); break; } return LastCheckStatus (ShapeExtend_DONE); } //======================================================================= //function : CheckConnected //purpose : //======================================================================= Standard_Boolean ShapeAnalysis_Wire::CheckConnected (const Standard_Integer num, const Standard_Real prec) { myStatus = ShapeExtend::EncodeStatus (ShapeExtend_OK); if ( ! IsLoaded() || NbEdges() < 1 ) return Standard_False; Standard_Integer n2 = ( num >0 ? num : NbEdges() ); Standard_Integer n1 = ( n2 >1 ? n2-1 : NbEdges() ); // if (n1 == n2) return 0; TopoDS_Edge E1 = WireData()->Edge ( n1 ); TopoDS_Edge E2 = WireData()->Edge ( n2 ); ShapeAnalysis_Edge sae; TopoDS_Vertex V1 = sae.LastVertex (E1); TopoDS_Vertex V2 = sae.FirstVertex (E2); if (V1.IsSame(V2)) return Standard_False; gp_Pnt p1 = BRep_Tool::Pnt (V1); gp_Pnt p2 = BRep_Tool::Pnt (V2); myMin3d = p1.Distance(p2); if ( myMin3d <= gp::Resolution() ) myStatus |= ShapeExtend::EncodeStatus (ShapeExtend_DONE1); else if ( myMin3d <= myPrecision ) myStatus |= ShapeExtend::EncodeStatus (ShapeExtend_DONE2); else if ( myMin3d <= prec ) myStatus |= ShapeExtend::EncodeStatus (ShapeExtend_DONE3); else { // et en inversant la derniere edge ? if ( n1 == n2 ) myStatus = ShapeExtend::EncodeStatus (ShapeExtend_FAIL1); else { V2 = sae.LastVertex (E2); p2 = BRep_Tool::Pnt (V2); Standard_Real dist = p1.Distance(p2); if ( dist > myPrecision ) myStatus = ShapeExtend::EncodeStatus (ShapeExtend_FAIL1); else { myMin3d = dist; myStatus = ShapeExtend::EncodeStatus (ShapeExtend_FAIL2); } } return Standard_False; } return Standard_True; } //======================================================================= //function : CheckSmall //purpose : //======================================================================= Standard_Boolean ShapeAnalysis_Wire::CheckSmall (const Standard_Integer num, const Standard_Real precsmall) { myStatus = ShapeExtend::EncodeStatus (ShapeExtend_OK); if ( ! IsLoaded() || NbEdges() <= 1 ) return Standard_False; //Standard_Integer n = ( num ? num : NbEdges() ); //szv#4:S4163:12Mar99 not needed TopoDS_Edge E = myWire->Edge ( num ? num : NbEdges() ); ShapeAnalysis_Edge sae; if ( BRep_Tool::Degenerated ( E ) ) { //:n2 abv 22 Jan 99: ma-test5.igs -> IGES (brep) -> read (pref3d): // degen edge with no pcurve should be removed if ( ! myFace.IsNull() && sae.HasPCurve ( E, Face() ) ) return Standard_False; myStatus = ShapeExtend::EncodeStatus (ShapeExtend_FAIL1); } TopoDS_Vertex V1 = sae.FirstVertex (E); TopoDS_Vertex V2 = sae.LastVertex (E); gp_Pnt p1 = BRep_Tool::Pnt (V1); gp_Pnt p2 = BRep_Tool::Pnt (V2); Standard_Real dist = p1.Distance(p2); Standard_Real prec = precsmall;//Min ( myPrecision, precsmall ); //Standard_Real prec = Min(BRep_Tool::Tolerance(V1),BRep_Tool::Tolerance(V2)); //skl if (dist > prec) return Standard_False; // pas nulle // La courbe 3D a present : est-elle FERMEE ou DE LONGUEUR NULLE ... ??? // Pour cela on prend le point milieu (y a-t-il mieux) // Si pas de C3D, on essaie la C2D ... gp_Pnt Pm; Standard_Real cf,cl; Handle(Geom_Curve) c3d; if ( sae.Curve3d (E,c3d,cf,cl,Standard_False) ) Pm = c3d->Value ( (cf+cl)/2. ); else { Handle(Geom2d_Curve) c2d; if ( ! myFace.IsNull() && sae.PCurve (E,myFace,c2d,cf,cl,Standard_False)) { gp_Pnt2d p2m = c2d->Value ( (cf+cl)/2. ); Pm = mySurf->Value (p2m); } else { myStatus = ShapeExtend::EncodeStatus (ShapeExtend_FAIL1); Pm = p1; //:n2 return Standard_False; } } if ( Pm.Distance(p1) > prec || Pm.Distance(p2) > prec ) return Standard_False; myStatus |= ShapeExtend::EncodeStatus ( V1.IsSame(V2) ? ShapeExtend_DONE1 : ShapeExtend_DONE2 ); return Standard_True; } //======================================================================= //function : CheckSeam //purpose : //======================================================================= Standard_Boolean ShapeAnalysis_Wire::CheckSeam(const Standard_Integer num, Handle(Geom2d_Curve)& C1, Handle(Geom2d_Curve)& C2, Standard_Real& cf, Standard_Real& cl) { myStatus = ShapeExtend::EncodeStatus (ShapeExtend_OK); if (!IsReady()) return Standard_False; Standard_Integer n = num; if (n == 0) n = NbEdges(); TopoDS_Edge E = myWire->Edge (n); if ( ! ShapeAnalysis_Edge().IsSeam ( E, myFace ) ) return Standard_False; // Extract the Two PCurves of the Seam TopoDS_Face ForwardFace = myFace; ForwardFace.Orientation (TopAbs_FORWARD); //szv#4:S4163:12Mar99 SGI warns TopoDS_Shape EF = E.Oriented(TopAbs_FORWARD); TopoDS_Shape ER = E.Oriented(TopAbs_REVERSED); C1 = BRep_Tool::CurveOnSurface(TopoDS::Edge(EF), ForwardFace, cf, cl); C2 = BRep_Tool::CurveOnSurface(TopoDS::Edge(ER), ForwardFace, cf, cl); if (C1.IsNull() || C2.IsNull()) return Standard_False; // SelectForward est destine a devenir un outil distinct Standard_Integer theCurveIndice = ShapeAnalysis_Curve().SelectForwardSeam (C1,C2); if ( theCurveIndice != 2 ) return Standard_False; myStatus = ShapeExtend::EncodeStatus (ShapeExtend_DONE1); return Standard_True; } //======================================================================= //function : CheckSeam //purpose : //======================================================================= Standard_Boolean ShapeAnalysis_Wire::CheckSeam(const Standard_Integer num) { Handle(Geom2d_Curve) C1, C2; Standard_Real cf, cl; return CheckSeam (num, C1, C2, cf, cl); } //======================================================================= //function : CheckDegenerated //purpose : //======================================================================= Standard_Boolean ShapeAnalysis_Wire::CheckDegenerated (const Standard_Integer num, gp_Pnt2d& p2d1, gp_Pnt2d& p2d2) { myStatus = ShapeExtend::EncodeStatus (ShapeExtend_OK); if ( ! IsReady() || NbEdges() < 1 ) return Standard_False; Standard_Integer n2 = (num > 0)? num : NbEdges(); Standard_Integer n1 = (n2 > 1)? n2-1 : NbEdges(); Standard_Integer n3 = (n2 < NbEdges())? n2+1 : 1; TopoDS_Edge E1 = myWire->Edge ( n1 ); TopoDS_Edge E2 = myWire->Edge ( n2 ); TopoDS_Edge E3 = myWire->Edge ( n3 ); ShapeAnalysis_Edge sae; // skip if edge is already marked as degenerated and has pcurve if ( BRep_Tool::Degenerated ( E2 ) && sae.HasPCurve ( E2, Face() ) ) { // skl 30.12.2004 for OCC7630 - we have to check pcurve if( sae.HasPCurve(E1,Face()) && sae.HasPCurve(E3,Face()) ) { Handle(Geom2d_Curve) c2d; Standard_Real fp,lp; sae.PCurve ( E2, myFace, c2d, fp, lp, Standard_True ); gp_Pnt2d p21 = c2d->Value(fp); gp_Pnt2d p22 = c2d->Value(lp); sae.PCurve ( E1, myFace, c2d, fp, lp, Standard_True ); gp_Pnt2d p12 = c2d->Value(lp); sae.PCurve ( E3, myFace, c2d, fp, lp, Standard_True ); gp_Pnt2d p31 = c2d->Value(fp); if( fabs(p12.Distance(p31)-p21.Distance(p22)) > 2*Precision::PConfusion() ) { // pcurve is bad => we can remove this edge in ShapeFix // if set needed status myStatus = ShapeExtend::EncodeStatus (ShapeExtend_FAIL2); } } return Standard_False; } //pdn allows to insert two sequences of degenerated edges (on separate bounds of surfaces) if ( n1 != n2 && BRep_Tool::Degenerated ( E1 ) && ! sae.HasPCurve ( E1, Face() ) ) { //:abv 13.05.02: OCC320 - fail (to remove edge) if two consequtive degenerated edges w/o pcurves if ( BRep_Tool::Degenerated ( E2 ) ) myStatus |= ShapeExtend::EncodeStatus (ShapeExtend_FAIL2); return Standard_False; } //:i8 if ( BRep_Tool::Degenerated ( E1 ) || //:i8 BRep_Tool::Degenerated ( E2 ) ) return Standard_False; // deja OK TopoDS_Vertex Vp = sae.FirstVertex (E1); //:i9 TopoDS_Vertex V0 = sae.LastVertex (E1); TopoDS_Vertex V1 = sae.FirstVertex (E2); TopoDS_Vertex V2 = sae.LastVertex (E2); gp_Pnt pp = BRep_Tool::Pnt (Vp); //:i9 gp_Pnt p0 = BRep_Tool::Pnt (V0); gp_Pnt p1 = BRep_Tool::Pnt (V1); gp_Pnt p2 = BRep_Tool::Pnt (V2); Standard_Real par1, par2; Standard_Boolean lack = Standard_False; Standard_Boolean dgnr = Standard_False; //pdn 12.03.99 minimal value processing first Standard_Real precFirst = Min(myPrecision,BRep_Tool::Tolerance(V1)); Standard_Real precFin = Max(myPrecision,BRep_Tool::Tolerance(V1)); Standard_Real precVtx = (myPrecision FWD/REV. Si LACK, toujours True Standard_Boolean forward = ( E2.Orientation() == TopAbs_FORWARD ); // FIX FEV 1998 : recompute singularity according precision if (p1.Distance(p2) <= precFirst) { // edge DGNR dgnr = mySurf->DegeneratedValues ( p1, precVtx, p2d1, p2d2, par1, par2, forward ); //smh#9 if ( dgnr ) { // abv 24 Feb 00: trj3_as1-ac-214.stp #6065: avoid making closed edge degenerated Standard_Real a, b; Handle(Geom_Curve) C3d = BRep_Tool::Curve ( E2, a, b ); if ( ! C3d.IsNull() ) { gp_Pnt p = C3d->Value ( 0.5 * ( a + b ) ); if ( p.SquareDistance ( p1 ) > precVtx * precVtx ) dgnr = Standard_False; } } } if ( ! dgnr ) { //:i9 abv 23 Sep 98: CTS20315-2 #63231: check that previous edge is not degenerated if ( n1 != n2 && p1.Distance(pp) <= precFirst && mySurf->IsDegenerated ( pp, precFirst ) && ! BRep_Tool::Degenerated ( E1 ) ) return Standard_False; //rln S4135 ShapeAnalysis_Surface new algorithms for singularities //:45 by abv 16 Dec 97: BUC60035 2659: precision increased to vertex tolerance //Standard_Real prec = Max ( myPrecision, BRep_Tool::Tolerance(V1) ); //:51 abv 22 Dec 97: recompute singularities if necessary //rln S4135 if ( prec > myPrecision ) mySurf->ComputeSingularities ( 2 * prec ); //:51 //:74 abv 15 Jan 97: *2 if ( p0.Distance ( p1 ) <= precFin ) {// ou DGNR manquante ? //rln S4135 singularity with precision = 2 * prec, but distance <= prec //lack = mySurf->DegeneratedValues ( p1, prec, p2d1, p2d2, par1, par2, forward); Standard_Real tmpPreci; gp_Pnt tmpP3d; Standard_Boolean tmpUIsoDeg; //#77 rln S4135: using singularity which has minimum gap between singular point and input 3D point Standard_Integer indMin = -1; Standard_Real gapMin2 = RealLast(); for (Standard_Integer i = 1; i <= mySurf->NbSingularities (precVtx); i++) { mySurf->Singularity (i, tmpPreci, tmpP3d, p2d1, p2d2, par1, par2, tmpUIsoDeg); Standard_Real gap2 = p1.SquareDistance (tmpP3d); if (gap2 <= precVtx * precVtx) if (gapMin2 > gap2) { gapMin2 = gap2; indMin = i; } } if (indMin >= 1) { mySurf->Singularity (indMin, tmpPreci, tmpP3d, p2d1, p2d2, par1, par2, tmpUIsoDeg); lack = Standard_True; } } //rln S4135 if ( prec > myPrecision ) mySurf->ComputeSingularities ( myPrecision ); //:51 } // voila, on a soit dgnr soit lack if ( ! lack && ! dgnr ) { //:abv 29.08.01: if singularity not detected but edge is marked // as degenerated, report fail if ( BRep_Tool::Degenerated ( E2 ) && ! sae.HasPCurve ( E2, Face() ) ) myStatus |= ShapeExtend::EncodeStatus (ShapeExtend_FAIL2); return Standard_False; } // OK, degenerated case detected; we will find its start and end in 2d if ( lack ) forward = Standard_True; //:24 by abv 28 Nov 97: // make degenerative pcurve parametrized exactly from end of pcurve of the // previous edge to the start of the next one if ( lack || n1 != n2 ) { //:i8 abv 18 Sep 98: ProSTEP TR9 r0501-ug.stp #182180: single degedge is a wire at apex of a cone Standard_Real a, b; Handle(Geom2d_Curve) c2d; if ( sae.PCurve ( E1, myFace, c2d, a, b, Standard_True ) ) { p2d1 = c2d->Value ( b ); //#84 rln gp_Pnt2d p2d = c2d->Value ( b ); //#84 rln par1 = ( p2d.XY() - aP2d.XY() ) * theDir2d.XY(); } else myStatus |= ShapeExtend::EncodeStatus (ShapeExtend_FAIL1); //pdn pcurves (fixing regression in f0 in degenerated case) if ( sae.PCurve ( ( dgnr ? E3 : E2 ), myFace, c2d, a, b, Standard_True ) ) { p2d2 = c2d->Value ( a ); //#84 rln gp_Pnt2d p2d = c2d->Value ( a ); //#84 rln par2 = ( p2d.XY() - aP2d.XY() ) * theDir2d.XY(); } else myStatus |= ShapeExtend::EncodeStatus (ShapeExtend_FAIL1); } /* if ( par2 < par1 ) { par1 = -par1; par2 = -par2; theDir2d.Reverse(); } */ //#84 rln 18.03.99 if pcurve is not degenerate anymore, the fix is postponned //to ShapeFix_Wire::FixLacking if ( ! mySurf->IsDegenerated ( p2d1, p2d2, precVtx, 10. ) ) { //:s1 abv 22 Apr 99: PRO7226 #489490 //smh#9 //:abv 24.05.02: OCC320 - fail (to remove edge) if two consequtive degenerated edges w/o pcurves if ( BRep_Tool::Degenerated ( E2 ) ) myStatus |= ShapeExtend::EncodeStatus (ShapeExtend_FAIL2); return Standard_False; } //added by rln 18/12/97 CSR# CTS18544 entity 13638 //the situation when degenerated edge already exists but flag is not set //(i.e. the parametric space is closed) GeomAdaptor_Surface& Ads = mySurf->Adaptor3d()->ChangeSurface(); Standard_Real max = Max ( Ads.UResolution(myPrecision), Ads.VResolution(myPrecision) ); if ( p2d1.Distance (p2d2) /*Abs (par1 - par2)*/ <= max + gp::Resolution() ) return Standard_False; //#84 rln p2d1 = aP2d.XY() + par1 * theDir2d.XY(); //#84 rln p2d2 = aP2d.XY() + par2 * theDir2d.XY(); myStatus = ShapeExtend::EncodeStatus ( dgnr ? ShapeExtend_DONE2 : ShapeExtend_DONE1 ); return Standard_True; } //======================================================================= //function : CheckDegenerated //purpose : //======================================================================= Standard_Boolean ShapeAnalysis_Wire::CheckDegenerated (const Standard_Integer num) { gp_Pnt2d p2d1, p2d2; return CheckDegenerated ( num, p2d1, p2d2 ); } //======================================================================= //function : CheckGap3d //purpose : //======================================================================= Standard_Boolean ShapeAnalysis_Wire::CheckGap3d(const Standard_Integer num) { myStatus = ShapeExtend::EncodeStatus (ShapeExtend_OK); //szv#4:S4163:12Mar99 optimized if ( !IsLoaded() || NbEdges() < 1 ) return Standard_False; //szvsh was nbedges < 2 Standard_Integer n2 = ( num >0 ? num : NbEdges() ); Standard_Integer n1 = ( n2 >1 ? n2-1 : NbEdges() ); TopoDS_Edge E1 = myWire->Edge(n1); TopoDS_Edge E2 = myWire->Edge(n2); Standard_Real uf1,ul1,uf2,ul2; Handle(Geom_Curve) C1,C2; ShapeAnalysis_Edge SAE; if (!SAE.Curve3d (E1,C1,uf1,ul1) || !SAE.Curve3d (E2,C2,uf2,ul2)) { myStatus = ShapeExtend::EncodeStatus (ShapeExtend_FAIL1); return Standard_False; } gp_Pnt p1 = C1->Value (ul1); gp_Pnt p2 = C2->Value (uf2); myMin3d = myMax3d = p1.Distance (p2); if (myMin3d > myPrecision) myStatus = ShapeExtend::EncodeStatus (ShapeExtend_DONE1); return LastCheckStatus (ShapeExtend_DONE); } //======================================================================= //function : CheckGap2d //purpose : //======================================================================= Standard_Boolean ShapeAnalysis_Wire::CheckGap2d(const Standard_Integer num) { myStatus = ShapeExtend::EncodeStatus (ShapeExtend_OK); //szv#4:S4163:12Mar99 optimized if ( !IsReady() || NbEdges() < 1 ) return Standard_False; //szvsh was nbedges < 2 Standard_Integer n2 = ( num >0 ? num : NbEdges() ); Standard_Integer n1 = ( n2 >1 ? n2-1 : NbEdges() ); TopoDS_Edge E1 = myWire->Edge(n1); TopoDS_Edge E2 = myWire->Edge(n2); Standard_Real uf1,ul1,uf2,ul2; Handle(Geom2d_Curve) C1,C2; ShapeAnalysis_Edge SAE; if (!SAE.PCurve (E1,myFace,C1,uf1,ul1) || !SAE.PCurve (E2,myFace,C2,uf2,ul2)) { myStatus = ShapeExtend::EncodeStatus (ShapeExtend_FAIL1); return Standard_False; } gp_Pnt2d p1 = C1->Value (ul1); gp_Pnt2d p2 = C2->Value (uf2); myMin2d = myMax2d = p1.Distance (p2); GeomAdaptor_Surface& SA = mySurf->Adaptor3d()->ChangeSurface(); if (myMin2d > (Max (SA.UResolution (myPrecision), SA.VResolution (myPrecision)) + Precision::PConfusion())) myStatus = ShapeExtend::EncodeStatus (ShapeExtend_DONE1); return LastCheckStatus (ShapeExtend_DONE); } //======================================================================= //function : CheckCurveGap //purpose : //======================================================================= Standard_Boolean ShapeAnalysis_Wire::CheckCurveGap(const Standard_Integer num) { myStatus = ShapeExtend::EncodeStatus (ShapeExtend_OK); if ( !IsLoaded() || NbEdges() < 1 ) return Standard_False; Standard_Integer n = ( num >0 ? num : NbEdges() ); TopoDS_Edge E = myWire->Edge(n); Standard_Real cuf,cul,pcuf,pcul; Handle(Geom_Curve) c; ShapeAnalysis_Edge SAE; if (!SAE.Curve3d (E,c,cuf,cul,Standard_False)) { myStatus = ShapeExtend::EncodeStatus (ShapeExtend_FAIL1); return Standard_False; } Handle(Geom2d_Curve) pc; if (!SAE.PCurve (E,myFace,pc,pcuf,pcul,Standard_False)) { myStatus = ShapeExtend::EncodeStatus (ShapeExtend_FAIL1); return Standard_False; } Handle(Geom2dAdaptor_HCurve) AC = new Geom2dAdaptor_HCurve(pc,pcuf,pcul); Handle(GeomAdaptor_HSurface) AS = new GeomAdaptor_HSurface(mySurf->Surface()); Adaptor3d_CurveOnSurface ACS(AC,AS); gp_Pnt cpnt, pcpnt; Standard_Integer nbp = 45; Standard_Real dist, maxdist=0.; for (Standard_Integer i=0; iValue(cuf + i*(cul-cuf)/(nbp-1)); pcpnt = ACS.Value(pcuf + i*(pcul-pcuf)/(nbp-1)); dist = cpnt.SquareDistance(pcpnt); if (maxdist myPrecision) myStatus = ShapeExtend::EncodeStatus (ShapeExtend_DONE1); return LastCheckStatus (ShapeExtend_DONE); } //======================================================================= //function : CheckSelfIntersectingEdge //purpose : //======================================================================= // auxiliary function //:h0 abv 29 May 98: PRO10105 1949: like in BRepCheck, point is to be taken // from 3d curve (but only if edge is SameParameter) static gp_Pnt GetPointOnEdge ( const TopoDS_Edge &edge, const Handle(ShapeAnalysis_Surface) &surf, const Handle(Geom2d_Curve) &Crv2d, const Standard_Real param ) { if ( BRep_Tool::SameParameter ( edge ) ) { Standard_Real f,l; TopLoc_Location L; const Handle(Geom_Curve) ConS = BRep_Tool::Curve ( edge, L, f, l ); if ( ! ConS.IsNull() ) return ConS->Value ( param ).Transformed ( L.Transformation() ); } return surf->Value ( Crv2d->Value ( param ) ); } //======================================================================= //function : CheckSelfIntersectingEdge //purpose : //======================================================================= Standard_Boolean ShapeAnalysis_Wire::CheckSelfIntersectingEdge (const Standard_Integer num, IntRes2d_SequenceOfIntersectionPoint& points2d, TColgp_SequenceOfPnt& points3d) { points2d.Clear(); points3d.Clear(); myStatus = ShapeExtend::EncodeStatus (ShapeExtend_OK); if ( ! IsReady() ) return Standard_False; TopoDS_Edge edge = WireData()->Edge ( num >0 ? num : NbEdges() ); ShapeAnalysis_Edge sae; Standard_Real a, b; Handle(Geom2d_Curve) Crv; if ( ! sae.PCurve ( edge, myFace, Crv, a, b, Standard_False ) ) { myStatus |= ShapeExtend::EncodeStatus (ShapeExtend_FAIL1); return Standard_False; } if ( Abs ( a - b ) <= ::Precision::PConfusion() ) return Standard_False; Standard_Real tolint = 1.0e-10; //szv#4:S4163:12Mar99 warning IntRes2d_Domain domain ( Crv->Value ( a ), a, tolint, Crv->Value ( b ), b, tolint ); Geom2dAdaptor_Curve AC ( Crv ); Geom2dInt_GInter Inter ( AC, domain, tolint, tolint ); if ( ! Inter.IsDone() ) return Standard_False; TopoDS_Vertex V1 = sae.FirstVertex ( edge ); TopoDS_Vertex V2 = sae.LastVertex ( edge ); if ( V1.IsNull() || V2.IsNull() ) { myStatus |= ShapeExtend::EncodeStatus (ShapeExtend_FAIL2); return Standard_False; } Standard_Real tol1 = BRep_Tool::Tolerance ( V1 ); Standard_Real tol2 = BRep_Tool::Tolerance ( V2 ); gp_Pnt pnt1 = BRep_Tool::Pnt ( V1 ); gp_Pnt pnt2 = BRep_Tool::Pnt ( V2 ); for ( Standard_Integer i=1; i <= Inter.NbPoints(); i++ ) { const IntRes2d_IntersectionPoint &IP = Inter.Point ( i ); const IntRes2d_Transition &Tr1 = IP.TransitionOfFirst(); const IntRes2d_Transition &Tr2 = IP.TransitionOfSecond(); if ( Tr1.PositionOnCurve() != IntRes2d_Middle && Tr2.PositionOnCurve() != IntRes2d_Middle ) continue; gp_Pnt pint = GetPointOnEdge ( edge, mySurf, Crv, IP.ParamOnFirst() ); Standard_Real dist21 = pnt1.SquareDistance ( pint ); Standard_Real dist22 = pnt2.SquareDistance ( pint ); if ( dist21 > tol1 * tol1 && dist22 > tol2 * tol2 ) { points2d.Append ( IP ); points3d.Append ( pint ); myStatus |= ShapeExtend::EncodeStatus (ShapeExtend_DONE1); } } return LastCheckStatus ( ShapeExtend_DONE ); } //======================================================================= //function : CheckSelfIntersectingEdge //purpose : //======================================================================= Standard_Boolean ShapeAnalysis_Wire::CheckSelfIntersectingEdge (const Standard_Integer num) { IntRes2d_SequenceOfIntersectionPoint points2d; TColgp_SequenceOfPnt points3d; return CheckSelfIntersectingEdge ( num, points2d, points3d ); } //======================================================================= //function : CheckIntersectingEdges //purpose : Test if two consequent edges are intersecting // It is made in accordance with the following check in BRepCheck: // - in BRepCheck_Wire::Orientation(), test for self-intersection //======================================================================= Standard_Boolean ShapeAnalysis_Wire::CheckIntersectingEdges (const Standard_Integer num, IntRes2d_SequenceOfIntersectionPoint& points2d, TColgp_SequenceOfPnt& points3d, TColStd_SequenceOfReal& errors) { points2d.Clear(); points3d.Clear(); errors.Clear(); myStatus = ShapeExtend::EncodeStatus (ShapeExtend_OK); if ( ! IsReady() || NbEdges() <2 ) return Standard_False; //szv#4:S4163:12Mar99 optimized Standard_Integer n2 = (num > 0)? num : NbEdges(); Standard_Integer n1 = (n2 > 1)? n2-1 : NbEdges(); TopoDS_Edge edge1 = myWire->Edge ( n1 ); TopoDS_Edge edge2 = myWire->Edge ( n2 ); ShapeAnalysis_Edge sae; TopoDS_Vertex V1 = sae.LastVertex ( edge1 ); TopoDS_Vertex V2 = sae.FirstVertex ( edge2 ); if ( V1.IsNull() || V2.IsNull() ) { myStatus |= ShapeExtend::EncodeStatus (ShapeExtend_FAIL1); return Standard_False; } if ( ! BRepTools::Compare ( V1, V2 ) ) { myStatus |= ShapeExtend::EncodeStatus (ShapeExtend_FAIL2); return Standard_False; } TopoDS_Vertex Vp = sae.FirstVertex ( edge1 ); TopoDS_Vertex Vn = sae.LastVertex ( edge2 ); Standard_Real a1, b1, a2, b2; Handle(Geom2d_Curve) Crv1, Crv2; if ( ! sae.PCurve ( edge1, myFace, Crv1, a1, b1, Standard_False ) ) { myStatus |= ShapeExtend::EncodeStatus (ShapeExtend_FAIL3); return Standard_False; } if ( ! sae.PCurve ( edge2, myFace, Crv2, a2, b2, Standard_False ) ) { myStatus |= ShapeExtend::EncodeStatus (ShapeExtend_FAIL3); return Standard_False; } if ( Abs ( a1 - b1 ) <= ::Precision::PConfusion() || Abs ( a2 - b2 ) <= ::Precision::PConfusion() ) return Standard_False; //:f7 abv 6 May 98: BUC50070 on #42276 Standard_Boolean isForward1 = ( edge1.Orientation() == TopAbs_FORWARD ); Standard_Boolean isForward2 = ( edge2.Orientation() == TopAbs_FORWARD ); Standard_Real tol0 = Max ( BRep_Tool::Tolerance ( V1 ), BRep_Tool::Tolerance ( V2 ) ); Standard_Real tol = tol0; gp_Pnt pnt = BRep_Tool::Pnt ( V1 ); // Standard_Boolean Status = Standard_False; Standard_Real tolint = 1.0e-10; //szv#4:S4163:12Mar99 warning IntRes2d_Domain d1 ( Crv1->Value ( a1 ), a1, tolint, Crv1->Value ( b1 ), b1, tolint ); IntRes2d_Domain d2 ( Crv2->Value ( a2 ), a2, tolint, Crv2->Value ( b2 ), b2, tolint ); Geom2dAdaptor_Curve C1 ( Crv1 ), C2 ( Crv2 ); //:64 abv 25 Dec 97: Attention! // Since Intersection algorithm is not symmetrical, for consistency with BRepCheck // edge with lower order number shoud be intersecting with edge with higher one // i.e., for intersection of last and first edges, they should go in reversed order // Example: entity #38285 from bug CSR #CTS17806 // NOTE: Tr1 and Tr2 are not reordered because they are used in the same manner Geom2dInt_GInter Inter; if ( num ==1 ) Inter.Perform ( C2, d2, C1, d1, tolint, tolint ); else Inter.Perform ( C1, d1, C2, d2, tolint, tolint ); if ( ! Inter.IsDone() ) return Standard_False; //:86 abv 22 Jan 98: fix self-intersection even if tolerance of vertex is enough // to annihilate it. This is done to prevent wrong effects if vertex tolerance // will be decreased (e.g., in FixLacking) Standard_Real tole = Max ( ( BRep_Tool::SameParameter ( edge1 ) ? BRep_Tool::Tolerance ( edge1 ) : tol0 ), ( BRep_Tool::SameParameter ( edge2 ) ? BRep_Tool::Tolerance ( edge2 ) : tol0 ) ); Standard_Real tolt = Min ( tol, Max ( tole, myPrecision ) ); //Standard_Real prevRange1 = RealLast(), prevRange2 = RealLast(); //SK Standard_Integer isLacking = -1; //:l0 abv: CATIA01 #1727: protect against adding lacking //#83 rln 19.03.99 sim2.igs, entity 4292 //processing also segments as in BRepCheck Standard_Integer NbPoints = Inter.NbPoints(), NbSegments = Inter.NbSegments(); for ( Standard_Integer i=1; i <= NbPoints + NbSegments; i++ ) { IntRes2d_IntersectionPoint IP; IntRes2d_Transition Tr1, Tr2; if (i <= NbPoints) IP = Inter.Point ( i ); else { const IntRes2d_IntersectionSegment &Seg = Inter.Segment ( i - NbPoints ); if (!Seg.HasFirstPoint() || !Seg.HasLastPoint()) continue; IP = Seg.FirstPoint(); Tr1 = IP.TransitionOfFirst(); Tr2 = IP.TransitionOfSecond(); if (Tr1.PositionOnCurve() == IntRes2d_Middle || Tr2.PositionOnCurve() == IntRes2d_Middle) IP = Seg.LastPoint(); } Tr1 = IP.TransitionOfFirst(); Tr2 = IP.TransitionOfSecond(); if ( Tr1.PositionOnCurve() != IntRes2d_Middle && Tr2.PositionOnCurve() != IntRes2d_Middle ) continue; Standard_Real param1, param2; param1 = ( num ==1 ? IP.ParamOnSecond() : IP.ParamOnFirst() ); param2 = ( num ==1 ? IP.ParamOnFirst() : IP.ParamOnSecond() ); //:r6 abv 8 Apr 99: r_47-sd.stp #173850: protect against working out of curve range if ( a1-param1 > ::Precision::PConfusion() || param1-b1 > ::Precision::PConfusion() || a2-param2 > ::Precision::PConfusion() || param2-b2 > ::Precision::PConfusion() ) continue; //:82 abv 21 Jan 98: point of intersection on Crv1 and Crv2 is different gp_Pnt pi1 = GetPointOnEdge ( edge1, mySurf, Crv1, param1 ); //:h0: thesurf.Value ( Crv1->Value ( param1 ) ); gp_Pnt pi2 = GetPointOnEdge ( edge2, mySurf, Crv2, param2 ); //:h0: thesurf.Value ( Crv2->Value ( param2 ) ); gp_Pnt pint = 0.5 * ( pi1.XYZ() + pi2.XYZ() ); Standard_Real di1 = pi1.SquareDistance ( pnt ); Standard_Real di2 = pi2.SquareDistance ( pnt ); Standard_Real dist2 = Max ( di1, di2 ); //rln 03/02/98: CSR#BUC50004 entity 56 (to avoid later inserting lacking edge) if ( isLacking <0 ) { //:l0 gp_Pnt2d end1 = Crv1->Value ( isForward1 ? b1 : a1 ); gp_Pnt2d end2 = Crv2->Value ( isForward2 ? a2 : b2 ); //:l0 Standard_Real distab2 = mySurf->Value ( end1 ).SquareDistance ( mySurf->Value ( end2 ) ); //:l0: test like in BRepCheck GeomAdaptor_Surface& Ads = mySurf->Adaptor3d()->ChangeSurface(); Standard_Real tol2d = 2 * Max ( Ads.UResolution(tol), Ads.VResolution(tol) ); isLacking = ( end1.SquareDistance(end2) >= tol2d * tol2d ); } if ( ( dist2 > tolt * tolt || //:86: tol -> tolt isLacking ) && //:l0 //:l0 distab2 > BRep_Tool::Tolerance ( edge1 ) + BRep_Tool::Tolerance ( edge2 ) ) && //rln ( ! BRepTools::Compare ( Vp, Vn ) || //:63 dist2 < pint.SquareDistance ( BRep_Tool::Pnt ( Vp ) ) ) ) { //:63 points2d.Append ( IP ); points3d.Append ( pint ); errors.Append ( 0.5 * pi1.Distance ( pi2 ) ); myStatus |= ShapeExtend::EncodeStatus (ShapeExtend_DONE1); } } return LastCheckStatus ( ShapeExtend_DONE ); } //======================================================================= //function : CheckIntersectingEdges //purpose : //======================================================================= Standard_Boolean ShapeAnalysis_Wire::CheckIntersectingEdges (const Standard_Integer num) { IntRes2d_SequenceOfIntersectionPoint points2d; TColgp_SequenceOfPnt points3d; TColStd_SequenceOfReal errors; return CheckIntersectingEdges ( num, points2d, points3d, errors ); } //======================================================================= //function : CheckIntersectingEdges //purpose : //======================================================================= Standard_Boolean ShapeAnalysis_Wire::CheckIntersectingEdges(const Standard_Integer num1, const Standard_Integer num2, IntRes2d_SequenceOfIntersectionPoint& points2d, TColgp_SequenceOfPnt& points3d, TColStd_SequenceOfReal& errors) { myStatus = ShapeExtend::EncodeStatus (ShapeExtend_OK); if ( ! IsReady() ) return Standard_False; Handle(ShapeExtend_WireData) sbwd = WireData(); Standard_Integer n2 = ( num2 >0 ? num2 : sbwd->NbEdges() ); Standard_Integer n1 = ( num1 >0 ? num1 : sbwd->NbEdges() ); TopoDS_Edge edge1 = sbwd->Edge ( n1 ); TopoDS_Edge edge2 = sbwd->Edge ( n2 ); ShapeAnalysis_Edge sae; Standard_Real a1, b1, a2, b2; Handle(Geom2d_Curve) Crv1, Crv2; if(!sae.PCurve ( edge1, myFace, Crv1, a1, b1, Standard_False )){ myStatus |= ShapeExtend::EncodeStatus (ShapeExtend_FAIL3); return Standard_False; } if(!sae.PCurve ( edge2, myFace, Crv2, a2, b2, Standard_False )){ myStatus |= ShapeExtend::EncodeStatus (ShapeExtend_FAIL3); return Standard_False; } if ( Abs ( a1 - b1 ) <= ::Precision::PConfusion() || Abs ( a2 - b2 ) <= ::Precision::PConfusion() ) return Standard_False; points2d.Clear(); points3d.Clear(); errors.Clear(); TColgp_Array1OfPnt vertexPoints(1,4); TColStd_Array1OfReal vertexTolers(1,4); vertexPoints(1) = BRep_Tool::Pnt(sae.FirstVertex(edge1)); vertexTolers(1) = BRep_Tool::Tolerance(sae.FirstVertex(edge1)); vertexPoints(2) = BRep_Tool::Pnt(sae.LastVertex(edge1)); vertexTolers(2) = BRep_Tool::Tolerance(sae.LastVertex(edge1)); vertexPoints(3) = BRep_Tool::Pnt(sae.FirstVertex(edge2)); vertexTolers(3) = BRep_Tool::Tolerance(sae.FirstVertex(edge2)); vertexPoints(4) = BRep_Tool::Pnt(sae.LastVertex(edge2)); vertexTolers(4) = BRep_Tool::Tolerance(sae.LastVertex(edge2)); Standard_Real tolint = 1.0e-10; IntRes2d_Domain d1 ( Crv1->Value ( a1 ), a1, tolint, Crv1->Value ( b1 ), b1, tolint ); IntRes2d_Domain d2 ( Crv2->Value ( a2 ), a2, tolint, Crv2->Value ( b2 ), b2, tolint ); Geom2dAdaptor_Curve C1 ( Crv1 ), C2 ( Crv2 ); Geom2dInt_GInter Inter; Inter.Perform ( C1, d1, C2, d2, tolint, tolint ); if ( ! Inter.IsDone() ) return Standard_False; //#83 rln 19.03.99 sim2.igs, entity 4292 //processing also segments as in BRepCheck Standard_Integer NbPoints = Inter.NbPoints(), NbSegments = Inter.NbSegments(); for ( Standard_Integer i=1; i <= NbPoints + NbSegments; i++ ) { IntRes2d_IntersectionPoint IP; IntRes2d_Transition Tr1, Tr2; if (i <= NbPoints) IP = Inter.Point ( i ); else { const IntRes2d_IntersectionSegment &Seg = Inter.Segment ( i - NbPoints ); if (!Seg.HasFirstPoint() || !Seg.HasLastPoint()) continue; IP = Seg.FirstPoint(); Tr1 = IP.TransitionOfFirst(); Tr2 = IP.TransitionOfSecond(); if (Tr1.PositionOnCurve() == IntRes2d_Middle || Tr2.PositionOnCurve() == IntRes2d_Middle) IP = Seg.LastPoint(); } Tr1 = IP.TransitionOfFirst(); Tr2 = IP.TransitionOfSecond(); if ( Tr1.PositionOnCurve() != IntRes2d_Middle && Tr2.PositionOnCurve() != IntRes2d_Middle ) continue; Standard_Real param1 = IP.ParamOnFirst(); Standard_Real param2 = IP.ParamOnSecond(); gp_Pnt pi1 = GetPointOnEdge ( edge1, mySurf, Crv1, param1 ); //:h0: thesurf.Value ( Crv1->Value ( param1 ) ); gp_Pnt pi2 = GetPointOnEdge ( edge2, mySurf, Crv2, param2 ); Standard_Boolean OK = Standard_False; for(Standard_Integer j=1; (j<=4)&&!OK; j++) { Standard_Real di1 = pi1.SquareDistance (vertexPoints(j)); Standard_Real di2 = pi2.SquareDistance (vertexPoints(j)); Standard_Real dist2 = Max ( di1, di2 ); if(dist2 < vertexTolers(j) * vertexTolers(j)) OK = Standard_True; } if(!OK) { gp_Pnt pint = 0.5 * ( pi1.XYZ() + pi2.XYZ() ); points2d.Append ( IP ); points3d.Append ( pint ); errors.Append ( 0.5 * pi1.Distance ( pi2 ) ); myStatus |= ShapeExtend::EncodeStatus (ShapeExtend_DONE1); } } return LastCheckStatus ( ShapeExtend_DONE ); } //======================================================================= //function : CheckIntersectingEdges //purpose : //======================================================================= Standard_Boolean ShapeAnalysis_Wire::CheckIntersectingEdges (const Standard_Integer num1, const Standard_Integer num2) { IntRes2d_SequenceOfIntersectionPoint points2d; TColgp_SequenceOfPnt points3d; TColStd_SequenceOfReal errors; return CheckIntersectingEdges(num1, num2, points2d, points3d, errors); } //======================================================================= //function : CheckLacking //purpose : Test if two edges are disconnected in 2d according to the // Adaptor_Surface::Resolution //======================================================================= Standard_Boolean ShapeAnalysis_Wire::CheckLacking (const Standard_Integer num, const Standard_Real Tolerance, gp_Pnt2d &p2d1, gp_Pnt2d &p2d2) { myStatus = ShapeExtend::EncodeStatus (ShapeExtend_OK); if ( ! IsReady() ) return Standard_False; //szv#4:S4163:12Mar99 optimized Standard_Integer n2 = (num > 0)? num : NbEdges(); Standard_Integer n1 = (n2 > 1)? n2-1 : NbEdges(); TopoDS_Edge E1 = myWire->Edge ( n1 ); TopoDS_Edge E2 = myWire->Edge ( n2 ); ShapeAnalysis_Edge sae; TopoDS_Vertex V1 = sae.LastVertex ( E1 ); TopoDS_Vertex V2 = sae.FirstVertex ( E2 ); // CKY 4 MAR 1998 : protection against null vertex if ( V1.IsNull() || V2.IsNull() ) { myStatus |= ShapeExtend::EncodeStatus (ShapeExtend_FAIL1); return Standard_False; } if ( ! BRepTools::Compare ( V1, V2 ) ) { myStatus |= ShapeExtend::EncodeStatus (ShapeExtend_FAIL2); return Standard_False; } Standard_Real a, b; gp_Vec2d v1, v2, v12; Handle(Geom2d_Curve) c2d; if ( ! sae.PCurve ( E1, myFace, c2d, a, b, Standard_True ) ) { myStatus |= ShapeExtend::EncodeStatus (ShapeExtend_FAIL3); return Standard_False; } c2d->D1 ( b, p2d1, v1 ); if ( E1.Orientation() == TopAbs_REVERSED ) v1.Reverse(); if ( ! sae.PCurve ( E2, myFace, c2d, a, b, Standard_True ) ) { myStatus |= ShapeExtend::EncodeStatus (ShapeExtend_FAIL3); return Standard_False; } c2d->D1 ( a, p2d2, v2 ); if ( E2.Orientation() == TopAbs_REVERSED ) v2.Reverse(); v12 = p2d2.XY() - p2d1.XY(); myMax2d = v12.SquareMagnitude(); // test like in BRepCheck Standard_Real tol = Max ( BRep_Tool::Tolerance ( V1 ), BRep_Tool::Tolerance ( V2 ) ); tol = ( Tolerance > gp::Resolution() && Tolerance < tol ? Tolerance : tol ); GeomAdaptor_Surface& Ads = mySurf->Adaptor3d()->ChangeSurface(); Standard_Real tol2d = 2 * Max ( Ads.UResolution(tol), Ads.VResolution(tol) ); if ( // tol2d < gp::Resolution() || //#2 smh 26.03.99 S4163 Zero divide myMax2d < tol2d * tol2d ) return Standard_False; myMax2d = Sqrt ( myMax2d ); myMax3d = tol * myMax2d / Max ( tol2d, gp::Resolution() ); myStatus |= ShapeExtend::EncodeStatus ( ShapeExtend_DONE1 ); if ( myMax2d < Precision::PConfusion() || //:abv 03.06.02 CTS21866.stp ( v1.SquareMagnitude() > gp::Resolution() && Abs ( v12.Angle ( v1 ) ) > 0.9 * PI ) || ( v2.SquareMagnitude() > gp::Resolution() && Abs ( v12.Angle ( v2 ) ) > 0.9 * PI ) ) myStatus |= ShapeExtend::EncodeStatus ( ShapeExtend_DONE2 ); return Standard_True; } //======================================================================= //function : CheckLacking //purpose : // //======================================================================= Standard_Boolean ShapeAnalysis_Wire::CheckLacking (const Standard_Integer num, const Standard_Real Tolerance) { gp_Pnt2d p1, p2; return CheckLacking (num, Tolerance, p1, p2); } //======================================================================= //function : CheckOuterBound //purpose : //======================================================================= Standard_Boolean ShapeAnalysis_Wire::CheckOuterBound(const Standard_Boolean APIMake) { myStatus = ShapeExtend::EncodeStatus (ShapeExtend_OK); if ( ! IsReady() ) return Standard_False; TopoDS_Wire wire; if (APIMake) wire = myWire->WireAPIMake(); else wire = myWire->Wire(); TopoDS_Shape sh = myFace.EmptyCopied(); //szv#4:S4163:12Mar99 SGI warns TopoDS_Face face = TopoDS::Face(sh); BRep_Builder B; B.Add (face, wire); if (ShapeAnalysis::IsOuterBound (face)) return Standard_False; myStatus = ShapeExtend::EncodeStatus (ShapeExtend_DONE1); return Standard_True; } //======================================================================= //function : CheckNotchedEdges //purpose : //======================================================================= static Standard_Real ProjectInside(const Adaptor3d_CurveOnSurface AD, const gp_Pnt pnt, const Standard_Real preci, gp_Pnt& proj, Standard_Real& param, const Standard_Boolean adjustToEnds = Standard_True) { ShapeAnalysis_Curve sac; Standard_Real dist = sac.Project(AD,pnt,preci,proj,param,adjustToEnds); Standard_Real uFirst = AD.FirstParameter(); Standard_Real uLast = AD.LastParameter(); if(paramuLast) { param = uLast; proj = AD.Value(uLast); return proj.Distance(pnt); } return dist; } Standard_Boolean ShapeAnalysis_Wire::CheckNotchedEdges(const Standard_Integer num, Standard_Integer& shortNum, Standard_Real& param, const Standard_Real Tolerance) { myStatus = ShapeExtend::EncodeStatus (ShapeExtend_OK); if ( ! IsReady() ) return Standard_False; Standard_Integer n2 = (num > 0)? num : NbEdges(); Standard_Integer n1 = (n2 > 1)? n2-1 : NbEdges(); TopoDS_Edge E1 = myWire->Edge ( n1 ); TopoDS_Edge E2 = myWire->Edge ( n2 ); if(BRep_Tool::Degenerated(E1)||BRep_Tool::Degenerated(E2)) return Standard_False; ShapeAnalysis_Edge sae; TopoDS_Vertex V1 = sae.LastVertex ( E1 ); TopoDS_Vertex V2 = sae.FirstVertex ( E2 ); if ( V1.IsNull() || V2.IsNull() ) { myStatus |= ShapeExtend::EncodeStatus (ShapeExtend_FAIL1); return Standard_False; } if ( ! BRepTools::Compare ( V1, V2 ) ) { myStatus |= ShapeExtend::EncodeStatus (ShapeExtend_FAIL2); return Standard_False; } Standard_Real a1, b1, a2, b2; gp_Pnt2d p2d1, p2d2; gp_Vec2d v1, v2; Handle(Geom2d_Curve) c2d1, c2d2; if ( ! sae.PCurve ( E1, myFace, c2d1, a1, b1, Standard_False ) ) { myStatus |= ShapeExtend::EncodeStatus (ShapeExtend_FAIL3); return Standard_False; } if(E1.Orientation()==TopAbs_REVERSED) c2d1->D1 ( a1, p2d1, v1 ); else { c2d1->D1 ( b1, p2d1, v1 ); v1.Reverse(); } if ( ! sae.PCurve ( E2, myFace, c2d2, a2, b2, Standard_False ) ) { myStatus |= ShapeExtend::EncodeStatus (ShapeExtend_FAIL3); return Standard_False; } if(E2.Orientation()==TopAbs_REVERSED) { c2d2->D1 ( b2, p2d2, v2 ); v2.Reverse(); } else c2d2->D1 ( a2, p2d2, v2 ); if ( v2.Magnitude() < gp::Resolution() || v1.Magnitude() < gp::Resolution()) return Standard_False; if ( Abs ( v2.Angle ( v1 ) ) > 0.1 || p2d1.Distance(p2d2) > Tolerance) return Standard_False; Handle(Geom2dAdaptor_HCurve) AC2d1 = new Geom2dAdaptor_HCurve(c2d1,a1,b1); Handle(GeomAdaptor_HSurface) AdS1 = new GeomAdaptor_HSurface(new Geom_Plane(gp_Pln())); Adaptor3d_CurveOnSurface Ad1(AC2d1,AdS1); Handle(Geom2dAdaptor_HCurve) AC2d2 = new Geom2dAdaptor_HCurve(c2d2,a2,b2); Handle(GeomAdaptor_HSurface) AdS2 = new GeomAdaptor_HSurface(new Geom_Plane(gp_Pln())); Adaptor3d_CurveOnSurface Ad2(AC2d2,AdS2); Adaptor3d_CurveOnSurface longAD, shortAD; Standard_Real lenP, firstP; ShapeAnalysis_Curve sac; gp_Pnt Proj1, Proj2; Standard_Real param1 = 0., param2 = 0.; p2d2=c2d2->Value(E2.Orientation()==TopAbs_FORWARD ? b2 : a2); p2d1=c2d1->Value(E1.Orientation()==TopAbs_FORWARD ? a1 : b1); Standard_Real dist1 = ProjectInside(Ad1,gp_Pnt(p2d2.X(),p2d2.Y(),0),Tolerance,Proj1,param1,Standard_False); Standard_Real dist2 = ProjectInside(Ad2,gp_Pnt(p2d1.X(),p2d1.Y(),0),Tolerance,Proj2,param2,Standard_False); if ( dist1 > Tolerance && dist2 > Tolerance) return Standard_False; if (dist1 < dist2 ) { shortAD = Ad2; longAD = Ad1; lenP = b2 - a2; firstP = a2; shortNum=n2; param=param1; } else { shortAD = Ad1; longAD = Ad2; lenP = b1 - a1; firstP = a1; shortNum=n1; param=param2; } Standard_Real step = lenP/23; for (Standard_Integer i = 1; i < 23; i++,firstP+=step) { Standard_Real d1 = sac.Project(longAD,shortAD.Value(firstP),Tolerance,Proj1,param1); if (d1 > Tolerance) { return Standard_False; } } return Standard_True; } //======================================================================= //function : CheckSmallArea //purpose : //======================================================================= Standard_Boolean ShapeAnalysis_Wire::CheckSmallArea(const Standard_Real prec2d) { myStatus = ShapeExtend::EncodeStatus (ShapeExtend_FAIL1); Standard_Integer NbEdges = myWire->NbEdges(); if ( !IsReady() || NbEdges <1 ) return Standard_False; myStatus = ShapeExtend::EncodeStatus (ShapeExtend_OK); Standard_Integer NbControl=23; Standard_Real area=0; gp_XY prev, cont; for (Standard_Integer nbe = 1; nbe <= NbEdges; nbe++) { Standard_Real First, Last; Handle(Geom2d_Curve) c2d; ShapeAnalysis_Edge sae; if (!sae.PCurve(myWire->Edge(nbe),myFace,c2d,First,Last)) { myStatus = ShapeExtend::EncodeStatus (ShapeExtend_FAIL2); return Standard_False; } Standard_Integer ibeg = 0; if( nbe == 1 ) { gp_Pnt2d pntIni = c2d->Value(First); prev = pntIni.XY(); cont = prev; ibeg = 1; } for ( Standard_Integer i = ibeg; i < NbControl; i++) { Standard_Real prm = ((NbControl-1-i)*First + i*Last)/(NbControl-1); gp_Pnt2d pntCurr = c2d->Value(prm); gp_XY curr = pntCurr.XY(); area += curr ^ prev; prev = curr; } } area += cont ^ prev; if ( Abs(area) < 2*prec2d*prec2d ) { myStatus = ShapeExtend::EncodeStatus (ShapeExtend_DONE1); return Standard_True; } return Standard_False; } //======================================================================= //function : CheckShapeConnect //purpose : //======================================================================= Standard_Boolean ShapeAnalysis_Wire::CheckShapeConnect(const TopoDS_Shape& shape,const Standard_Real prec) { Standard_Real tailhead, tailtail, headhead, headtail; return CheckShapeConnect (tailhead, tailtail, headtail, headhead, shape, prec); } //======================================================================= //function : CheckShapeConnect //purpose : //======================================================================= Standard_Boolean ShapeAnalysis_Wire::CheckShapeConnect(Standard_Real& tailhead, Standard_Real& tailtail, Standard_Real& headtail, Standard_Real& headhead, const TopoDS_Shape& shape, const Standard_Real prec) { myStatus = ShapeExtend::EncodeStatus (ShapeExtend_FAIL1); if (!IsLoaded () || shape.IsNull()) return Standard_False; TopoDS_Vertex V1,V2; TopoDS_Edge E; TopoDS_Wire W; ShapeAnalysis_Edge SAE; if (shape.ShapeType() == TopAbs_EDGE) { E = TopoDS::Edge (shape); V1 = SAE.FirstVertex (E); V2 = SAE.LastVertex (E); } else if (shape.ShapeType() == TopAbs_WIRE) { W = TopoDS::Wire (shape); ShapeAnalysis::FindBounds (W,V1,V2); } else return Standard_False; myStatus = ShapeExtend::EncodeStatus (ShapeExtend_DONE1); // on va comparer les points avec ceux de thevfirst et thevlast gp_Pnt p1 = BRep_Tool::Pnt(V1); gp_Pnt p2 = BRep_Tool::Pnt(V2); TopoDS_Vertex vfirst = SAE.FirstVertex (myWire->Edge (1)), vlast = SAE.LastVertex (myWire->Edge (NbEdges())); gp_Pnt pf = BRep_Tool::Pnt(vfirst); gp_Pnt pl = BRep_Tool::Pnt(vlast); tailhead = p1.Distance(pl); tailtail = p2.Distance(pl); headhead = p1.Distance(pf); headtail = p2.Distance(pf); Standard_Real dm1 = tailhead, dm2 = headtail; Standard_Integer res1 = 0, res2 = 0; if (tailhead > tailtail) {res1 = 1; dm1 = tailtail;} if (headtail > headhead) {res2 = 1; dm2 = headhead;} Standard_Integer result = res1; myMin3d = Min (dm1, dm2); myMax3d = Max (dm1, dm2); if (dm1 > dm2) {dm1 = dm2; result = res2 + 2;} switch (result) { case 1: myStatus = ShapeExtend::EncodeStatus (ShapeExtend_DONE2); break; case 2: myStatus = ShapeExtend::EncodeStatus (ShapeExtend_DONE3); break; case 3: myStatus = ShapeExtend::EncodeStatus (ShapeExtend_DONE4); break; } if (!res1) myStatus |= ShapeExtend::EncodeStatus (ShapeExtend_DONE5); if (!res2) myStatus |= ShapeExtend::EncodeStatus (ShapeExtend_DONE6); if (myMin3d > Max (myPrecision, prec)) myStatus = ShapeExtend::EncodeStatus (ShapeExtend_FAIL2); return LastCheckStatus (ShapeExtend_DONE); } //======================================================================= //function : CheckLoop //purpose : //======================================================================= Standard_Boolean isMultiVertex(const TopTools_ListOfShape& alshape, const TopTools_MapOfShape& aMapSmallEdges, const TopTools_MapOfShape& aMapSeemEdges) { TopTools_ListIteratorOfListOfShape lIt1(alshape); Standard_Integer nbNotAccount =0; for( ; lIt1.More() ; lIt1.Next()) { if(aMapSmallEdges.Contains(lIt1.Value())) nbNotAccount++; else if(aMapSeemEdges.Contains(lIt1.Value())) nbNotAccount++; } return ((alshape.Extent() -nbNotAccount) >2); } Standard_Boolean ShapeAnalysis_Wire::CheckLoop(TopTools_IndexedMapOfShape& aMapLoopVertices, TopTools_DataMapOfShapeListOfShape& aMapVertexEdges, TopTools_MapOfShape& aMapSmallEdges, TopTools_MapOfShape& aMapSeemEdges) { myStatus = ShapeExtend::EncodeStatus(ShapeExtend_OK); if (!IsLoaded() || NbEdges() < 2) return Standard_False; Standard_Real aSavPreci = Precision(); SetPrecision(Precision::Infinite()); Standard_Integer i =1; for( ; i <= myWire->NbEdges(); i++) { TopoDS_Edge aedge = myWire->Edge(i); TopoDS_Vertex aV1,aV2; TopExp::Vertices(aedge,aV1,aV2); Standard_Boolean isSame = aV1.IsSame(aV2); if(myWire->IsSeam(i)) aMapSeemEdges.Add(aedge); ///continue; else if(BRep_Tool::Degenerated(aedge)) aMapSmallEdges.Add(aedge); else if(isSame && CheckSmall(i,BRep_Tool::Tolerance(aV1))) aMapSmallEdges.Add(aedge); if(!aMapVertexEdges.IsBound(aV1)) { TopTools_ListOfShape alshape; aMapVertexEdges.Bind(aV1,alshape); } if(!aMapVertexEdges.IsBound(aV2)) { TopTools_ListOfShape alshape; aMapVertexEdges.Bind(aV2,alshape); } if(isSame) { TopTools_ListOfShape& alshape = aMapVertexEdges.ChangeFind(aV1); alshape.Append(aedge); alshape.Append(aedge); if(alshape.Extent() >2 && isMultiVertex( alshape,aMapSmallEdges,aMapSeemEdges)) aMapLoopVertices.Add(aV1); } else { TopTools_ListOfShape& alshape = aMapVertexEdges.ChangeFind(aV1); alshape.Append(aedge); if(alshape.Extent() >2 && isMultiVertex( alshape,aMapSmallEdges,aMapSeemEdges)) aMapLoopVertices.Add(aV1); TopTools_ListOfShape& alshape2 = aMapVertexEdges.ChangeFind(aV2); alshape2.Append(aedge); if(alshape2.Extent() >2 && isMultiVertex( alshape2,aMapSmallEdges,aMapSeemEdges)) aMapLoopVertices.Add(aV2); } } SetPrecision(aSavPreci); if(aMapLoopVertices.Extent()) { myStatus = ShapeExtend::EncodeStatus (ShapeExtend_DONE1); myStatusLoop |= myStatus; return Standard_True; } return Standard_False; }