0031035: Coding - uninitialized class fields reported by Visual Studio Code Analysis
[occt.git] / src / AdvApp2Var / AdvApp2Var_Context.cxx
1 // Created on: 1996-07-02
2 // Created by: Joelle CHAUVET
3 // Copyright (c) 1996-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 <AdvApp2Var_ApproxF2var.hxx>
19 #include <AdvApp2Var_Context.hxx>
20 #include <Standard_ConstructionError.hxx>
21
22 // Calculaton of parameters
23 static Standard_Boolean lesparam(const Standard_Integer iordre,
24                                  const Standard_Integer ncflim, 
25                                  const Standard_Integer icodeo, 
26                                  Standard_Integer& nbpnts, 
27                                  Standard_Integer& ndgjac)
28 {
29   // jacobi degree
30   ndgjac = ncflim; // it always keeps a reserve coefficient
31   if (icodeo< 0) return Standard_False;
32   if (icodeo > 0) {
33     ndgjac += (9 - (iordre+1)); //iordre rescales the frequences upwards
34     ndgjac += (icodeo-1)*10;
35   }
36   // ---> Min Number of required pointss.
37   if (ndgjac < 8) { nbpnts = 8;  } 
38   else if (ndgjac < 10) { nbpnts = 10; }
39   //  else if (ndgjac < 15) { nbpnt = 15; } Bug Uneven number
40   else if (ndgjac < 20) { nbpnts = 20;}
41   //  else if (ndgjac < 25) { nbpnt = 25; } Bug Uneven number
42   else if (ndgjac < 30) { nbpnts = 30;}
43   else if (ndgjac < 40) { nbpnts = 40;}
44   else if (ndgjac < 50) { nbpnts = 50;}
45   //  else if (*ndgjac < 61) { nbpnt = 61;} Bug Uneven number
46   else {
47     nbpnts = 50;
48 #ifdef OCCT_DEBUG
49     std::cout << "F(U, V) : Not enough points of discretization" << std::endl; 
50 #endif
51   }
52
53   // If constraints are on borders, this adds 2 points
54   if (iordre>-1) {  nbpnts += 2;}
55
56   return Standard_True;
57 }
58 //============================================================================
59 //function : AdvApp2Var_Context
60 //purpose  :
61 //============================================================================
62
63  AdvApp2Var_Context::
64 AdvApp2Var_Context()
65 : myFav(0),
66   myOrdU(0),
67   myOrdV(0),
68   myLimU(0),
69   myLimV(0),
70   myNb1DSS(0),
71   myNb2DSS(0),
72   myNb3DSS(0),
73   myNbURoot(0),
74   myNbVRoot(0),
75   myJDegU(0),
76   myJDegV(0)
77 {
78 }
79
80 //============================================================================
81 //function : AdvApp2Var_Context
82 //purpose  :
83 //============================================================================
84
85  AdvApp2Var_Context::
86 AdvApp2Var_Context(const Standard_Integer ifav,
87                    const Standard_Integer iu,
88                    const Standard_Integer iv, 
89                    const Standard_Integer nlimu,
90                    const Standard_Integer nlimv, 
91                    const Standard_Integer iprecis,
92                    const Standard_Integer nb1Dss,
93                    const Standard_Integer nb2Dss,
94                    const Standard_Integer nb3Dss,
95                    const Handle(TColStd_HArray1OfReal)& tol1D,
96                    const Handle(TColStd_HArray1OfReal)& tol2D,
97                    const Handle(TColStd_HArray1OfReal)& tol3D,
98                    const Handle(TColStd_HArray2OfReal)& tof1D,
99                    const Handle(TColStd_HArray2OfReal)& tof2D,
100                    const Handle(TColStd_HArray2OfReal)& tof3D) :
101 myFav(ifav),
102 myOrdU(iu),
103 myOrdV(iv),
104 myLimU(nlimu),
105 myLimV(nlimv),
106 myNb1DSS(nb1Dss),
107 myNb2DSS(nb2Dss),
108 myNb3DSS(nb3Dss)
109 {
110 Standard_Integer ErrorCode=0,NbPntU=0,JDegU=0,NbPntV=0,JDegV=0;
111 Standard_Integer ncfl;
112
113 // myNbURoot,myJDegU
114 ncfl = nlimu;
115 if (ncfl<2*iu+2) ncfl = 2*iu+2;
116 if (!lesparam(iu,ncfl,iprecis,NbPntU,JDegU) )
117   { throw Standard_ConstructionError("AdvApp2Var_Context");}
118 myNbURoot = NbPntU;
119 myJDegU = JDegU;
120 if (iu>-1) NbPntU = myNbURoot - 2;
121
122 // myJMaxU
123 Standard_Integer i,j,size = JDegU-2*iu-1;
124 Handle (TColStd_HArray1OfReal) JMaxU =
125   new TColStd_HArray1OfReal(1,size);
126 Standard_Real *JU_array =
127   (Standard_Real *) &JMaxU->ChangeArray1()(JMaxU->Lower());
128 AdvApp2Var_ApproxF2var::mma2jmx_(&JDegU,(integer *)&iu,JU_array);
129 myJMaxU = JMaxU;
130
131 // myNbVRoot,myJDegV
132 ncfl = nlimv;
133 if (ncfl<2*iv+2) ncfl = 2*iv+2;
134 //Ma1nbp(&iv,&ncfl,&iprec,&NbPntV,&JDegV,&ErrorCode);
135 if (!lesparam(iv, ncfl, iprecis, NbPntV, JDegV) )
136   { throw Standard_ConstructionError("AdvApp2Var_Context");}
137 myNbVRoot = NbPntV;
138 myJDegV = JDegV;
139 if (iv>-1) NbPntV = myNbVRoot - 2;
140
141 // myJMaxV
142 size = JDegV-2*iv-1;
143 Handle (TColStd_HArray1OfReal) JMaxV =
144   new TColStd_HArray1OfReal(1,size);
145 Standard_Real *JV_array =
146   (Standard_Real *) &JMaxV->ChangeArray1()(JMaxV->Lower());
147 AdvApp2Var_ApproxF2var::mma2jmx_(&JDegV,(integer *)&iv,JV_array);
148 myJMaxV = JMaxV;
149
150 // myURoots, myVRoots
151 Handle (TColStd_HArray1OfReal) URoots =
152   new TColStd_HArray1OfReal(1,myNbURoot);
153 Standard_Real *U_array =
154   (Standard_Real *) &URoots->ChangeArray1()(URoots->Lower());
155 Handle (TColStd_HArray1OfReal) VRoots =
156   new TColStd_HArray1OfReal(1,myNbVRoot);
157 Standard_Real *V_array =
158   (Standard_Real *) &VRoots->ChangeArray1()(VRoots->Lower());
159 AdvApp2Var_ApproxF2var::mma2roo_(&NbPntU,&NbPntV,U_array,V_array);
160 myURoots = URoots;
161 myVRoots = VRoots;
162
163 // myUGauss
164 size = (NbPntU/2+1)*(myJDegU-2*iu-1);
165 Handle (TColStd_HArray1OfReal) UGauss =
166   new TColStd_HArray1OfReal(1,size);
167 Standard_Real *UG_array =
168   (Standard_Real *) &UGauss->ChangeArray1()(UGauss->Lower());
169 AdvApp2Var_ApproxF2var::mmapptt_(&JDegU,&NbPntU,&iu,UG_array,&ErrorCode);
170 if (ErrorCode != 0 ) {
171   throw Standard_ConstructionError("AdvApp2Var_Context : Error in FORTRAN");
172 }
173 myUGauss = UGauss;
174
175 // myVGauss
176 size = (NbPntV/2+1)*(myJDegV-2*iv-1);
177 Handle (TColStd_HArray1OfReal) VGauss =
178   new TColStd_HArray1OfReal(1,size);
179 Standard_Real *VG_array =
180   (Standard_Real *) &VGauss->ChangeArray1()(VGauss->Lower());
181 AdvApp2Var_ApproxF2var::mmapptt_(&JDegV,&NbPntV,&iv,VG_array,&ErrorCode);
182 if (ErrorCode != 0 ) {
183   throw Standard_ConstructionError("AdvApp2Var_Context : Error in FORTRAN");
184 }
185 myVGauss = VGauss;
186
187 //myInternalTol, myFrontierTol, myCuttingTol
188 Standard_Integer nbss = nb1Dss + nb2Dss + nb3Dss;
189 Handle (TColStd_HArray1OfReal) ITol =
190   new TColStd_HArray1OfReal(1,nbss);
191 for (i=1;i<=nb1Dss;i++) {
192   ITol->SetValue(i,tol1D->Value(i));
193 }
194 for (i=1;i<=nb2Dss;i++) {
195   ITol->SetValue(i+nb1Dss,tol2D->Value(i));
196 }
197 for (i=1;i<=nb3Dss;i++) {
198   ITol->SetValue(i+nb1Dss+nb2Dss,tol3D->Value(i));
199 }
200 if (iu>-1||iv>-1) {
201   for (i=1;i<=nbss;i++) { ITol->SetValue(i,ITol->Value(i)/2); }
202 }
203 Handle (TColStd_HArray2OfReal) FTol =
204   new TColStd_HArray2OfReal(1,nbss,1,4);
205 Handle (TColStd_HArray2OfReal) CTol =
206   new TColStd_HArray2OfReal(1,nbss,1,4);
207 for (i=1;i<=nb1Dss;i++) {
208 for (j=1;j<=4;j++) {
209   FTol->SetValue(i,j,tof1D->Value(i,j));
210   CTol->SetValue(i,j,0);
211 }
212 }
213 for (i=1;i<=nb2Dss;i++) {
214 for (j=1;j<=4;j++) {
215   FTol->SetValue(nb1Dss+i,j,tof2D->Value(i,j));
216   CTol->SetValue(nb1Dss+i,j,0);
217 }
218 }
219 for (i=1;i<=nb3Dss;i++) {
220 for (j=1;j<=4;j++) {
221   FTol->SetValue(nb1Dss+nb2Dss+i,j,tof3D->Value(i,j));
222   CTol->SetValue(nb1Dss+nb2Dss+i,j,0);
223 }
224 }
225 if (iu>-1||iv>-1) {
226   Standard_Real tolmin, poids, hmax[4];
227   hmax[0] = 0;
228   hmax[1] = 1;
229   hmax[2] = 1.5;
230   hmax[3] = 1.75;
231   poids = hmax[iu+1]*hmax[iv+1] + hmax[iu+1] + hmax[iv+1];
232   for (i=1;i<=nbss;i++) {
233     for (j=1;j<=4;j++) {
234       tolmin = (ITol->Value(i))/poids;
235       if (tolmin<FTol->Value(i,j)) FTol->SetValue(i,j,tolmin);
236       CTol->SetValue(i,j,tolmin);
237     }
238   }
239 }
240 myInternalTol = ITol;
241 myFrontierTol = FTol;
242 myCuttingTol = CTol;
243 }
244
245 //============================================================================
246 //function : TotalDimension
247 //purpose  : 
248 //============================================================================
249
250 Standard_Integer AdvApp2Var_Context::TotalDimension() const 
251 {
252   return myNb1DSS + 2*myNb2DSS + 3*myNb3DSS; 
253 }
254
255 //============================================================================
256 //function : TotalNumberSSP
257 //purpose  :
258 //============================================================================
259
260 Standard_Integer AdvApp2Var_Context::TotalNumberSSP() const 
261 {
262   return myNb1DSS + myNb2DSS + myNb3DSS; 
263 }
264
265 //============================================================================
266 //function : FavorIso
267 //purpose  : return 1 for IsoU, 2 for IsoV, 2 by default
268 //============================================================================
269
270 Standard_Integer AdvApp2Var_Context::FavorIso() const 
271 {
272   return myFav; 
273 }
274
275 //============================================================================
276 //function : UOrder
277 //purpose  : return the order of continuity requested in U
278 //============================================================================
279
280 Standard_Integer AdvApp2Var_Context::UOrder() const 
281 {
282   return myOrdU; 
283 }
284
285 //============================================================================
286 //function : VOrder
287 //purpose  : return the order of continuity requested in V
288 //============================================================================
289
290 Standard_Integer AdvApp2Var_Context::VOrder() const 
291 {
292   return myOrdV; 
293 }
294
295 //============================================================================
296 //function : ULimit
297 //purpose  : return the max number of coeff. in U of the polynomial approx.
298 //============================================================================
299
300 Standard_Integer AdvApp2Var_Context::ULimit() const 
301 {
302   return myLimU; 
303 }
304
305 //============================================================================
306 //function : VLimit
307 //purpose  : return the max number of coeff. in V of the polynomial approx.
308 //============================================================================
309
310 Standard_Integer AdvApp2Var_Context::VLimit() const 
311 {
312   return myLimV;
313 }
314
315 //============================================================================
316 //function : UJacDeg
317 //purpose  : return the max degree of the Jacobi functions for U parameter
318 //============================================================================
319
320 Standard_Integer AdvApp2Var_Context::UJacDeg() const 
321 {
322   return myJDegU; 
323 }
324
325 //============================================================================
326 //function : VJacDeg
327 //purpose  : return the max degree of the Jacobi functions for V parameter
328 //============================================================================
329
330 Standard_Integer AdvApp2Var_Context::VJacDeg() const 
331 {
332   return myJDegV; 
333 }
334
335 //============================================================================
336 //function : UJacMax
337 //purpose  : return the max value of the Jacobi functions for U parameter
338 //============================================================================
339
340 Handle(TColStd_HArray1OfReal) AdvApp2Var_Context::UJacMax() const 
341 {
342   return myJMaxU;
343 }
344
345 //============================================================================
346 //function : VJacMax
347 //purpose  : return the max value of the Jacobi functions for V parameter
348 //============================================================================
349
350 Handle(TColStd_HArray1OfReal) AdvApp2Var_Context::VJacMax() const 
351 {
352   return myJMaxV;
353 }
354
355 //============================================================================
356 //function : URoots
357 //purpose  : return Legendre roots for U parameter
358 //============================================================================
359
360 Handle(TColStd_HArray1OfReal) AdvApp2Var_Context::URoots() const 
361 {
362   return myURoots;
363 }
364
365 //============================================================================
366 //function : VRoots
367 //purpose  : return Legendre roots for V parameter
368 //============================================================================
369
370 Handle(TColStd_HArray1OfReal) AdvApp2Var_Context::VRoots() const 
371 {
372   return myVRoots;
373 }
374
375 //============================================================================
376 //function : UGauss
377 //purpose  : return Gauss roots for U parameter
378 //============================================================================
379
380 Handle(TColStd_HArray1OfReal) AdvApp2Var_Context::UGauss() const 
381 {
382   return myUGauss;
383 }
384
385 //============================================================================
386 //function : VGauss
387 //purpose  : return Gauss roots for V parameter
388 //============================================================================
389
390 Handle(TColStd_HArray1OfReal) AdvApp2Var_Context::VGauss() const 
391 {
392   return myVGauss;
393 }
394
395 //============================================================================
396 //function : IToler
397 //purpose  : return tolerances for the approximation of patches
398 //============================================================================
399
400 Handle(TColStd_HArray1OfReal) AdvApp2Var_Context::IToler() const 
401 {
402   return myInternalTol;
403 }
404
405 //============================================================================
406 //function : FToler
407 //purpose  : return tolerances for the approximation of frontiers
408 //============================================================================
409
410 Handle(TColStd_HArray2OfReal) AdvApp2Var_Context::FToler() const 
411 {
412   return myFrontierTol;
413 }
414
415 //============================================================================
416 //function : CToler
417 //purpose  : return tolerances for the approximation of cutting lines
418 //============================================================================
419
420 Handle(TColStd_HArray2OfReal) AdvApp2Var_Context::CToler() const 
421 {
422   return myCuttingTol;
423 }
424