src/BRepFill/BRepFill_Pipe.cxx

   1 // Created on: 1994-06-07
   2 // Created by: Bruno DUMORTIER
   3 // Copyright (c) 1994-1999 Matra Datavision
   4 // Copyright (c) 1999-2014 OPEN CASCADE SAS
   5 //
   6 // This file is part of Open CASCADE Technology software library.
   7 //
   8 // This library is free software; you can redistribute it and/or modify it under
   9 // the terms of the GNU Lesser General Public License version 2.1 as published
  10 // by the Free Software Foundation, with special exception defined in the file
  11 // OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
  12 // distribution for complete text of the license and disclaimer of any warranty.
  13 //
  14 // Alternatively, this file may be used under the terms of Open CASCADE
  15 // commercial license or contractual agreement.
  16
  17 #include <BRepFill_Pipe.ixx>
  18
  19 #include <Standard_ErrorHandler.hxx>
  20
  21 #include <BRep_Tool.hxx>
  22 #include <BRep_Builder.hxx>
  23 #include <BRepClass3d_SolidClassifier.hxx>
  24 #include <BRepLib_MakeVertex.hxx>
  25 #include <BRepTools_Substitution.hxx>
  26
  27 #include <GeomFill_CorrectedFrenet.hxx>
  28 #include <GeomFill_Frenet.hxx>
  29 #include <GeomFill_DiscreteTrihedron.hxx>
  30 #include <GeomFill_CurveAndTrihedron.hxx>
  31
  32 #include <BRepFill_SectionPlacement.hxx>
  33 #include <BRepFill_ShapeLaw.hxx>
  34 #include <BRepFill_Edge3DLaw.hxx>
  35 #include <BRepFill_Sweep.hxx>
  36
  37 #include <GeomAbs_Shape.hxx>
  38 #include <TopExp.hxx>
  39 #include <TopAbs_ShapeEnum.hxx>
  40 #include <TopoDS.hxx>
  41 #include <TopoDS_Shell.hxx>
  42 #include <TopoDS_Solid.hxx>
  43 #include <TopoDS_Compound.hxx>
  44 #include <TopoDS_Iterator.hxx>
  45 #include <TopTools_DataMapOfShapeInteger.hxx>
  46 #include <TColStd_DataMapOfIntegerInteger.hxx>
  47 #include <TColStd_DataMapIteratorOfDataMapOfIntegerInteger.hxx>
  48
  49 #include <Precision.hxx>
  50 #include <Standard_NotImplemented.hxx>
  51
  52 #include <Geom_TrimmedCurve.hxx>
  53 #include <Geom_OffsetCurve.hxx>
  54 #include <Geom_BSplineCurve.hxx>
  55 #include <BRepBuilderAPI_Transform.hxx>
  56 #include <BRepBuilderAPI_Copy.hxx>
  57 #include <TopTools_SequenceOfShape.hxx>
  58 #include <TopTools_ListIteratorOfListOfShape.hxx>
  59 #include <BRepLib.hxx>
  60
  61 #include <Geom2dAdaptor_HCurve.hxx>
  62 #include <GeomAdaptor_HSurface.hxx>
  63 #include <Adaptor3d_CurveOnSurface.hxx>
  64
  65 #include <ShapeUpgrade_RemoveLocations.hxx>
  66
  67 #ifdef DRAW
  68 #include <DBRep.hxx>
  69 static Standard_Boolean Affich = 0;
  70 #endif
  71
  72 static void ReverseModifiedEdges(TopoDS_Shape& aShape,
  73                                  TopTools_MapOfShape& Emap)
  74 {
  75   TopExp_Explorer Explo(aShape, TopAbs_FACE);
  76   BRep_Builder BB;
  77
  78   for (; Explo.More(); Explo.Next())
  79   {
  80     TopoDS_Shape aFace = Explo.Current();
  81     TopoDS_Iterator itf(aFace);
  82     for (; itf.More(); itf.Next())
  83     {
  84       TopoDS_Shape aWire = itf.Value();
  85       TopTools_ListOfShape Ledges;
  86       TopoDS_Iterator itw(aWire);
  87       for (; itw.More(); itw.Next())
  88         Ledges.Append(itw.Value());
  89
  90       aWire.Free(Standard_True);
  91       TopTools_ListIteratorOfListOfShape itl(Ledges);
  92       for (; itl.More(); itl.Next())
  93         BB.Remove(aWire, itl.Value());
  94
  95       for (itl.Initialize(Ledges); itl.More(); itl.Next())
  96       {
  97         TopoDS_Shape anEdge = itl.Value();
  98         if (Emap.Contains(anEdge))
  99           anEdge.Reverse();
 100         BB.Add(aWire, anEdge);
 101       }
 102     }
 103   }
 104 }
 105
 106 static void UpdateTolFromTopOrBottomPCurve(const TopoDS_Face& aFace,
 107                                            TopoDS_Edge& anEdge)
 108 {
 109   Standard_Real fpar, lpar;
 110   Handle(Geom2d_Curve) aPCurve = BRep_Tool::CurveOnSurface(anEdge, aFace, fpar, lpar);
 111   if (aPCurve.IsNull())
 112     return;
 113
 114   Handle(Geom_Curve) aCurve = BRep_Tool::Curve(anEdge, fpar, lpar);
 115   if (aCurve.IsNull())
 116     return;
 117
 118   Handle(Geom2dAdaptor_HCurve) GAHC2d = new Geom2dAdaptor_HCurve(aPCurve, fpar, lpar);
 119   Handle(Geom_Surface) aSurf = BRep_Tool::Surface(aFace);
 120   Handle(GeomAdaptor_HSurface) GAHS = new GeomAdaptor_HSurface(aSurf);
 121   Adaptor3d_CurveOnSurface ConS(GAHC2d, GAHS);
 122
 123   Standard_Real Tol = BRep_Tool::Tolerance(anEdge);
 124   Standard_Real InitTol = Tol;
 125   Standard_Real TolTol = Tol*Tol;
 126   const Standard_Integer NCONTROL = 22;
 127   Standard_Real delta = (lpar - fpar)/NCONTROL;
 128   for (Standard_Integer i = 0; i <= NCONTROL; i++)
 129   {
 130     Standard_Real par = fpar + i*delta;
 131     gp_Pnt pnt = aCurve->Value(par);
 132     gp_Pnt prj = ConS.Value(par);
 133     Standard_Real sqdist = pnt.SquareDistance(prj);
 134     if (sqdist > TolTol)
 135       TolTol = sqdist;
 136   }
 137   Tol = 1.00005 * Sqrt(TolTol);
 138   if (Tol >= InitTol)
 139   {
 140     BRep_Builder BB;
 141     BB.UpdateEdge(anEdge, Tol);
 142     TopoDS_Iterator itv(anEdge);
 143     for (; itv.More(); itv.Next())
 144     {
 145       TopoDS_Vertex aVertex = TopoDS::Vertex(itv.Value());
 146       BB.UpdateVertex(aVertex, Tol);
 147     }
 148   }
 149 }
 150
 151 //=======================================================================
 152 //function : BRepFill_Pipe
 153 //purpose  :
 154 //=======================================================================
 155
 156 BRepFill_Pipe::BRepFill_Pipe()
 157 {
 158   myDegmax = 11;
 159   mySegmax = 100;
 160   myContinuity = GeomAbs_C2;
 161   myMode = GeomFill_IsCorrectedFrenet;
 162   myForceApproxC1 = Standard_False;
 163
 164   myCurIndexOfSectionEdge = 1;
 165 }
 166
 167
 168 //=======================================================================
 169 //function : BRepFill_Pipe
 170 //purpose  :
 171 //=======================================================================
 172
 173 BRepFill_Pipe::BRepFill_Pipe(const TopoDS_Wire&  Spine,
 174                              const TopoDS_Shape& Profile,
 175                              const GeomFill_Trihedron aMode,
 176                              const Standard_Boolean ForceApproxC1,
 177                              const Standard_Boolean KPart)
 178
 179 {
 180   myDegmax = 11;
 181   mySegmax = 100;
 182
 183   myMode = GeomFill_IsCorrectedFrenet;
 184   if (aMode == GeomFill_IsFrenet ||
 185       aMode == GeomFill_IsCorrectedFrenet ||
 186       aMode == GeomFill_IsDiscreteTrihedron)
 187     myMode = aMode;
 188
 189   myContinuity = GeomAbs_C2;
 190   if (myMode == GeomFill_IsDiscreteTrihedron)
 191     myContinuity = GeomAbs_C0;
 192
 193   myForceApproxC1 = ForceApproxC1;
 194
 195   myCurIndexOfSectionEdge = 1;
 196
 197   Perform(Spine, Profile, KPart);
 198 }
 199
 200
 201 //=======================================================================
 202 //function : Perform
 203 //purpose  :
 204 //=======================================================================
 205
 206 void BRepFill_Pipe::Perform(const TopoDS_Wire&  Spine,
 207                             const TopoDS_Shape& Profile,
 208                             const Standard_Boolean /*KPart*/)
 209
 210 {
 211   mySections.Nullify();
 212   myFaces.Nullify();
 213   myEdges.Nullify();
 214
 215   mySpine   = Spine;
 216   myProfile = Profile;
 217
 218   DefineRealSegmax();
 219
 220   BRepTools_WireExplorer wexp;
 221   TopoDS_Shape TheProf;
 222
 223   Handle(GeomFill_TrihedronLaw) TLaw;
 224   switch (myMode)
 225   {
 226   case GeomFill_IsFrenet:
 227     TLaw = new GeomFill_Frenet();
 228     break;
 229   case GeomFill_IsCorrectedFrenet:
 230     TLaw = new GeomFill_CorrectedFrenet();
 231     break;
 232   case GeomFill_IsDiscreteTrihedron:
 233     TLaw = new GeomFill_DiscreteTrihedron();
 234     break;
 235   default:
 236     break;
 237   }
 238   Handle(GeomFill_CurveAndTrihedron) Loc =
 239     new (GeomFill_CurveAndTrihedron) (TLaw);
 240   myLoc = new (BRepFill_Edge3DLaw) (mySpine, Loc);
 241   if (myLoc->NbLaw() == 0) {
 242     return; // Degenerated case
 243   }
 244   myLoc->TransformInG0Law(); // Set into continuity
 245
 246   BRepFill_SectionPlacement Place(myLoc, Profile);
 247   myTrsf = Place.Transformation();
 248
 249   TopLoc_Location Loc2(myTrsf), Loc1;
 250   Loc1 = Profile.Location();
 251   TopoDS_Shape aux;
 252   TheProf =  myProfile;
 253   TheProf.Location(Loc2.Multiplied(Loc1));
 254
 255   // Construct First && Last Shape
 256   Handle(GeomFill_LocationLaw) law;
 257
 258   gp_Mat M;
 259   gp_Vec V;
 260   gp_Trsf fila;
 261   Standard_Real first, last;
 262   myLoc->Law(1)->GetDomain(first, last);
 263   myLoc->Law(1)->D0(first, M, V);
 264     fila.SetValues(M(1,1), M(1,2), M(1,3), V.X(),
 265                    M(2,1), M(2,2), M(2,3), V.Y(),
 266                    M(3,1), M(3,2), M(3,3), V.Z(),
 267                    1.e-12, 1.e-14);
 268
 269   fila.Multiply(myTrsf);
 270   TopLoc_Location LocFirst(fila);
 271   myFirst = myProfile;
 272   if ( ! LocFirst.IsIdentity()) {
 273     //myFirst.Location( LocFirst.Multiplied(myProfile.Location()) );
 274     myFirst = BRepBuilderAPI_Transform(myProfile, fila, Standard_True); //copy
 275   }
 276
 277   ShapeUpgrade_RemoveLocations RemLoc;
 278   RemLoc.Remove(myFirst);
 279   myFirst = RemLoc.GetResult();
 280   TopLoc_Location theLoc = myFirst.Location();
 281   if (!theLoc.IsIdentity())
 282   {
 283     TopoDS_Shape NewMyFirst = BRepBuilderAPI_Copy(myFirst);
 284     RemLoc.Remove(NewMyFirst);
 285     NewMyFirst = RemLoc.GetResult();
 286     TopLoc_Location theIdentity;
 287     NewMyFirst.Location(theIdentity);
 288     myFirst = BRepBuilderAPI_Transform(NewMyFirst, theLoc.Transformation(), Standard_True);
 289   }
 290
 291   myLoc->Law(myLoc->NbLaw())->GetDomain(first, last);
 292   myLoc->Law(myLoc->NbLaw())->D0(last,M, V);
 293 //    try { // Not good, but there are no other means to test SetValues
 294   fila.SetValues(M(1,1), M(1,2), M(1,3), V.X(),
 295                  M(2,1), M(2,2), M(2,3), V.Y(),
 296                  M(3,1), M(3,2), M(3,3), V.Z(),
 297                  1.e-12, 1.e-14);
 298   fila.Multiply(myTrsf);
 299   TopLoc_Location LocLast(fila);
 300   if (! myLoc->IsClosed() || LocFirst != LocLast) {
 301     myLast = myProfile;
 302     if ( ! LocLast.IsIdentity()) {
 303       //myLast.Location(LocLast.Multiplied(myProfile.Location()) );
 304       myLast = BRepBuilderAPI_Transform(myProfile, fila, Standard_True); //copy
 305     }
 306   }
 307   else {
 308     myLast = myFirst;
 309   }
 310
 311   RemLoc.Remove(myLast);
 312   myLast = RemLoc.GetResult();
 313   theLoc = myLast.Location();
 314   if (!theLoc.IsIdentity())
 315   {
 316     TopoDS_Shape NewMyLast = BRepBuilderAPI_Copy(myLast);
 317     RemLoc.Remove(NewMyLast);
 318     NewMyLast = RemLoc.GetResult();
 319     TopLoc_Location theIdentity;
 320     NewMyLast.Location(theIdentity);
 321     myLast = BRepBuilderAPI_Transform(NewMyLast, theLoc.Transformation(), Standard_True);
 322   }
 323
 324 #if DRAW
 325   if (Affich) {
 326     DBRep::Set("theprof",  TheProf);
 327     DBRep::Set("thefirst", myFirst);
 328     DBRep::Set("thelast" , myLast);
 329  }
 330 #endif
 331
 332   myShape = MakeShape(TheProf, myFirst, myLast);
 333 }
 334
 335
 336 //=======================================================================
 337 //function : Spine
 338 //purpose  :
 339 //=======================================================================
 340
 341 const TopoDS_Shape& BRepFill_Pipe::Spine() const
 342 {
 343   return mySpine;
 344 }
 345
 346 //=======================================================================
 347 //function : Profile
 348 //purpose  :
 349 //=======================================================================
 350
 351 const TopoDS_Shape& BRepFill_Pipe::Profile() const
 352 {
 353   return myProfile;
 354 }
 355
 356 //=======================================================================
 357 //function : Shape
 358 //purpose  :
 359 //=======================================================================
 360
 361 const TopoDS_Shape& BRepFill_Pipe::Shape() const
 362 {
 363   return myShape;
 364 }
 365
 366
 367 //=======================================================================
 368 //function : FirstShape
 369 //purpose  :
 370 //=======================================================================
 371
 372 const TopoDS_Shape& BRepFill_Pipe::FirstShape() const
 373 {
 374   return myFirst;
 375 }
 376
 377
 378 //=======================================================================
 379 //function : LastShape
 380 //purpose  :
 381 //=======================================================================
 382
 383 const TopoDS_Shape& BRepFill_Pipe::LastShape() const
 384 {
 385   return myLast;
 386 }
 387
 388
 389 //=======================================================================
 390 //function : Face
 391 //purpose  :
 392 //=======================================================================
 393
 394 TopoDS_Face BRepFill_Pipe::Face(const TopoDS_Edge& ESpine,
 395                                 const TopoDS_Edge& EProfile)
 396 {
 397   TopoDS_Face theFace;
 398
 399   if ( BRep_Tool::Degenerated(EProfile))
 400     return theFace;
 401
 402   Standard_Integer ii, ispin = 0, iprof = 0, count = 0;
 403
 404   // *************************************************
 405   // Search if EProfile is an edge of myProfile
 406   // *************************************************
 407   iprof = FindEdge(myProfile, EProfile, count);
 408
 409   if (!iprof) Standard_DomainError::Raise(
 410               "BRepFill_Pipe::Face : Edge not in the Profile");
 411
 412
 413   // *************************************************
 414   // Search if ESpine  is an edge of mySpine and find
 415   // the index of the corresponding Filler
 416   // *************************************************
 417  for (ii=1; ii<=myLoc->NbLaw() && (!ispin); ii++)
 418     if  (ESpine.IsSame(myLoc->Edge(ii))) ispin = ii;
 419
 420   if (!ispin) Standard_DomainError::Raise(
 421     "BRepFill_Pipe::Edge  : Edge not in the Spine");
 422
 423   theFace = TopoDS::Face(myFaces->Value(iprof, ispin));
 424   return theFace;
 425
 426 }
 427
 428 //=======================================================================
 429 //function : Edge
 430 //purpose  :
 431 //=======================================================================
 432 TopoDS_Edge BRepFill_Pipe::Edge(const TopoDS_Edge&   ESpine,
 433                                 const TopoDS_Vertex& VProfile)
 434 {
 435   Standard_Integer ii, ispin = 0, iprof = 0, count = 0;;
 436
 437   // *************************************************
 438   // Search if VProfile is a Vertex of myProfile
 439   // *************************************************
 440   iprof = FindVertex(myProfile, VProfile, count);
 441   if (!iprof) Standard_DomainError::Raise(
 442         "BRepFill_Pipe::Edge : Vertex not in the Profile");
 443
 444
 445   // *************************************************
 446   // Search if ESpine  is an edge of mySpine and find
 447   // the index of the corresponding Filler
 448   // *************************************************
 449
 450   for (ii=1; ii<=myLoc->NbLaw() && (!ispin); ii++)
 451     if  (ESpine.IsSame(myLoc->Edge(ii))) ispin = ii;
 452
 453   if (!ispin) Standard_DomainError::Raise(
 454     "BRepFill_Pipe::Edge  : Edge not in the Spine");
 455
 456
 457   // *************************************************
 458   // Generate the corresponding Shape
 459   // *************************************************
 460   TopoDS_Edge theEdge;
 461   theEdge = TopoDS::Edge(myEdges->Value(iprof, ispin));
 462
 463   return theEdge;
 464
 465 }
 466
 467
 468 //=======================================================================
 469 //function : Section
 470 //purpose  :
 471 //=======================================================================
 472
 473 TopoDS_Shape BRepFill_Pipe::Section(const TopoDS_Vertex& VSpine) const
 474 {
 475   TopoDS_Iterator it, itv;
 476
 477   Standard_Integer ii, ispin = 0;
 478
 479   TopoDS_Shape curSect = myProfile;
 480
 481   // *************************************************
 482   // Search if ESpine  is an edge of mySpine and find
 483   // the index of the corresponding Filler
 484   // *************************************************
 485
 486   // iterate on all the edges of mySpine
 487  for (ii=1; ii<=myLoc->NbLaw()+1 && (!ispin); ii++)
 488     if  (VSpine.IsSame(myLoc->Vertex(ii))) ispin = ii;
 489
 490   if (!ispin) Standard_DomainError::Raise(
 491     "BRepFill_Pipe::Section  : Vertex not in the Spine");
 492
 493   BRep_Builder B;
 494   TopoDS_Compound Comp;
 495   B.MakeCompound(Comp);
 496   for (ii=1; ii<=mySections->ColLength(); ii++)
 497     B.Add(Comp, mySections->Value(ii, ispin));
 498
 499   return Comp;
 500 }
 501
 502 //=======================================================================
 503 //function : PipeLine
 504 //purpose  : Construct a wire by sweeping of a point
 505 //=======================================================================
 506
 507 TopoDS_Wire BRepFill_Pipe::PipeLine(const gp_Pnt& Point)
 508 {
 509  // Postioning
 510  gp_Pnt P;
 511  P = Point;
 512  P.Transform(myTrsf);
 513
 514  BRepLib_MakeVertex MkV(P);
 515  Handle(BRepFill_ShapeLaw) Section =
 516         new (BRepFill_ShapeLaw) (MkV.Vertex());
 517
 518  // Sweeping
 519  BRepFill_Sweep MkSw(Section, myLoc, Standard_True);
 520  MkSw.SetForceApproxC1(myForceApproxC1);
 521  MkSw.Build( myReversedEdges, myTapes,
 522              BRepFill_Modified, myContinuity, GeomFill_Location, myDegmax, mySegmax );
 523  TopoDS_Shape aLocalShape = MkSw.Shape();
 524  return TopoDS::Wire(aLocalShape);
 525 // return TopoDS::Wire(MkSw.Shape());
 526 }
 527
 528 //=======================================================================
 529 //function : MakeShape
 530 //purpose  :
 531 //=======================================================================
 532
 533 TopoDS_Shape BRepFill_Pipe::MakeShape(const TopoDS_Shape& S,
 534                                       const TopoDS_Shape& FirstShape,
 535                                       const TopoDS_Shape& LastShape)
 536 {
 537   TopoDS_Shape result;
 538   BRep_Builder B;
 539   Standard_Boolean explode = Standard_False;
 540   TopoDS_Shape TheS, TheFirst, TheLast;
 541   Standard_Integer InitialLength = 0;
 542   TheS = S;
 543   TheFirst = FirstShape;
 544   TheLast = LastShape;
 545   if (! myFaces.IsNull()) InitialLength = myFaces->ColLength();
 546
 547   // there are two kinds of generation
 548   //  1. generate with S from each Filler (Vertex, Edge)
 549   //  2. call MakeShape recursively on the subshapes of S
 550   //
 551   // explode is True in the second case
 552
 553   // create the result empty
 554
 555   switch (S.ShapeType()) {
 556
 557   case TopAbs_VERTEX :
 558     {
 559       B.MakeWire(TopoDS::Wire(result));
 560       break;
 561     }
 562
 563   case TopAbs_EDGE :
 564     {
 565       TopoDS_Wire W;
 566       B.MakeShell(TopoDS::Shell(result));
 567       B.MakeWire(W);
 568       B.Add(W, S);
 569       W.Closed(S.Closed());
 570       TheS = W;
 571       if (!FirstShape.IsNull()) {
 572         B.MakeWire(W);
 573         B.Add(W, FirstShape);
 574         W.Closed(FirstShape.Closed());
 575         TheFirst = W;
 576       }
 577       if (!LastShape.IsNull()) {
 578         B.MakeWire(W);
 579         B.Add(W, LastShape);
 580         W.Closed(LastShape.Closed());
 581         TheLast = W;
 582       }
 583       break;
 584     }
 585
 586   case TopAbs_WIRE :
 587     B.MakeShell(TopoDS::Shell(result));
 588     break;
 589
 590   case TopAbs_FACE :
 591     {
 592       B.MakeShell(TopoDS::Shell(result));
 593       explode = Standard_True;
 594       if ( !mySpine.Closed() && !TheFirst.IsNull()) {
 595          B.Add(result, TheFirst.Reversed());
 596       }
 597       break;
 598     }
 599
 600   case TopAbs_SHELL :
 601     {
 602       B.MakeCompSolid(TopoDS::CompSolid(result));
 603       explode = Standard_True;
 604       break;
 605     }
 606
 607   case TopAbs_SOLID :
 608   case TopAbs_COMPSOLID :
 609     Standard_DomainError::Raise("BRepFill_Pipe::profile contains solids");
 610     break;
 611
 612   case TopAbs_COMPOUND :
 613     {
 614       B.MakeCompound(TopoDS::Compound(result));
 615       explode = Standard_True;
 616       break;
 617     }
 618   default:
 619     break;
 620   }
 621
 622   if (explode) {
 623     // add the subshapes
 624     TopoDS_Iterator itFirst, itLast;
 625     TopoDS_Shape first, last;
 626     if (!TheFirst.IsNull()) itFirst.Initialize(TheFirst);
 627     if (!TheLast.IsNull()) itLast.Initialize(TheLast);
 628
 629     for (TopoDS_Iterator it(S); it.More(); it.Next()) {
 630       if (!TheFirst.IsNull()) first = itFirst.Value();
 631       if (!TheLast.IsNull())  last = itLast.Value();
 632       if (TheS.ShapeType() == TopAbs_FACE )
 633         MakeShape(it.Value(), first, last);
 634       else
 635         B.Add(result,MakeShape(it.Value(), first, last));
 636
 637       if (!TheFirst.IsNull()) itFirst.Next();
 638       if (!TheLast.IsNull())  itLast.Next();
 639     }
 640   }
 641
 642   else {
 643     if (TheS.ShapeType() == TopAbs_VERTEX ) {
 644      Handle(BRepFill_ShapeLaw) Section =
 645         new (BRepFill_ShapeLaw) (TopoDS::Vertex(TheS));
 646       BRepFill_Sweep MkSw(Section, myLoc, Standard_True);
 647       MkSw.SetForceApproxC1(myForceApproxC1);
 648       MkSw.Build( myReversedEdges, myTapes,
 649                   BRepFill_Modified, myContinuity, GeomFill_Location, myDegmax, mySegmax );
 650       result = MkSw.Shape();
 651
 652       Handle(TopTools_HArray2OfShape) aSections = MkSw.Sections();
 653
 654       if (aSections.IsNull() == Standard_False) {
 655         const Standard_Integer aVLast = aSections->UpperCol();
 656
 657         myFirst = aSections->Value(1, 1);
 658         myLast  = aSections->Value(1, aVLast);
 659       }
 660     }
 661
 662     if (TheS.ShapeType() == TopAbs_WIRE ) {
 663       Handle(BRepFill_ShapeLaw) Section =
 664         new (BRepFill_ShapeLaw) (TopoDS::Wire(TheS));
 665       BRepFill_Sweep MkSw(Section, myLoc, Standard_True);
 666       MkSw.SetBounds(TopoDS::Wire(TheFirst),
 667                      TopoDS::Wire(TheLast));
 668       MkSw.SetForceApproxC1(myForceApproxC1);
 669       MkSw.Build( myReversedEdges, myTapes,
 670                   BRepFill_Modified, myContinuity, GeomFill_Location, myDegmax, mySegmax );
 671       result = MkSw.Shape();
 672       //Correct <myFirst> and <myLast>
 673       ReverseModifiedEdges(myFirst, myReversedEdges);
 674       ReverseModifiedEdges(myLast, myReversedEdges);
 675
 676       // Labeling of elements
 677       if (mySections.IsNull()) {
 678         myFaces    = MkSw.SubShape();
 679         mySections = MkSw.Sections();
 680         myEdges    = MkSw.InterFaces();
 681       }
 682       else {
 683         Handle(TopTools_HArray2OfShape) Aux, Somme;
 684         Standard_Integer length;
 685         Standard_Integer ii, jj, kk;
 686
 687         Aux = MkSw.SubShape();
 688         length = Aux->ColLength() + myFaces->ColLength();
 689         Somme = new (TopTools_HArray2OfShape) (1, length, 1,
 690                                                Aux->RowLength());
 691         for (jj=1; jj<=myFaces->RowLength(); jj++) {
 692           for (ii=1; ii<=myFaces->ColLength(); ii++)
 693             Somme->SetValue(ii, jj, myFaces->Value(ii, jj));
 694
 695           for (kk=1, ii=myFaces->ColLength()+1;
 696                kk <=Aux->ColLength(); kk++, ii++)
 697             Somme->SetValue(ii, jj, Aux->Value(kk, jj));
 698         }
 699         myFaces = Somme;
 700
 701         Aux = MkSw.Sections();
 702         length = Aux->ColLength() + mySections->ColLength();
 703         Somme = new (TopTools_HArray2OfShape) (1, length, 1,
 704                                                Aux->RowLength());
 705         for (jj=1; jj<=mySections->RowLength(); jj++) {
 706           for (ii=1; ii<=mySections->ColLength(); ii++)
 707             Somme->SetValue(ii, jj, mySections->Value(ii, jj));
 708
 709           myCurIndexOfSectionEdge = mySections->ColLength()+1;
 710
 711           for (kk=1, ii=mySections->ColLength()+1;
 712                kk <=Aux->ColLength(); kk++, ii++)
 713             Somme->SetValue(ii, jj, Aux->Value(kk, jj));
 714         }
 715         mySections = Somme;
 716
 717         Aux = MkSw.InterFaces();
 718         length = Aux->ColLength() + myEdges->ColLength();
 719         Somme = new (TopTools_HArray2OfShape) (1, length, 1,
 720                                                Aux->RowLength());
 721         for (jj=1; jj<=myEdges->RowLength(); jj++) {
 722           for (ii=1; ii<=myEdges->ColLength(); ii++)
 723             Somme->SetValue(ii, jj, myEdges->Value(ii, jj));
 724
 725           for (kk=1, ii=myEdges->ColLength()+1;
 726                kk <=Aux->ColLength(); kk++, ii++)
 727             Somme->SetValue(ii, jj, Aux->Value(kk, jj));
 728         }
 729
 730         myEdges = Somme;
 731       }
 732     }
 733   }
 734
 735   if ( TheS.ShapeType() == TopAbs_FACE ) {
 736     Standard_Integer ii, jj;
 737     //jgv
 738     TopExp_Explorer Explo(result, TopAbs_FACE);
 739     for (; Explo.More(); Explo.Next())
 740     {
 741       TopoDS_Shape aFace = Explo.Current();
 742       RebuildTopOrBottomFace(aFace.Reversed(), Standard_True); //top face was reversed
 743     }
 744     /////
 745     TopoDS_Face F;
 746     for (ii=InitialLength+1; ii<=myFaces->ColLength(); ii++) {
 747       for (jj=1; jj<=myFaces->RowLength(); jj++) {
 748         F = TopoDS::Face(myFaces->Value(ii, jj));
 749         if (!F.IsNull()) B.Add(result, F);
 750       }
 751     }
 752
 753     if ( !mySpine.Closed()) {
 754       // if Spine is not closed
 755       // add the last face of the solid
 756
 757       //jgv
 758       RebuildTopOrBottomFace(TheLast, Standard_False); //bottom face
 759       /////
 760       B.Add(result, TopoDS::Face(TheLast));
 761     }
 762
 763     TopoDS_Solid solid;
 764     BRep_Builder BS;
 765     BS.MakeSolid(solid);
 766
 767     result.Closed(Standard_True);
 768     BS.Add(solid,TopoDS::Shell(result));
 769
 770     BRepClass3d_SolidClassifier SC(solid);
 771     SC.PerformInfinitePoint(Precision::Confusion());
 772     if ( SC.State() == TopAbs_IN) {
 773       BS.MakeSolid(solid);
 774       TopoDS_Shape aLocalShape = result.Reversed();
 775       BS.Add(solid,TopoDS::Shell(aLocalShape));
 776 //      BS.Add(solid,TopoDS::Shell(result.Reversed()));
 777     }
 778     return solid;
 779   }
 780   else {
 781     return result;
 782   }
 783 }
 784
 785 //============================================================================
 786 //function : FindEdge
 787 //purpose  : Find the number of edge corresponding to the edge of the profile.
 788 //============================================================================
 789
 790 Standard_Integer BRepFill_Pipe::FindEdge(const TopoDS_Shape& S,
 791                                          const TopoDS_Edge& E,
 792                                          Standard_Integer& InitialLength) const
 793 {
 794   Standard_Integer result = 0;
 795
 796   switch (S.ShapeType()) {
 797
 798   case TopAbs_EDGE :
 799     {
 800       InitialLength++;
 801       if (S.IsSame(E)) result = InitialLength;
 802       break;
 803     }
 804
 805   case TopAbs_WIRE :
 806     {
 807       Standard_Integer ii = InitialLength+1;
 808       Handle(BRepFill_ShapeLaw) Section =
 809         new (BRepFill_ShapeLaw) (TopoDS::Wire(S), Standard_False);
 810       InitialLength += Section->NbLaw();
 811
 812       for (; (ii<=InitialLength) && (!result); ii++) {
 813         if (E.IsSame(Section->Edge(ii)) ) result = ii;
 814       }
 815       break;
 816     }
 817
 818   case TopAbs_FACE :
 819   case TopAbs_SHELL :
 820   case TopAbs_COMPOUND :
 821     {
 822      for (TopoDS_Iterator it(S); it.More() && (!result); it.Next())
 823        result = FindEdge(it.Value(), E, InitialLength );
 824      break;
 825     }
 826
 827   case TopAbs_SOLID :
 828   case TopAbs_COMPSOLID :
 829     Standard_DomainError::Raise("BRepFill_Pipe::SOLID or COMPSOLID");
 830     break;
 831   default:
 832     break;
 833   }
 834
 835   return result;
 836 }
 837
 838 //=======================================================================
 839 //function : FindVertex
 840 //purpose  : Find the number of edge corresponding to an edge of the profile.
 841 //=======================================================================
 842
 843 Standard_Integer BRepFill_Pipe::FindVertex(const TopoDS_Shape& S,
 844                                            const TopoDS_Vertex& V,
 845                                            Standard_Integer& InitialLength) const
 846 {
 847   Standard_Integer result = 0;
 848
 849   switch (S.ShapeType()) {
 850   case TopAbs_VERTEX :
 851     {
 852       InitialLength++;
 853       if (S.IsSame(V)) result = InitialLength;
 854       break;
 855     }
 856
 857   case TopAbs_EDGE :
 858     {
 859       TopoDS_Vertex VF, VL;
 860       TopExp::Vertices(TopoDS::Edge(S), VF, VL);
 861       if (S.Orientation() == TopAbs_REVERSED) {
 862         TopoDS_Vertex aux;
 863         aux = VF; VF = VL; VL = aux;
 864       }
 865       if (VF.IsSame(V)) result = InitialLength+1;
 866       else if (VL.IsSame(V)) result = InitialLength+2;
 867       InitialLength += 2;
 868       break;
 869     }
 870
 871   case TopAbs_WIRE :
 872     {
 873       Standard_Integer ii = InitialLength+1;
 874       Handle(BRepFill_ShapeLaw) Section =
 875         new (BRepFill_ShapeLaw) (TopoDS::Wire(S), Standard_False);
 876       InitialLength += Section->NbLaw()+1;
 877
 878       for (; (ii<=InitialLength) && (!result); ii++) {
 879         if (V.IsSame(Section->Vertex(ii, 0.)) ) result = ii;
 880       }
 881       break;
 882     }
 883
 884   case TopAbs_FACE :
 885   case TopAbs_SHELL :
 886   case TopAbs_COMPOUND :
 887     {
 888      for (TopoDS_Iterator it(S); it.More() && (!result); it.Next())
 889        result = FindVertex(it.Value(), V, InitialLength);
 890      break;
 891     }
 892
 893   case TopAbs_SOLID :
 894   case TopAbs_COMPSOLID :
 895     Standard_DomainError::Raise("BRepFill_Pipe::SOLID or COMPSOLID");
 896     break;
 897   default:
 898     break;
 899   }
 900
 901   return result;
 902 }
 903
 904 //=======================================================================
 905 //function : DefineRealSegmax
 906 //purpose  : Defines the real number of segments
 907 //           required in the case of bspline spine
 908 //=======================================================================
 909
 910 void BRepFill_Pipe::DefineRealSegmax()
 911 {
 912   Standard_Integer RealSegmax = 0;
 913
 914   TopoDS_Iterator iter(mySpine);
 915   for (; iter.More(); iter.Next())
 916     {
 917       TopoDS_Edge E = TopoDS::Edge(iter.Value());
 918       Standard_Real first, last;
 919       Handle(Geom_Curve) C = BRep_Tool::Curve( E, first, last );
 920       if (C.IsNull())
 921         continue;
 922       while (C->DynamicType() == STANDARD_TYPE(Geom_TrimmedCurve) ||
 923              C->DynamicType() == STANDARD_TYPE(Geom_OffsetCurve))
 924         {
 925           if (C->DynamicType() == STANDARD_TYPE(Geom_TrimmedCurve))
 926             C = (*((Handle(Geom_TrimmedCurve)*)&C))->BasisCurve();
 927           if (C->DynamicType() == STANDARD_TYPE(Geom_OffsetCurve))
 928             C = (*((Handle(Geom_OffsetCurve)*)&C))->BasisCurve();
 929         }
 930       if (C->DynamicType() == STANDARD_TYPE(Geom_BSplineCurve))
 931         {
 932           const Handle(Geom_BSplineCurve)& BC = *((Handle(Geom_BSplineCurve)*)&C);
 933           Standard_Integer NbKnots = BC->NbKnots();
 934           Standard_Integer RealNbKnots = NbKnots;
 935           if (first > BC->FirstParameter())
 936             {
 937               Standard_Integer I1, I2;
 938               BC->LocateU( first, Precision::PConfusion(), I1, I2 );
 939               RealNbKnots -= I1-1;
 940             }
 941           if (last < BC->LastParameter())
 942             {
 943               Standard_Integer I1, I2;
 944               BC->LocateU( last, Precision::PConfusion(), I1, I2 );
 945               RealNbKnots -= NbKnots-I2;
 946             }
 947           RealSegmax += RealNbKnots-1;
 948         }
 949     }
 950
 951   if (mySegmax < RealSegmax)
 952     mySegmax = RealSegmax;
 953 }
 954
 955 //=======================================================================
 956 //function : RebuildTopOrBottomFace
 957 //purpose  : Correct orientation of v-iso edges
 958 //           according to new 3d and 2d curves taken from swept surfaces
 959 //=======================================================================
 960
 961 void BRepFill_Pipe::RebuildTopOrBottomFace(const TopoDS_Shape& aFace,
 962                                            const Standard_Boolean IsTop) const
 963 {
 964   Standard_Integer IndexOfSection =
 965     (IsTop)? 1 : mySections->RowLength();
 966
 967   Standard_Integer ii;
 968   BRep_Builder BB;
 969   TopoDS_Iterator itf(aFace);
 970   for (; itf.More(); itf.Next())
 971   {
 972     TopoDS_Shape aWire = itf.Value();
 973     TopTools_SequenceOfShape InitEdges;
 974     TopTools_SequenceOfShape ResEdges;
 975     TopoDS_Iterator itw(aWire);
 976     for (; itw.More(); itw.Next())
 977     {
 978       TopoDS_Shape anEdge = itw.Value();
 979       for (ii = myCurIndexOfSectionEdge; ii <= mySections->ColLength(); ii++)
 980       {
 981         TopoDS_Shape aVisoEdge = mySections->Value(ii, IndexOfSection);
 982         if (anEdge.IsSame(aVisoEdge))
 983         {
 984           InitEdges.Append(anEdge);
 985           ResEdges.Append(aVisoEdge);
 986           break;
 987         }
 988       }
 989     }
 990     aWire.Free(Standard_True);
 991     for (ii = 1; ii <= InitEdges.Length(); ii++)
 992     {
 993       BB.Remove(aWire, InitEdges(ii));
 994       UpdateTolFromTopOrBottomPCurve(TopoDS::Face(aFace), TopoDS::Edge(ResEdges(ii)));
 995       BB.Add(aWire, ResEdges(ii));
 996     }
 997   }
 998 }