0031035: Coding - uninitialized class fields reported by Visual Studio Code Analysis
[occt.git] / src / ProjLib / ProjLib_ProjectOnSurface.cxx
1 // Created on: 1994-09-15
2 // Created by: Bruno DUMORTIER
3 // Copyright (c) 1994-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 <Adaptor3d_HCurve.hxx>
19 #include <Adaptor3d_HSurface.hxx>
20 #include <AppCont_Function.hxx>
21 #include <AppParCurves_MultiCurve.hxx>
22 #include <Approx_FitAndDivide.hxx>
23 #include <BSplCLib.hxx>
24 #include <Extrema_ExtPS.hxx>
25 #include <Extrema_POnSurf.hxx>
26 #include <Geom_BSplineCurve.hxx>
27 #include <PLib.hxx>
28 #include <Precision.hxx>
29 #include <ProjLib_ProjectOnSurface.hxx>
30 #include <Standard_NoSuchObject.hxx>
31 #include <TColgp_Array1OfPnt.hxx>
32 #include <TColStd_Array1OfInteger.hxx>
33 #include <TColStd_Array1OfReal.hxx>
34
35 //=======================================================================
36 //function : OnSurface_Value
37 //purpose  : Evaluate current point of the projected curve
38 //=======================================================================
39 static gp_Pnt OnSurface_Value(const Standard_Real U,
40                               const Handle(Adaptor3d_HCurve)& myCurve,
41                               Extrema_ExtPS * myExtPS)
42 {
43   // on essaie de rendre le point solution le plus proche.
44   myExtPS->Perform(myCurve->Value(U));
45   
46   Standard_Real Dist2Min = RealLast();
47   Standard_Integer Index = 0;
48   
49   for ( Standard_Integer i = 1; i <= myExtPS->NbExt(); i++) {
50     if ( myExtPS->SquareDistance(i) < Dist2Min) {
51       Index = i;
52       Dist2Min = myExtPS->SquareDistance(Index);
53     }
54   }
55   if ( Index == 0 ) {
56     std::cout << " Extrema non trouve pour U = " << U << std::endl;
57     return gp_Pnt(0.,0.,0.);
58   }
59   else {
60     return (myExtPS->Point(Index)).Value();
61   }
62 }
63
64 //=======================================================================
65 //function : OnSurface_D1
66 //purpose  : 
67 //=======================================================================
68
69 static Standard_Boolean OnSurface_D1(const Standard_Real , // U,
70                                      gp_Pnt& ,       // P,
71                                      gp_Vec&       , // V,
72                                      const  Handle(Adaptor3d_HCurve)& , //  myCurve,
73                                      Extrema_ExtPS *) // myExtPS)
74 {
75   return Standard_False;
76 }
77
78
79 //=======================================================================
80 //  class  : ProjLib_OnSurface
81 //purpose  : Use to approximate the projection on a plane
82 //=======================================================================
83
84 class ProjLib_OnSurface : public AppCont_Function
85
86 {
87 public:
88
89   ProjLib_OnSurface(const Handle(Adaptor3d_HCurve)   & C, 
90                     const Handle(Adaptor3d_HSurface) & S)
91  : myCurve(C)
92   {
93     myNbPnt = 1;
94     myNbPnt2d = 0;
95     Standard_Real U = myCurve->FirstParameter();
96     gp_Pnt P = myCurve->Value(U);
97     Standard_Real Tol = Precision::PConfusion();
98     myExtPS = new Extrema_ExtPS(P,S->Surface(),Tol,Tol);
99   }
100
101   ~ProjLib_OnSurface() { delete myExtPS; }
102
103   Standard_Real FirstParameter() const
104     {return myCurve->FirstParameter();}
105
106   Standard_Real LastParameter() const
107     {return myCurve->LastParameter();}
108
109   Standard_Boolean Value(const Standard_Real   theT,
110                          NCollection_Array1<gp_Pnt2d>& /*thePnt2d*/,
111                          NCollection_Array1<gp_Pnt>&   thePnt) const
112   {
113       thePnt(1) = OnSurface_Value(theT, myCurve, myExtPS);
114       return Standard_True;
115   }
116
117   Standard_Boolean D1(const Standard_Real   theT,
118                       NCollection_Array1<gp_Vec2d>& /*theVec2d*/,
119                       NCollection_Array1<gp_Vec>&   theVec) const
120   {
121     gp_Pnt aPnt;
122     return OnSurface_D1(theT, aPnt, theVec(1), myCurve, myExtPS);
123   }
124
125 private:
126   ProjLib_OnSurface (const ProjLib_OnSurface&);
127   ProjLib_OnSurface& operator= (const ProjLib_OnSurface&);
128
129 private:
130   Handle(Adaptor3d_HCurve)       myCurve;
131   Extrema_ExtPS*                 myExtPS;
132 };
133
134
135 //=====================================================================//
136 //                                                                     //
137 //  D E S C R I P T I O N   O F   T H E   C L A S S  :                 // 
138 //                                                                     //
139 //         P r o j L i b _ A p p r o x P r o j e c t O n P l a n e     //
140 //                                                                     //
141 //=====================================================================//
142
143
144 //=======================================================================
145 //function : ProjLib_ProjectOnSurface
146 //purpose  : 
147 //=======================================================================
148
149 ProjLib_ProjectOnSurface::ProjLib_ProjectOnSurface() :
150 myTolerance(0.0),
151 myIsDone(Standard_False)
152 {
153 }
154
155 //=======================================================================
156 //function : ProjLib_ProjectOnSurface
157 //purpose  : 
158 //=======================================================================
159
160 ProjLib_ProjectOnSurface::ProjLib_ProjectOnSurface
161 (const Handle(Adaptor3d_HSurface)& S ) :
162 myTolerance(0.0),
163 myIsDone(Standard_False)
164 {
165   mySurface = S;
166 }
167
168 void ProjLib_ProjectOnSurface::Load(const Handle(Adaptor3d_HCurve)& C,
169                                     const Standard_Real  Tolerance) 
170 {
171   myTolerance = Tolerance ;
172   myCurve = C;
173   myIsDone = Standard_False ; 
174   if (!mySurface.IsNull()) { 
175       
176     ProjLib_OnSurface F(myCurve, mySurface);
177
178     Standard_Integer Deg1, Deg2;
179     Deg1 = 8; Deg2 = 8;
180     
181     Approx_FitAndDivide Fit(F,Deg1,Deg2,Precision::Approximation(),
182                             Precision::PApproximation(),Standard_True);
183     Standard_Integer i;
184     Standard_Integer NbCurves = Fit.NbMultiCurves();
185     Standard_Integer MaxDeg = 0;
186     
187     // Pour transformer la MultiCurve en BSpline, il faut que toutes 
188     // les Bezier la constituant aient le meme degre -> Calcul de MaxDeg
189     Standard_Integer NbPoles  = 1;
190     for (i = 1; i <= NbCurves; i++) {
191       Standard_Integer Deg = Fit.Value(i).Degree();
192       MaxDeg = Max ( MaxDeg, Deg);
193     }
194     NbPoles = MaxDeg * NbCurves + 1;               //Poles sur la BSpline
195     TColgp_Array1OfPnt  Poles( 1, NbPoles);
196     
197     TColgp_Array1OfPnt TempPoles( 1, MaxDeg + 1);  //pour augmentation du degre
198     
199     TColStd_Array1OfReal Knots( 1, NbCurves + 1);  //Noeuds de la BSpline
200     
201     Standard_Integer Compt = 1;
202     for (i = 1; i <= Fit.NbMultiCurves(); i++) {
203       Fit.Parameters(i, Knots(i), Knots(i+1)); 
204       
205       AppParCurves_MultiCurve MC = Fit.Value( i);   //Charge la Ieme Curve
206       TColgp_Array1OfPnt LocalPoles( 1, MC.Degree() + 1);//Recupere les poles
207       MC.Curve(1, Poles);
208       
209       //Augmentation eventuelle du degre
210       Standard_Integer Inc = MaxDeg - MC.Degree();
211       if ( Inc > 0) {
212         BSplCLib::IncreaseDegree( Inc, LocalPoles, BSplCLib::NoWeights(),
213                                  TempPoles, BSplCLib::NoWeights());
214           //mise a jour des poles de la PCurve
215           for (Standard_Integer j = 1 ; j <= MaxDeg + 1; j++) {
216             Poles.SetValue( Compt, TempPoles( j));
217             Compt++;
218           }
219         }
220       else {
221         //mise a jour des poles de la PCurve
222         for (Standard_Integer j = 1 ; j <= MaxDeg + 1; j++) {
223           Poles.SetValue( Compt, LocalPoles( j));
224           Compt++;
225         }
226       } 
227       
228       Compt--;
229     }
230     
231     //mise a jour des fields de ProjLib_Approx
232     
233     Standard_Integer NbKnots = NbCurves + 1;
234     
235     TColStd_Array1OfInteger  Mults( 1, NbKnots);
236     Mults.SetValue( 1, MaxDeg + 1);
237     for ( i = 2; i <= NbCurves; i++) {
238       Mults.SetValue( i, MaxDeg);
239     }
240     Mults.SetValue(NbKnots, MaxDeg + 1);
241     myResult = 
242       new Geom_BSplineCurve(Poles,
243                             Knots,
244                             Mults,
245                             MaxDeg,
246                             Standard_False) ;
247     myIsDone = Standard_True ;
248   }
249 }
250
251 //=======================================================================
252 //function : ~ProjLib_ProjectOnSurface
253 //purpose  : 
254 //=======================================================================
255
256 ProjLib_ProjectOnSurface::~ProjLib_ProjectOnSurface()
257 {}
258
259 //=======================================================================
260 //function : BSpline
261 //purpose  : 
262 //=======================================================================
263
264 Handle(Geom_BSplineCurve) ProjLib_ProjectOnSurface::BSpline() const 
265 {
266   Standard_NoSuchObject_Raise_if
267     (!myIsDone,
268      "ProjLib_ProjectOnSurface:BSpline");
269   return myResult ;
270 }
271
272 //=======================================================================
273 //function : IsDone
274 //purpose  : 
275 //=======================================================================
276
277 Standard_Boolean ProjLib_ProjectOnSurface::IsDone() const 
278 {
279   return myIsDone;
280 }