0022922: Clean up warnings on uninitialized / unused variables
[occt.git] / src / Extrema / Extrema_ExtCS.cxx
1 // File:        Extrema_ExtCS.cxx
2 // Created:     Tue Jul 18 10:33:04 1995
3 // Author:      Modelistation
4 //              <model@metrox>
5
6 //  Modified by skv - Thu Jul  7 12:29:34 2005 OCC9134
7
8 #include <Extrema_ExtCS.ixx>
9 #include <Extrema_GenExtCS.hxx>
10 #include <StdFail_NotDone.hxx>
11 #include <Standard_NotImplemented.hxx>
12 #include <StdFail_InfiniteSolutions.hxx>
13 #include <Precision.hxx>
14 #include <GeomAbs_CurveType.hxx>
15 #include <gp_Pnt.hxx>
16 #include <gp_Pln.hxx>
17 #include <gp_Cylinder.hxx>
18 #include <gp_Cone.hxx>
19 #include <gp_Sphere.hxx>
20 #include <gp_Torus.hxx>
21 #include <gp_Lin.hxx>
22
23 #include <ElCLib.hxx>
24 #include <Extrema_ExtPElC.hxx>
25 #include <Bnd_Box.hxx>
26 #include <BndLib_AddSurface.hxx>
27 #include <Extrema_ExtPElS.hxx>
28
29 Extrema_ExtCS::Extrema_ExtCS() 
30 {
31   myDone = Standard_False;
32 }
33
34 Extrema_ExtCS::Extrema_ExtCS(const Adaptor3d_Curve&   C,
35                              const Adaptor3d_Surface& S,
36                              const Standard_Real    TolC,
37                              const Standard_Real    TolS)
38
39 {
40   Initialize(S, S.FirstUParameter(), S.LastUParameter(), 
41              S.FirstVParameter(), S.LastVParameter(), 
42              TolC, TolS);
43   Perform(C, C.FirstParameter(), C.LastParameter());
44 }
45
46 Extrema_ExtCS::Extrema_ExtCS(const Adaptor3d_Curve&   C,
47                              const Adaptor3d_Surface& S,
48                              const Standard_Real    UCinf,
49                              const Standard_Real    UCsup,
50                              const Standard_Real    Uinf,       
51                              const Standard_Real    Usup,
52                              const Standard_Real    Vinf,       
53                              const Standard_Real    Vsup,
54                              const Standard_Real    TolC,
55                              const Standard_Real    TolS)
56
57 {
58   Initialize(S, Uinf, Usup, Vinf, Vsup, TolC, TolS);
59   Perform(C, UCinf, UCsup);
60 }
61
62
63 void Extrema_ExtCS::Initialize(const Adaptor3d_Surface& S,
64                                const Standard_Real    Uinf,     
65                                const Standard_Real    Usup,
66                                const Standard_Real    Vinf,     
67                                const Standard_Real    Vsup,
68                                const Standard_Real    TolC,
69                                const Standard_Real    TolS)
70 {
71   myS = (Adaptor3d_SurfacePtr)&S;
72   myIsPar = Standard_False;
73   myuinf  = Uinf;
74   myusup  = Usup;
75   myvinf  = Vinf;
76   myvsup  = Vsup;
77   mytolC  = TolC;
78   mytolS  = TolS;
79   myStype  = myS->GetType();
80 }
81
82                                 
83 void Extrema_ExtCS::Perform(const Adaptor3d_Curve& C,
84                             const Standard_Real  Uinf,
85                             const Standard_Real  Usup)
86 {
87   myucinf = Uinf;
88   myucsup = Usup;
89   myPOnS.Clear();
90   myPOnC.Clear();
91   mySqDist.Clear();
92   Standard_Integer i;
93   Standard_Integer NbT, NbU, NbV;
94   NbT = NbU = NbV = 10;
95   GeomAbs_CurveType myCtype  = C.GetType();
96
97
98   switch(myCtype) {
99
100   case GeomAbs_Line: 
101     {
102       
103       switch(myStype) {
104       case GeomAbs_Sphere:
105         myExtElCS.Perform(C.Line(), myS->Sphere());
106         break;   
107       case GeomAbs_Cylinder:
108         myExtElCS.Perform(C.Line(), myS->Cylinder());
109         break;
110       case GeomAbs_Plane:
111         myExtElCS.Perform(C.Line(), myS->Plane());
112         if (myExtElCS.IsParallel())   break;
113
114       case GeomAbs_Torus:
115       case GeomAbs_Cone:
116       case GeomAbs_BezierSurface:
117       case GeomAbs_BSplineSurface:
118       case GeomAbs_SurfaceOfRevolution:
119       case GeomAbs_SurfaceOfExtrusion:
120       case GeomAbs_OtherSurface:
121         {
122           Standard_Real cfirst = myucinf, clast = myucsup;
123           Standard_Real ufirst = myS->FirstUParameter(), ulast = myS->LastUParameter(), 
124                         vfirst = myS->FirstVParameter(), vlast = myS->LastVParameter();
125
126           if(Precision::IsInfinite(Abs(cfirst)) || Precision::IsInfinite(Abs(clast))) {
127
128             Bnd_Box aSurfBox;
129       BndLib_AddSurface::Add(*myS, ufirst, ulast, vfirst, vlast, Precision::Confusion(), aSurfBox);
130             Standard_Real xmin, ymin, zmin, xmax, ymax, zmax;
131             aSurfBox.Get(xmin, ymin, zmin, xmax, ymax, zmax);
132             Standard_Real tmin = Precision::Infinite(), tmax = -tmin;
133             gp_Lin aLin = C.Line();
134             
135
136             if(!( Precision::IsInfinite(Abs(xmin)) || Precision::IsInfinite(Abs(xmax)) || 
137                   Precision::IsInfinite(Abs(ymin)) || Precision::IsInfinite(Abs(ymax)) || 
138                   Precision::IsInfinite(Abs(zmin)) || Precision::IsInfinite(Abs(zmax)))  ) {
139               
140               Extrema_ExtPElC anExt;
141               Extrema_POnCurv aPntOnLin;
142               Standard_Real aParOnLin;
143               Standard_Real lim = Precision::Infinite();
144               gp_Pnt aLimPntArray[8];
145               
146               aLimPntArray[0].SetCoord(xmin, ymin, zmin);
147               aLimPntArray[1].SetCoord(xmax, ymin, zmin);
148               aLimPntArray[2].SetCoord(xmin, ymax, zmin);
149               aLimPntArray[3].SetCoord(xmax, ymax, zmin);
150               aLimPntArray[4].SetCoord(xmin, ymin, zmax);
151               aLimPntArray[5].SetCoord(xmax, ymin, zmax);
152               aLimPntArray[6].SetCoord(xmin, ymax, zmax);
153               aLimPntArray[7].SetCoord(xmax, ymax, zmax);
154
155               for(i = 0; i <= 7; i++) {
156                 anExt.Perform(aLimPntArray[i], aLin, Precision::Confusion(), -lim, lim);
157                 aPntOnLin = anExt.Point(1);
158                 aParOnLin = aPntOnLin.Parameter();
159                 tmin = Min(aParOnLin, tmin);
160                 tmax = Max(aParOnLin, tmax);
161               }
162               
163             }
164             else {
165               tmin = -1.e+50;
166               tmax =  1.e+50;
167             }
168
169
170             cfirst = Max(cfirst, tmin);
171             clast  = Min(clast,  tmax);
172
173           }
174
175
176             
177           Extrema_GenExtCS Ext(C, *myS, NbT, NbU, NbV, cfirst, clast, ufirst, ulast,
178                                vfirst, vlast, mytolC, mytolS);
179
180           myDone = Ext.IsDone();
181           if (myDone) {
182             Standard_Integer NbExt = Ext.NbExt();
183             Standard_Real T,U,V;
184             Extrema_POnCurv PC;
185             Extrema_POnSurf PS;
186             for (i = 1; i <= NbExt; i++) {
187               PC = Ext.PointOnCurve(i);
188               PS = Ext.PointOnSurface(i);
189               T = PC.Parameter();
190               PS.Parameter(U, V);
191               if (myS->IsUPeriodic())
192                 U = ElCLib::InPeriod(U, myuinf, myuinf+myS->UPeriod());
193               if (myS->IsVPeriodic())
194                 V = ElCLib::InPeriod(V, myvinf, myvinf+myS->VPeriod());
195
196               if ((myucinf-T) <= mytolC && (T-myucsup) <= mytolC &&
197                   (myuinf-U) <= mytolS && (U-myusup) <= mytolS &&
198                   (myvinf-V) <= mytolS && (V-myvsup) <= mytolS) {
199                 mySqDist.Append(Ext.SquareDistance(i));
200                 myPOnC.Append(PC);
201                 myPOnS.Append(Extrema_POnSurf(U, V, PS.Value()));
202
203               
204               }
205             }
206           }
207           return;
208           
209         }
210 #ifndef DEB
211       default:
212 #endif
213         break;
214       }
215       break;
216     }
217 //  Modified by skv - Thu Jul  7 12:29:34 2005 OCC9134 Begin
218   case GeomAbs_Circle:
219     {
220       if(myStype == GeomAbs_Cylinder) {
221         myExtElCS.Perform(C.Circle(), myS->Cylinder());
222         break;
223       }
224     }
225   case GeomAbs_Hyperbola: 
226     {
227       if(myCtype == GeomAbs_Hyperbola && myStype == GeomAbs_Plane) {
228 //  Modified by skv - Thu Jul  7 12:29:34 2005 OCC9134 End
229         myExtElCS.Perform(C.Hyperbola(), myS->Plane());
230         break;
231       }
232     }
233   default:
234     {
235       Extrema_GenExtCS Ext;
236       Ext.Initialize(*myS, NbU, NbV, mytolS);
237       if(myCtype == GeomAbs_Hyperbola) {
238         Standard_Real tmin = Max(-20., C.FirstParameter());
239         Standard_Real tmax = Min(20., C.LastParameter());
240         Ext.Perform(C, NbT, tmin, tmax, mytolC); // to avoid overflow
241       }
242       else {
243         if(myCtype == GeomAbs_Circle && NbT < 13) {
244           NbT = 13;
245         }
246         Ext.Perform(C, NbT, mytolC);
247       }
248         
249       myDone = Ext.IsDone();
250       if (myDone) {
251         Standard_Integer NbExt = Ext.NbExt();
252         Standard_Real T,U,V;
253         Extrema_POnCurv PC;
254         Extrema_POnSurf PS;
255         for (i = 1; i <= NbExt; i++) {
256           PC = Ext.PointOnCurve(i);
257           PS = Ext.PointOnSurface(i);
258           T = PC.Parameter();
259           PS.Parameter(U, V);
260           if (C.IsPeriodic())
261             T = ElCLib::InPeriod(T, myucinf, myucinf+C.Period());
262           if (myS->IsUPeriodic())
263             U = ElCLib::InPeriod(U, myuinf, myuinf+myS->UPeriod());
264           if (myS->IsVPeriodic())
265             V = ElCLib::InPeriod(V, myvinf, myvinf+myS->VPeriod());
266           
267           if ((myucinf-T) <= mytolC && (T-myucsup) <= mytolC &&
268               (myuinf-U) <= mytolS && (U-myusup) <= mytolS &&
269               (myvinf-V) <= mytolS && (V-myvsup) <= mytolS) {
270             mySqDist.Append(Ext.SquareDistance(i));
271             PC.SetValues(T, PC.Value());
272             myPOnC.Append(PC);
273             myPOnS.Append(Extrema_POnSurf(U, V, PS.Value()));
274           }
275         }
276       }
277       return;
278     }
279     break;
280   }
281   
282   myDone = myExtElCS.IsDone();
283   if (myDone) {
284     myIsPar = myExtElCS.IsParallel();
285     if (myIsPar) {
286       mySqDist.Append(myExtElCS.SquareDistance(1));
287     }
288     else {
289       Standard_Integer NbExt = myExtElCS.NbExt();
290       Standard_Real U, V;
291       for (i = 1; i <= NbExt; i++) {
292         Extrema_POnCurv PC;
293         Extrema_POnSurf PS;
294         myExtElCS.Points(i, PC, PS);
295         Standard_Real Ucurve = PC.Parameter();
296         PS.Parameter(U, V);
297
298         if((myStype == GeomAbs_Sphere) || (myStype == GeomAbs_Cylinder)) {
299           U = ElCLib::InPeriod(U, myuinf, myuinf+2.*M_PI);
300         }
301
302         if ((myuinf-U) <= mytolS && (U-myusup) <= mytolS &&
303             (myvinf-V) <= mytolS && (V-myvsup) <= mytolS &&
304             (myucinf-Ucurve) <= mytolC && (Ucurve-myucsup) <= mytolC) {
305           mySqDist.Append(myExtElCS.SquareDistance(i));
306           myPOnS.Append(Extrema_POnSurf(U, V, PS.Value()));
307           myPOnC.Append(PC);
308         }
309       }
310     }
311   }
312   
313 }
314
315
316 Standard_Boolean Extrema_ExtCS::IsDone() const
317 {
318   return myDone;
319 }
320
321 Standard_Boolean Extrema_ExtCS::IsParallel() const
322 {
323   return myIsPar;
324 }
325
326
327 Standard_Real Extrema_ExtCS::SquareDistance(const Standard_Integer N) const
328 {
329   if(!myDone) StdFail_NotDone::Raise();
330   if (myIsPar && N != 1) StdFail_InfiniteSolutions::Raise();
331   if ((N < 1) || (N > mySqDist.Length())) Standard_OutOfRange::Raise();
332   return mySqDist.Value(N);
333 }
334
335
336 Standard_Integer Extrema_ExtCS::NbExt() const
337 {
338   if(!myDone) StdFail_NotDone::Raise();
339   return mySqDist.Length();
340 }
341
342
343
344 void Extrema_ExtCS::Points(const Standard_Integer N,
345                             Extrema_POnCurv&       P1,
346                             Extrema_POnSurf&       P2) const
347 {
348   if(!myDone) StdFail_NotDone::Raise();
349   P1 = myPOnC.Value(N);
350   P2 = myPOnS.Value(N);
351 }