0031035: Coding - uninitialized class fields reported by Visual Studio Code Analysis
[occt.git] / src / Extrema / Extrema_FuncPSNorm.cxx
1 // Created on: 1995-07-18
2 // Created by: Modelistation
3 // Copyright (c) 1995-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 <Extrema_FuncPSNorm.hxx>
19 #include <Adaptor3d_Surface.hxx>
20 #include <Extrema_POnSurf.hxx>
21 #include <GeomAbs_IsoType.hxx>
22 #include <gp_Pnt.hxx>
23 #include <gp_Vec.hxx>
24 #include <math_Matrix.hxx>
25 #include <Precision.hxx>
26 #include <Standard_TypeMismatch.hxx>
27
28 Extrema_FuncPSNorm::Extrema_FuncPSNorm ()
29 : myS(NULL),
30   myU(0.0),
31   myV(0.0)
32 {
33   myPinit = Standard_False;
34   mySinit = Standard_False;
35 }
36
37 //=============================================================================
38 Extrema_FuncPSNorm::Extrema_FuncPSNorm (const gp_Pnt& P,
39                                       const Adaptor3d_Surface& S)
40 : myU(0.0),
41   myV(0.0)
42 {
43   myP = P;
44   myS = (Adaptor3d_SurfacePtr)&S;
45   myPinit = Standard_True;
46   mySinit = Standard_True;
47 }
48
49 //=============================================================================
50 void Extrema_FuncPSNorm::Initialize(const Adaptor3d_Surface& S)
51 {
52   myS = (Adaptor3d_SurfacePtr)&S;
53   mySinit = Standard_True;
54   myPoint.Clear();
55   mySqDist.Clear();
56 }
57
58 //=============================================================================
59
60 void Extrema_FuncPSNorm::SetPoint(const gp_Pnt& P)
61 {
62   myP = P;
63   myPinit = Standard_True;
64   myPoint.Clear();
65   mySqDist.Clear();
66 }
67
68 //=============================================================================
69
70 //=============================================================================
71
72 Standard_Integer Extrema_FuncPSNorm::NbVariables () const { return 2;}
73 //=============================================================================
74
75 Standard_Integer Extrema_FuncPSNorm::NbEquations () const { return 2;}
76 //=============================================================================
77
78 Standard_Boolean Extrema_FuncPSNorm::Value (const math_Vector& UV, 
79                                            math_Vector& F)
80 {
81   if (!myPinit || !mySinit) throw Standard_TypeMismatch();
82   myU = UV(1);
83   myV = UV(2);
84   gp_Vec Dus, Dvs;
85   myS->D1(myU,myV,myPs,Dus,Dvs);
86
87   gp_Vec PPs (myP,myPs);
88
89   F(1) = PPs.Dot(Dus);
90   F(2) = PPs.Dot(Dvs);
91
92   return Standard_True;
93 }
94 //=============================================================================
95
96 Standard_Boolean Extrema_FuncPSNorm::Derivatives (const math_Vector& UV, 
97                                                  math_Matrix& Df)
98 {
99   math_Vector F(1,2);
100   return Values(UV,F,Df);
101 }
102 //=============================================================================
103
104 Standard_Boolean Extrema_FuncPSNorm::Values (const math_Vector& UV, 
105                                             math_Vector& F,
106                                             math_Matrix& Df)
107 {
108   if (!myPinit || !mySinit) throw Standard_TypeMismatch();
109   myU = UV(1);
110   myV = UV(2);
111   gp_Vec Dus, Dvs, Duus, Dvvs, Duvs;
112   myS->D2(myU,myV,myPs,Dus,Dvs,Duus,Dvvs,Duvs);
113
114   gp_Vec PPs (myP,myPs);
115
116   Df(1,1) = Dus.SquareMagnitude() + PPs.Dot(Duus);
117   Df(1,2) = Dvs.Dot(Dus)          + PPs.Dot(Duvs);
118   Df(2,1) = Df(1,2);
119   Df(2,2) = Dvs.SquareMagnitude() + PPs.Dot(Dvvs);
120
121   // 3. Value
122   F(1) = PPs.Dot(Dus);
123   F(2) = PPs.Dot(Dvs);
124
125   return Standard_True;
126 }
127 //=============================================================================
128
129 Standard_Integer Extrema_FuncPSNorm::GetStateNumber ()
130 {
131   if (!myPinit || !mySinit) throw Standard_TypeMismatch();
132   //comparison of solution with previous solutions
133   Standard_Integer i = 1, nbSol = mySqDist.Length();
134   Standard_Real tol2d = Precision::PConfusion() * Precision::PConfusion();
135    
136   for( ; i <=  nbSol; i++)
137   {
138     Standard_Real aU, aV;
139         myPoint(i).Parameter(aU, aV);
140         if( ((myU - aU ) * (myU - aU ) + (myV - aV ) * (myV - aV )) <= tol2d )
141       break;
142   }
143   if( i <= nbSol)
144           return 0;
145   mySqDist.Append(myPs.SquareDistance(myP));
146   myPoint.Append(Extrema_POnSurf(myU,myV,myPs));
147   return 0;
148 }
149 //=============================================================================
150
151 Standard_Integer Extrema_FuncPSNorm::NbExt () const
152 {
153   return mySqDist.Length();
154 }
155 //=============================================================================
156
157 Standard_Real Extrema_FuncPSNorm::SquareDistance (const Standard_Integer N) const
158 {
159   if (!myPinit || !mySinit) throw Standard_TypeMismatch();
160   return mySqDist.Value(N);
161 }
162 //=============================================================================
163
164 const Extrema_POnSurf& Extrema_FuncPSNorm::Point (const Standard_Integer N) const
165 {
166   if (!myPinit || !mySinit) throw Standard_TypeMismatch();
167   return myPoint.Value(N);
168 }