src/PARTITION/Partition_Inter3d.cxx

   1 //  GEOM PARTITION : partition algorithm
   2 //
   3 //  Copyright (C) 2003  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
   4 //  CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
   5 //
   6 //  This library is free software; you can redistribute it and/or
   7 //  modify it under the terms of the GNU Lesser General Public
   8 //  License as published by the Free Software Foundation; either
   9 //  version 2.1 of the License.
  10 //
  11 //  This library is distributed in the hope that it will be useful,
  12 //  but WITHOUT ANY WARRANTY; without even the implied warranty of
  13 //  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  14 //  Lesser General Public License for more details.
  15 //
  16 //  You should have received a copy of the GNU Lesser General Public
  17 //  License along with this library; if not, write to the Free Software
  18 //  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  19 //
  20 //  See http://www.opencascade.org/SALOME/ or email : webmaster.salome@opencascade.org
  21 //
  22 //
  23 //
  24 //  File   : Partition_Inter3d.cxx
  25 //  Author : Benedicte MARTIN
  26 //  Module : GEOM
  27 //  $Header$
  28
  29 using namespace std;
  30 #include "Partition_Inter3d.ixx"
  31 #include "Partition_Inter2d.hxx"
  32 #include "utilities.h"
  33
  34 #include <BRepOffset_Tool.hxx>
  35 #include <BRep_Builder.hxx>
  36 #include <BRep_Tool.hxx>
  37 #include <BRepAlgo_AsDes.hxx>
  38 #include <BRepAlgo_Image.hxx>
  39 #include <BRepLib.hxx>
  40
  41 #include <TopExp.hxx>
  42 #include <TopExp_Explorer.hxx>
  43
  44 #include <TopoDS.hxx>
  45 #include <TopoDS_Vertex.hxx>
  46 #include <TopoDS_Edge.hxx>
  47 #include <TopoDS_Face.hxx>
  48 #include <TopoDS_Compound.hxx>
  49 #include <TopTools_ListOfShape.hxx>
  50 #include <TopTools_ListIteratorOfListOfShape.hxx>
  51 #include <TopTools_DataMapIteratorOfDataMapOfShapeListOfShape.hxx>
  52 #include <TopOpeBRepTool_BoxSort.hxx>
  53
  54 #ifdef DEB
  55 #include <DBRep.hxx>
  56 #endif
  57
  58 #include <stdio.h>
  59 #include <TopOpeBRepDS_HDataStructure.hxx>
  60 #include <TopOpeBRep_DSFiller.hxx>
  61 #include <TopOpeBRepTool_GeomTool.hxx>
  62 #include <TopOpeBRepTool_OutCurveType.hxx>
  63 #include <TopOpeBRepDS_BuildTool.hxx>
  64 #include <TopOpeBRepBuild_Builder.hxx>
  65 #include <TopOpeBRepDS_CurveExplorer.hxx>
  66 #include <Geom2d_Curve.hxx>
  67 #include <TopOpeBRepDS_PointIterator.hxx>
  68 #include <TopOpeBRepDS_Transition.hxx>
  69 #include <Geom_Curve.hxx>
  70 #include <TopOpeBRepTool_CurveTool.hxx>
  71 #include <TopOpeBRepDS_Interference.hxx>
  72 #include <TopTools_DataMapIteratorOfDataMapOfShapeShape.hxx>
  73 #include <BRepLib_MakeVertex.hxx>
  74 #include <Precision.hxx>
  75 #include <TColStd_MapOfInteger.hxx>
  76 #include <BRepTools.hxx>
  77 #include <Geom_RectangularTrimmedSurface.hxx>
  78 #include <Geom_Surface.hxx>
  79 #include <Geom_TrimmedCurve.hxx>
  80 #include <Geom2dAPI_ProjectPointOnCurve.hxx>
  81 #include <GeomAPI_ProjectPointOnCurve.hxx>
  82
  83 //=======================================================================
  84 //function : Partition_Inter3d
  85 //purpose  :
  86 //=======================================================================
  87
  88 Partition_Inter3d::Partition_Inter3d()
  89 {
  90 }
  91 //=======================================================================
  92 //function : Partition_Inter3d
  93 //purpose  :
  94 //=======================================================================
  95
  96 Partition_Inter3d::Partition_Inter3d(const Handle(BRepAlgo_AsDes)& AsDes)
  97   :myAsDes(AsDes)
  98 {
  99   mySectionEdgesAD = new BRepAlgo_AsDes;
 100 }
 101
 102 //=======================================================================
 103 //function : CompletPart3d
 104 //purpose  : FaceShapeMap is just to know the shape a face belongs to
 105 //=======================================================================
 106
 107 void Partition_Inter3d::CompletPart3d(const TopTools_ListOfShape& SetOfFaces1,
 108                                       const TopTools_DataMapOfShapeShape& FaceShapeMap)
 109 {
 110   if (myAsDes.IsNull())
 111     myAsDes = new BRepAlgo_AsDes;
 112
 113   TopTools_ListIteratorOfListOfShape it;
 114
 115   //---------------------------------------------------------------
 116   // Construction of bounding boxes.
 117   //---------------------------------------------------------------
 118
 119   BRep_Builder B;
 120   TopoDS_Compound CompOS;
 121   B.MakeCompound(CompOS);
 122   for (it.Initialize(SetOfFaces1); it.More(); it.Next())
 123     B.Add(CompOS, it.Value());
 124
 125   TopOpeBRepTool_BoxSort BOS;
 126   BOS.AddBoxesMakeCOB(CompOS,TopAbs_FACE);
 127
 128   for (it.Initialize(SetOfFaces1); it.More(); it.Next()) {
 129     TopoDS_Face F1 = TopoDS::Face(it.Value());
 130
 131     // avoid intersecting faces of one shape
 132     TopoDS_Shape S1;
 133     if (FaceShapeMap.IsBound(F1)) S1 = FaceShapeMap.Find(F1);
 134     // avoid intersecting faces sharing vertices, suppose they belong to
 135     // shapes sharing same faces
 136     TopTools_IndexedMapOfShape VM;
 137     TopExp::MapShapes( F1, TopAbs_VERTEX, VM);
 138
 139     TColStd_ListIteratorOfListOfInteger itLI = BOS.Compare(F1);
 140     for (; itLI.More(); itLI.Next()) {
 141       TopoDS_Face F2 = TopoDS::Face(BOS.TouchedShape(itLI));
 142       if (F1.IsSame(F2) || IsDone(F1,F2))
 143         continue;
 144
 145       TopoDS_Shape S2;
 146       if (FaceShapeMap.IsBound(F2)) S2 = FaceShapeMap.Find(F2);
 147       if (!S1.IsNull() && S1.IsSame(S2))
 148         continue; // descendants of one shape
 149
 150       TopExp_Explorer expV (F2, TopAbs_VERTEX);
 151       for ( ; expV.More(); expV.Next())
 152         if (VM.Contains( expV.Current() ))
 153           break;
 154       if (expV.More())
 155         continue; // faces have a common edge
 156
 157       F1.Orientation(TopAbs_FORWARD);
 158       F2.Orientation(TopAbs_FORWARD);
 159       FacesPartition(F1,F2);
 160     }
 161
 162     // mark as modified a face which has at least one new edge
 163     if (!myAsDes->HasDescendant( F1 ))
 164       continue;
 165     TopTools_ListIteratorOfListOfShape itE (myAsDes->Descendant( F1 ));
 166     for ( ; itE.More(); itE.Next()) {
 167       if (myNewEdges.Contains( itE.Value())) {
 168         myTouched.Add( F1 );
 169         break;
 170       }
 171     }
 172   }
 173 }
 174
 175 //=======================================================================
 176 //function : PutInBounds
 177 //purpose  :
 178 //=======================================================================
 179
 180 static void PutInBounds (const TopoDS_Face&          F,
 181                          const TopoDS_Edge&          E,
 182                          Handle(Geom2d_Curve)&       C2d)
 183 {
 184   Standard_Real   umin,umax,vmin,vmax;
 185   Standard_Real   f,l;
 186   BRep_Tool::Range(E,f,l);
 187
 188   TopLoc_Location L; // Recup S avec la location pour eviter la copie.
 189   Handle (Geom_Surface) S   = BRep_Tool::Surface(F,L);
 190
 191   if (S->IsKind(STANDARD_TYPE(Geom_RectangularTrimmedSurface))) {
 192     S = (*(Handle_Geom_RectangularTrimmedSurface*)&S)->BasisSurface();
 193   }
 194   //---------------
 195   // Recadre en U.
 196   //---------------
 197   if (!S->IsUPeriodic() && !S->IsVPeriodic()) return;
 198
 199   BRepTools::UVBounds(F,umin,umax,vmin,vmax);
 200
 201   if (S->IsUPeriodic()) {
 202     Standard_Real period  = S->UPeriod();
 203     Standard_Real eps     = period*1.e-6;
 204     gp_Pnt2d      Pf      = C2d->Value(f);
 205     gp_Pnt2d      Pl      = C2d->Value(l);
 206     gp_Pnt2d      Pm      = C2d->Value(0.34*f + 0.66*l);
 207     Standard_Real minC    = Min(Pf.X(),Pl.X()); minC = Min(minC,Pm.X());
 208     Standard_Real maxC    = Max(Pf.X(),Pl.X()); maxC = Max(maxC,Pm.X());
 209     Standard_Real du = 0.;
 210     if (minC< umin - eps) {
 211       du = (int((umin - minC)/period) + 1)*period;
 212     }
 213     if (minC > umax + eps) {
 214       du = -(int((minC - umax)/period) + 1)*period;
 215     }
 216     if (du != 0) {
 217       gp_Vec2d T1(du,0.);
 218       C2d->Translate(T1);
 219       minC += du; maxC += du;
 220     }
 221     // Ajuste au mieux la courbe dans le domaine.
 222     if (maxC > umax +100*eps) {
 223       Standard_Real d1 = maxC - umax;
 224       Standard_Real d2 = umin - minC + period;
 225       if (d2 < d1) du =-period;
 226       if ( du != 0.) {
 227         gp_Vec2d T2(du,0.);
 228         C2d->Translate(T2);
 229       }
 230     }
 231   }
 232   //------------------
 233   // Recadre en V.
 234   //------------------
 235   if (S->IsVPeriodic()) {
 236     Standard_Real period  = S->VPeriod();
 237     Standard_Real eps     = period*1.e-6;
 238     gp_Pnt2d      Pf      = C2d->Value(f);
 239     gp_Pnt2d      Pl      = C2d->Value(l);
 240     gp_Pnt2d      Pm      = C2d->Value(0.34*f + 0.66*l);
 241     Standard_Real minC    = Min(Pf.Y(),Pl.Y()); minC = Min(minC,Pm.Y());
 242     Standard_Real maxC    = Max(Pf.Y(),Pl.Y()); maxC = Max(maxC,Pm.Y());
 243     Standard_Real dv = 0.;
 244     if (minC< vmin - eps) {
 245       dv = (int((vmin - minC)/period) + 1)*period;
 246     }
 247     if (minC > vmax + eps) {
 248       dv = -(int((minC - vmax)/period) + 1)*period;
 249     }
 250     if (dv != 0) {
 251       gp_Vec2d T1(0.,dv);
 252       C2d->Translate(T1);
 253       minC += dv; maxC += dv;
 254     }
 255     // Ajuste au mieux la courbe dans le domaine.
 256     if (maxC > vmax +100*eps) {
 257       Standard_Real d1 = maxC - vmax;
 258       Standard_Real d2 = vmin - minC + period;
 259       if (d2 < d1) dv =-period;
 260       if ( dv != 0.) {
 261         gp_Vec2d T2(0.,dv);
 262         C2d->Translate(T2);
 263       }
 264     }
 265   }
 266 }
 267
 268 //=======================================================================
 269 //function : Inter3D
 270 //purpose  :
 271 //=======================================================================
 272
 273 void Partition_Inter3d::Inter3D(const TopoDS_Face& F1,
 274                                 const TopoDS_Face& F2,
 275                                 TopTools_ListOfShape& L)
 276 {
 277   BRep_Builder B;
 278
 279   // fill the data Structure
 280   Handle(TopOpeBRepDS_HDataStructure) DatStr = new TopOpeBRepDS_HDataStructure();
 281   TopOpeBRep_DSFiller DSFiller;
 282   DSFiller.Insert(F1,F2,DatStr);
 283
 284   // define the GeomTool used by the DSFiller :
 285   // compute BSpline of degree 1 on intersection curves.
 286   Standard_Real tol3dAPPROX = 1e-7;
 287   Standard_Real tol2dAPPROX = 1e-7;
 288   TopOpeBRepTool_GeomTool GT2 (TopOpeBRepTool_APPROX);
 289   GT2.SetTolerances(tol3dAPPROX,tol2dAPPROX);
 290   TopOpeBRepDS_BuildTool  BT(GT2);
 291
 292   // Perform Section
 293   TopOpeBRepBuild_Builder TopB(BT);
 294   TopB.Perform(DatStr);
 295
 296   // ===============
 297   // Store new edges
 298   // ===============
 299
 300   L.Clear();
 301   TopOpeBRepDS_CurveExplorer cex(DatStr->DS());
 302   for (; cex.More(); cex.Next()) {
 303     const TopOpeBRepDS_Curve& CDS = cex.Curve();
 304     Standard_Integer ic = cex.Index();
 305     Handle(Geom2d_Curve) pc1 = CDS.Curve1();
 306     Handle(Geom2d_Curve) pc2 = CDS.Curve2();
 307
 308     TopTools_ListIteratorOfListOfShape itLE = TopB.NewEdges(ic);
 309     while (itLE.More()) {
 310       TopoDS_Edge E = TopoDS::Edge(itLE.Value());
 311
 312 //       Standard_Real f,l;
 313 //       BRep_Tool::Range(E,f,l);
 314       PutInBounds (F1,E,pc1);
 315       PutInBounds (F2,E,pc2);
 316
 317       B.UpdateEdge (E,pc1,F1,0.);
 318       B.UpdateEdge (E,pc2,F2,0.);
 319
 320       L.Append (E);
 321
 322       itLE.Next();
 323       if (itLE.More()) {
 324         pc1 = Handle(Geom2d_Curve)::DownCast(pc1->Copy());
 325         pc2 = Handle(Geom2d_Curve)::DownCast(pc2->Copy());
 326       }
 327     }
 328   }
 329
 330   // ===================================================
 331   // Store section edges, same domain faces and verives
 332   // ===================================================
 333
 334   TopTools_ListOfShape empty, LSP, LSE;
 335
 336   if ( DatStr->HasSameDomain( F1 )) { // same domain faces
 337     if (!mySameDomainFM.IsBound(F1))
 338       mySameDomainFM.Bind(F1,empty);
 339     if (!mySameDomainFM.IsBound(F2))
 340       mySameDomainFM.Bind(F2,empty);
 341     mySameDomainFM(F1).Append(F2);
 342     mySameDomainFM(F2).Append(F1);
 343   }
 344
 345   const TopOpeBRepDS_DataStructure& DS = DatStr->DS();
 346   Standard_Integer j,i,nes = DS.NbSectionEdges();
 347   if (!nes) return;
 348
 349
 350   TopoDS_Vertex V, sdeV1, sdeV2;
 351   TopTools_MapOfShape MV;
 352
 353   // put vertices on section edges
 354   for (i=1;i<=nes;i++) {
 355
 356     TopoDS_Edge se, sde, oe; // section, same domain, other edge
 357     se = DS.SectionEdge(i);
 358     if (! TopB.IsSplit(se,TopAbs_ON))
 359       continue;
 360
 361     if (DatStr->HasSameDomain(se)) {
 362       sde = TopoDS::Edge( DatStr->SameDomain(se).Value() );
 363       TopExp::Vertices( sde, sdeV1, sdeV2);
 364     }
 365
 366     TColStd_MapOfInteger MIV;
 367     TopOpeBRepDS_PointIterator itP (DS.ShapeInterferences( se ));
 368     itP.SupportKind( TopOpeBRepDS_EDGE );
 369     for (; itP.More(); itP.Next()) {
 370       oe = TopoDS::Edge( DS.Shape( itP.Support()));
 371       if (itP.IsVertex()) {
 372         if ( !MIV.Add( itP.Current() ))
 373           continue;
 374         V = TopoDS::Vertex( DS.Shape( itP.Current()));
 375         if ( !sde.IsNull() && (V.IsSame(sdeV1) || V.IsSame(sdeV2)) )
 376           oe = sde;
 377         V = ReplaceSameDomainV( V , oe );
 378         V.Orientation( TopAbs_INTERNAL);
 379         B.UpdateVertex( V, itP.Parameter(), se, 0.);
 380       }
 381       else {
 382         const TopOpeBRepDS_Point& DSP = DS.Point( itP.Current());
 383         V = BRepLib_MakeVertex( DSP.Point() );
 384         V.Orientation( TopAbs_INTERNAL);
 385         B.UpdateVertex( V, itP.Parameter(), se, DSP.Tolerance());
 386         // make V be on the other edge
 387         TopOpeBRepDS_PointIterator itOP (DS.ShapeInterferences( oe ));
 388         for (; itOP.More(); itOP.Next()) {
 389           const TopOpeBRepDS_Point& ODSP = DS.Point( itOP.Current());
 390           if ( DSP.IsEqual (ODSP)) {
 391             B.UpdateVertex( V, itOP.Parameter(), TopoDS::Edge(oe), ODSP.Tolerance());
 392             break;
 393           }
 394         }
 395       }
 396       TopoDS_Vertex addedV = Partition_Inter2d::AddVonE( V,se,oe,myAsDes);
 397       if (!addedV.IsSame( V ))
 398         mySameDomainVM.Bind (V, addedV);
 399       MV.Add( addedV );
 400     }
 401   }
 402
 403   TopB.SplitSectionEdges();
 404
 405   TopTools_DataMapOfShapeShape SEM; // map split - section edge
 406   TopTools_IndexedMapOfShape ME[2];
 407   TopExp::MapShapes( F1, TopAbs_EDGE, ME[1]);
 408   TopExp::MapShapes( F2, TopAbs_EDGE, ME[0]);
 409
 410   // add section edge to the face it intersects and find
 411   // splits ON that do not have same domain pair
 412
 413   for (i=1;i<=nes;i++) {
 414
 415     const TopoDS_Edge& se = DS.SectionEdge(i);
 416     if (! TopB.IsSplit(se,TopAbs_ON))
 417       continue;
 418
 419     Standard_Integer ancRank = DS.AncestorRank(se);
 420     if (ME[ancRank-1].Contains( se ))
 421       continue; // se is an edge of face it intersects
 422
 423     const TopoDS_Face& F = (ancRank == 1) ? F2 : F1;
 424
 425     // add se to face but dont add twice
 426     TopTools_ListIteratorOfListOfShape itE;
 427     if (myAsDes->HasDescendant( F )) {
 428       for (itE.Initialize( (myAsDes->Descendant( F )) ); itE.More(); itE.Next())
 429         if (se.IsSame( itE.Value() ))
 430           break;
 431     }
 432     if (!itE.More()) {
 433       myAsDes->Add( F, se );
 434       Standard_Real tol, f,l, umin=1e100, umax=-1e100;
 435       Handle(Geom2d_Curve) pc = BRep_Tool::CurveOnSurface( se, F, f,l);
 436       if (pc.IsNull()) {
 437         TopTools_ListIteratorOfListOfShape it( TopB.Splits(se,TopAbs_ON) );
 438         for ( ;it.More();it.Next()) {
 439           const TopoDS_Edge& E = TopoDS::Edge ( it.Value());
 440           BRep_Tool::Range(E, f, l);
 441           umin = Min( umin, f);
 442           umax = Max( umax, l);
 443         }
 444         Handle(Geom_Curve) C3d = BRep_Tool::Curve( se, f, l);
 445         if (umin < umax) // sometimes umin == umax for closed edge
 446           C3d = new Geom_TrimmedCurve( C3d, umin, umax);
 447         pc = TopOpeBRepTool_CurveTool::MakePCurveOnFace (F,C3d,tol);
 448         if (pc.IsNull()) {
 449           MESSAGE (" CANT BUILD PCURVE ");
 450         }
 451         B.UpdateEdge( se, pc, F, tol);
 452       }
 453     }
 454
 455     // to detect splits that do not have same domain pair
 456     TopTools_ListIteratorOfListOfShape it( TopB.Splits(se,TopAbs_ON) );
 457     for ( ;it.More();it.Next()) {
 458       const TopoDS_Edge& S = TopoDS::Edge ( it.Value());
 459       if (SEM.IsBound( S ))
 460         SEM.UnBind( S );
 461       else
 462         SEM.Bind( S, se);
 463     }
 464   }
 465
 466   // store vertices of ON splits and bind section edges to faces
 467   for (i=1;i<=nes;i++) {
 468
 469     const TopoDS_Edge& se = DS.SectionEdge(i);
 470     if (! TopB.IsSplit(se,TopAbs_ON))
 471       continue;
 472
 473     Standard_Integer ancRank = DS.AncestorRank(se);
 474     if (ME[ancRank-1].Contains( se ))
 475       continue; // se is an edge of face it intersects
 476
 477     TopoDS_Face F = (ancRank == 1) ? F2 : F1;
 478
 479     // add vertices of splits
 480     Standard_Boolean added = Standard_False;
 481     TopTools_ListIteratorOfListOfShape it( TopB.Splits(se,TopAbs_ON) );
 482     for ( ;it.More();it.Next()) {
 483       const TopoDS_Edge& S = TopoDS::Edge ( it.Value());
 484       if (!SEM.IsBound( S ))
 485         continue;
 486
 487       added = Standard_True;
 488       mySectionEdgesAD->Add( F, se );
 489
 490       TopoDS_Vertex VS[2];
 491       TopExp::Vertices (S, VS[0], VS[1]);
 492       for (j=0; j<2; ++j) {
 493         if (mySameDomainVM.IsBound( VS[j] ))
 494           VS[j] = TopoDS::Vertex( mySameDomainVM( VS[j] ));
 495         if ( !MV.Contains( VS[j] )) {
 496           // find equal vertex on se - point interference
 497           gp_Pnt P1 = BRep_Tool::Pnt( VS[j] );
 498           TopTools_ListIteratorOfListOfShape itV( myAsDes->Descendant(se) );
 499           for (; itV.More(); itV.Next()) {
 500             V = TopoDS::Vertex( itV.Value() );
 501             gp_Pnt P2 = BRep_Tool::Pnt( V );
 502             if (P1.IsEqual( P2, Precision::Confusion())) {
 503               mySameDomainVM.Bind (VS[j], V);
 504               VS[j] = V;
 505               break;
 506             }
 507           }
 508           if (!itV.More())  // no interferences with edges
 509             myAsDes->Add( se, VS[j]);
 510         }
 511         mySectionEdgesAD->Add( F, VS[j]);
 512       }
 513       mySectionEdgesAD->Add( F, S );
 514     }
 515     if (!added)
 516       mySectionEdgesAD->Add( F, se );
 517
 518     myNewEdges.Add( se );
 519   }
 520 }
 521
 522 //=======================================================================
 523 //function : FacesPartition
 524 //purpose  :
 525 //=======================================================================
 526
 527 void Partition_Inter3d::FacesPartition(const TopoDS_Face& F1,
 528                                        const TopoDS_Face& F2)
 529      //(const TopTools_DataMapOfShapeListOfShape& /*SetOfFaces2*/)
 530 {
 531   TopTools_ListOfShape LInt;
 532
 533   Inter3D (F1,F2,LInt);
 534
 535   StorePart3d (F1,F2,LInt);
 536 }
 537
 538 //=======================================================================
 539 //function : SetDone
 540 //purpose  :
 541 //=======================================================================
 542
 543 void Partition_Inter3d::SetDone(const TopoDS_Face& F1,
 544                                 const TopoDS_Face& F2)
 545 {
 546   if (!myDone.IsBound(F1)) {
 547     TopTools_ListOfShape empty;
 548     myDone.Bind(F1,empty);
 549   }
 550   myDone(F1).Append(F2);
 551   if (!myDone.IsBound(F2)) {
 552     TopTools_ListOfShape empty;
 553     myDone.Bind(F2,empty);
 554   }
 555   myDone(F2).Append(F1);
 556 }
 557
 558 //=======================================================================
 559 //function : IsDone
 560 //purpose  :
 561 //=======================================================================
 562
 563 Standard_Boolean Partition_Inter3d::IsDone(const TopoDS_Face& F1,
 564                                            const TopoDS_Face& F2)
 565
 566   const
 567 {
 568   if (myDone.IsBound(F1)) {
 569     TopTools_ListIteratorOfListOfShape it (myDone(F1));
 570     for (; it.More(); it.Next()) {
 571       if (it.Value().IsSame(F2)) return Standard_True;
 572     }
 573   }
 574   return Standard_False;
 575 }
 576
 577 //=======================================================================
 578 //function : StorePart3d
 579 //purpose  :
 580 //=======================================================================
 581
 582 void Partition_Inter3d::StorePart3d(const TopoDS_Face& F1,
 583                                     const TopoDS_Face& F2,
 584                                     const TopTools_ListOfShape& LInt)
 585 {
 586
 587   if (!LInt.IsEmpty()) {
 588     myAsDes->Add( F1,LInt);
 589     myAsDes->Add( F2,LInt);
 590
 591     TopTools_ListIteratorOfListOfShape it(LInt);
 592     for (; it.More(); it.Next()) {
 593
 594       TopoDS_Edge E = TopoDS::Edge(it.Value());
 595
 596       BRep_Builder B;
 597       B.SameParameter(E,Standard_False);
 598       BRepLib::SameParameter(E,1.0e-7);
 599
 600       myNewEdges.Add(E);
 601     }
 602   }
 603   SetDone(F1,F2);
 604 }
 605
 606 //=======================================================================
 607 //function : TouchedFaces
 608 //purpose  :
 609 //=======================================================================
 610
 611 TopTools_MapOfShape& Partition_Inter3d::TouchedFaces()
 612 {
 613   return myTouched;
 614 }
 615
 616 //=======================================================================
 617 //function : AsDes
 618 //purpose  :
 619 //=======================================================================
 620
 621 Handle(BRepAlgo_AsDes) Partition_Inter3d::AsDes() const
 622 {
 623   return myAsDes;
 624 }
 625
 626 //=======================================================================
 627 //function : NewEdges
 628 //purpose  :
 629 //=======================================================================
 630
 631 TopTools_MapOfShape& Partition_Inter3d::NewEdges()
 632 {
 633   return myNewEdges;
 634 }
 635
 636 //=======================================================================
 637 //function : Affiche
 638 //purpose  :
 639 //=======================================================================
 640
 641 void Partition_Inter3d::Affiche(const TopTools_ListOfShape& SetOfFaces) const
 642 {
 643 #ifdef DEB
 644   char PSection[1024];
 645   char *section=PSection;
 646   Standard_Integer i = 0;
 647   Standard_Real j=1;
 648   TopTools_ListOfShape aList;
 649   TopTools_ListIteratorOfListOfShape it;
 650   for (it.Initialize(SetOfFaces); it.More(); it.Next()) {
 651     const TopoDS_Shape& OS = it.Value();
 652     aList=myAsDes->Descendant(OS);
 653     MESSAGE ( " the number of items stored in the list " << j << " :  " << aList.Extent() )
 654     j++;
 655     TopTools_ListIteratorOfListOfShape itaList;
 656     for (itaList.Initialize(aList); itaList.More(); itaList.Next()) {
 657       const TopoDS_Shape& SS = itaList.Value();
 658       i++;
 659       sprintf(PSection,"section_%d",i);
 660       DBRep::Set(section,SS);
 661     }
 662   }
 663 #endif
 664 }
 665
 666 //=======================================================================
 667 //function : SameDomain
 668 //purpose  :
 669 //=======================================================================
 670
 671 const TopTools_ListOfShape& Partition_Inter3d::SameDomain(const TopoDS_Face& F) const
 672 {
 673   if (mySameDomainFM.IsBound( F ))
 674     return mySameDomainFM (F);
 675
 676   static TopTools_ListOfShape empty;
 677   return empty;
 678 }
 679
 680 //=======================================================================
 681 //function : HasSameDomainF
 682 //purpose  : Return true if F has same domain faces
 683 //=======================================================================
 684
 685 Standard_Boolean Partition_Inter3d::HasSameDomainF(const TopoDS_Shape& F) const
 686 {
 687   return mySameDomainFM.IsBound( F );
 688 }
 689
 690 //=======================================================================
 691 //function : IsSameDomain
 692 //purpose  : Return true if F1 and F2 are same domain faces
 693 //=======================================================================
 694
 695 Standard_Boolean Partition_Inter3d::IsSameDomainF(const TopoDS_Shape& F1,
 696                                                  const TopoDS_Shape& F2) const
 697 {
 698   if (mySameDomainFM.IsBound( F1 )) {
 699     TopTools_ListIteratorOfListOfShape it (mySameDomainFM( F1 ));
 700     for (; it.More(); it.Next())
 701       if (F2.IsSame( it.Value()))
 702         return Standard_True;
 703   }
 704   return F1.IsSame( F2 );
 705 }
 706
 707 //=======================================================================
 708 //function : ReplaceSameDomainV
 709 //purpose  : return same domain vertex of  V if it was replaced
 710 //           and make this vertex to be on E too, else return V
 711 //=======================================================================
 712
 713 TopoDS_Vertex Partition_Inter3d::ReplaceSameDomainV(const TopoDS_Vertex& V,
 714                                                     const TopoDS_Edge&   E) const
 715 {
 716   TopoDS_Vertex SDV = V;
 717   if (mySameDomainVM.IsBound( V )) {
 718
 719     TopoDS_Vertex V1,V2;
 720     TopExp::Vertices(E,V1,V2);
 721     Standard_Boolean isClosed = V1.IsSame( V2 ) && V.IsSame(V1);
 722
 723     SDV = TopoDS::Vertex( mySameDomainVM(V) );
 724     Standard_Real tol = BRep_Tool::Tolerance( V );
 725     BRep_Builder B;
 726     SDV.Orientation( V.Orientation());
 727
 728     if (isClosed) {
 729       Standard_Real f, l;
 730       BRep_Tool::Range (E, f, l);
 731       Standard_Boolean isFirst = IsEqual( BRep_Tool::Parameter(V,E), f );
 732       B.UpdateVertex(SDV, (isFirst ? f : l), E, tol);
 733       SDV.Reverse();
 734       B.UpdateVertex(SDV, (isFirst ? l : f), E, tol);
 735     }
 736     else
 737       B.UpdateVertex (SDV, BRep_Tool::Parameter(V,E), E, tol);
 738
 739   }
 740   return SDV;
 741 }
 742
 743 //=======================================================================
 744 //function : SectionEdgesAD
 745 //purpose  :
 746 //=======================================================================
 747
 748 Handle(BRepAlgo_AsDes) Partition_Inter3d::SectionEdgesAD() const
 749 {
 750   return mySectionEdgesAD;
 751 }
 752
 753 //=======================================================================
 754 //function : IsSectionEdge
 755 //purpose  : return True if  E  is  an  edge  of  a face and it
 756 //           intersects an other face
 757 //=======================================================================
 758
 759 Standard_Boolean
 760   Partition_Inter3d::IsSectionEdge(const TopoDS_Edge& E) const
 761 {
 762   return mySectionEdgesAD->HasAscendant(E);
 763 }
 764
 765 //=======================================================================
 766 //function : HasSectionEdge
 767 //purpose  : return True if an  edge  of  F intersects an other
 768 //           face or F is intersected by edge of an other face
 769 //=======================================================================
 770
 771 Standard_Boolean
 772   Partition_Inter3d::HasSectionEdge(const TopoDS_Face& F) const
 773 {
 774   return mySectionEdgesAD->HasDescendant(F);
 775 }
 776
 777 //=======================================================================
 778 //function : IsSplitOn
 779 //purpose  : return True if NewE is split of OldE on F
 780 //=======================================================================
 781
 782 Standard_Boolean
 783   Partition_Inter3d::IsSplitOn(const TopoDS_Edge& NewE,
 784                                const TopoDS_Edge& OldE,
 785                                const TopoDS_Face& F) const
 786 {
 787   if (! mySectionEdgesAD->HasDescendant(F))
 788     return Standard_False;
 789
 790   TopTools_ListIteratorOfListOfShape itE ( mySectionEdgesAD->Descendant(F) );
 791   for ( ; itE.More(); itE.Next()) {
 792     if ( itE.Value().ShapeType() != TopAbs_EDGE ||
 793         ! OldE.IsSame ( itE.Value() ))
 794       continue;
 795     // an edge encountered, its vertices and a split come next
 796     itE.Next();
 797     if (!itE.More()) break;
 798     const TopoDS_Shape& V3 = itE.Value();
 799     if (V3.ShapeType() != TopAbs_VERTEX) continue;
 800     itE.Next();
 801     if (!itE.More()) break;
 802     const TopoDS_Shape& V4 = itE.Value();
 803     if (V4.ShapeType() != TopAbs_VERTEX) continue;
 804
 805     TopoDS_Vertex V1, V2;
 806     TopExp::Vertices( OldE, V1, V2);
 807
 808     if ( V1.IsSame(V2) &&
 809         (V1.IsSame(V3) || V1.IsSame(V4)) ) {
 810       // closed old edge; use the split for the test
 811       itE.Next();
 812       if (!itE.More()) break;
 813       const TopoDS_Edge& split = TopoDS::Edge( itE.Value() );
 814       // check distance at middle point of NewE
 815       Standard_Real f1,l1, f2,l2;
 816       Handle(Geom2d_Curve) PC1 = BRep_Tool::CurveOnSurface( split, F ,f1,l1);
 817       if (!PC1.IsNull()) {
 818         Handle(Geom2d_Curve) PC2 = BRep_Tool::CurveOnSurface(NewE, F ,f2,l2);
 819         gp_Pnt2d P = PC2->Value( 0.5*(f2+l2) );
 820         Geom2dAPI_ProjectPointOnCurve proj (P, PC1, f1, l1);
 821         if (proj.NbPoints() &&
 822             proj.LowerDistance() <= Precision::Confusion())
 823           return Standard_True;
 824       }
 825       else {
 826         Handle(Geom_Curve) C1 = BRep_Tool::Curve( split ,f1,l1);
 827         Handle(Geom_Curve) C2 = BRep_Tool::Curve( NewE  ,f2,l2);
 828         gp_Pnt P = C2->Value( 0.5*(f2+l2) );
 829         GeomAPI_ProjectPointOnCurve proj (P, C1, f1, l1);
 830         if (proj.NbPoints() &&
 831             proj.LowerDistance() <= Precision::Confusion())
 832           return Standard_True;
 833       }
 834     }
 835     else {
 836       Standard_Real u3 = BRep_Tool::Parameter( TopoDS::Vertex(V3), OldE);
 837       Standard_Real u4 = BRep_Tool::Parameter( TopoDS::Vertex(V4), OldE);
 838
 839       Standard_Real f,l, u;
 840       BRep_Tool::Range( NewE, f,l);
 841       u = 0.5*(f+l);
 842       f = Min(u3,u4);
 843       l = Max(u3,u4);
 844
 845       if (u <= l && u >= f)
 846         return Standard_True;
 847     }
 848   }
 849   return Standard_False;
 850 }
 851
 852 //=======================================================================
 853 //function : SectionEdgeFaces
 854 //purpose  : return faces cut by section edge
 855 //=======================================================================
 856
 857 const TopTools_ListOfShape&
 858   Partition_Inter3d::SectionEdgeFaces(const TopoDS_Edge& SecE) const
 859 {
 860   return mySectionEdgesAD->Ascendant( SecE );
 861 }