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