1 // Created on: 1996-04-23
2 // Created by: Jacques GOUSSARD
3 // Copyright (c) 1996-1999 Matra Datavision
4 // Copyright (c) 1999-2012 OPEN CASCADE SAS
6 // The content of this file is subject to the Open CASCADE Technology Public
7 // License Version 6.5 (the "License"). You may not use the content of this file
8 // except in compliance with the License. Please obtain a copy of the License
9 // at http://www.opencascade.org and read it completely before using this file.
11 // The Initial Developer of the Original Code is Open CASCADE S.A.S., having its
12 // main offices at: 1, place des Freres Montgolfier, 78280 Guyancourt, France.
14 // The Original Code and all software distributed under the License is
15 // distributed on an "AS IS" basis, without warranty of any kind, and the
16 // Initial Developer hereby disclaims all such warranties, including without
17 // limitation, any warranties of merchantability, fitness for a particular
18 // purpose or non-infringement. Please see the License for the specific terms
19 // and conditions governing the rights and limitations under the License.
23 #include <BRepFeat.ixx>
25 //#include <BRepAlgo_Cut.hxx>
27 #include <BRepBuilderAPI.hxx>
28 #include <BRepAdaptor_Curve.hxx>
29 #include <TopExp_Explorer.hxx>
30 #include <TopTools_MapOfShape.hxx>
31 #include <BRep_Tool.hxx>
32 #include <BRep_Builder.hxx>
33 #include <Geom_Curve.hxx>
34 #include <Geom_TrimmedCurve.hxx>
35 #include <Geom2d_TrimmedCurve.hxx>
36 #include <Extrema_ExtPC.hxx>
37 #include <GeomAdaptor_Curve.hxx>
38 #include <BRepLProp.hxx>
40 #include <TopoDS_Vertex.hxx>
41 #include <TopoDS_Face.hxx>
42 #include <TopoDS_Edge.hxx>
43 #include <TopoDS_Shell.hxx>
44 #include <TopoDS_Solid.hxx>
45 #include <Precision.hxx>
46 #include <GCPnts_QuasiUniformDeflection.hxx>
47 #include <BRepTopAdaptor_FClass2d.hxx>
48 #include <Geom2dAdaptor_Curve.hxx>
49 #include <GeomProjLib.hxx>
50 #include <gp_Vec2d.hxx>
51 #include <BRepTools.hxx>
52 #include <Geom_Surface.hxx>
53 #include <Bnd_Box.hxx>
54 #include <BRepBndLib.hxx>
55 #include <BRepLib_MakeFace.hxx>
56 #include <Geom_RectangularTrimmedSurface.hxx>
57 #include <Geom_Plane.hxx>
58 #include <Geom_CylindricalSurface.hxx>
59 #include <Geom_ConicalSurface.hxx>
60 #include <LocOpe_CSIntersector.hxx>
61 #include <LocOpe_PntFace.hxx>
62 #include <LocOpe_BuildShape.hxx>
64 #include <TColGeom_SequenceOfCurve.hxx>
69 #define NECHANTBARYC 11
71 //=======================================================================
72 //function : SampleEdges
74 //=======================================================================
76 void BRepFeat::SampleEdges(const TopoDS_Shape& theShape,
77 TColgp_SequenceOfPnt& theSeq)
79 LocOpe::SampleEdges(theShape,theSeq);
84 //=======================================================================
85 //function : Barycenter
86 //purpose : Calcul du barycentre des edges d'un shape
87 //=======================================================================
89 void BRepFeat::Barycenter(const TopoDS_Shape& S,
92 TopTools_MapOfShape theMap;
93 TopExp_Explorer exp(S,TopAbs_EDGE);
96 Standard_Real f,l,prm;
98 Standard_Integer i, nbp= 0;
100 for (; exp.More(); exp.Next()) {
101 const TopoDS_Edge& edg = TopoDS::Edge(exp.Current());
102 if (!theMap.Add(edg)) {
105 if (!BRep_Tool::Degenerated(edg)) {
106 C = BRep_Tool::Curve(edg,Loc,f,l);
107 C = Handle(Geom_Curve)::DownCast(C->Transformed(Loc.Transformation()));
108 for (i=1;i<NECHANTBARYC; i++) {
109 prm = ((NECHANTBARYC-i)*f + i*l)/NECHANTBARYC;
110 Bar += C->Value(prm).XYZ();
116 for (exp.Init(S,TopAbs_VERTEX); exp.More(); exp.Next()) {
117 if (theMap.Add(exp.Current())) {
118 Bar += (BRep_Tool::Pnt(TopoDS::Vertex(exp.Current()))).XYZ();
123 Bar.Divide((Standard_Real)nbp);
128 //=======================================================================
129 //function : ParametricBarycenter
130 //purpose : Calcul du barycentre "parametrique" shape sur une courbe
131 //=======================================================================
133 Standard_Real BRepFeat::ParametricBarycenter(const TopoDS_Shape& S,
134 const Handle(Geom_Curve)& CC)
136 TopTools_MapOfShape theMap;
137 TopExp_Explorer exp(S,TopAbs_EDGE);
139 Handle(Geom_Curve) C;
140 Standard_Real f,l,prm;
141 Standard_Integer i, nbp= 0;
142 GeomAdaptor_Curve TheCurve(CC);
144 extpc.Initialize(TheCurve,CC->FirstParameter(),CC->LastParameter());
145 Standard_Real parbar = 0;
147 for (; exp.More(); exp.Next()) {
148 const TopoDS_Edge& edg = TopoDS::Edge(exp.Current());
149 if (!theMap.Add(edg)) {
152 if (!BRep_Tool::Degenerated(edg)) {
153 C = BRep_Tool::Curve(edg,Loc,f,l);
154 C = Handle(Geom_Curve)::DownCast(C->Transformed(Loc.Transformation()));
155 for (i=1;i<NECHANTBARYC; i++) {
156 prm = ((NECHANTBARYC-i)*f + i*l)/NECHANTBARYC;
157 gp_Pnt pone = C->Value(prm);
158 // On projette sur CC
160 if (extpc.IsDone() && extpc.NbExt() >= 1) {
161 Standard_Real Dist2Min = extpc.SquareDistance(1);
162 Standard_Integer kmin = 1;
163 for (Standard_Integer k=2; k<=extpc.NbExt(); k++) {
164 Standard_Real Dist2 = extpc.SquareDistance(k);
165 if (Dist2 < Dist2Min) {
171 Standard_Real prmp = extpc.Point(kmin).Parameter();
178 for (exp.Init(S,TopAbs_VERTEX); exp.More(); exp.Next()) {
179 if (theMap.Add(exp.Current())) {
180 gp_Pnt pone = BRep_Tool::Pnt(TopoDS::Vertex(exp.Current()));
181 // On projette sur CC
183 if (extpc.IsDone() && extpc.NbExt() >= 1) {
184 Standard_Real Dist2Min = extpc.SquareDistance(1);
185 Standard_Integer kmin = 1;
186 for (Standard_Integer k=2; k<=extpc.NbExt(); k++) {
187 Standard_Real Dist2 = extpc.SquareDistance(k);
188 if (Dist2 < Dist2Min) {
195 Standard_Real prmp = extpc.Point(kmin).Parameter();
197 extpc.Point(kmin).Parameter();
203 parbar /=((Standard_Real)nbp);
208 //=======================================================================
209 //function : ParametricBarycenter
210 //purpose : Calcul du barycentre "parametrique" shape sur une courbe
211 //=======================================================================
213 void BRepFeat::ParametricMinMax(const TopoDS_Shape& S,
214 const Handle(Geom_Curve)& CC,
215 Standard_Real& prmin,
216 Standard_Real& prmax,
217 Standard_Real& prbmin,
218 Standard_Real& prbmax,
219 Standard_Boolean& flag,
220 const Standard_Boolean Ori)
222 LocOpe_CSIntersector ASI(S);
223 TColGeom_SequenceOfCurve scur;
226 if(ASI.IsDone() && ASI.NbPoints(1) >=1) {
228 prmin = Min(ASI.Point(1,1).Parameter(),
229 ASI.Point(1, ASI.NbPoints(1)).Parameter());
230 prmax = Max(ASI.Point(1,1).Parameter(),
231 ASI.Point(1, ASI.NbPoints(1)).Parameter());
234 TopAbs_Orientation Ori = ASI.Point(1,1).Orientation();
235 if (Ori == TopAbs_FORWARD) {
236 prmin = ASI.Point(1,1).Parameter();
237 prmax = ASI.Point(1, ASI.NbPoints(1)).Parameter();
240 prmax = ASI.Point(1,1).Parameter();
241 prmin = ASI.Point(1, ASI.NbPoints(1)).Parameter();
244 flag = Standard_True;
249 flag = Standard_False;
252 TopTools_MapOfShape theMap;
253 TopExp_Explorer exp(S,TopAbs_EDGE);
255 Handle(Geom_Curve) C;
256 Standard_Real f,l,prm;
257 // Standard_Integer i, nbp= 0;
259 GeomAdaptor_Curve TheCurve(CC);
261 extpc.Initialize(TheCurve,CC->FirstParameter(),CC->LastParameter());
263 prbmax = RealFirst();
264 for (; exp.More(); exp.Next()) {
265 const TopoDS_Edge& edg = TopoDS::Edge(exp.Current());
266 if (!theMap.Add(edg)) {
269 if (!BRep_Tool::Degenerated(edg)) {
270 C = BRep_Tool::Curve(edg,Loc,f,l);
271 C = Handle(Geom_Curve)::DownCast(C->Transformed(Loc.Transformation()));
272 for (i=1;i<NECHANTBARYC; i++) {
273 prm = ((NECHANTBARYC-i)*f + i*l)/NECHANTBARYC;
274 gp_Pnt pone = C->Value(prm);
275 // On projette sur CC
277 if (extpc.IsDone() && extpc.NbExt() >= 1) {
278 Standard_Real Dist2Min = extpc.SquareDistance(1);
279 Standard_Integer kmin = 1;
280 for (Standard_Integer k=2; k<=extpc.NbExt(); k++) {
281 Standard_Real Dist2 = extpc.SquareDistance(k);
282 if (Dist2 < Dist2Min) {
287 Standard_Real prmp = extpc.Point(kmin).Parameter();
288 if (prmp <= prbmin) {
291 if (prmp >= prbmax) {
299 for (exp.Init(S,TopAbs_VERTEX); exp.More(); exp.Next()) {
300 if (theMap.Add(exp.Current())) {
301 gp_Pnt pone = BRep_Tool::Pnt(TopoDS::Vertex(exp.Current()));
302 // On projette sur CC
304 if (extpc.IsDone() && extpc.NbExt() >= 1) {
305 Standard_Real Dist2Min = extpc.SquareDistance(1);
306 Standard_Integer kmin = 1;
307 for (Standard_Integer k=2; k<=extpc.NbExt(); k++) {
308 Standard_Real Dist2 = extpc.SquareDistance(k);
309 if (Dist2 < Dist2Min) {
314 Standard_Real prmp = extpc.Point(kmin).Parameter();
315 if (prmp <= prbmin) {
318 if (prmp >= prbmax) {
329 //=======================================================================
332 //=======================================================================
334 static Standard_Boolean IsIn (BRepTopAdaptor_FClass2d& FC,
335 Geom2dAdaptor_Curve AC)
337 Standard_Real Def = 100*Precision::Confusion();
338 GCPnts_QuasiUniformDeflection QU(AC,Def);
340 for (Standard_Integer i = 1; i <= QU.NbPoints(); i++) {
341 gp_Pnt2d P = AC.Value(QU.Parameter(i));
342 if (FC.Perform(P, Standard_False) == TopAbs_OUT) {
343 return Standard_False;
347 return Standard_True;
351 //=======================================================================
352 //function : PutInBoundsU
353 //purpose : Recadre la courbe 2d dans les bounds de la face
354 //=======================================================================
359 static void PutInBoundsU (Standard_Real umin,
362 Standard_Real period,
365 Handle(Geom2d_Curve)& C2d)
367 gp_Pnt2d Pf = C2d->Value(f);
368 gp_Pnt2d Pl = C2d->Value(l);
369 gp_Pnt2d Pm = C2d->Value(0.34*f + 0.66*l);
370 Standard_Real minC = Min(Pf.X(),Pl.X()); minC = Min(minC,Pm.X());
371 Standard_Real maxC = Max(Pf.X(),Pl.X()); maxC = Max(maxC,Pm.X());
372 Standard_Real du = 0.;
373 if (minC< umin - eps) {
374 du = (int((umin - minC)/period) + 1)*period;
376 if (minC > umax + eps) {
377 du = -(int((minC - umax)/period) + 1)*period;
382 minC += du; maxC += du;
384 // Ajuste au mieux la courbe dans le domaine.
385 if (maxC > umax +100*eps) {
386 Standard_Real d1 = maxC - umax;
387 Standard_Real d2 = umin - minC + period;
388 if (d2 < d1) du =-period;
397 //=======================================================================
398 //function : PutInBoundsU
399 //purpose : Recadre la courbe 2d dans les bounds de la face
400 //=======================================================================
405 static void PutInBoundsV (Standard_Real vmin,
408 Standard_Real period,
411 Handle(Geom2d_Curve)& C2d)
413 gp_Pnt2d Pf = C2d->Value(f);
414 gp_Pnt2d Pl = C2d->Value(l);
415 gp_Pnt2d Pm = C2d->Value(0.34*f + 0.66*l);
416 Standard_Real minC = Min(Pf.Y(),Pl.Y()); minC = Min(minC,Pm.Y());
417 Standard_Real maxC = Max(Pf.Y(),Pl.Y()); maxC = Max(maxC,Pm.Y());
418 Standard_Real dv = 0.;
419 if (minC< vmin - eps) {
420 dv = (int((vmin - minC)/period) + 1)*period;
422 if (minC > vmax + eps) {
423 dv = -(int((minC - vmax)/period) + 1)*period;
428 minC += dv; maxC += dv;
430 // Ajuste au mieux la courbe dans le domaine.
431 if (maxC > vmax +100*eps) {
432 Standard_Real d1 = maxC - vmax;
433 Standard_Real d2 = vmin - minC + period;
434 if (d2 < d1) dv =-period;
443 //=======================================================================
444 //function : IsInside
446 //=======================================================================
449 Standard_Boolean BRepFeat::IsInside(const TopoDS_Face& F1,
450 const TopoDS_Face& F2)
453 exp.Init(F1, TopAbs_EDGE);
455 Standard_Real umin,umax,vmin,vmax, uperiod=0, vperiod=0;
456 Standard_Integer flagu = 0, flagv = 0;
457 TopLoc_Location L; // Recup S avec la location pour eviter la copie.
458 Handle (Geom_Surface) S = BRep_Tool::Surface(F2);
459 // Standard_Real periodu, periodv;
460 BRepTools::UVBounds(F2,umin,umax,vmin,vmax);
462 if (S->IsUPeriodic()) {
464 uperiod = S->UPeriod();
467 if (S->IsVPeriodic()) {
469 vperiod = S->VPeriod();
471 TopoDS_Shape aLocalShape = F2.Oriented(TopAbs_FORWARD);
472 BRepTopAdaptor_FClass2d FC (TopoDS::Face(aLocalShape),Precision::Confusion());
473 // BRepTopAdaptor_FClass2d FC (TopoDS::Face(F2.Oriented(TopAbs_FORWARD)),
474 // Precision::Confusion());
475 for(; exp.More(); exp.Next()) {
477 Handle(Geom_Curve) C0 = BRep_Tool::Curve(TopoDS::Edge(exp.Current()),f1,l1);
478 Handle(Geom2d_Curve) C = GeomProjLib::Curve2d(C0,f1,l1,S);
479 TopoDS_Edge E = TopoDS::Edge(exp.Current());
480 if(flagu == 1 || flagv == 1) {
481 Standard_Real eps = BRep_Tool::Tolerance(E);
482 BRep_Tool::Range(E,f1,l1);
483 if(flagu == 1) PutInBoundsU(umin, umax, eps, uperiod, f1, l1, C);
484 if(flagv == 1) PutInBoundsV(vmin, vmax, eps, vperiod, f1, l1, C);
486 Geom2dAdaptor_Curve AC(C,f1,l1);
488 return Standard_False;
492 return Standard_True;
497 //=======================================================================
498 //function : FaceUntil
500 //=======================================================================
503 void BRepFeat::FaceUntil(const TopoDS_Shape& Sbase,
507 BRepBndLib::Add(Sbase,B);
508 Standard_Real c[6], bnd;
509 B.Get(c[0],c[2],c[4],c[1],c[3],c[5]);
511 for(Standard_Integer i = 1 ; i < 6; i++) {
512 if(c[i] > bnd) bnd = c[i];
517 Handle(Geom_Surface) s = BRep_Tool::Surface(FUntil);
518 Handle(Standard_Type) styp = s->DynamicType();
519 if (styp == STANDARD_TYPE(Geom_RectangularTrimmedSurface)) {
520 s = Handle(Geom_RectangularTrimmedSurface)::DownCast(s)->BasisSurface();
521 styp = s->DynamicType();
523 Handle(Geom_RectangularTrimmedSurface) str;
524 if (styp == STANDARD_TYPE(Geom_Plane)) {
525 str = new Geom_RectangularTrimmedSurface
526 (s, bnd, -bnd, bnd, -bnd, Standard_True, Standard_True);
528 else if (styp == STANDARD_TYPE(Geom_CylindricalSurface)) {
529 str = new Geom_RectangularTrimmedSurface
530 (s, 0., 2.*M_PI, bnd, -bnd, Standard_True, Standard_True);
532 else if (styp == STANDARD_TYPE(Geom_ConicalSurface)) {
533 str = new Geom_RectangularTrimmedSurface
534 (s, 0., 2.*M_PI, bnd, -bnd, Standard_True, Standard_True);
541 FUntil = BRepLib_MakeFace(str, Precision::Confusion());
546 //=======================================================================
549 //=======================================================================
551 TopoDS_Solid BRepFeat::Tool(const TopoDS_Shape& SRef,
552 const TopoDS_Face& Fac,
553 const TopAbs_Orientation Orf)
555 TopTools_ListOfShape lfaces;
556 // for (TopExp_Explorer exp(SRef,TopAbs_FACE); exp.More(); exp.Next()) {
557 TopExp_Explorer exp(SRef,TopAbs_FACE) ;
558 for ( ; exp.More(); exp.Next()) {
559 if (exp.Current().ShapeType() == TopAbs_FACE) {
560 lfaces.Append(exp.Current());
564 LocOpe_BuildShape bs(lfaces);
565 const TopoDS_Shape& Res = bs.Shape();
567 if (Res.ShapeType() == TopAbs_SHELL) {
569 Sh = TopoDS::Shell(Res);
571 else if (Res.ShapeType() == TopAbs_SOLID) {
572 exp.Init(Res,TopAbs_SHELL);
573 Sh = TopoDS::Shell(exp.Current());
586 Sh.Orientation(TopAbs_FORWARD);
588 TopAbs_Orientation orient;
590 TopAbs_Orientation orient = TopAbs_FORWARD;
592 for (exp.Init(Sh,TopAbs_FACE); exp.More(); exp.Next()) {
593 if (exp.Current().IsSame(Fac)) {
594 orient = exp.Current().Orientation();
599 Standard_Boolean reverse = Standard_False;
600 if ((orient == Fac.Orientation() && Orf == TopAbs_REVERSED) ||
601 (orient != Fac.Orientation() && Orf == TopAbs_FORWARD)) {
602 reverse = Standard_True;
617 //=======================================================================
619 //purpose : Print the error Description of a StatusError on a stream.
620 //=======================================================================
622 Standard_OStream& BRepFeat::Print(const BRepFeat_StatusError se,
629 case BRepFeat_BadDirect :
630 s << "Directions must be opposite";
632 case BRepFeat_BadIntersect :
633 s << "Intersection failure";
635 case BRepFeat_EmptyBaryCurve :
636 s << "Empty BaryCurve";
638 case BRepFeat_EmptyCutResult :
639 s << "Failure in Cut : Empty resulting shape";
641 case BRepFeat_FalseSide :
642 s << "Verify plane and wire orientation";
644 case BRepFeat_IncDirection :
645 s << "Incoherent Direction for shapes From and Until";
647 case BRepFeat_IncSlidFace :
648 s << "Sliding face not in Base shape";
650 case BRepFeat_IncParameter :
651 s << "Incoherent Parameter : shape Until before shape From";
653 case BRepFeat_IncTypes :
654 s << "Invalid option for faces From and Until : 1 Support and 1 not";
656 case BRepFeat_IntervalOverlap :
657 s << "Shapes From and Until overlap";
659 case BRepFeat_InvFirstShape :
660 s << "Invalid First shape : more than 1 face";
662 case BRepFeat_InvOption :
663 s << "Invalid option";
665 case BRepFeat_InvShape :
666 s << "Invalid shape";
668 case BRepFeat_LocOpeNotDone :
669 s << "Local Operation not done";
671 case BRepFeat_LocOpeInvNotDone :
672 s << "Local Operation : intersection line conflict";
674 case BRepFeat_NoExtFace :
675 s << "No Extreme faces";
677 case BRepFeat_NoFaceProf :
678 s << "No Face Profile";
680 case BRepFeat_NoGluer :
681 s << "Gluer Failure";
683 case BRepFeat_NoIntersectF :
684 s << "No intersection between Feature and shape From";
686 case BRepFeat_NoIntersectU :
687 s << "No intersection between Feature and shape Until";
689 case BRepFeat_NoParts :
690 s << "No parts of tool kept";
692 case BRepFeat_NoProjPt :
693 s << "No projection points";
695 case BRepFeat_NotInitialized :
696 s << "Fields not initialized";
698 case BRepFeat_NotYetImplemented :
699 s << "Not yet implemented";
701 case BRepFeat_NullRealTool :
702 s << "Real Tool : Null DPrism";
704 case BRepFeat_NullToolF :
705 s << "Null Tool : Invalid type for shape Form";
707 case BRepFeat_NullToolU :
708 s << "Null Tool : Invalid type for shape Until";