Integration of OCCT 6.5.0 from SVN
[occt.git] / src / BRepLib / BRepLib_MakeShell.cxx
1 // File:        BRepLib_MakeShell.cxx
2 // Created:     Wed Jan  4 13:37:46 1995
3 // Author:      Bruno DUMORTIER
4 //              <dub@fuegox>
5
6
7 #include <BRepLib_MakeShell.ixx>
8
9 #include <BRepLib.hxx>
10 #include <Precision.hxx>
11 #include <TColStd_Array1OfReal.hxx>
12
13 #include <Geom_RectangularTrimmedSurface.hxx>
14 #include <Geom_BSplineSurface.hxx>
15 #include <Geom2d_Line.hxx>
16 #include <GeomAdaptor_Surface.hxx>
17 #include <TColGeom2d_Array1OfCurve.hxx>
18 #include <gp_Pnt.hxx>
19
20 #include <TopoDS.hxx>
21 #include <TopoDS_Face.hxx>
22 #include <TopoDS_Wire.hxx>
23 #include <TopoDS_Edge.hxx>
24 #include <TopoDS_Vertex.hxx>
25 #include <TopTools_Array1OfShape.hxx>
26 #include <BRep_Builder.hxx>
27
28
29 //=======================================================================
30 //function : BRepLib_MakeShell
31 //purpose  : 
32 //=======================================================================
33
34 BRepLib_MakeShell::BRepLib_MakeShell() :
35        myError(BRepLib_EmptyShell)
36 {
37 }
38
39
40 //=======================================================================
41 //function : BRepLib_MakeShell
42 //purpose  : 
43 //=======================================================================
44
45 BRepLib_MakeShell::BRepLib_MakeShell(const Handle(Geom_Surface)& S,
46                                      const Standard_Boolean Segment)
47 {
48   Standard_Real UMin,UMax,VMin,VMax;
49   S->Bounds(UMin,UMax,VMin,VMax);
50   Init(S,UMin,UMax,VMin,VMax,Segment);
51 }
52
53
54 //=======================================================================
55 //function : BRepLib_MakeShell
56 //purpose  : 
57 //=======================================================================
58
59 BRepLib_MakeShell::BRepLib_MakeShell(const Handle(Geom_Surface)& S, 
60                                      const Standard_Real UMin,
61                                      const Standard_Real UMax, 
62                                      const Standard_Real VMin, 
63                                      const Standard_Real VMax,
64                                      const Standard_Boolean Segment)
65 {
66   Init(S,UMin,UMax,VMin,VMax,Segment);
67 }
68
69
70 //=======================================================================
71 //function : Init
72 //purpose  : 
73 //=======================================================================
74
75 void BRepLib_MakeShell::Init(const Handle(Geom_Surface)& S, 
76                              const Standard_Real UMin, 
77                              const Standard_Real UMax, 
78                              const Standard_Real VMin, 
79                              const Standard_Real VMax,
80                              const Standard_Boolean Segment)
81 {
82   Handle(Geom_Surface) BS = S;
83   if ( S->DynamicType() == STANDARD_TYPE(Geom_RectangularTrimmedSurface)) {
84     Handle(Geom_RectangularTrimmedSurface) RTS = 
85       Handle(Geom_RectangularTrimmedSurface)::DownCast(S);
86     BS = RTS->BasisSurface();
87   }
88   myError = BRepLib_EmptyShell;
89   Standard_Real tol = Precision::Confusion();
90
91   // Make a shell from a surface
92   GeomAdaptor_Surface GS(BS,UMin,UMax,VMin,VMax);
93
94   Standard_Integer nu = GS.NbUIntervals(GeomAbs_C2);
95   Standard_Integer nv = GS.NbVIntervals(GeomAbs_C2);
96
97   Standard_Boolean uperiodic = GS.IsUPeriodic();
98   Standard_Boolean vperiodic = GS.IsVPeriodic();
99
100   if (nu == 0 || nv == 0) return;
101
102   // arrays of parameters and pcurves
103   TColStd_Array1OfReal upars(1,nu+1);
104   TColStd_Array1OfReal vpars(1,nv+1);
105   TColGeom2d_Array1OfCurve uisos(1,nu+1);
106   TColGeom2d_Array1OfCurve visos(1,nv+1);
107
108   Standard_Integer iu,iv;
109
110   GS.UIntervals(upars,GeomAbs_C2);
111   gp_Dir2d dv(0,1);
112   for (iu = 1; iu <= nu+1; iu++) {
113     Standard_Real u = upars(iu);
114     if (!Precision::IsInfinite(u))
115       uisos(iu)  = new Geom2d_Line(gp_Pnt2d(u,0.),dv);
116   }
117   
118   GS.VIntervals(vpars,GeomAbs_C2);
119   gp_Dir2d du(1,0);
120   for (iv = 1; iv <= nv+1 ; iv++) {
121     Standard_Real v = vpars(iv);
122     if (!Precision::IsInfinite(v))
123       visos(iv)  = new Geom2d_Line(gp_Pnt2d(0.,v),du);
124   }
125     
126   // create row by row
127
128   // create the shell
129   BRep_Builder B;
130   B.MakeShell(TopoDS::Shell(myShape));
131
132   // arrays of edges and vertices for each row
133   TopTools_Array1OfShape botedges(1,nu);
134   TopTools_Array1OfShape botvertices(1,nu+1);
135
136   // copies of the first ones for periodic case
137   TopTools_Array1OfShape fbotedges(1,nu);
138   TopTools_Array1OfShape fbotvertices(1,nu+1);
139
140   TopoDS_Face F;
141   TopoDS_Wire W;
142   TopoDS_Edge eleft,eright,etop,ebot,feleft;
143   TopoDS_Vertex vlb,vlt,vrb,vrt,fvlt;
144
145
146   // init the botedges and botvertices
147   if (!Precision::IsInfinite(vpars(1))) {
148     if (!Precision::IsInfinite(upars(1)))
149       B.MakeVertex(vrt,S->Value(upars(1),vpars(1)),tol);
150     fbotvertices(1) = botvertices(1) = vrt;
151
152     for (iu = 1; iu <= nu; iu++) {
153       vlt = vrt;
154
155       if (uperiodic && iu == nu)
156         vrt = TopoDS::Vertex(botvertices(1));
157       else if (!Precision::IsInfinite(upars(iu+1)))
158         B.MakeVertex(vrt,S->Value(upars(iu+1),vpars(1)),tol);
159
160       fbotvertices(iu+1) = botvertices(iu+1) = vrt;
161       B.MakeEdge(etop);
162       if (!vlt.IsNull()) {
163         vlt.Orientation(TopAbs_FORWARD);
164         B.Add(etop,vlt);
165       }
166       if (!vrt.IsNull()) {
167         vrt.Orientation(TopAbs_REVERSED);
168         B.Add(etop,vrt);
169       }
170       fbotedges(iu) = botedges(iu) = etop;
171     }
172   }
173
174   for (iv = 1; iv <= nv; iv++) {
175
176     // compute the first edge and vertices of the line
177     vrb = TopoDS::Vertex(botvertices(1));
178
179     if (vperiodic && iv == nv) {
180       vrt = TopoDS::Vertex(fbotvertices(1));
181     }
182     else {
183       vrt.Nullify();
184       if (!Precision::IsInfinite(vpars(iv+1))) {
185         if (!Precision::IsInfinite(upars(1)))
186           B.MakeVertex(vrt,S->Value(upars(1),vpars(iv+1)),tol);
187       }
188     }
189
190     eright.Nullify();
191     if (!Precision::IsInfinite(upars(1))) {
192       B.MakeEdge(eright);
193       if (!vrb.IsNull()) {
194         vrb.Orientation(TopAbs_FORWARD);
195         B.Add(eright,vrb);
196       }
197       if (!vrt.IsNull()) {
198         vrt.Orientation(TopAbs_REVERSED);
199         B.Add(eright,vrt);
200       }
201     }
202
203     fvlt  = vrt;
204     feleft = eright;
205
206
207     // make the row of faces
208     
209     for (iu = 1; iu <= nu; iu++) {
210       
211       // create the face at iu, iv
212
213       // the surface
214       Handle(Geom_Surface) SS = Handle(Geom_Surface)::DownCast(BS->Copy());
215       if (GS.GetType() == GeomAbs_BSplineSurface && Segment) {
216         Handle(Geom_BSplineSurface)::DownCast(SS)
217           ->Segment(upars(iu),upars(iu+1),
218                     vpars(iv),vpars(iv+1) );
219       }
220       B.MakeFace(F,SS,tol);
221
222       // the wire
223
224       B.MakeWire(W);
225
226       // the vertices
227
228       vlb = vrb;
229       vrb = TopoDS::Vertex(botvertices(iu+1));
230       vlt = vrt;
231
232       if (uperiodic && iu == nu)
233         vrt = fvlt;
234       else {
235         vrt.Nullify();
236         if (!Precision::IsInfinite(vpars(iv+1))) {
237           if (!Precision::IsInfinite(upars(iu+1)))
238             B.MakeVertex(vrt,S->Value(upars(iu+1),vpars(iv+1)),tol);
239         }
240       }
241
242       botvertices(iu)   = vlt;
243       botvertices(iu+1) = vrt;
244
245       // the edges
246
247       eleft = eright;
248
249       if  (uperiodic && iu == nu)
250         eright = feleft;
251       else {
252         eright.Nullify();
253         if (!Precision::IsInfinite(upars(iu+1))) {
254           B.MakeEdge(eright);
255           if (!vrb.IsNull()) {
256             vrb.Orientation(TopAbs_FORWARD);
257             B.Add(eright,vrb);
258           }
259           if (!vrt.IsNull()) {
260             vrt.Orientation(TopAbs_REVERSED);
261             B.Add(eright,vrt);
262           }
263         }
264       }
265
266       if ( uperiodic && nu == 1) {
267         if (!eleft.IsNull() && !eright.IsNull()) {
268           B.UpdateEdge(eleft,uisos(2),uisos(1),F,tol);
269           B.Range(eleft,F,vpars(iv),vpars(iv+1));
270         }
271       }
272       else {
273         if (!eleft.IsNull()) {
274           B.UpdateEdge(eleft,uisos(iu),F,tol);
275           B.Range(eleft,F,vpars(iv),vpars(iv+1));
276         }
277         if (!eright.IsNull()) {
278           B.UpdateEdge(eright,uisos(iu+1),F,tol);
279           B.Range(eright,F,vpars(iv),vpars(iv+1));
280         }
281       }
282       
283       ebot = TopoDS::Edge(botedges(iu));
284
285       if (vperiodic && iv == nv) 
286         etop = TopoDS::Edge(fbotedges(iu));
287       else {
288         etop.Nullify();
289         if (!Precision::IsInfinite(vpars(iv+1))) {
290           B.MakeEdge(etop);
291           if (!vlt.IsNull()) {
292             vlt.Orientation(TopAbs_FORWARD);
293             B.Add(etop,vlt);
294           }
295           if (!vrt.IsNull()) {
296             vrt.Orientation(TopAbs_REVERSED);
297             B.Add(etop,vrt);
298           }
299         }
300       }
301
302       if ( vperiodic && nv == 1) {
303         if (!ebot.IsNull() && !etop.IsNull()) {
304           B.UpdateEdge(ebot,visos(1),visos(2),F,tol);
305           B.Range(ebot,F,vpars(iv),vpars(iv+1));
306         }
307       }
308       else {
309         if (!ebot.IsNull()) {
310           B.UpdateEdge(ebot,visos(iv),F,tol);
311           B.Range(ebot,F,upars(iu),upars(iu+1));
312         }
313         if (!etop.IsNull()) {
314           B.UpdateEdge(etop,visos(iv+1),F,tol);
315           B.Range(etop,F,upars(iu),upars(iu+1));
316         }
317       }
318       
319       botedges(iu) = etop;
320
321       if (!eleft.IsNull()) {
322         eleft.Orientation(TopAbs_REVERSED);
323         B.Add(W,eleft);
324       }
325       if (!ebot.IsNull()) {
326         ebot.Orientation(TopAbs_FORWARD);
327         B.Add(W,ebot);
328       }
329       if (!eright.IsNull()) {
330         eright.Orientation(TopAbs_FORWARD);
331         B.Add(W,eright);
332       }
333       if (!etop.IsNull()) {
334         etop.Orientation(TopAbs_REVERSED);
335         B.Add(W,etop);
336       }
337       
338       B.Add(F,W);
339       B.Add(myShape,F);
340     }
341   }
342   
343   // codage des courbes 3d et regularites.
344   BRepLib::BuildCurves3d(myShape,tol);
345   BRepLib::EncodeRegularity(myShape);
346   
347   myError = BRepLib_ShellDone;
348   Done();
349 }
350
351
352 //=======================================================================
353 //function : Error
354 //purpose  : 
355 //=======================================================================
356
357 BRepLib_ShellError BRepLib_MakeShell::Error() const 
358 {
359   return myError;
360 }
361
362
363 //=======================================================================
364 //function : TopoDS_Shell&
365 //purpose  : 
366 //=======================================================================
367
368 const TopoDS_Shell& BRepLib_MakeShell::Shell() const 
369 {
370   return TopoDS::Shell(myShape);
371 }
372
373
374
375 //=======================================================================
376 //function : TopoDS_Shell
377 //purpose  : 
378 //=======================================================================
379
380 BRepLib_MakeShell::operator TopoDS_Shell() const
381 {
382   return Shell();
383 }
384
385