0031007: Coding - eliminate warnings issued while compiling with -pedantic flag
[occt.git] / src / BRepLib / BRepLib_MakeFace.cxx
1 // Created on: 1993-07-23
2 // Created by: Remi LEQUETTE
3 // Copyright (c) 1993-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 <BRep_Builder.hxx>
19 #include <BRepLib.hxx>
20 #include <BRepLib_FindSurface.hxx>
21 #include <BRepLib_MakeFace.hxx>
22 #include <BRepTopAdaptor_FClass2d.hxx>
23 #include <ElCLib.hxx>
24 #include <Geom2d_Line.hxx>
25 #include <Geom_BezierCurve.hxx>
26 #include <Geom_BSplineCurve.hxx>
27 #include <Geom_ConicalSurface.hxx>
28 #include <Geom_Curve.hxx>
29 #include <Geom_CylindricalSurface.hxx>
30 #include <Geom_OffsetSurface.hxx>
31 #include <Geom_Plane.hxx>
32 #include <Geom_RectangularTrimmedSurface.hxx>
33 #include <Geom_SphericalSurface.hxx>
34 #include <Geom_Surface.hxx>
35 #include <Geom_SurfaceOfLinearExtrusion.hxx>
36 #include <Geom_SurfaceOfRevolution.hxx>
37 #include <Geom_ToroidalSurface.hxx>
38 #include <GeomAbs_CurveType.hxx>
39 #include <GeomAdaptor_Curve.hxx>
40 #include <gp_Cone.hxx>
41 #include <gp_Cylinder.hxx>
42 #include <gp_Pln.hxx>
43 #include <gp_Sphere.hxx>
44 #include <gp_Torus.hxx>
45 #include <Precision.hxx>
46 #include <StdFail_NotDone.hxx>
47 #include <TopAbs_State.hxx>
48 #include <TopoDS.hxx>
49 #include <TopoDS_Edge.hxx>
50 #include <TopoDS_Face.hxx>
51 #include <TopoDS_Iterator.hxx>
52 #include <TopoDS_Vertex.hxx>
53 #include <TopoDS_Wire.hxx>
54
55 //=======================================================================
56 //function : BRepLib_MakeFace
57 //purpose  : 
58 //=======================================================================
59 BRepLib_MakeFace::BRepLib_MakeFace() :
60        myError(BRepLib_NoFace)
61 {
62 }
63
64
65 //=======================================================================
66 //function : BRepLib_MakeFace
67 //purpose  : 
68 //=======================================================================
69
70 BRepLib_MakeFace::BRepLib_MakeFace(const TopoDS_Face& F)
71 {
72   Init(F);
73 }
74
75
76 //=======================================================================
77 //function : BRepLib_MakeFace
78 //purpose  : 
79 //=======================================================================
80
81 BRepLib_MakeFace::BRepLib_MakeFace(const gp_Pln& P)
82 {
83   Handle(Geom_Plane) GP = new Geom_Plane(P);
84   Init(GP, Standard_True, Precision::Confusion());
85 }
86
87
88 //=======================================================================
89 //function : BRepLib_MakeFace
90 //purpose  : 
91 //=======================================================================
92
93 BRepLib_MakeFace::BRepLib_MakeFace(const gp_Cylinder& C)
94 {
95   Handle(Geom_CylindricalSurface) GC = new Geom_CylindricalSurface(C);
96   Init(GC, Standard_True, Precision::Confusion());
97 }
98
99
100 //=======================================================================
101 //function : BRepLib_MakeFace
102 //purpose  : 
103 //=======================================================================
104
105 BRepLib_MakeFace::BRepLib_MakeFace(const gp_Cone& C)
106 {
107   Handle(Geom_ConicalSurface) GC = new Geom_ConicalSurface(C);
108   Init(GC, Standard_True, Precision::Confusion());
109 }
110
111
112 //=======================================================================
113 //function : BRepLib_MakeFace
114 //purpose  : 
115 //=======================================================================
116
117 BRepLib_MakeFace::BRepLib_MakeFace(const gp_Sphere& S)
118 {
119   Handle(Geom_SphericalSurface) GS = new Geom_SphericalSurface(S);
120   Init(GS, Standard_True, Precision::Confusion());
121 }
122
123
124 //=======================================================================
125 //function : BRepLib_MakeFace
126 //purpose  : 
127 //=======================================================================
128
129 BRepLib_MakeFace::BRepLib_MakeFace(const gp_Torus& T)
130 {
131   Handle(Geom_ToroidalSurface) GT = new Geom_ToroidalSurface(T);
132   Init(GT, Standard_True, Precision::Confusion());
133 }
134
135
136 //=======================================================================
137 //function : BRepLib_MakeFace
138 //purpose  : 
139 //=======================================================================
140
141 BRepLib_MakeFace::BRepLib_MakeFace(const Handle(Geom_Surface)& S,
142                                    const Standard_Real TolDegen)
143 {
144   Init(S, Standard_True, TolDegen);
145 }
146
147
148 //=======================================================================
149 //function : BRepLib_MakeFace
150 //purpose  : 
151 //=======================================================================
152
153 BRepLib_MakeFace::BRepLib_MakeFace(const gp_Pln& P,
154                                    const Standard_Real UMin,
155                                    const Standard_Real UMax,
156                                    const Standard_Real VMin,
157                                    const Standard_Real VMax)
158 {
159   Handle(Geom_Plane) GP = new Geom_Plane(P);
160   Init(GP, UMin, UMax, VMin, VMax, Precision::Confusion());
161 }
162
163
164 //=======================================================================
165 //function : BRepLib_MakeFace
166 //purpose  : 
167 //=======================================================================
168
169 BRepLib_MakeFace::BRepLib_MakeFace(const gp_Cylinder& C,
170                                    const Standard_Real UMin,
171                                    const Standard_Real UMax,
172                                    const Standard_Real VMin,
173                                    const Standard_Real VMax)
174 {
175   Handle(Geom_CylindricalSurface) GC = new Geom_CylindricalSurface(C);
176   Init(GC, UMin, UMax, VMin, VMax, Precision::Confusion());
177 }
178
179
180 //=======================================================================
181 //function : BRepLib_MakeFace
182 //purpose  : 
183 //=======================================================================
184
185 BRepLib_MakeFace::BRepLib_MakeFace(const gp_Cone& C,
186                                    const Standard_Real UMin,
187                                    const Standard_Real UMax,
188                                    const Standard_Real VMin,
189                                    const Standard_Real VMax)
190 {
191   Handle(Geom_ConicalSurface) GC = new Geom_ConicalSurface(C);
192   Init(GC, UMin, UMax, VMin, VMax, Precision::Confusion());
193 }
194
195
196 //=======================================================================
197 //function : BRepLib_MakeFace
198 //purpose  : 
199 //=======================================================================
200
201 BRepLib_MakeFace::BRepLib_MakeFace(const gp_Sphere& S,
202                                    const Standard_Real UMin,
203                                    const Standard_Real UMax,
204                                    const Standard_Real VMin,
205                                    const Standard_Real VMax)
206 {
207   Handle(Geom_SphericalSurface) GS = new Geom_SphericalSurface(S);
208   Init(GS, UMin, UMax, VMin, VMax, Precision::Confusion());
209 }
210
211
212 //=======================================================================
213 //function : BRepLib_MakeFace
214 //purpose  : 
215 //=======================================================================
216
217 BRepLib_MakeFace::BRepLib_MakeFace(const gp_Torus& T,
218                                    const Standard_Real UMin,
219                                    const Standard_Real UMax,
220                                    const Standard_Real VMin,
221                                    const Standard_Real VMax)
222 {
223   Handle(Geom_ToroidalSurface) GT = new Geom_ToroidalSurface(T);
224   Init(GT, UMin, UMax, VMin, VMax, Precision::Confusion());
225 }
226
227
228 //=======================================================================
229 //function : BRepLib_MakeFace
230 //purpose  : 
231 //=======================================================================
232
233 BRepLib_MakeFace::BRepLib_MakeFace(const Handle(Geom_Surface)& S,
234                                    const Standard_Real UMin,
235                                    const Standard_Real UMax,
236                                    const Standard_Real VMin,
237                                    const Standard_Real VMax,
238                                    const Standard_Real TolDegen)
239 {
240   Init(S, UMin, UMax, VMin, VMax, TolDegen);
241 }
242
243
244 //=======================================================================
245 //function : BRepLib_MakeFace
246 //purpose  : 
247 //=======================================================================
248
249 BRepLib_MakeFace::BRepLib_MakeFace(const TopoDS_Wire& W,
250                                    const Standard_Boolean OnlyPlane)
251
252 {
253   // Find a surface through the wire
254   BRepLib_FindSurface FS(W, -1, OnlyPlane, Standard_True);
255   if (!FS.Found()) {
256     myError = BRepLib_NotPlanar;
257     return;
258   }
259  
260   // build the face and add the wire
261   BRep_Builder B;
262   myError = BRepLib_FaceDone;
263
264   Standard_Real tol = Max(1.2*FS.ToleranceReached(), FS.Tolerance());
265
266   B.MakeFace(TopoDS::Face(myShape),FS.Surface(),FS.Location(),tol);
267   Add(W);
268   //
269   BRepLib::UpdateTolerances(myShape);
270   //
271   BRepLib::SameParameter(myShape, tol, Standard_True);
272   //
273   if (BRep_Tool::IsClosed(W))
274     CheckInside();
275 }
276
277
278 //=======================================================================
279 //function : BRepLib_MakeFace
280 //purpose  : 
281 //=======================================================================
282
283 BRepLib_MakeFace::BRepLib_MakeFace(const gp_Pln& P,
284                                    const TopoDS_Wire& W,
285                                    const Standard_Boolean Inside)
286 {
287   Handle(Geom_Plane) Pl = new Geom_Plane(P);
288   Init(Pl, Standard_False, Precision::Confusion());
289   Add(W);
290   if (Inside && BRep_Tool::IsClosed(W)) CheckInside();
291 }
292
293
294 //=======================================================================
295 //function : BRepLib_MakeFace
296 //purpose  : 
297 //=======================================================================
298
299 BRepLib_MakeFace::BRepLib_MakeFace(const gp_Cylinder& C,
300                                    const TopoDS_Wire& W,
301                                    const Standard_Boolean Inside)
302 {
303   Handle(Geom_CylindricalSurface) GC = new Geom_CylindricalSurface(C);
304   Init(GC, Standard_False, Precision::Confusion());
305   Add(W);
306   if (Inside && BRep_Tool::IsClosed(W)) CheckInside();
307 }
308
309
310 //=======================================================================
311 //function : BRepLib_MakeFace
312 //purpose  : 
313 //=======================================================================
314
315 BRepLib_MakeFace::BRepLib_MakeFace(const gp_Cone& C,
316                                    const TopoDS_Wire& W,
317                                    const Standard_Boolean Inside)
318 {
319   Handle(Geom_ConicalSurface) GC = new Geom_ConicalSurface(C);
320   Init(GC, Standard_False, Precision::Confusion());
321   Add(W);
322   if (Inside && BRep_Tool::IsClosed(W)) CheckInside();
323 }
324
325
326 //=======================================================================
327 //function : BRepLib_MakeFace
328 //purpose  : 
329 //=======================================================================
330
331 BRepLib_MakeFace::BRepLib_MakeFace(const gp_Sphere& S,
332                                    const TopoDS_Wire& W,
333                                    const Standard_Boolean Inside)
334 {
335   Handle(Geom_SphericalSurface) GS = new Geom_SphericalSurface(S);
336   Init(GS, Standard_False, Precision::Confusion());
337   Add(W);
338   if (Inside && BRep_Tool::IsClosed(W)) CheckInside();
339 }
340
341
342 //=======================================================================
343 //function : BRepLib_MakeFace
344 //purpose  : 
345 //=======================================================================
346
347 BRepLib_MakeFace::BRepLib_MakeFace(const gp_Torus& T,
348                                    const TopoDS_Wire& W,
349                                    const Standard_Boolean Inside)
350 {
351   Handle(Geom_ToroidalSurface) GT = new Geom_ToroidalSurface(T);
352   Init(GT, Standard_False, Precision::Confusion());
353   Add(W);
354   if (Inside && BRep_Tool::IsClosed(W)) CheckInside();
355 }
356
357
358 //=======================================================================
359 //function : BRepLib_MakeFace
360 //purpose  : 
361 //=======================================================================
362
363 BRepLib_MakeFace::BRepLib_MakeFace(const Handle(Geom_Surface)& S,
364                                    const TopoDS_Wire& W,
365                                    const Standard_Boolean Inside)
366 {
367   Init(S, Standard_False, Precision::Confusion());
368   Add(W);
369   if (Inside && BRep_Tool::IsClosed(W)) CheckInside();
370 }
371
372
373 //=======================================================================
374 //function : BRepLib_MakeFace
375 //purpose  : 
376 //=======================================================================
377
378 BRepLib_MakeFace::BRepLib_MakeFace(const TopoDS_Face& F,
379                                    const TopoDS_Wire& W)
380 {
381   Init(F);
382   Add(W);
383 }
384
385 //=======================================================================
386 //function : Init
387 //purpose  : 
388 //=======================================================================
389
390 void BRepLib_MakeFace::Init(const TopoDS_Face& F)
391 {
392   // copy the face
393   myShape = F.EmptyCopied();
394   myError = BRepLib_FaceDone;
395   
396   BRep_Builder B;
397   TopoDS_Iterator It(F);
398   while (It.More()) {
399     B.Add(myShape,It.Value());
400     It.Next();
401   }
402 }
403
404 //=======================================================================
405 //function : Init
406 //purpose  : 
407 //=======================================================================
408
409 void  BRepLib_MakeFace::Init(const Handle(Geom_Surface)& S,
410                              const Standard_Boolean Bound,
411                              const Standard_Real TolDegen)
412 {
413   myError = BRepLib_FaceDone;
414   if (Bound) {
415     Standard_Real UMin,UMax,VMin,VMax;
416     S->Bounds(UMin,UMax,VMin,VMax);
417     Init(S,UMin,UMax,VMin,VMax,TolDegen);
418   }
419   else {
420     BRep_Builder B;
421     B.MakeFace(TopoDS::Face(myShape),S,Precision::Confusion());
422   }
423   BRep_Builder B;
424   B.NaturalRestriction(TopoDS::Face(myShape),Standard_True);
425 }
426
427
428 //=======================================================================
429 //function : IsDegenerated
430 //purpose  : Checks whether the passed curve is degenerated with the
431 //           passed tolerance value
432 //=======================================================================
433
434 Standard_Boolean BRepLib_MakeFace::IsDegenerated(
435     const Handle(Geom_Curve)& theCurve,
436     const Standard_Real       theMaxTol,
437           Standard_Real&      theActTol)
438 {
439   GeomAdaptor_Curve AC(theCurve);
440   Standard_Real aConfusion = Precision::Confusion();
441   theActTol = aConfusion;
442   GeomAbs_CurveType Type = AC.GetType();
443
444   if (Type == GeomAbs_Circle) {
445     gp_Circ Circ = AC.Circle();
446     if ( Circ.Radius() > theMaxTol )
447       return Standard_False;
448     theActTol = Max(Circ.Radius(), aConfusion);
449     return Standard_True;
450   }
451   else if (Type == GeomAbs_BSplineCurve) {
452     Handle(Geom_BSplineCurve) BS = AC.BSpline();
453     Standard_Integer NbPoles = BS->NbPoles();
454     Standard_Real aMaxPoleDist2 = 0.0, aMaxTol2 = theMaxTol*theMaxTol;
455     gp_Pnt P1,P2;
456     P1 = BS->Pole(1);
457     for (Standard_Integer i = 2; i <= NbPoles; i++) {
458       P2 = BS->Pole(i);
459       Standard_Real aPoleDist2 = P1.SquareDistance(P2);
460       if (aPoleDist2 > aMaxTol2)
461         return Standard_False;
462       if (aPoleDist2 > aMaxPoleDist2)
463         aMaxPoleDist2 = aPoleDist2;
464     }
465     theActTol = Max(1.000001*Sqrt(aMaxPoleDist2), aConfusion);
466     return Standard_True;
467   }
468   else if (Type == GeomAbs_BezierCurve) {
469     Handle(Geom_BezierCurve) BZ = AC.Bezier();
470     Standard_Integer NbPoles = BZ->NbPoles();
471     Standard_Real aMaxPoleDist2 = 0.0, aMaxTol2 = theMaxTol*theMaxTol;
472     gp_Pnt P1,P2;
473     P1 = BZ->Pole(1);
474     for (Standard_Integer i = 2; i <= NbPoles; i++) {
475       P2 = BZ->Pole(i);
476       Standard_Real aPoleDist2 = P1.SquareDistance(P2);
477       if (aPoleDist2 > aMaxTol2)
478         return Standard_False;
479       if (aPoleDist2 > aMaxPoleDist2)
480         aMaxPoleDist2 = aPoleDist2;
481     }
482     theActTol = Max(1.000001*Sqrt(aMaxPoleDist2), aConfusion);
483     return Standard_True;
484   }
485   
486   return Standard_False;
487 }
488
489 //=======================================================================
490 //function : Init
491 //purpose  : 
492 //=======================================================================
493
494 void  BRepLib_MakeFace::Init(const Handle(Geom_Surface)& SS,
495                              const Standard_Real Um,
496                              const Standard_Real UM,
497                              const Standard_Real Vm,
498                              const Standard_Real VM,
499                              const Standard_Real TolDegen)
500 {
501   myError = BRepLib_FaceDone;
502
503   Standard_Real UMin = Um;
504   Standard_Real UMax = UM;
505   Standard_Real VMin = Vm;
506   Standard_Real VMax = VM;
507
508   Standard_Real umin,umax,vmin,vmax,T;
509
510   Handle(Geom_Surface) S = SS, BS = SS;
511   Handle(Geom_RectangularTrimmedSurface) RS = 
512     Handle(Geom_RectangularTrimmedSurface)::DownCast(S);
513   if (!RS.IsNull())
514     BS = RS->BasisSurface();
515
516   Standard_Boolean OffsetSurface = 
517     (BS->DynamicType() == STANDARD_TYPE(Geom_OffsetSurface));
518      
519   // adjust periodical surface or reordonate
520   // check if the values are in the natural range
521   Standard_Real epsilon = Precision::PConfusion();
522   
523   BS->Bounds(umin,umax,vmin,vmax);
524
525   if (OffsetSurface) {
526     Handle(Geom_OffsetSurface) OS = Handle(Geom_OffsetSurface)::DownCast(BS);
527     Handle(Geom_Surface) Base = OS->BasisSurface();
528
529     if (Base->DynamicType() == STANDARD_TYPE(Geom_SurfaceOfLinearExtrusion)) {
530       if (Precision::IsInfinite(umin) || Precision::IsInfinite(umax))
531         S = new Geom_RectangularTrimmedSurface(OS, UMin, UMax, VMin, VMax);
532       else
533         S = new Geom_RectangularTrimmedSurface(OS, VMin, VMax, Standard_False);
534     }  else if (Base->DynamicType() == STANDARD_TYPE(Geom_SurfaceOfRevolution)) {
535       if (Precision::IsInfinite(vmin) || Precision::IsInfinite(vmax))
536         S = new Geom_RectangularTrimmedSurface(OS, VMin, VMax, Standard_False);
537     }
538   }
539
540   if (S->IsUPeriodic()) {
541     ElCLib::AdjustPeriodic(umin,umax,epsilon,UMin,UMax);
542   }
543   else if (UMin > UMax) {
544     T = UMin;
545     UMin = UMax;
546     UMax = T;
547     if ((umin - UMin > epsilon) || (UMax - umax > epsilon)) {
548       myError = BRepLib_ParametersOutOfRange;
549       return;
550     }
551   }
552
553   if (S->IsVPeriodic()) {
554     ElCLib::AdjustPeriodic(vmin,vmax,epsilon,VMin,VMax);
555   }
556   else if (VMin > VMax) {
557     T = VMin;
558     VMin = VMax;
559     VMax = T;
560     if ((vmin - VMin > epsilon) || (VMax - vmax > epsilon)) {
561       myError = BRepLib_ParametersOutOfRange;
562       return;
563     }
564   }
565   
566
567   // compute infinite flags
568   Standard_Boolean umininf = Precision::IsNegativeInfinite(UMin);
569   Standard_Boolean umaxinf = Precision::IsPositiveInfinite(UMax);
570   Standard_Boolean vmininf = Precision::IsNegativeInfinite(VMin);
571   Standard_Boolean vmaxinf = Precision::IsPositiveInfinite(VMax);
572   
573   // closed flag
574   Standard_Boolean uclosed = 
575     S->IsUClosed() && 
576       Abs(UMin - umin) < epsilon && 
577         Abs(UMax - umax) < epsilon;
578   
579   Standard_Boolean vclosed = 
580     S->IsVClosed() && 
581       Abs(VMin - vmin) < epsilon && 
582         Abs(VMax - vmax) < epsilon;
583   
584   
585   // compute 3d curves and degenerate flag
586   Standard_Real maxTol = TolDegen;
587   Handle(Geom_Curve) Cumin,Cumax,Cvmin,Cvmax;
588   Standard_Boolean Dumin,Dumax,Dvmin,Dvmax;
589   Dumin = Dumax = Dvmin = Dvmax = Standard_False;
590   Standard_Real uminTol = Precision::Confusion(), 
591                 umaxTol = Precision::Confusion(), 
592                 vminTol = Precision::Confusion(), 
593                 vmaxTol = Precision::Confusion();
594
595   if (!umininf) {
596     Cumin = S->UIso(UMin);
597     Dumin = IsDegenerated(Cumin, maxTol, uminTol);
598   }
599   if (!umaxinf) {
600     Cumax = S->UIso(UMax);
601     Dumax = IsDegenerated(Cumax, maxTol, umaxTol);
602   }
603   if (!vmininf) {
604     Cvmin = S->VIso(VMin);
605     Dvmin = IsDegenerated(Cvmin, maxTol, vminTol);
606   }
607   if (!vmaxinf) {
608     Cvmax = S->VIso(VMax);
609     Dvmax = IsDegenerated(Cvmax, maxTol, vmaxTol);
610   }
611
612   // compute vertices
613   BRep_Builder B;
614   
615   TopoDS_Vertex V00,V10,V11,V01;
616   
617   if (!umininf) {
618     if (!vmininf) B.MakeVertex(V00,S->Value(UMin,VMin), Max(uminTol, vminTol));
619     if (!vmaxinf) B.MakeVertex(V01,S->Value(UMin,VMax), Max(uminTol, vmaxTol));
620   }
621   if (!umaxinf) {
622     if (!vmininf) B.MakeVertex(V10,S->Value(UMax,VMin), Max(umaxTol, vminTol));
623     if (!vmaxinf) B.MakeVertex(V11,S->Value(UMax,VMax), Max(umaxTol, vmaxTol));
624   }
625
626   if (uclosed) {
627     V10 = V00;
628     V11 = V01;
629   }
630
631   if (vclosed) {
632     V01 = V00;
633     V11 = V10;
634   }
635
636   if (Dumin) V00 = V01;
637   if (Dumax) V10 = V11;
638   if (Dvmin) V00 = V10;
639   if (Dvmax) V01 = V11;
640
641   // make the lines
642   Handle(Geom2d_Line) Lumin,Lumax,Lvmin,Lvmax;
643   if (!umininf)
644     Lumin = new Geom2d_Line(gp_Pnt2d(UMin,0),gp_Dir2d(0,1));
645   if (!umaxinf)
646     Lumax = new Geom2d_Line(gp_Pnt2d(UMax,0),gp_Dir2d(0,1));
647   if (!vmininf)
648     Lvmin = new Geom2d_Line(gp_Pnt2d(0,VMin),gp_Dir2d(1,0));
649   if (!vmaxinf)
650     Lvmax = new Geom2d_Line(gp_Pnt2d(0,VMax),gp_Dir2d(1,0));
651   
652   // make the face
653   TopoDS_Face& F = TopoDS::Face(myShape);
654   B.MakeFace(F,S,Precision::Confusion());
655
656   // make the edges
657   TopoDS_Edge eumin,eumax,evmin,evmax;
658
659   if (!umininf) {
660     if (!Dumin)
661       B.MakeEdge(eumin,Cumin,uminTol);
662     else
663       B.MakeEdge(eumin);
664     if (uclosed) 
665       B.UpdateEdge(eumin,Lumax,Lumin,F,Max(uminTol, umaxTol));
666     else
667       B.UpdateEdge(eumin,Lumin,F,uminTol);
668     B.Degenerated(eumin,Dumin);
669     if (!vmininf) {
670       V00.Orientation(TopAbs_FORWARD);
671       B.Add(eumin,V00);
672     }
673     if (!vmaxinf) {
674       V01.Orientation(TopAbs_REVERSED);
675       B.Add(eumin,V01);
676     }
677     B.Range(eumin,VMin,VMax);
678   }
679
680   if (!umaxinf) {
681     if (uclosed)
682       eumax = eumin;
683     else {
684       if (!Dumax)
685         B.MakeEdge(eumax,Cumax,umaxTol);
686       else
687         B.MakeEdge(eumax);
688       B.UpdateEdge(eumax,Lumax,F,umaxTol);
689       B.Degenerated(eumax,Dumax);
690       if (!vmininf) {
691         V10.Orientation(TopAbs_FORWARD);
692         B.Add(eumax,V10);
693       }
694       if (!vmaxinf) {
695         V11.Orientation(TopAbs_REVERSED);
696         B.Add(eumax,V11);
697       }
698       B.Range(eumax,VMin,VMax);
699     }
700   }
701
702   if (!vmininf) {
703     if (!Dvmin)
704       B.MakeEdge(evmin,Cvmin,vminTol);
705     else
706       B.MakeEdge(evmin);
707     if (vclosed)
708       B.UpdateEdge(evmin,Lvmin,Lvmax,F,Max(vminTol, vmaxTol));
709     else
710       B.UpdateEdge(evmin,Lvmin,F,vminTol);
711     B.Degenerated(evmin,Dvmin);
712     if (!umininf) {
713       V00.Orientation(TopAbs_FORWARD);
714       B.Add(evmin,V00);
715     }
716     if (!umaxinf) {
717       V10.Orientation(TopAbs_REVERSED);
718       B.Add(evmin,V10);
719     }
720     B.Range(evmin,UMin,UMax);
721   }
722
723   if (!vmaxinf) {
724     if (vclosed)
725       evmax = evmin;
726     else {
727       if (!Dvmax)
728         B.MakeEdge(evmax,Cvmax,vmaxTol);
729       else
730         B.MakeEdge(evmax);
731       B.UpdateEdge(evmax,Lvmax,F,vmaxTol);
732       B.Degenerated(evmax,Dvmax);
733       if (!umininf) {
734         V01.Orientation(TopAbs_FORWARD);
735         B.Add(evmax,V01);
736       }
737       if (!umaxinf) {
738         V11.Orientation(TopAbs_REVERSED);
739         B.Add(evmax,V11);
740       }
741       B.Range(evmax,UMin,UMax);
742     }
743   }
744
745   // make the wires and add them to the face
746   eumin.Orientation(TopAbs_REVERSED);
747   evmax.Orientation(TopAbs_REVERSED);
748   
749   TopoDS_Wire W;
750
751   if (!umininf && !umaxinf && vmininf && vmaxinf) {
752     // two wires in u
753     B.MakeWire(W);
754     B.Add(W,eumin);
755     B.Add(F,W);
756     B.MakeWire(W);
757     B.Add(W,eumax);
758     B.Add(F,W);
759     F.Closed(uclosed);
760   }
761     
762   else if (umininf && umaxinf && !vmininf && !vmaxinf) {
763     // two wires in v
764     B.MakeWire(W);
765     B.Add(W,evmin);
766     B.Add(F,W);
767     B.MakeWire(W);
768     B.Add(W,evmax);
769     B.Add(F,W);
770     F.Closed(vclosed);
771   }
772     
773   else if (!umininf || !umaxinf || !vmininf || !vmaxinf) {
774     // one wire
775     B.MakeWire(W);
776     if (!umininf) B.Add(W,eumin);
777     if (!vmininf) B.Add(W,evmin);
778     if (!umaxinf) B.Add(W,eumax);
779     if (!vmaxinf) B.Add(W,evmax);
780     B.Add(F,W);
781     W.Closed(!umininf && !umaxinf && !vmininf && !vmaxinf);
782     F.Closed(uclosed && vclosed);
783   }
784
785   if (OffsetSurface) { 
786     // Les Isos sont Approximees a Precision::Approximation()
787     // et on code Precision::Confusion() dans l'arete.
788     // ==> Un petit passage dans SamePrameter pour regler les tolerances.
789     BRepLib::SameParameter(F, Precision::Confusion(), Standard_True);
790   }
791
792   Done();
793 }
794
795
796 //=======================================================================
797 //function : Add
798 //purpose  : 
799 //=======================================================================
800
801 void  BRepLib_MakeFace::Add(const TopoDS_Wire& W)
802 {
803   BRep_Builder B;
804   B.Add(myShape,W);
805   B.NaturalRestriction(TopoDS::Face(myShape),Standard_False);
806   Done();
807 }
808
809
810 //=======================================================================
811 //function : Face
812 //purpose  : 
813 //=======================================================================
814
815 const TopoDS_Face&  BRepLib_MakeFace::Face()const 
816 {
817   return TopoDS::Face(myShape);
818 }
819
820
821
822 //=======================================================================
823 //function : operator
824 //purpose  : 
825 //=======================================================================
826
827 BRepLib_MakeFace::operator TopoDS_Face() const
828 {
829   return Face();
830 }
831
832 //=======================================================================
833 //function : Error
834 //purpose  : 
835 //=======================================================================
836
837 BRepLib_FaceError BRepLib_MakeFace::Error() const
838 {
839   return myError;
840 }
841
842
843 //=======================================================================
844 //function : CheckInside
845 //purpose  : Reverses the current face if not a bounded area
846 //=======================================================================
847
848 void BRepLib_MakeFace::CheckInside()
849 {
850   // compute the area and return the face if the area is negative
851   TopoDS_Face F = TopoDS::Face(myShape);
852   BRepTopAdaptor_FClass2d FClass(F,0.);
853   if ( FClass.PerformInfinitePoint() == TopAbs_IN) {
854     BRep_Builder B;
855     TopoDS_Shape S = myShape.EmptyCopied();
856     TopoDS_Iterator it(myShape);
857     while (it.More()) {
858       B.Add(S,it.Value().Reversed());
859       it.Next();
860     }
861     myShape = S;
862   }
863 }