MAGiC  V5.0
Mailleurs Automatiques de Géometries intégrés à la Cao
sld_surface.cpp
Aller à la documentation de ce fichier.
1 //####//------------------------------------------------------------
2 //####//------------------------------------------------------------
3 //####// MAGiC
4 //####// Jean Christophe Cuilliere et Vincent FRANCOIS
5 //####// Departement de Genie Mecanique - UQTR
6 //####//------------------------------------------------------------
7 //####// MAGIC est un projet de recherche de l equipe ERICCA
8 //####// du departement de genie mecanique de l Universite du Quebec a Trois Rivieres
9 //####// http://www.uqtr.ca/ericca
10 //####// http://www.uqtr.ca/
11 //####//------------------------------------------------------------
12 //####//------------------------------------------------------------
13 //####//
14 //####// sld_surface.cpp
15 //####//
16 //####//------------------------------------------------------------
17 //####//------------------------------------------------------------
18 //####// COPYRIGHT 2000-2024
19 //####// jeu 13 jun 2024 11:58:54 EDT
20 //####//------------------------------------------------------------
21 //####//------------------------------------------------------------
22 
23 
24 #include "gestionversion.h"
25 
26 #ifdef BREP_SLD
27 
28 #include "mg_surface.h"
29 #include "sld_surface.h"
30 #include <system.hpp>
31 #include <atl\atlmod.h>
32 #include "smartvars.h"
33 #include "sld_fonction.h"
34 #include "stplane.h"
35 #include "stcylindrical.h"
36 #include "stconical.h"
37 #include "stspherical.h"
38 #include "sttoroidal.h"
39 #include "stbsplines.h"
40 
41 
42 
43 SLD_SURFACE::SLD_SURFACE(unsigned long num,string idface,SLD_FONCTION& fonc):MG_SURFACE(num),fonction(fonc)
44 {
45  initialisation(idface);
46 }
47 
48 
49 SLD_SURFACE::SLD_SURFACE(string idface,SLD_FONCTION& fonc):MG_SURFACE(),fonction(fonc)
50 {
51  initialisation(idface);
52 }
53 
54 SLD_SURFACE::SLD_SURFACE(string idface, SLD_FONCTION& fonc, CComPtr < IFace2 > & p_swFace) : MG_SURFACE(), swFace(p_swFace), fonction(fonc)
55 {
56  swFace->IGetSurface(&swSurface);
57  initialisation();
58 }
59 
60 void SLD_SURFACE::initialisation(string idface)
61 {
62  idoriginal=idface;
63  fonction.GetParID((char*)idface.c_str(),swFace);
64  swFace->IGetSurface(&swSurface);
65  initialisation();
66 }
67 
68 void SLD_SURFACE::initialisation()
69 {
70  CComVariant vparam;
71  swSurface->Parameterization( &vparam);
72  SafeDoubleArray saparam(vparam);
73  u_min=saparam[0];
74  u_max=saparam[1];
75  v_min=saparam[2];
76  v_max=saparam[3];
77  stsurface=NULL;
78  VARIANT_BOOL retval;
79  swSurface->IsPlane(&retval);
80  if (retval==1)
81  {
82  CComVariant vRetval;
83  vRetval=swSurface->get_PlaneParams();
84  SafeDoubleArray sda(vRetval);
85  double origine[3];
86  double dir[3];
87  origine[0]=sda[3];
88  origine[1]=sda[4];
89  origine[2]=sda[5];
90  dir[0]=sda[0];
91  dir[1]=sda[1];
92  dir[2]=sda[2];
93  stsurface=new ST_PLANE(origine,dir);
94  return;
95  }
96  swSurface->IsCylinder(&retval);
97  if (retval==1)
98  {
99  CComVariant vRetval;
100  vRetval=swSurface->get_CylinderParams();
101  SafeDoubleArray sda(vRetval);
102  double origine[3];
103  double dir[3];
104  origine[0]=sda[0];
105  origine[1]=sda[1];
106  origine[2]=sda[2];
107  dir[0]=sda[3];
108  dir[1]=sda[4];
109  dir[2]=sda[5];
110  double rayon=sda[6];
111  stsurface=new ST_CYLINDRICAL(origine,dir,NULL,rayon);
112  return;
113  }
114  swSurface->IsCone(&retval);
115  if (retval==1)
116  {
117  CComVariant vRetval;
118  vRetval=swSurface->get_ConeParams();
119  SafeDoubleArray sda(vRetval);
120  double origine[3];
121  double dir[3];
122  origine[0]=sda[0];
123  origine[1]=sda[1];
124  origine[2]=sda[2];
125  dir[0]=-sda[3];
126  dir[1]=-sda[4];
127  dir[2]=-sda[5];
128  double rayon=sda[6];
129  double angle=sda[7];
130  stsurface=new ST_CONICAL(origine,dir,NULL,rayon,angle);
131  return;
132  }
133  swSurface->IsSphere(&retval);
134  if (retval==1)
135  {
136  CComVariant vRetval;
137  vRetval=swSurface->get_SphereParams() ;
138  SafeDoubleArray sda(vRetval);
139  double origine[3];
140  origine[0]=sda[0];
141  origine[1]=sda[1];
142  origine[2]=sda[2];
143  double rayon=sda[3];
144  stsurface=new ST_SPHERICAL(origine,NULL,NULL,rayon);
145  return;
146  }
147  swSurface->IsTorus(&retval) ;
148  if (retval==1)
149  {
150  CComVariant vRetval;
151  vRetval=swSurface->get_TorusParams();
152  SafeDoubleArray sda(vRetval);
153  double origine[3];
154  double dir[3];
155  origine[0]=sda[0];
156  origine[1]=sda[1];
157  origine[2]=sda[2];
158  dir[0]=sda[3];
159  dir[1]=sda[4];
160  dir[2]=sda[5];
161  double grayon=sda[6];
162  double prayon=sda[7];
163  stsurface=new ST_TOROIDAL(origine,dir,grayon,prayon);
164  return;
165  }
166  VARIANT vRetval;
167  VARIANT_BOOL sense;
168  VARIANT vp0;
169  vp0=swFace->GetUVBounds();
170  double tolerance = fonction.GetValeurTolerance() ;
171  swSurface->GetBSurfParams2(false,false,vp0,tolerance,&sense,&vRetval);
172  SafeDoubleArray sda(vRetval);
173  int uordre,vordre;
174  int numcpt,numlpt;
175  int uperiodic,vperiodic;
176  int dim,dummy;
177  OPERATEUR::doubleto2int(sda[0],uordre,vordre);
178  OPERATEUR::doubleto2int(sda[1],numcpt,numlpt);
179  OPERATEUR::doubleto2int(sda[2],uperiodic,vperiodic);
180  OPERATEUR::doubleto2int(sda[3],dim,dummy);
181  vector<double> knots_u;
182  vector<double> knots_v;
183  vector<double> poids;
184  vector<double> ptsctr;
185  for (int i=0;i<uordre+numcpt;i++)
186  knots_u.insert(knots_u.end(),sda[4+i]);
187  for (int i=0;i<vordre+numlpt;i++)
188  knots_v.insert(knots_v.end(),sda[4+uordre+numcpt+i]);
189  for (int i=0;i<numcpt;i++)
190  for (int j=0;j<numlpt;j++)
191  {
192  int num=j*numcpt+i;
193  ptsctr.insert(ptsctr.end(),sda[4+uordre+numcpt+vordre+numlpt+dim*num]);
194  ptsctr.insert(ptsctr.end(),sda[4+uordre+numcpt+vordre+numlpt+dim*num+1]);
195  ptsctr.insert(ptsctr.end(),sda[4+uordre+numcpt+vordre+numlpt+dim*num+2]);
196  if (dim==4) poids.insert(poids.end(),sda[4+uordre+numcpt+vordre+numlpt+dim*num+3]);
197  else poids.insert(poids.end(),1.);
198  }
199  int sens;
200  sens=1;
201  stsurface=new ST_B_SPLINE_SURF(uordre-1,vordre-1,knots_u,knots_v,ptsctr,poids,sens);
202  ((ST_B_SPLINE_SURF*) stsurface)->set_periodique_u( uperiodic );
203  ((ST_B_SPLINE_SURF*) stsurface)->set_periodique_v( vperiodic );
204 };
205 
206 
207 SLD_SURFACE::~SLD_SURFACE()
208 {
209  if (swFace.p != NULL) swFace.p = NULL;
210  if (swSurface.p != NULL) swSurface.p = NULL;
211  if (stsurface!=NULL) delete stsurface;
212 }
213 
214 void SLD_SURFACE::evaluer(double *uv,double *xyz)
215 {
216  if (stsurface!=NULL)
217  {
218  stsurface->evaluer(uv,xyz);
219  return;
220  }
221  double UParam=uv[0];
222  double VParam=uv[1];
223  long NumUDeriv=0;
224  long NumVDeriv=0;
225  CComVariant vretval;
226 
227  swSurface->Evaluate ( UParam, VParam, NumUDeriv, NumVDeriv, &vretval);
228  SafeDoubleArray saretval(vretval);
229  xyz[0]=saretval[0];
230  xyz[1]=saretval[1];
231  xyz[2]=saretval[2];
232 }
233 
234 void SLD_SURFACE::deriver(double *uv,double *xyzdu, double *xyzdv)
235 {
236  if (stsurface!=NULL)
237  {
238  stsurface->deriver(uv,xyzdu,xyzdv);
239  return;
240  }
241  double UParam=uv[0];
242  double VParam=uv[1];
243  long NumUDeriv=1;
244  long NumVDeriv=1;
245  CComVariant vretval;
246 
247  swSurface->Evaluate ( UParam, VParam, NumUDeriv, NumVDeriv, &vretval);
248  SafeDoubleArray saretval(vretval);
249  xyzdu[0]=saretval[3];
250  xyzdu[1]=saretval[4];
251  xyzdu[2]=saretval[5];
252  xyzdv[0]=saretval[6];
253  xyzdv[1]=saretval[7];
254  xyzdv[2]=saretval[8];
255 }
256 
257 void SLD_SURFACE::deriver_seconde(double *uv,double* xyzduu,double* xyzduv,double* xyzdvv,double *xyz, double *xyzdu, double *xyzdv)
258 {
259  if (stsurface!=NULL)
260  {
261  stsurface->deriver_seconde(uv,xyzduu,xyzduv,xyzdvv,xyz,xyzdu,xyzdv);
262  return;
263  }
264  double UParam=uv[0];
265  double VParam=uv[1];
266  long NumUDeriv=2;
267  long NumVDeriv=2;
268  CComVariant vretval;
269 
270  swSurface->Evaluate ( UParam, VParam, NumUDeriv, NumVDeriv, &vretval);
271  SafeDoubleArray saretval(vretval);
272  if (xyz!=NULL) //or xyz~=NULL
273  {
274  xyz[0]=saretval[0];
275  xyz[1]=saretval[1];
276  xyz[2]=saretval[2];
277  }
278  if (xyzdu!= NULL)
279  {
280  xyzdu[0]=saretval[3];
281  xyzdu[1]=saretval[4];
282  xyzdu[2]=saretval[5];
283  }
284  if (xyzdv!=NULL)
285  {
286  xyzdv[0]=saretval[9];
287  xyzdv[1]=saretval[10];
288  xyzdv[2]=saretval[11];
289  }
290  xyzduu[0]=saretval[6];
291  xyzduu[1]=saretval[7];
292  xyzduu[2]=saretval[8];
293  xyzduv[0]=saretval[12];
294  xyzduv[1]=saretval[13];
295  xyzduv[2]=saretval[14];
296  xyzdvv[0]=saretval[18];
297  xyzdvv[1]=saretval[19];
298  xyzdvv[2]=saretval[20];
299 }
300 
301 void SLD_SURFACE::inverser(double* uv,double *xyz,double precision)
302 {//Param�tre pr�cision
303  if (stsurface!=NULL)
304  {
305  stsurface->inverser(uv,xyz,precision);
306  return;
307  }
308  CComVariant vUV;
309  swSurface->ReverseEvaluate(xyz[0], xyz[1], xyz[2], &vUV);
310  SafeDoubleArray saUV(vUV);
311  uv[0] = saUV[0];
312  uv[1] = saUV[1];
313 }
314 
315 int SLD_SURFACE::est_periodique_u(void)
316 {
317  return stsurface->est_periodique_u();
318 }
319 
320 int SLD_SURFACE::est_periodique_v(void)
321 {
322  return stsurface->est_periodique_v();
323 }
324 
325 double SLD_SURFACE::get_periode_u(void)
326 {
327  return stsurface->get_periode_u();
328 }
329 
330 double SLD_SURFACE::get_periode_v(void)
331 {
332  return stsurface->get_periode_v();
333 }
334 
335 double SLD_SURFACE::get_umin()
336 {
337  if (stsurface != NULL)
338  return stsurface->get_umin();
339  return u_min;
340 }
341 
342 double SLD_SURFACE::get_umax()
343 {
344  if (stsurface != NULL)
345  return stsurface->get_umax();
346  return u_max;
347 }
348 
349 double SLD_SURFACE::get_vmin()
350 {
351  if (stsurface != NULL)
352  return stsurface->get_vmin();
353  return v_min;
354 }
355 
356 double SLD_SURFACE::get_vmax()
357 {
358  if (stsurface != NULL)
359  return stsurface->get_vmax();
360  return v_max;
361 }
362 
363 void SLD_SURFACE::enregistrer(ostream& o,double version)
364 {
365  o << "%" << get_id() << "=SURFACE_SLD("<< idoriginal.c_str() << ");" << endl;
366 }
367 
368 
369 
370 int SLD_SURFACE::get_type_geometrique(TPL_LISTE_ENTITE<double> &param)
371 {
372  if (stsurface!=NULL)
373  return stsurface->get_type_geometrique(param);;
374  return 0;
375 }
376 
377 
378 
379 void SLD_SURFACE::get_param_NURBS(int& indx_premier_ptctr,TPL_LISTE_ENTITE<double> &param)
380 {
381  stsurface->get_param_NURBS(indx_premier_ptctr,param);
382 }
383 
384 
385 #endif
386 
387 
388 
389 
SLD_SURFACE::deriver_seconde
virtual void deriver_seconde(double *uv, double *xyzduu, double *xyzduv, double *xyzdvv, double *xyz=NULL, double *xyzdu=NULL, double *xyzdv=NULL)
Definition: sld_surface.cpp:257
ST_SURFACE::get_periode_u
virtual double get_periode_u(void)=0
sld_surface.h
SLD_SURFACE::get_umin
virtual double get_umin()
Definition: sld_surface.cpp:335
SLD_SURFACE::get_type_geometrique
virtual int get_type_geometrique(TPL_LISTE_ENTITE< double > &param)
Definition: sld_surface.cpp:370
stspherical.h
gestionversion.h
MG_SURFACE
Definition: mg_surface.h:31
stconical.h
ST_SURFACE::est_periodique_v
virtual int est_periodique_v(void)=0
ST_SURFACE::get_type_geometrique
virtual int get_type_geometrique(TPL_LISTE_ENTITE< double > &param)=0
mg_surface.h
SLD_SURFACE::fonction
SLD_FONCTION & fonction
Definition: sld_surface.h:66
SLD_FONCTION::GetValeurTolerance
double GetValeurTolerance()
Definition: sld_fonction.cpp:812
MG_IDENTIFICATEUR::get_id
unsigned long get_id()
Definition: mg_identificateur.cpp:53
SLD_SURFACE::get_umax
virtual double get_umax()
Definition: sld_surface.cpp:342
ST_SURFACE::evaluer
virtual void evaluer(double *uv, double *xyz)=0
SLD_FONCTION
Definition: sld_fonction.h:37
SLD_SURFACE::inverser
virtual void inverser(double *uv, double *xyz, double precision=1e-6)
Definition: sld_surface.cpp:301
ST_B_SPLINE_SURF
Definition: stbsplines.h:34
SLD_SURFACE::get_periode_u
virtual double get_periode_u(void)
Definition: sld_surface.cpp:325
ST_SURFACE::get_param_NURBS
virtual void get_param_NURBS(int &indx_premier_ptctr, TPL_LISTE_ENTITE< double > &param)=0
SLD_SURFACE::get_vmin
virtual double get_vmin()
Definition: sld_surface.cpp:349
SLD_SURFACE::evaluer
virtual void evaluer(double *uv, double *xyz)
Definition: sld_surface.cpp:214
ST_SURFACE::get_vmax
virtual double get_vmax(void)=0
sld_fonction.h
ST_SURFACE::deriver_seconde
virtual void deriver_seconde(double *uv, double *xyzduu, double *xyzduv, double *xyzdvv, double *xyz=NULL, double *xyzdu=NULL, double *xyzdv=NULL)=0
MG_SURFACE::u_max
double u_max
Definition: mg_surface.h:64
MG_SURFACE::v_min
double v_min
Definition: mg_surface.h:65
SLD_SURFACE::~SLD_SURFACE
virtual ~SLD_SURFACE()
Definition: sld_surface.cpp:207
MG_SURFACE::u_min
double u_min
Definition: mg_surface.h:63
SLD_SURFACE::stsurface
class ST_SURFACE * stsurface
Definition: sld_surface.h:67
ST_CONICAL
Definition: stconical.h:33
ST_SURFACE::get_vmin
virtual double get_vmin(void)=0
SLD_SURFACE::est_periodique_u
virtual int est_periodique_u(void)
Definition: sld_surface.cpp:315
ST_SURFACE::get_umax
virtual double get_umax(void)=0
stcylindrical.h
ST_SURFACE::get_umin
virtual double get_umin(void)=0
SLD_SURFACE::initialisation
virtual void initialisation()
Definition: sld_surface.cpp:68
SLD_SURFACE::swSurface
CComPtr< ISurface > swSurface
Definition: sld_surface.h:63
SLD_SURFACE::enregistrer
virtual void enregistrer(std::ostream &o, double version)
Definition: sld_surface.cpp:363
SLD_SURFACE::swFace
CComPtr< IFace2 > swFace
Definition: sld_surface.h:62
SLD_SURFACE::SLD_SURFACE
SLD_SURFACE(unsigned long num, std::string idface, class SLD_FONCTION &fonc)
MG_SURFACE::v_max
double v_max
Definition: mg_surface.h:66
OPERATEUR::doubleto2int
static void doubleto2int(double val, int &val1, int &val2)
Definition: ot_mathematique.cpp:1666
ST_CYLINDRICAL
Definition: stcylindrical.h:33
ST_PLANE
Definition: stplane.h:33
ST_SURFACE::inverser
virtual void inverser(double *uv, double *xyz, double precision=1e-6)=0
SLD_SURFACE::get_param_NURBS
virtual void get_param_NURBS(int &indx_premier_ptctr, TPL_LISTE_ENTITE< double > &param)
Definition: sld_surface.cpp:379
SLD_SURFACE::get_vmax
virtual double get_vmax()
Definition: sld_surface.cpp:356
SLD_SURFACE::est_periodique_v
virtual int est_periodique_v(void)
Definition: sld_surface.cpp:320
ST_SPHERICAL
Definition: stspherical.h:33
sttoroidal.h
TPL_LISTE_ENTITE< double >
ST_SURFACE::est_periodique_u
virtual int est_periodique_u(void)=0
SLD_FONCTION::GetParID
BOOL GetParID(char *ID, CComPtr< IFace2 > &)
Definition: sld_fonction.cpp:381
stplane.h
ST_SURFACE::get_periode_v
virtual double get_periode_v(void)=0
SLD_SURFACE::deriver
virtual void deriver(double *uv, double *xyzdu, double *xyzdv)
Definition: sld_surface.cpp:234
SLD_SURFACE::get_periode_v
virtual double get_periode_v(void)
Definition: sld_surface.cpp:330
ST_SURFACE::deriver
virtual void deriver(double *uv, double *xyzdu, double *xyzdv)=0
SLD_SURFACE::idoriginal
std::string idoriginal
Definition: sld_surface.h:68
stbsplines.h
ST_TOROIDAL
Definition: sttoroidal.h:33