projects / occt.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
0025656: Specification of semantic of Closed flag of an edge
[occt.git] / src / BRepFill / BRepFill_Filling.cxx
CommitLineData
b311480e 1// Created on: 1998-08-26
2// Created by: Julia GERASIMOVA
3// Copyright (c) 1998-1999 Matra Datavision
973c2be1 4// Copyright (c) 1999-2014 OPEN CASCADE SAS
b311480e 5//
973c2be1 6// This file is part of Open CASCADE Technology software library.
b311480e 7//
d5f74e42 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
973c2be1 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.
b311480e 13//
973c2be1 14// Alternatively, this file may be used under the terms of Open CASCADE
15// commercial license or contractual agreement.
7fd59977 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>
01697018 56#include <TopTools_ListIteratorOfListOfShape.hxx>
7fd59977 57#include <BRepTools_WireExplorer.hxx>
58#include <TColgp_Array1OfPnt2d.hxx>
01697018 59#include <TColGeom2d_HArray1OfCurve.hxx>
7fd59977 60
61#include <Geom2d_BezierCurve.hxx>
62#include <Geom2d_TrimmedCurve.hxx>
63#include <Geom_Surface.hxx>
64
65#include <TColgp_SequenceOfXY.hxx>
66#include <GeomPlate_PlateG0Criterion.hxx>
67#include <BRep_TVertex.hxx>
68
69static 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
01697018
J		 90static 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//======================================================================
7fd59977 123BRepFill_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
01697018
J		 152//=======================================================================
153//function : SetConstrParam
154//purpose :
155//======================================================================
7fd59977 156void 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
01697018
J		 167//=======================================================================
168//function : SetResolParam
169//purpose :
170//======================================================================
7fd59977 171void 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
01697018
J		 182//=======================================================================
183//function : SetApproxParam
184//purpose :
185//======================================================================
7fd59977 186void BRepFill_Filling::SetApproxParam( const Standard_Integer MaxDeg,
187 const Standard_Integer MaxSegments )
188{
189 myMaxDeg = MaxDeg;
190 myMaxSegments = MaxSegments;
191}
192
01697018
J		 193//=======================================================================
194//function : LoadInitSurface
195//purpose :
196//======================================================================
7fd59977 197void BRepFill_Filling::LoadInitSurface( const TopoDS_Face& aFace )
198{
199 myInitFace = aFace;
200 myIsInitFaceGiven = Standard_True;
201}
202
01697018
J		 203//=======================================================================
204//function : Add
205//purpose : adds an edge as a constraint
206//======================================================================
7fd59977 207Standard_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 );
01697018
J		 216 TopTools_ListOfShape EmptyList;
217 myOldNewMap.Bind(anEdge, EmptyList);
7fd59977 218 return myBoundary.Length();
219 }
220 else
221 {
222 myConstraints.Append( EdgeFaceAndOrder );
223 return (myBoundary.Length() + myFreeConstraints.Length() + myConstraints.Length());
224 }
225}
226
01697018
J		 227//=======================================================================
228//function : Add
229//purpose : adds an edge with supporting face as a constraint
230//======================================================================
7fd59977 231Standard_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 );
01697018
J		 240 TopTools_ListOfShape EmptyList;
241 myOldNewMap.Bind(anEdge, EmptyList);
7fd59977 242 return myBoundary.Length();
243 }
244 else
245 {
246 myConstraints.Append( EdgeFaceAndOrder );
247 return (myBoundary.Length() + myFreeConstraints.Length() + myConstraints.Length());
248 }
249}
250
01697018
J		 251//=======================================================================
252//function : Add
253//purpose : adds a "free constraint": face without edge
254//======================================================================
7fd59977 255Standard_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
01697018
J		 263//=======================================================================
264//function : Add
265//purpose : adds a point constraint
266//======================================================================
7fd59977 267Standard_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
01697018
J		 274//=======================================================================
275//function : Add
276//purpose : adds a point constraint on a face
277//======================================================================
7fd59977 278Standard_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
01697018
J		 293//=======================================================================
294//function : AddConstraints
295//purpose :
296//======================================================================
7fd59977 297void 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;
01697018
J		 304 Standard_Integer i;
305 for (i = 1; i <= SeqOfConstraints.Length(); i++)
7fd59977 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;
7fd59977 325 TopLoc_Location loc;
326 Standard_Real f, l;
7fd59977 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
01697018
J		 382//=======================================================================
383//function : BuildWires
384//purpose :
385//======================================================================
386void BRepFill_Filling::BuildWires( TopTools_ListOfShape& EdgeList, TopTools_ListOfShape& WireList )
7fd59977 387{
388 TopoDS_Wire CurWire;
389 TopoDS_Edge CurEdge;
01697018
J		 390 TopTools_ListIteratorOfListOfShape Itl;
391 Standard_Integer i, j;
7fd59977 392
01697018 393 while (! EdgeList.IsEmpty())
7fd59977 394 {
395 BRepLib_MakeWire MW;
01697018
J		 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())
7fd59977 464}
465
01697018
J		 466//=======================================================================
467//function : FindExtremitiesOfHoles
468//purpose :
469//======================================================================
470void BRepFill_Filling::FindExtremitiesOfHoles(const TopTools_ListOfShape& WireList,
471 TopTools_SequenceOfShape& VerSeq ) const
7fd59977 472{
01697018
J		 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);
7fd59977 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 );
01697018
J		 488 TopoDS_Vertex theVertex = Vlast;
489 VerSeq.Append( Vlast );
7fd59977 490
01697018 491 while (! WireSeq.IsEmpty())
7fd59977 492 {
01697018 493 TopoDS_Vertex MinVtx;
1d47d8d0 494 Standard_Integer i, MinInd = 1;
7fd59977 495 Standard_Boolean IsLast = Standard_False;
c6541a0c 496 Standard_Real MinAngle = M_PI;
7fd59977 497
01697018 498 for (i = 1; i <= WireSeq.Length(); i++)
7fd59977 499 {
01697018 500 const TopoDS_Wire& CurWire = TopoDS::Wire( WireSeq(i) );
7fd59977 501 TopExp::Vertices( CurWire, Vfirst, Vlast );
502
01697018
J		 503 Standard_Real angle =
504 FinVec.Angle(gp_Vec(BRep_Tool::Pnt(theVertex), BRep_Tool::Pnt(Vfirst)));
7fd59977 505 if (angle < MinAngle)
506 {
507 MinAngle = angle;
01697018
J		 508 MinVtx = Vfirst;
509 MinInd = i;
7fd59977 510 IsLast = Standard_True;
511 }
01697018 512 angle = FinVec.Angle(gp_Vec(BRep_Tool::Pnt(theVertex), BRep_Tool::Pnt(Vlast)));
7fd59977 513 if (angle < MinAngle)
514 {
515 MinAngle = angle;
01697018
J		 516 MinVtx = Vlast;
517 MinInd = i;
7fd59977 518 IsLast = Standard_False;
519 }
520 }
01697018
J		 521 VerSeq.Append( MinVtx );
522 const TopoDS_Wire& MinWire = TopoDS::Wire(WireSeq(MinInd));
7fd59977 523 TopExp::Vertices( MinWire, Vfirst, Vlast );
524 if (IsLast)
525 {
526 FinVec = MakeFinVec( MinWire, Vlast );
01697018 527 theVertex = Vlast;
7fd59977 528 }
529 else
530 {
531 FinVec = MakeFinVec( MinWire, Vfirst );
01697018 532 theVertex = Vfirst;
7fd59977 533 }
01697018
J		 534 VerSeq.Append( theVertex );
535 WireSeq.Remove(MinInd);
7fd59977 536 }
537 TopExp::Vertices( theWire, Vfirst, Vlast );
01697018 538 VerSeq.Append( Vfirst );
7fd59977 539}
540
01697018
J		 541//=======================================================================
542//function : Build
543//purpose : builds the resulting face
544//======================================================================
7fd59977 545void BRepFill_Filling::Build()
546{
547 GeomPlate_BuildPlateSurface thebuild( myDegree, myNbPtsOnCur, myNbIter,
01697018 548 myTol2d, myTol3d, myTolAng, myTolCurv, myAnisotropie );
7fd59977 549
550 myBuilder = thebuild;
551 TopoDS_Edge CurEdge;
552 TopoDS_Face CurFace;
553 Standard_Integer i, j;
554
01697018
J		 555 //Creating array of vertices: extremities of wires
556 TopTools_SequenceOfShape VerSeq;
7fd59977 557
558 //Building missing bounds
01697018 559 TopTools_ListOfShape EdgeList, WireList;
7fd59977 560 for (i = 1; i <= myBoundary.Length(); i++)
01697018 561 EdgeList.Append( myBoundary(i).myEdge );
7fd59977 562
01697018
J		 563 BuildWires( EdgeList, WireList );
564 FindExtremitiesOfHoles( WireList, VerSeq );
7fd59977 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
01697018 578 for (i = 1; i <= VerSeq.Length(); i += 2)
7fd59977 579 {
01697018
J		 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 );
7fd59977 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 */
01697018
J		 595
596 gp_Pnt LastPnt = BRep_Tool::Pnt(LastVtx);
597 Projector.Init( LastPnt, CurSurface );
7fd59977 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 );
01697018 614 TopoDS_Edge E = BRepLib_MakeEdge( Line2d, CurSurface, FirstVtx, LastVtx );
7fd59977 615 Add( E, CurFace, myFreeConstraints(j).myOrder );
01697018 616 VerSeq.Remove( i, i+1 );
7fd59977 617 break;
01697018 618 } //for (i = 1; i <= VerSeq.Length(); i += 2)
7fd59977 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;
7fd59977 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 }
01697018
J		 670
671 //Build the final wire and final face
672 TopTools_ListOfShape FinalEdges;
673 Handle(TColGeom2d_HArray1OfCurve) CurvesOnPlate = myBuilder.Curves2d();
674 BRep_Builder BB;
7fd59977 675 for (i = 1; i <= myBoundary.Length(); i++)
01697018
J		 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());
da72a17c 684
01697018
J		 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 }
7fd59977 699
01697018
J		 700 TopoDS_Wire FinalWire = WireFromList(FinalEdges);
701 if (!(FinalWire.Closed()))
7fd59977 702 Standard_Failure::Raise("Wire is not closed");
01697018
J		 703
704 myFace = BRepLib_MakeFace( Surface, FinalWire );
7fd59977 705}
706
01697018
J		 707//=======================================================================
708//function : IsDone
709//purpose :
710//======================================================================
7fd59977 711Standard_Boolean BRepFill_Filling::IsDone() const
712{
713 return myIsDone;
714}
715
01697018
J		 716//=======================================================================
717//function : Face
718//purpose : returns the result
719//======================================================================
7fd59977 720TopoDS_Face BRepFill_Filling::Face() const
721{
722 return myFace;
723}
724
01697018
J		 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//==========================================================================
7fd59977 743Standard_Real BRepFill_Filling::G0Error() const
744{
745 return myBuilder.G0Error();
746}
747
01697018
J		 748//=======================================================================
749//function : G1Error
750//purpose : returns maximum angle between the resulting surface
751// and constraint surfaces at boundaries
752//======================================================================
7fd59977 753Standard_Real BRepFill_Filling::G1Error() const
754{
755 return myBuilder.G1Error();
756}
757
01697018
J		 758//=======================================================================
759//function : G2Error
760//purpose : returns maximum difference of curvature between
761// the resulting surface and constraint surfaces at boundaries
762//======================================================================
7fd59977 763Standard_Real BRepFill_Filling::G2Error() const
764{
765 return myBuilder.G2Error();
766}
767
01697018
J		 768//==========================================================================
769//function : G0Error
770//purpose : returns maximum distance between the constraint number Index
771// and the resulting surface
772//==========================================================================
7fd59977 773Standard_Real BRepFill_Filling::G0Error( const Standard_Integer Index )
774{
775 return myBuilder.G0Error( Index );
776}
777
01697018
J		 778//==========================================================================
779//function : G1Error
780//purpose : returns maximum angle between the constraint number Index
781// and the resulting surface
782//==========================================================================
7fd59977 783Standard_Real BRepFill_Filling::G1Error( const Standard_Integer Index )
784{
785 return myBuilder.G1Error( Index );
786}
787
01697018
J		 788//==========================================================================
789//function : G2Error
790//purpose : returns maximum difference of curvature between
791// the constraint number Index and the resulting surface
792//==========================================================================
7fd59977 793Standard_Real BRepFill_Filling::G2Error( const Standard_Integer Index )
794{
795 return myBuilder.G2Error( Index );
796}
Open CASCADE Technology repository