1 // Copyright (C) 2005 OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
2 // CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
4 // This library is free software; you can redistribute it and/or
5 // modify it under the terms of the GNU Lesser General Public
6 // License as published by the Free Software Foundation; either
7 // version 2.1 of the License.
9 // This library is distributed in the hope that it will be useful
10 // but WITHOUT ANY WARRANTY; without even the implied warranty of
11 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
12 // Lesser General Public License for more details.
14 // You should have received a copy of the GNU Lesser General Public
15 // License along with this library; if not, write to the Free Software
16 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
18 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
20 // File: NMTTools_Tools.cxx
21 // Created: Mon Dec 8 10:35:15 2003
22 // Author: Peter KURNEV
26 #include <NMTTools_Tools.ixx>
28 #include <TColStd_IndexedMapOfInteger.hxx>
32 #include <gp_Pnt2d.hxx>
34 #include <Geom_Surface.hxx>
35 #include <GeomAPI_ProjectPointOnSurf.hxx>
38 #include <TopoDS_Vertex.hxx>
39 #include <TopoDS_Shape.hxx>
40 #include <TopoDS_Edge.hxx>
44 #include <TopTools_ListIteratorOfListOfShape.hxx>
45 #include <TopTools_IndexedMapOfShape.hxx>
47 #include <BRep_Tool.hxx>
48 #include <BRep_Builder.hxx>
49 #include <BRepTools.hxx>
51 #include <BOPTColStd_IndexedDataMapOfIntegerIndexedMapOfInteger.hxx>
52 #include <BOPTColStd_IndexedDataMapOfIntegerIndexedMapOfInteger.hxx>
54 #include <BOPTools_VVInterference.hxx>
55 #include <BOPTools_SSInterference.hxx>
57 #include <BOPTools_Tools2D.hxx>
58 #include <BOPTools_Tools.hxx>
59 #include <NMTTools_ListIteratorOfListOfCoupleOfShape.hxx>
60 #include <NMTTools_IndexedDataMapOfShapeIndexedMapOfShape.hxx>
61 #include <NMTTools_CoupleOfShape.hxx>
62 #include <TopTools_IndexedMapOfShape.hxx>
63 #include <Geom2d_Curve.hxx>
64 #include <Geom_Curve.hxx>
65 #include <Geom_TrimmedCurve.hxx>
66 #include <BOPTools_Tools2D.hxx>
67 #include <BRepLib.hxx>
68 #include <BOPTools_Tools3D.hxx>
69 #include <TopExp_Explorer.hxx>
72 void ProcessBlock(const Standard_Integer iV,
73 const BOPTColStd_IndexedDataMapOfIntegerIndexedMapOfInteger& aMCV,
74 TColStd_IndexedMapOfInteger& aProcessed,
75 TColStd_IndexedMapOfInteger& aChain);
77 void ProcessBlock(const TopoDS_Shape& aF,
78 const NMTTools_IndexedDataMapOfShapeIndexedMapOfShape& aMCV,
79 TopTools_IndexedMapOfShape& aProcessed,
80 TopTools_IndexedMapOfShape& aChain);
82 //=======================================================================
83 // function: MakePCurve
85 //=======================================================================
86 void NMTTools_Tools::MakePCurve(const TopoDS_Edge& aE,
87 const TopoDS_Face& aF,
88 const Handle(Geom2d_Curve)& aC2Dx,
89 const Standard_Real aTolR2D)
91 Standard_Integer k, aNbV;
92 Standard_Real aTolEdge, aTolFact, aTolV, aTolVmax;
93 Standard_Real aTFirst, aTLast, aOutFirst, aOutLast, aOutTol;
95 TopTools_IndexedMapOfShape aVMap;
99 aFFWD.Orientation(TopAbs_FORWARD);
101 aTolEdge=BRep_Tool::Tolerance(aE);
102 aTolFact=Max(aTolEdge, aTolR2D);
104 TopExp::MapShapes(aE, TopAbs_VERTEX, aVMap);
108 for (k=1; k<=aNbV; ++k) {
109 const TopoDS_Vertex& aV=TopoDS::Vertex(aVMap(k));
110 aTolV=BRep_Tool::Tolerance(aV);
111 if (aTolV>aTolVmax) {
116 if (aTolFact>aTolVmax) {
120 const Handle(Geom_Curve)& aC3DE=BRep_Tool::Curve(aE, aTFirst, aTLast);
121 Handle(Geom_TrimmedCurve)aC3DETrim=new Geom_TrimmedCurve(aC3DE, aTFirst, aTLast);
123 Handle(Geom2d_Curve) aC2D, aC2DA;
127 BOPTools_Tools2D::BuildPCurveForEdgeOnFace(aE, aFFWD);
128 BOPTools_Tools2D::CurveOnSurface(aE, aFFWD, aC2D, aOutFirst, aOutLast, aOutTol, Standard_True);
130 if (aC3DE->IsPeriodic()) {
131 BOPTools_Tools2D::AdjustPCurveOnFace(aFFWD, aTFirst, aTLast, aC2D, aC2DA);
134 BOPTools_Tools2D::AdjustPCurveOnFace(aFFWD, aC3DETrim, aC2D, aC2DA);
137 aBB.UpdateEdge(aE, aC2DA, aFFWD, aTolFact);
138 BRepLib::SameParameter(aE);
141 //=======================================================================
142 // function: IsSplitInOnFace
144 //=======================================================================
145 Standard_Boolean NMTTools_Tools::IsSplitInOnFace(const TopoDS_Edge& aE,
146 const TopoDS_Face& aF,
147 IntTools_Context& aContext)
149 Standard_Boolean bFlag;
150 Standard_Real aT, aTolE, aTolF, aTol, aDist, aU, aV;
154 aTolE=BRep_Tool::Tolerance(aE);
155 aTolF=BRep_Tool::Tolerance(aF);
158 GeomAPI_ProjectPointOnSurf& aProjector=aContext.ProjPS(aF);
160 aT=BOPTools_Tools2D::IntermediatePoint(aE);
161 BOPTools_Tools::PointOnEdge(aE, aT, aP);
163 aProjector.Perform(aP);
164 bFlag=aProjector.IsDone();
169 aDist=aProjector.LowerDistance();
170 bFlag=(aDist <= aTol);
175 aProjector.LowerDistanceParameters(aU, aV);
176 aP2D.SetCoord(aU, aV);
177 bFlag=aContext.IsPointInOnFace (aF, aP2D);
180 //=======================================================================
181 // function: NMTTools_Tools::MakeNewVertex
183 //=======================================================================
184 void NMTTools_Tools::MakeNewVertex(const TopTools_ListOfShape& aLVs,
185 TopoDS_Vertex& aNewVertex)
187 Standard_Integer aNb;
188 Standard_Real aTi, aDi, aDmax=-1.e5;
190 gp_XYZ aXYZ(0.,0.,0.), aXYZi;
191 TopTools_ListIteratorOfListOfShape anIt;
198 anIt.Initialize(aLVs);
199 for (; anIt.More(); anIt.Next()) {
200 TopoDS_Vertex aVi=TopoDS::Vertex(anIt.Value());
201 aPi=BRep_Tool::Pnt(aVi);
206 aXYZ.Divide((Standard_Real)aNb);
209 anIt.Initialize(aLVs);
210 for (; anIt.More(); anIt.Next()) {
211 TopoDS_Vertex aVi=TopoDS::Vertex(anIt.Value());
212 aPi=BRep_Tool::Pnt(aVi);
213 aTi=BRep_Tool::Tolerance(aVi);
214 aDi=aP.Distance(aPi);
221 aBB.MakeVertex (aNewVertex, aP, aDmax);
223 //=======================================================================
224 // function: FindChains
226 //=======================================================================
227 void NMTTools_Tools::FindChains(const BOPTools_CArray1OfSSInterference& FFs,
228 BOPTColStd_IndexedDataMapOfIntegerIndexedMapOfInteger& aMapChains)
230 Standard_Boolean bIsTangentFaces;
231 Standard_Integer j, aNb, anIndex1, anIndex2;
232 BOPTColStd_IndexedDataMapOfIntegerIndexedMapOfInteger aMCV;
235 for (j=1; j<=aNb; ++j) {
236 const BOPTools_SSInterference& aFF=FFs(j);
238 bIsTangentFaces=aFF.IsTangentFaces();
239 if (!bIsTangentFaces) {
243 aFF.Indices(anIndex1, anIndex2);
245 if (aMCV.Contains(anIndex1)) {
246 TColStd_IndexedMapOfInteger& aMV=aMCV.ChangeFromKey(anIndex1);
251 TColStd_IndexedMapOfInteger aMV;
254 aMCV.Add(anIndex1, aMV);
257 if (aMCV.Contains(anIndex2)) {
258 TColStd_IndexedMapOfInteger& aMV=aMCV.ChangeFromKey(anIndex2);
263 TColStd_IndexedMapOfInteger aMV;
266 aMCV.Add(anIndex2, aMV);
269 NMTTools_Tools::FindChains(aMCV, aMapChains);
271 //=======================================================================
272 // function: FindChains
274 //=======================================================================
275 void NMTTools_Tools::FindChains(const BOPTools_CArray1OfVVInterference& VVs,
276 BOPTColStd_IndexedDataMapOfIntegerIndexedMapOfInteger& aMapChains)
278 Standard_Integer j, aNb, anIndex1, anIndex2;
279 BOPTColStd_IndexedDataMapOfIntegerIndexedMapOfInteger aMCV;
282 for (j=1; j<=aNb; ++j) {
283 const BOPTools_VVInterference& VV=VVs(j);
284 VV.Indices(anIndex1, anIndex2);
286 if (aMCV.Contains(anIndex1)) {
287 TColStd_IndexedMapOfInteger& aMV=aMCV.ChangeFromKey(anIndex1);
292 TColStd_IndexedMapOfInteger aMV;
295 aMCV.Add(anIndex1, aMV);
298 if (aMCV.Contains(anIndex2)) {
299 TColStd_IndexedMapOfInteger& aMV=aMCV.ChangeFromKey(anIndex2);
304 TColStd_IndexedMapOfInteger aMV;
307 aMCV.Add(anIndex2, aMV);
310 NMTTools_Tools::FindChains(aMCV, aMapChains);
313 //=======================================================================
314 // function: FindChains
316 //=======================================================================
317 void NMTTools_Tools::FindChains(const BOPTColStd_IndexedDataMapOfIntegerIndexedMapOfInteger& aMCV,
318 BOPTColStd_IndexedDataMapOfIntegerIndexedMapOfInteger& aMapChains)
320 Standard_Integer i, j, aNbCV, aNbV, iV, iVx;
321 TColStd_IndexedMapOfInteger aProcessed, aChain;
324 for (i=1; i<=aNbCV; ++i) {
326 if (aProcessed.Contains(iV)) {
333 const TColStd_IndexedMapOfInteger& aMV=aMCV(i);
335 for (j=1; j<=aNbV; ++j) {
337 ProcessBlock(iVx, aMCV, aProcessed, aChain);
339 aMapChains.Add(i, aChain);
343 //=======================================================================
344 // function: ProcessBlock
346 //=======================================================================
347 void ProcessBlock(const Standard_Integer iV,
348 const BOPTColStd_IndexedDataMapOfIntegerIndexedMapOfInteger& aMCV,
349 TColStd_IndexedMapOfInteger& aProcessed,
350 TColStd_IndexedMapOfInteger& aChain)
352 Standard_Integer j, aNbV, iVx;
354 if (aProcessed.Contains(iV)) {
360 const TColStd_IndexedMapOfInteger& aMV=aMCV.FindFromKey(iV);
362 for (j=1; j<=aNbV; ++j) {
364 ProcessBlock(iVx, aMCV, aProcessed, aChain);
367 //=======================================================================
368 // function: AreFacesSameDomain
370 //=======================================================================
371 Standard_Boolean NMTTools_Tools::AreFacesSameDomain(const TopoDS_Face& aF1x,
372 const TopoDS_Face& aF2y,
373 IntTools_Context& aCtx)
375 Standard_Boolean bFlag;
376 Standard_Integer i, aNbE1, aNbE2;
378 TopTools_IndexedMapOfShape aME1, aME2;
379 TopoDS_Edge aEF1, aEF2;
380 TopoDS_Face aF1, aF2;
383 aF1.Orientation(TopAbs_FORWARD);
385 aF2.Orientation(TopAbs_FORWARD);
387 TopExp::MapShapes(aF1, TopAbs_EDGE, aME1);
388 TopExp::MapShapes(aF2, TopAbs_EDGE, aME2);
390 bFlag=Standard_False;
395 if(!aNbE1 || !aNbE2){
403 for (i=1; i<=aNbE1; ++i) {
404 const TopoDS_Edge& aE1=TopoDS::Edge(aME1(i));
405 if (BRep_Tool::Degenerated(aE1)) {
406 // not try to compare deg edges because it
407 // can not have same TShape on different faces at all
410 if (!aME2.Contains(aE1)) {
415 Standard_Real aTolF1, aTolF2, aTol;
418 TopExp_Explorer anExp;
420 aTolF1=BRep_Tool::Tolerance(aF1);
421 aTolF2=BRep_Tool::Tolerance(aF2);
424 anExp.Init(aF1, TopAbs_EDGE);
425 for (; anExp.More(); anExp.Next()) {
426 const TopoDS_Edge& aE1=TopoDS::Edge(anExp.Current());
427 if (!BRep_Tool::Degenerated(aE1)) {
428 BOPTools_Tools3D::PointNearEdge(aE1, aF1, aP2D, aP);
429 bFlag=aCtx.IsValidPointForFace(aP, aF2, aTol);
436 Standard_Real aU1, aU2, aV1, aV2;
437 Standard_Real dU, dV, aU, aV;
438 Standard_Integer aNbP=5, aNbP1, j;
442 aTolF1=BRep_Tool::Tolerance(aF1);
443 aTolF2=BRep_Tool::Tolerance(aF2);
446 BRepTools::UVBounds(aF1, aU1, aU2, aV1, aV2);
447 Handle(Geom_Surface) aS1=BRep_Tool::Surface(aF1);
453 for (i=1; i<=aNbP; ++i) {
455 for (j=1; j<=aNbP; ++j) {
457 aP2D.SetCoord(aU, aV);
459 if(aCtx.IsPointInFace(aF1, aP2D)) {
460 aP=aS1->Value(aU, aV);
461 bFlag=aCtx.IsValidPointForFace(aP, aF2, aTol);
472 //=======================================================================
473 // function: FindChains
475 //=======================================================================
476 void NMTTools_Tools::FindChains(const NMTTools_ListOfCoupleOfShape& aLCS,
477 NMTTools_IndexedDataMapOfShapeIndexedMapOfShape& aMapChains)
479 NMTTools_ListIteratorOfListOfCoupleOfShape aItCS;
480 NMTTools_IndexedDataMapOfShapeIndexedMapOfShape aMCV;
482 aItCS.Initialize(aLCS);
483 for (; aItCS.More(); aItCS.Next()) {
484 const NMTTools_CoupleOfShape& aCS=aItCS.Value();
486 const TopoDS_Shape& aF1=aCS.Shape1();
487 const TopoDS_Shape& aF2=aCS.Shape2();
490 if (aMCV.Contains(aF1)) {
491 TopTools_IndexedMapOfShape& aMV=aMCV.ChangeFromKey(aF1);
496 TopTools_IndexedMapOfShape aMV;
502 if (aMCV.Contains(aF2)) {
503 TopTools_IndexedMapOfShape& aMV=aMCV.ChangeFromKey(aF2);
508 TopTools_IndexedMapOfShape aMV;
514 NMTTools_Tools::FindChains(aMCV, aMapChains);
516 //=======================================================================
517 // function: FindChains
519 //=======================================================================
520 void NMTTools_Tools::FindChains(const NMTTools_IndexedDataMapOfShapeIndexedMapOfShape& aMCV,
521 NMTTools_IndexedDataMapOfShapeIndexedMapOfShape& aMapChains)
523 Standard_Integer i, j, aNbCV, aNbV;
524 TopTools_IndexedMapOfShape aProcessed, aChain;
527 for (i=1; i<=aNbCV; ++i) {
528 const TopoDS_Shape& aF=aMCV.FindKey(i);
529 if (aProcessed.Contains(aF)) {
536 const TopTools_IndexedMapOfShape& aMV=aMCV(i);
538 for (j=1; j<=aNbV; ++j) {
539 const TopoDS_Shape& aFx=aMV(j);
540 ProcessBlock(aFx, aMCV, aProcessed, aChain);
542 aMapChains.Add(aF, aChain);
546 //=======================================================================
547 // function: ProcessBlock
549 //=======================================================================
550 void ProcessBlock(const TopoDS_Shape& aF,
551 const NMTTools_IndexedDataMapOfShapeIndexedMapOfShape& aMCV,
552 TopTools_IndexedMapOfShape& aProcessed,
553 TopTools_IndexedMapOfShape& aChain)
555 Standard_Integer j, aNbV;
557 if (aProcessed.Contains(aF)) {
563 const TopTools_IndexedMapOfShape& aMV=aMCV.FindFromKey(aF);
565 for (j=1; j<=aNbV; ++j) {
566 const TopoDS_Shape& aFx=aMV(j);
567 ProcessBlock(aFx, aMCV, aProcessed, aChain);