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