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