src/BRepFill/BRepFill_Filling.cxx

   1 // Created on: 1998-08-26
   2 // Created by: Julia GERASIMOVA
   3 // Copyright (c) 1998-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_Filling.ixx>
  18 #include <BRepFill_EdgeFaceAndOrder.hxx>
  19 #include <BRepFill_FaceAndOrder.hxx>
  20 #include <BRepAdaptor_HSurface.hxx>
  21 #include <BRepAdaptor_HCurve.hxx>
  22 #include <BRepAdaptor_HCurve2d.hxx>
  23 #include <BRepFill_CurveConstraint.hxx>
  24
  25 #include <Geom2dAdaptor_HCurve.hxx>
  26 #include <GeomAdaptor_HSurface.hxx>
  27 #include <Adaptor3d_CurveOnSurface.hxx>
  28 #include <Adaptor3d_HCurveOnSurface.hxx>
  29
  30 #include <GeomPlate_MakeApprox.hxx>
  31 #include <GeomPlate_CurveConstraint.hxx>
  32 #include <GeomPlate_PointConstraint.hxx>
  33
  34 #include <BRepLib_MakeEdge.hxx>
  35 #include <BRepLib_MakeEdge2d.hxx>
  36 #include <BRepLib_MakeWire.hxx>
  37 #include <TopTools_ListIteratorOfListOfShape.hxx>
  38 #include <TopoDS.hxx>
  39 #include <TopoDS_Wire.hxx>
  40 #include <BRepLib_MakeFace.hxx>
  41 #include <TColStd_HArray1OfReal.hxx>
  42 #include <BRep_ListIteratorOfListOfCurveRepresentation.hxx>
  43 #include <BRep_TEdge.hxx>
  44 #include <BRep_CurveRepresentation.hxx>
  45 #include <TopExp.hxx>
  46 #include <BRep_Tool.hxx>
  47 #include <GeomAPI_ProjectPointOnSurf.hxx>
  48 #include <TopoDS_Vertex.hxx>
  49
  50 #include <Precision.hxx>
  51 #include <GeomProjLib.hxx>
  52 #include <BRepTools.hxx>
  53 #include <BRep_Builder.hxx>
  54 #include <BRepLib.hxx>
  55 #include <TColgp_SequenceOfPnt.hxx>
  56 #include <TopTools_ListIteratorOfListOfShape.hxx>
  57 #include <BRepTools_WireExplorer.hxx>
  58 #include <TColgp_Array1OfPnt2d.hxx>
  59 #include <TColGeom2d_HArray1OfCurve.hxx>
  60
  61 #include <Geom2d_BezierCurve.hxx>
  62 #include <Geom2d_TrimmedCurve.hxx>
  63 #include <Geom_BSplineSurface.hxx>
  64
  65 #include <TColgp_SequenceOfXY.hxx>
  66 #include <GeomPlate_PlateG0Criterion.hxx>
  67 #include <BRep_TVertex.hxx>
  68
  69 static gp_Vec MakeFinVec( const TopoDS_Wire aWire, const TopoDS_Vertex aVertex )
  70 {
  71   TopoDS_Vertex Vfirst, Vlast, Origin;
  72   BRepTools_WireExplorer Explo( aWire );
  73   for (; Explo.More(); Explo.Next())
  74     {
  75       TopExp::Vertices( Explo.Current(), Vfirst, Vlast );
  76       if (Vfirst.IsSame( aVertex ))
  77         {
  78           Origin = Vlast;
  79           break;
  80         }
  81       if (Vlast.IsSame( aVertex ))
  82         {
  83           Origin = Vfirst;
  84           break;
  85         }
  86     }
  87   return gp_Vec( BRep_Tool::Pnt( Origin ), BRep_Tool::Pnt( aVertex ) );
  88 }
  89
  90 static TopoDS_Wire WireFromList(TopTools_ListOfShape& Edges)
  91 {
  92   BRepLib_MakeWire MW;
  93   TopoDS_Edge anEdge = TopoDS::Edge(Edges.First());
  94   MW.Add(anEdge);
  95   Edges.RemoveFirst();
  96
  97   while (!Edges.IsEmpty())
  98   {
  99     TopoDS_Wire CurWire = MW.Wire();
 100     TopoDS_Vertex V1, V2;
 101     TopExp::Vertices(CurWire, V1, V2);
 102     TopTools_ListIteratorOfListOfShape itl(Edges);
 103     for (; itl.More(); itl.Next())
 104     {
 105       anEdge = TopoDS::Edge(itl.Value());
 106       TopoDS_Vertex V3, V4;
 107       TopExp::Vertices(anEdge, V3, V4);
 108       if (V1.IsSame(V3) || V1.IsSame(V4) ||
 109           V2.IsSame(V3) || V2.IsSame(V4))
 110         break;
 111     }
 112     MW.Add(anEdge);
 113     Edges.Remove(itl);
 114   }
 115
 116   return (MW.Wire());
 117 }
 118
 119 //=======================================================================
 120 //function : Constructor
 121 //purpose  :
 122 //======================================================================
 123 BRepFill_Filling::BRepFill_Filling( const Standard_Integer Degree,
 124                                     const Standard_Integer NbPtsOnCur,
 125                                     const Standard_Integer NbIter,
 126                                     const Standard_Boolean Anisotropie,
 127                                     const Standard_Real Tol2d,
 128                                     const Standard_Real Tol3d,
 129                                     const Standard_Real TolAng,
 130                                     const Standard_Real TolCurv,
 131                                     const Standard_Integer MaxDeg,
 132                                     const Standard_Integer MaxSegments )
 133 {
 134   myDegree = Degree;
 135   myNbPtsOnCur = NbPtsOnCur;
 136   myNbIter = NbIter;
 137   myAnisotropie = Anisotropie;
 138
 139   myTol2d = Tol2d;
 140   myTol3d = Tol3d;
 141   myTolAng = TolAng;
 142   myTolCurv = TolCurv;
 143
 144   myMaxDeg = MaxDeg;
 145   myMaxSegments = MaxSegments;
 146
 147   myIsInitFaceGiven = Standard_False;
 148
 149   myIsDone = Standard_False;
 150 }
 151
 152 //=======================================================================
 153 //function : SetConstrParam
 154 //purpose  :
 155 //======================================================================
 156 void BRepFill_Filling::SetConstrParam( const Standard_Real Tol2d,
 157                                        const Standard_Real Tol3d,
 158                                        const Standard_Real TolAng,
 159                                        const Standard_Real TolCurv )
 160 {
 161   myTol2d = Tol2d;
 162   myTol3d = Tol3d;
 163   myTolAng = TolAng;
 164   myTolCurv = TolCurv;
 165 }
 166
 167 //=======================================================================
 168 //function : SetResolParam
 169 //purpose  :
 170 //======================================================================
 171 void BRepFill_Filling::SetResolParam( const Standard_Integer Degree,
 172                                       const Standard_Integer NbPtsOnCur,
 173                                       const Standard_Integer NbIter,
 174                                       const Standard_Boolean Anisotropie )
 175 {
 176   myDegree = Degree;
 177   myNbPtsOnCur = NbPtsOnCur;
 178   myNbIter = NbIter;
 179   myAnisotropie = Anisotropie;
 180 }
 181
 182 //=======================================================================
 183 //function : SetApproxParam
 184 //purpose  :
 185 //======================================================================
 186 void BRepFill_Filling::SetApproxParam( const Standard_Integer MaxDeg,
 187                                        const Standard_Integer MaxSegments )
 188 {
 189   myMaxDeg      = MaxDeg;
 190   myMaxSegments = MaxSegments;
 191 }
 192
 193 //=======================================================================
 194 //function : LoadInitSurface
 195 //purpose  :
 196 //======================================================================
 197 void BRepFill_Filling::LoadInitSurface( const TopoDS_Face& aFace )
 198 {
 199   myInitFace = aFace;
 200   myIsInitFaceGiven = Standard_True;
 201 }
 202
 203 //=======================================================================
 204 //function : Add
 205 //purpose  : adds an edge as a constraint
 206 //======================================================================
 207 Standard_Integer BRepFill_Filling::Add( const TopoDS_Edge& anEdge,
 208                                         const GeomAbs_Shape Order,
 209                                         const Standard_Boolean IsBound )
 210 {
 211   TopoDS_Face NullFace;
 212   BRepFill_EdgeFaceAndOrder EdgeFaceAndOrder( anEdge, NullFace, Order );
 213   if (IsBound)
 214     {
 215       myBoundary.Append( EdgeFaceAndOrder );
 216       TopTools_ListOfShape EmptyList;
 217       myOldNewMap.Bind(anEdge, EmptyList);
 218       return myBoundary.Length();
 219     }
 220   else
 221     {
 222       myConstraints.Append( EdgeFaceAndOrder );
 223       return (myBoundary.Length() + myFreeConstraints.Length() + myConstraints.Length());
 224     }
 225 }
 226
 227 //=======================================================================
 228 //function : Add
 229 //purpose  : adds an edge with supporting face as a constraint
 230 //======================================================================
 231 Standard_Integer BRepFill_Filling::Add( const TopoDS_Edge& anEdge,
 232                                         const TopoDS_Face& Support,
 233                                         const GeomAbs_Shape Order,
 234                                         const Standard_Boolean IsBound )
 235 {
 236   BRepFill_EdgeFaceAndOrder EdgeFaceAndOrder( anEdge, Support, Order );
 237   if (IsBound)
 238     {
 239       myBoundary.Append( EdgeFaceAndOrder );
 240       TopTools_ListOfShape EmptyList;
 241       myOldNewMap.Bind(anEdge, EmptyList);
 242       return myBoundary.Length();
 243     }
 244   else
 245     {
 246       myConstraints.Append( EdgeFaceAndOrder );
 247       return (myBoundary.Length() + myFreeConstraints.Length() + myConstraints.Length());
 248     }
 249 }
 250
 251 //=======================================================================
 252 //function : Add
 253 //purpose  : adds a "free constraint": face without edge
 254 //======================================================================
 255 Standard_Integer BRepFill_Filling::Add( const TopoDS_Face& Support,
 256                                         const GeomAbs_Shape Order )
 257 {
 258   BRepFill_FaceAndOrder FaceAndOrder( Support, Order );
 259   myFreeConstraints.Append( FaceAndOrder );
 260   return (myBoundary.Length() + myFreeConstraints.Length());
 261 }
 262
 263 //=======================================================================
 264 //function : Add
 265 //purpose  : adds a point constraint
 266 //======================================================================
 267 Standard_Integer BRepFill_Filling::Add( const gp_Pnt& Point )
 268 {
 269   Handle( GeomPlate_PointConstraint ) aPC = new GeomPlate_PointConstraint( Point, GeomAbs_C0, myTol3d );
 270   myPoints.Append( aPC );
 271   return (myBoundary.Length() + myFreeConstraints.Length() + myConstraints.Length() + myPoints.Length());
 272 }
 273
 274 //=======================================================================
 275 //function : Add
 276 //purpose  : adds a point constraint on a face
 277 //======================================================================
 278 Standard_Integer BRepFill_Filling::Add( const Standard_Real U,
 279                                         const Standard_Real V,
 280                                         const TopoDS_Face& Support,
 281                                         const GeomAbs_Shape Order )
 282 {
 283   Handle( BRepAdaptor_HSurface ) HSurf = new BRepAdaptor_HSurface();
 284   HSurf->ChangeSurface().Initialize( Support );
 285   Handle( GeomPlate_PointConstraint ) aPC =
 286     new GeomPlate_PointConstraint( U, V, BRep_Tool::Surface( HSurf->ChangeSurface().Face() ), Order,
 287                                    myTol3d, myTolAng, myTolCurv );
 288   myPoints.Append( aPC );
 289   return (myBoundary.Length() + myFreeConstraints.Length() + myConstraints.Length() + myPoints.Length());
 290 }
 291
 292
 293 //=======================================================================
 294 //function : AddConstraints
 295 //purpose  :
 296 //======================================================================
 297 void BRepFill_Filling::AddConstraints( const BRepFill_SequenceOfEdgeFaceAndOrder& SeqOfConstraints )
 298 {
 299   TopoDS_Edge CurEdge;
 300   TopoDS_Face CurFace;
 301   GeomAbs_Shape CurOrder;
 302
 303   Handle(GeomPlate_CurveConstraint) Constr;
 304   Standard_Integer i;
 305   for (i = 1; i <= SeqOfConstraints.Length(); i++)
 306     {
 307       CurEdge = SeqOfConstraints(i).myEdge;
 308       CurFace = SeqOfConstraints(i).myFace;
 309       CurOrder = SeqOfConstraints(i).myOrder;
 310
 311       if (CurFace.IsNull()) {
 312         if (CurOrder == GeomAbs_C0) {
 313           Handle( BRepAdaptor_HCurve ) HCurve = new BRepAdaptor_HCurve();
 314           HCurve->ChangeCurve().Initialize( CurEdge );
 315
 316           Constr = new BRepFill_CurveConstraint(HCurve,
 317                                                 CurOrder,
 318                                                 myNbPtsOnCur,
 319                                                 myTol3d );
 320         }
 321         else { // Pas de representation Topologique
 322           // On prend une representation Geometrique : au pif !
 323           Handle( Geom_Surface ) Surface;
 324           Handle( Geom2d_Curve ) C2d;
 325           TopLoc_Location loc;
 326           Standard_Real f, l;
 327           BRep_Tool::CurveOnSurface( CurEdge, C2d, Surface, loc, f, l);
 328           if (Surface.IsNull()) {
 329             Standard_Failure::Raise( "Add" );
 330             return;
 331           }
 332           Surface = Handle(Geom_Surface)::DownCast(Surface->Copy());
 333           Surface->Transform(loc.Transformation());
 334           Handle( GeomAdaptor_HSurface ) Surf = new GeomAdaptor_HSurface(Surface);
 335           Handle( Geom2dAdaptor_HCurve ) Curve2d = new Geom2dAdaptor_HCurve(C2d);
 336
 337           Adaptor3d_CurveOnSurface CurvOnSurf( Curve2d, Surf );
 338           Handle (Adaptor3d_HCurveOnSurface) HCurvOnSurf = new Adaptor3d_HCurveOnSurface( CurvOnSurf );
 339
 340           Constr = new GeomPlate_CurveConstraint(HCurvOnSurf,
 341                                                  CurOrder,
 342                                                  myNbPtsOnCur,
 343                                                  myTol3d,
 344                                                  myTolAng,
 345                                                  myTolCurv );
 346         }
 347       }
 348       else
 349         {
 350           Handle( BRepAdaptor_HSurface ) Surf = new BRepAdaptor_HSurface();
 351           Surf->ChangeSurface().Initialize( CurFace );
 352           Handle( BRepAdaptor_HCurve2d ) Curve2d = new BRepAdaptor_HCurve2d();
 353           Curve2d->ChangeCurve2d().Initialize( CurEdge, CurFace );
 354           // If CurEdge has no 2d representation on CurFace,
 355           // there will be exception "Attempt to access to null object"
 356           // in this initialization (null pcurve).
 357           Adaptor3d_CurveOnSurface CurvOnSurf( Curve2d, Surf );
 358           Handle (Adaptor3d_HCurveOnSurface) HCurvOnSurf = new Adaptor3d_HCurveOnSurface( CurvOnSurf );
 359
 360           Constr = new BRepFill_CurveConstraint( HCurvOnSurf,
 361                                                  CurOrder,
 362                                                  myNbPtsOnCur,
 363                                                  myTol3d,
 364                                                  myTolAng,
 365                                                  myTolCurv );
 366         }
 367       if (myIsInitFaceGiven)
 368         {
 369           Handle( Geom2d_Curve ) Curve2d;
 370           Standard_Real FirstPar, LastPar;
 371           Curve2d = BRep_Tool::CurveOnSurface( CurEdge, myInitFace, FirstPar, LastPar );
 372           if (! Curve2d.IsNull())
 373             {
 374               Curve2d = new Geom2d_TrimmedCurve( Curve2d, FirstPar, LastPar );
 375               Constr->SetCurve2dOnSurf( Curve2d );
 376             }
 377         }
 378       myBuilder.Add( Constr );
 379     }
 380 }
 381
 382 //=======================================================================
 383 //function : BuildWires
 384 //purpose  :
 385 //======================================================================
 386 void BRepFill_Filling::BuildWires( TopTools_ListOfShape& EdgeList, TopTools_ListOfShape& WireList )
 387 {
 388   TopoDS_Wire CurWire;
 389   TopoDS_Edge CurEdge;
 390   TopTools_ListIteratorOfListOfShape Itl;
 391   Standard_Integer i, j;
 392
 393   while (! EdgeList.IsEmpty())
 394     {
 395       BRepLib_MakeWire MW;
 396       TopoDS_Edge FirstEdge = TopoDS::Edge(EdgeList.First());
 397       MW.Add(FirstEdge);
 398       EdgeList.RemoveFirst();
 399       TopoDS_Vertex V_wire[2], V_edge[2];
 400
 401       for (;;)
 402       {
 403         TopoDS_Wire CurWire = MW.Wire();
 404         TopExp::Vertices(CurWire, V_wire[0], V_wire[1]);
 405         Standard_Boolean found = Standard_False;
 406         for (Itl.Initialize( EdgeList ); Itl.More(); Itl.Next())
 407         {
 408           TopoDS_Edge CurEdge = TopoDS::Edge(Itl.Value());
 409           TopExp::Vertices(CurEdge, V_edge[0], V_edge[1]);
 410           for (i = 0; i < 2; i++)
 411           {
 412             for (j = 0; j < 2; j++)
 413               if (V_wire[i].IsSame(V_edge[j]))
 414               {
 415                 MW.Add(CurEdge);
 416                 EdgeList.Remove(Itl);
 417                 found = Standard_True;
 418                 break;
 419               }
 420             if (found)
 421               break;
 422           }
 423           if (found)
 424             break;
 425         }
 426         if (!found) //try to find geometric coincidence
 427         {
 428           gp_Pnt P_wire[2];
 429           P_wire[0] = BRep_Tool::Pnt(V_wire[0]);
 430           P_wire[1] = BRep_Tool::Pnt(V_wire[1]);
 431           for (Itl.Initialize( EdgeList ); Itl.More(); Itl.Next())
 432           {
 433             TopoDS_Edge CurEdge = TopoDS::Edge(Itl.Value());
 434             TopExp::Vertices(CurEdge, V_edge[0], V_edge[1]);
 435             for (i = 0; i < 2; i++)
 436             {
 437               for (j = 0; j < 2; j++)
 438               {
 439                 Standard_Real aDist = P_wire[i].Distance(BRep_Tool::Pnt(V_edge[j]));
 440                 if (aDist < BRep_Tool::Tolerance(V_wire[i]) &&
 441                     aDist < BRep_Tool::Tolerance(V_edge[j]))
 442                 {
 443                   MW.Add(CurEdge);
 444                   myOldNewMap(CurEdge).Append(MW.Edge());
 445                   EdgeList.Remove(Itl);
 446                   found = Standard_True;
 447                   break;
 448                 }
 449               }
 450               if (found)
 451                 break;
 452             }
 453             if (found)
 454               break;
 455           }
 456         }
 457         if (!found) //end of current wire, begin next wire
 458         {
 459           WireList.Append( MW.Wire() );
 460           break;
 461         }
 462       } //end of for (;;)
 463     } //end of while (! EdgeList.IsEmpty())
 464 }
 465
 466 //=======================================================================
 467 //function : FindExtremitiesOfHoles
 468 //purpose  :
 469 //======================================================================
 470 void BRepFill_Filling::FindExtremitiesOfHoles(const TopTools_ListOfShape& WireList,
 471                                               TopTools_SequenceOfShape& VerSeq ) const
 472 {
 473   TopTools_SequenceOfShape WireSeq;
 474   TopTools_ListIteratorOfListOfShape Itl(WireList);
 475   for (; Itl.More(); Itl.Next())
 476     WireSeq.Append(Itl.Value());
 477
 478   TopoDS_Wire theWire;
 479   theWire = TopoDS::Wire(WireSeq(1));
 480   WireSeq.Remove(1);
 481
 482   if (theWire.Closed())
 483     return;
 484
 485   TopoDS_Vertex Vfirst, Vlast;
 486   TopExp::Vertices( theWire, Vfirst, Vlast );
 487   gp_Vec FinVec = MakeFinVec( theWire, Vlast );
 488   TopoDS_Vertex theVertex = Vlast;
 489   VerSeq.Append( Vlast );
 490
 491   while (! WireSeq.IsEmpty())
 492     {
 493       TopoDS_Vertex MinVtx;
 494       Standard_Integer i, MinInd = 1;
 495       Standard_Boolean IsLast = Standard_False;
 496       Standard_Real MinAngle = M_PI;
 497
 498       for (i = 1; i <= WireSeq.Length(); i++)
 499         {
 500           const TopoDS_Wire& CurWire = TopoDS::Wire( WireSeq(i) );
 501           TopExp::Vertices( CurWire, Vfirst, Vlast );
 502
 503           Standard_Real angle =
 504             FinVec.Angle(gp_Vec(BRep_Tool::Pnt(theVertex), BRep_Tool::Pnt(Vfirst)));
 505           if (angle < MinAngle)
 506             {
 507               MinAngle = angle;
 508               MinVtx = Vfirst;
 509               MinInd = i;
 510               IsLast = Standard_True;
 511             }
 512           angle = FinVec.Angle(gp_Vec(BRep_Tool::Pnt(theVertex), BRep_Tool::Pnt(Vlast)));
 513           if (angle < MinAngle)
 514             {
 515               MinAngle = angle;
 516               MinVtx = Vlast;
 517               MinInd = i;
 518               IsLast = Standard_False;
 519             }
 520         }
 521       VerSeq.Append( MinVtx );
 522       const TopoDS_Wire& MinWire = TopoDS::Wire(WireSeq(MinInd));
 523       TopExp::Vertices( MinWire, Vfirst, Vlast );
 524       if (IsLast)
 525         {
 526           FinVec = MakeFinVec( MinWire, Vlast );
 527           theVertex = Vlast;
 528         }
 529       else
 530         {
 531           FinVec = MakeFinVec( MinWire, Vfirst );
 532           theVertex = Vfirst;
 533         }
 534       VerSeq.Append( theVertex );
 535       WireSeq.Remove(MinInd);
 536     }
 537   TopExp::Vertices( theWire, Vfirst, Vlast );
 538   VerSeq.Append( Vfirst );
 539 }
 540
 541 //=======================================================================
 542 //function : Build
 543 //purpose  : builds the resulting face
 544 //======================================================================
 545 void BRepFill_Filling::Build()
 546 {
 547   GeomPlate_BuildPlateSurface thebuild( myDegree, myNbPtsOnCur, myNbIter,
 548                                         myTol2d, myTol3d, myTolAng, myTolCurv, myAnisotropie );
 549
 550   myBuilder = thebuild;
 551   TopoDS_Edge CurEdge;
 552   TopoDS_Face CurFace;
 553   Standard_Integer i, j;
 554
 555   //Creating array of vertices: extremities of wires
 556   TopTools_SequenceOfShape VerSeq;
 557
 558   //Building missing bounds
 559   TopTools_ListOfShape EdgeList, WireList;
 560   for (i = 1; i <= myBoundary.Length(); i++)
 561     EdgeList.Append( myBoundary(i).myEdge );
 562
 563   BuildWires( EdgeList, WireList );
 564   FindExtremitiesOfHoles( WireList, VerSeq );
 565
 566   //Searching for surfaces for missing bounds
 567   for (j = 1; j <= myFreeConstraints.Length(); j++)
 568     {
 569       GeomAPI_ProjectPointOnSurf Projector;
 570       Quantity_Parameter U1, V1, U2, V2;
 571
 572       CurFace = myFreeConstraints(j).myFace;
 573       Handle( BRepAdaptor_HSurface ) HSurf = new BRepAdaptor_HSurface();
 574       HSurf->ChangeSurface().Initialize( CurFace );
 575       Handle( Geom_Surface ) CurSurface = BRep_Tool::Surface( HSurf->ChangeSurface().Face() );
 576       //BRepTopAdaptor_FClass2d Classifier( CurFace, Precision::Confusion() );
 577
 578       for (i = 1; i <= VerSeq.Length(); i += 2)
 579         {
 580           const TopoDS_Vertex& FirstVtx = TopoDS::Vertex(VerSeq(i));
 581           const TopoDS_Vertex& LastVtx  = TopoDS::Vertex(VerSeq(i+1));
 582
 583           gp_Pnt FirstPnt = BRep_Tool::Pnt(FirstVtx);
 584           Projector.Init( FirstPnt, CurSurface );
 585           if (Projector.LowerDistance() > Precision::Confusion())
 586             continue;
 587           Projector.LowerDistanceParameters( U1, V1 );
 588
 589           /*
 590           State = Classifier.Perform(gp_Pnt2d( U1, V1 ),
 591                                      ((HSurf->IsUPeriodic() || HSurf->IsVPeriodic())? Standard_True : Standard_False));
 592           if (State == TopAbs_OUT || State == TopAbs_UNKNOWN)
 593             continue;
 594           */
 595
 596           gp_Pnt LastPnt = BRep_Tool::Pnt(LastVtx);
 597           Projector.Init( LastPnt, CurSurface );
 598           if (Projector.LowerDistance() > Precision::Confusion())
 599             continue;
 600           Projector.LowerDistanceParameters( U2, V2 );
 601
 602           /*
 603           State = Classifier.Perform(gp_Pnt2d( U2, V2 ),
 604                                      ((HSurf->IsUPeriodic() || HSurf->IsVPeriodic())? Standard_True : Standard_False));
 605           if (State == TopAbs_OUT || State == TopAbs_UNKNOWN)
 606             continue;
 607           */
 608
 609           //Making the constraint
 610           TColgp_Array1OfPnt2d Points( 1, 2 );
 611           Points(1) = gp_Pnt2d( U1, V1 );
 612           Points(2) = gp_Pnt2d( U2, V2 );
 613           Handle( Geom2d_BezierCurve ) Line2d = new Geom2d_BezierCurve( Points );
 614           TopoDS_Edge E = BRepLib_MakeEdge( Line2d, CurSurface, FirstVtx, LastVtx );
 615           Add( E, CurFace, myFreeConstraints(j).myOrder );
 616           VerSeq.Remove( i, i+1 );
 617           break;
 618         } //for (i = 1; i <= VerSeq.Length(); i += 2)
 619     } //for (j = 1; j <= myFreeConstraints.Length(); j++)
 620
 621   //Load initial surface to myBuilder if it is given
 622   if (myIsInitFaceGiven)
 623     {
 624       Handle( BRepAdaptor_HSurface ) HSurfInit = new BRepAdaptor_HSurface();
 625       HSurfInit->ChangeSurface().Initialize( myInitFace );
 626       myBuilder.LoadInitSurface( BRep_Tool::Surface( HSurfInit->ChangeSurface().Face() ) );
 627     }
 628
 629   //Adding constraints to myBuilder
 630   AddConstraints( myBoundary );
 631   myBuilder.SetNbBounds( myBoundary.Length() );
 632   AddConstraints( myConstraints );
 633   for (i = 1; i <= myPoints.Length(); i++)
 634     myBuilder.Add( myPoints(i) );
 635
 636   myBuilder.Perform();
 637   if (myBuilder.IsDone())
 638     myIsDone = Standard_True;
 639   else
 640     {
 641       myIsDone = Standard_False;
 642       return;
 643     }
 644
 645   Handle( GeomPlate_Surface ) GPlate = myBuilder.Surface();
 646   Handle( Geom_BSplineSurface ) Surface;
 647   // Approximation
 648   Standard_Real dmax = 1.1 * myBuilder.G0Error(); //???????????
 649   //Standard_Real dmax = myTol3d;
 650   if (! myIsInitFaceGiven)
 651     {
 652      Standard_Real seuil; //?????
 653
 654      TColgp_SequenceOfXY S2d;
 655      TColgp_SequenceOfXYZ S3d;
 656      myBuilder.Disc2dContour(4,S2d);
 657      myBuilder.Disc3dContour(4,0,S3d);
 658      seuil = Max( myTol3d, 10*myBuilder.G0Error() ); //????????
 659      GeomPlate_PlateG0Criterion Criterion( S2d, S3d, seuil );
 660      GeomPlate_MakeApprox Approx( GPlate, Criterion, myTol3d, myMaxSegments, myMaxDeg );
 661      Surface = Approx.Surface();
 662    }
 663   else
 664     {
 665       GeomPlate_MakeApprox Approx( GPlate, myTol3d, myMaxSegments, myMaxDeg, dmax, 0 ); //?????????????
 666       //GeomConvert_ApproxSurface Approx( GPlate, myTol3d, GeomAbs_C1, GeomAbs_C1, myMaxDeg, myMaxDeg, myMaxSegments, 1 );
 667       //Approx.Dump( cout );
 668       Surface = Approx.Surface();
 669     }
 670
 671   //Build the final wire and final face
 672   TopTools_ListOfShape FinalEdges;
 673   Handle(TColGeom2d_HArray1OfCurve) CurvesOnPlate = myBuilder.Curves2d();
 674   BRep_Builder BB;
 675   for (i = 1; i <= myBoundary.Length(); i++)
 676   {
 677     const TopoDS_Edge& InitEdge = myBoundary(i).myEdge;
 678     TopoDS_Edge anEdge = InitEdge;
 679     if (!myOldNewMap(anEdge).IsEmpty())
 680       anEdge = TopoDS::Edge( myOldNewMap(anEdge).First() );
 681     Handle(Geom2d_Curve) aCurveOnPlate = CurvesOnPlate->Value(i);
 682
 683     TopoDS_Edge NewEdge = TopoDS::Edge(anEdge.EmptyCopied());
 684
 685     TopoDS_Vertex V1, V2;
 686     TopExp::Vertices(anEdge, V1, V2, Standard_True); //with orientation
 687     BB.UpdateVertex(V1, dmax);
 688     BB.UpdateVertex(V2, dmax);
 689     BB.Add(NewEdge, V1);
 690     BB.Add(NewEdge, V2);
 691     TopLoc_Location Loc;
 692     BB.UpdateEdge(NewEdge, aCurveOnPlate, Surface, Loc, dmax);
 693     //BRepLib::SameRange(NewEdge);
 694     BRepLib::SameParameter(NewEdge, dmax, Standard_True);
 695     FinalEdges.Append(NewEdge);
 696     myOldNewMap(InitEdge).Clear();
 697     myOldNewMap(InitEdge).Append(NewEdge);
 698   }
 699
 700   TopoDS_Wire FinalWire = WireFromList(FinalEdges);
 701   if (!(FinalWire.Closed()))
 702     Standard_Failure::Raise("Wire is not closed");
 703
 704   myFace = BRepLib_MakeFace( Surface, FinalWire );
 705 }
 706
 707 //=======================================================================
 708 //function : IsDone
 709 //purpose  :
 710 //======================================================================
 711 Standard_Boolean BRepFill_Filling::IsDone() const
 712 {
 713   return myIsDone;
 714 }
 715
 716 //=======================================================================
 717 //function : Face
 718 //purpose  : returns the result
 719 //======================================================================
 720 TopoDS_Face BRepFill_Filling::Face() const
 721 {
 722   return myFace;
 723 }
 724
 725 //=======================================================================
 726 //function : Generated
 727 //purpose  : returns the new edge (first in list) made from old edge "S"
 728 //=======================================================================
 729  const TopTools_ListOfShape& BRepFill_Filling::Generated(const TopoDS_Shape& S)
 730 {
 731   myGenerated.Clear();
 732
 733   if (myOldNewMap.IsBound(S))
 734     myGenerated.Append(myOldNewMap(S));
 735
 736   return myGenerated;
 737 }
 738
 739 //==========================================================================
 740 //function : G0Error
 741 //purpose  : returns maximum distance from boundary to the resulting surface
 742 //==========================================================================
 743 Standard_Real BRepFill_Filling::G0Error() const
 744 {
 745   return myBuilder.G0Error();
 746 }
 747
 748 //=======================================================================
 749 //function : G1Error
 750 //purpose  : returns maximum angle between the resulting surface
 751 //           and constraint surfaces at boundaries
 752 //======================================================================
 753 Standard_Real BRepFill_Filling::G1Error() const
 754 {
 755   return myBuilder.G1Error();
 756 }
 757
 758 //=======================================================================
 759 //function : G2Error
 760 //purpose  : returns maximum difference of curvature between
 761 //           the resulting surface and constraint surfaces at boundaries
 762 //======================================================================
 763 Standard_Real BRepFill_Filling::G2Error() const
 764 {
 765   return myBuilder.G2Error();
 766 }
 767
 768 //==========================================================================
 769 //function : G0Error
 770 //purpose  : returns maximum distance between the constraint number Index
 771 //           and the resulting surface
 772 //==========================================================================
 773 Standard_Real BRepFill_Filling::G0Error( const Standard_Integer Index )
 774 {
 775   return myBuilder.G0Error( Index );
 776 }
 777
 778 //==========================================================================
 779 //function : G1Error
 780 //purpose  : returns maximum angle between the constraint number Index
 781 //           and the resulting surface
 782 //==========================================================================
 783 Standard_Real BRepFill_Filling::G1Error( const Standard_Integer Index )
 784 {
 785   return myBuilder.G1Error( Index );
 786 }
 787
 788 //==========================================================================
 789 //function : G2Error
 790 //purpose  : returns maximum difference of curvature between
 791 //           the constraint number Index and the resulting surface
 792 //==========================================================================
 793 Standard_Real BRepFill_Filling::G2Error( const Standard_Integer Index )
 794 {
 795   return myBuilder.G2Error( Index );
 796 }