ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/REPOS_ERICCA/magic/lib/microstructure/src/mstruct_generateur.cpp
Revision: 951
Committed: Fri Aug 10 15:17:17 2018 UTC (6 years, 9 months ago) by couturad
File size: 44719 byte(s)
Log Message: 
-> Ajout de Project Chrono (voir CMakeLists.txt).
-> Ajout d'un générateur de microstructure basé sur la dynamique des corps rigides (MSTRUCT_GENERATEUR_DCR).
-> Ajout d'un opérateur de décallage de la topologie (MG_CG_OP_TRANSF_DECALLAGE).
-> Retrait de «using namespace std»  (conflit avec namespace chrono) et modification des fichiers affectés.
-> Modification de mailleur2d.cpp afin d'enregistrer un fichier MAGiC (void.magic) lorsque le nombre d'itération dépasse la valeur maximale.

File Contents

# User Rev Content
1 couturad 919 #include "mstruct_generateur.h"
2     #include "mg_cg_modele.h"
3     #include "mg_gestionnaire.h"
4     #include "mg_cg.h"
5 couturad 951 #include "ot_geometrie.h"
6     #include "mg_cg_forme_volume.h"
7     #include "mg_cg_forme_multi_volume.h"
8     #include "mg_cg_op_inclusion_matrice.h"
9     #include "mg_cg_op_bool_difference.h"
10     #include "mg_cg_op_bool_intersection.h"
11     #include "mg_cg_op_transf_decallage.h"
12     #include <stack>
13 couturad 919
14 couturad 951
15     MSTRUCT_GENERATEUR::MSTRUCT_GENERATEUR(MG_GESTIONNAIRE* mg_gestionnaire,char* nom_mgcg_modele)
16 couturad 919 {
17     m_mg_gestionnaire=mg_gestionnaire;
18     m_mgcg_modele=MG_CG::creer_MG_CG_MODELE(m_mg_gestionnaire,nom_mgcg_modele);
19     mg_gestionnaire->ajouter_mgcg_modele(m_mgcg_modele);
20     m_mgcg_modele->change_importer_triangulation(true);
21     m_mgcg_modele->change_fusionner_entite_similaire(false);
22     m_mgcg_modele->change_epsilon_triangulation(0.01);
23     }
24    
25 couturad 951 MSTRUCT_GENERATEUR::~MSTRUCT_GENERATEUR(void)
26 couturad 919 {
27     }
28    
29 couturad 951 MG_GESTIONNAIRE* MSTRUCT_GENERATEUR::get_mg_gestionnaire(void)
30 couturad 919 {
31     return m_mg_gestionnaire;
32     }
33    
34 couturad 951 MG_CG_MODELE* MSTRUCT_GENERATEUR::get_mgcg_modele(void)
35 couturad 919 {
36     return m_mgcg_modele;
37     }
38    
39 couturad 951 void MSTRUCT_GENERATEUR::change_mgcg_modele(MG_CG_MODELE* mgcg_modele)
40 couturad 919 {
41     m_mgcg_modele=mgcg_modele;
42     }
43    
44 couturad 951 void MSTRUCT_GENERATEUR::change_boite3d_ves(BOITE_3D boite3d_ves)
45     {
46     m_boite3d_ves=boite3d_ves;
47     }
48 couturad 919
49 couturad 951 void MSTRUCT_GENERATEUR::change_boite3d_distribution(BOITE_3D boite3d_distribution)
50 couturad 929 {
51 couturad 951 m_boite3d_distribution=boite3d_distribution;
52     }
53    
54     void MSTRUCT_GENERATEUR::change_intersection_bords_ves(bool intersection_bords_ves)
55     {
56     m_intersection_bords_ves=intersection_bords_ves;
57     }
58    
59     void MSTRUCT_GENERATEUR::change_angle_min(double angle_min)
60     {
61     m_angle_min=angle_min;
62     }
63    
64     void MSTRUCT_GENERATEUR::change_longueur_min(double longueur_min)
65     {
66     m_longueur_min=longueur_min;
67     }
68    
69     void MSTRUCT_GENERATEUR::change_aire_min(double aire_min)
70     {
71     m_aire_min=aire_min;
72     }
73    
74     void MSTRUCT_GENERATEUR::change_volume_min(double volume_min)
75     {
76     m_volume_min=volume_min;
77     }
78    
79     void MSTRUCT_GENERATEUR::change_distance_inter_volume_min(double distance_min)
80     {
81     m_distance_min=distance_min;
82     }
83    
84     void MSTRUCT_GENERATEUR::change_nb_pas_grille(int nb_pas_grille)
85     {
86     m_nb_pas_grille=nb_pas_grille;
87     }
88    
89     void MSTRUCT_GENERATEUR::change_nb_pas_echantillonage(int nb_pas_echantillonnage)
90     {
91     m_nb_pas_echantillonnage=nb_pas_echantillonnage;
92     }
93    
94     void MSTRUCT_GENERATEUR::change_precision(double precision)
95     {
96     m_precision=precision;
97     }
98    
99     void MSTRUCT_GENERATEUR::initialiser_grille(void)
100     {
101     if(!m_grille_init)
102     m_grille_volume_inclusion.initialiser(m_boite3d_distribution.get_xmin(),m_boite3d_distribution.get_ymin(),m_boite3d_distribution.get_zmin(),
103     m_boite3d_distribution.get_xmax(),m_boite3d_distribution.get_ymax(),m_boite3d_distribution.get_zmax(),
104     m_nb_pas_echantillonnage,m_nb_pas_echantillonnage,m_nb_pas_echantillonnage);
105     }
106    
107     int MSTRUCT_GENERATEUR::position_relative_inclusion_ves(MG_VOLUME* inclusion,bool &plan_x0yz,bool &plan_x1yz,bool &plan_xy0z,bool &plan_xy1z,bool &plan_xyz0,bool &plan_xyz1)
108     {
109     plan_x0yz=false;plan_x1yz=false;plan_xy0z=false;plan_xy1z=false;plan_xyz0=false;plan_xyz1=false;
110     BOITE_3D boite_inclusion=inclusion->get_boite_3D(m_nb_pas_echantillonnage);
111     int position_relative;
112     if(m_boite3d_ves.contient(boite_inclusion.get_xmin(),boite_inclusion.get_ymin(),boite_inclusion.get_zmin()))
113     if(m_boite3d_ves.contient(boite_inclusion.get_xmin(),boite_inclusion.get_ymin(),boite_inclusion.get_zmax()))
114     if(m_boite3d_ves.contient(boite_inclusion.get_xmin(),boite_inclusion.get_ymax(),boite_inclusion.get_zmin()))
115     if(m_boite3d_ves.contient(boite_inclusion.get_xmin(),boite_inclusion.get_ymax(),boite_inclusion.get_zmax()))
116     if(m_boite3d_ves.contient(boite_inclusion.get_xmax(),boite_inclusion.get_ymin(),boite_inclusion.get_zmin()))
117     if(m_boite3d_ves.contient(boite_inclusion.get_xmax(),boite_inclusion.get_ymin(),boite_inclusion.get_zmax()))
118     if(m_boite3d_ves.contient(boite_inclusion.get_xmax(),boite_inclusion.get_ymax(),boite_inclusion.get_zmin()))
119     if(m_boite3d_ves.contient(boite_inclusion.get_xmax(),boite_inclusion.get_ymax(),boite_inclusion.get_zmax())) return POSITION_RELATIVE::INTERIEUR;
120    
121     if(!m_boite3d_ves.contient(boite_inclusion.get_xmin(),boite_inclusion.get_ymin(),boite_inclusion.get_zmin()))
122     if(!m_boite3d_ves.contient(boite_inclusion.get_xmin(),boite_inclusion.get_ymin(),boite_inclusion.get_zmax()))
123     if(!m_boite3d_ves.contient(boite_inclusion.get_xmin(),boite_inclusion.get_ymax(),boite_inclusion.get_zmin()))
124     if(!m_boite3d_ves.contient(boite_inclusion.get_xmin(),boite_inclusion.get_ymax(),boite_inclusion.get_zmax()))
125     if(!m_boite3d_ves.contient(boite_inclusion.get_xmax(),boite_inclusion.get_ymin(),boite_inclusion.get_zmin()))
126     if(!m_boite3d_ves.contient(boite_inclusion.get_xmax(),boite_inclusion.get_ymin(),boite_inclusion.get_zmax()))
127     if(!m_boite3d_ves.contient(boite_inclusion.get_xmax(),boite_inclusion.get_ymax(),boite_inclusion.get_zmin()))
128     if(!m_boite3d_ves.contient(boite_inclusion.get_xmax(),boite_inclusion.get_ymax(),boite_inclusion.get_zmax())) return POSITION_RELATIVE::EXTERIEUR;
129    
130     if(boite_inclusion.get_xmin()<m_boite3d_ves.get_xmin()) plan_x0yz=true;
131     if(boite_inclusion.get_xmax()>m_boite3d_ves.get_xmax()) plan_x1yz=true;
132     if(boite_inclusion.get_ymin()<m_boite3d_ves.get_ymin()) plan_xy0z=true;
133     if(boite_inclusion.get_ymax()>m_boite3d_ves.get_ymax()) plan_xy0z=true;
134     if(boite_inclusion.get_zmin()<m_boite3d_ves.get_zmin()) plan_xyz0=true;
135     if(boite_inclusion.get_zmax()>m_boite3d_ves.get_zmax()) plan_xyz1=true;
136    
137     return POSITION_RELATIVE::AU_BORD;
138     }
139    
140     bool MSTRUCT_GENERATEUR::angle_entre_face_acceptable(MG_ARETE* mgarete)
141     {
142     long nb_mg_coarete = mgarete->get_nb_mg_coarete();
143     for(int i=0;i<nb_mg_coarete-1;i++)
144     {
145     MG_FACE* face1 = mgarete->get_mg_coarete(i)->get_boucle()->get_mg_face();
146     for(int j=i+1;j<nb_mg_coarete;j++)
147     {
148     MG_FACE* face2 = mgarete->get_mg_coarete(j)->get_boucle()->get_mg_face();
149     double tmin = mgarete->get_tmin();
150     double tmax = mgarete->get_tmax();
151     if(mgarete->get_courbe()->est_periodique())
152     {
153     tmin=0;
154     tmax=mgarete->get_courbe()->get_periode();
155     }
156     double pas_t = (tmin-tmax)/m_nb_pas_echantillonnage;
157     bool planaire=true;
158     for(int k=0;k<m_nb_pas_echantillonnage;k++)
159     {
160     double t = tmin+k*pas_t;
161     double xyz[3];
162     mgarete->evaluer(t,xyz);
163     double uv1[2];
164     face1->inverser(uv1,xyz);
165     double uv2[2];
166     face2->inverser(uv2,xyz);
167     double normale1[3];
168     face1->calcul_normale_unitaire(uv1,normale1);
169     double normale2[3];
170     face2->calcul_normale_unitaire(uv2,normale2);
171     OT_VECTEUR_3D vect_normale1(normale1);
172     OT_VECTEUR_3D vect_normale2(normale2);
173     double angle = fabs(acos((vect_normale1*vect_normale2)/(vect_normale1.get_longueur()*vect_normale2.get_longueur())));
174     if(planaire)
175     if(angle<m_angle_min)
176     if(!OPERATEUR::egal(angle,0.0,m_precision))planaire=false;
177     if(!planaire)
178     {
179     if(angle<m_angle_min) return false;
180     }
181     }
182     }
183     }
184     return true;
185     }
186    
187     bool MSTRUCT_GENERATEUR::longueur_acceptable(MG_ARETE* mgarete)
188     {
189     if(OT_GEOMETRIE::get_longueur(mgarete)<m_longueur_min) return false;
190     return true;
191     }
192    
193     bool MSTRUCT_GENERATEUR::aire_acceptable(MG_FACE* mgface)
194     {
195     if(OT_GEOMETRIE::get_aire(mgface,m_precision)<m_aire_min) return false;
196     return true;
197     }
198    
199     bool MSTRUCT_GENERATEUR::volume_acceptable(MG_VOLUME* mgvolume)
200     {
201     if(OT_GEOMETRIE::get_volume(mgvolume,m_precision)<m_volume_min) return false;
202     return true;
203     }
204    
205     bool MSTRUCT_GENERATEUR::volume_intersection_ves_acceptable(MG_VOLUME* mgvolume)
206     {
207     if(m_matrice->get_type_forme()==MG_CG_FORME::TYPE_FORME::VOLUME)
208     {
209     MG_CG_FORME_VOLUME* forme_volume_matrice = (MG_CG_FORME_VOLUME*)m_matrice;
210     double volume_inclusion = OT_GEOMETRIE::get_volume_intersection(mgvolume,forme_volume_matrice->get_mg_volume());
211     if(volume_inclusion>m_volume_min) return true;
212     return false;
213     }
214     else if(m_matrice->get_type_forme()==MG_CG_FORME::TYPE_FORME::MULTI_VOLUME)
215     {
216     MG_CG_FORME_MULTI_VOLUME* forme_multi_volume_matrice = (MG_CG_FORME_MULTI_VOLUME*)m_matrice;
217     double volume_inclusion=0;
218     std::map<long,MG_VOLUME*>::iterator it_volume;
219     for(MG_VOLUME* volume_matrice=forme_multi_volume_matrice->get_premier_mg_volume(it_volume);volume_matrice!=NULL;volume_matrice=forme_multi_volume_matrice->get_suivant_mg_volume(it_volume))
220     {
221     volume_inclusion+=OT_GEOMETRIE::get_volume_intersection(mgvolume,volume_matrice);
222     }
223     if(volume_inclusion>m_volume_min) return true;
224     return false;
225     }
226     }
227    
228     bool MSTRUCT_GENERATEUR::distance_acceptable_bord_ves(MG_VOLUME* mgvolume)
229     {
230     BOITE_3D boite_volume=mgvolume->get_boite_3D(m_nb_pas_echantillonnage);
231     if(boite_volume.get_xmin()<(m_boite3d_ves.get_xmin()+m_distance_min)) return false;
232     if(boite_volume.get_ymin()<(m_boite3d_ves.get_ymin()+m_distance_min)) return false;
233     if(boite_volume.get_zmin()<(m_boite3d_ves.get_zmin()+m_distance_min)) return false;
234     if(boite_volume.get_xmax()>(m_boite3d_ves.get_xmax()-m_distance_min)) return false;
235     if(boite_volume.get_ymax()>(m_boite3d_ves.get_ymax()-m_distance_min)) return false;
236     if(boite_volume.get_zmax()>(m_boite3d_ves.get_zmax()-m_distance_min)) return false;
237     return true;
238     }
239    
240    
241     bool MSTRUCT_GENERATEUR::configuration_acceptable(MG_CG_FORME_VOLUME* forme_volume_inclusion,MG_CG_FORME* forme_matrice)
242     {
243     TPL_MAP_ENTITE<MG_FACE*> tpl_map_face;
244     TPL_MAP_ENTITE<MG_ARETE*> tpl_map_arete;
245     TPL_MAP_ENTITE<MG_SOMMET*> tpl_map_sommet;
246     TPL_MAP_ENTITE<MG_ELEMENT_TOPOLOGIQUE*> tpl_map_ele_topo_matrice;
247     TPL_GRILLE<MG_ELEMENT_TOPOLOGIQUE*> tpl_grille_ele_topo;
248     tpl_grille_ele_topo.initialiser(m_boite3d_ves.get_xmin()-0.1,m_boite3d_ves.get_ymin()-0.1,m_boite3d_ves.get_zmin()-0.1,
249     m_boite3d_ves.get_xmax()+0.1,m_boite3d_ves.get_ymax()+0.1,m_boite3d_ves.get_zmax()+0.1,
250     m_nb_pas_grille,m_nb_pas_grille,m_nb_pas_grille);
251     if(OT_GEOMETRIE::get_volume(forme_volume_inclusion->get_mg_volume())<m_volume_min) return FAIL;
252     OT_GEOMETRIE::get_map_mg_ele_topo_sous_jacent(forme_volume_inclusion->get_mg_volume(),tpl_map_sommet,tpl_map_arete,tpl_map_face);
253     if(forme_matrice!=NULL)
254     {
255     if(forme_matrice->get_type_forme()==MG_CG_FORME::TYPE_FORME::VOLUME)
256     {
257     MG_CG_FORME_VOLUME* forme_volume_matrice=(MG_CG_FORME_VOLUME*)forme_matrice;
258     OT_GEOMETRIE::get_map_mg_ele_topo_sous_jacent(forme_volume_matrice->get_mg_volume(),tpl_map_sommet,tpl_map_arete,tpl_map_face);
259     }
260     else if(forme_matrice->get_type_forme()==MG_CG_FORME::TYPE_FORME::MULTI_VOLUME)
261     {
262     MG_CG_FORME_MULTI_VOLUME* forme_multi_volume_matrice=(MG_CG_FORME_MULTI_VOLUME*)forme_matrice;
263     std::map<long,MG_VOLUME*>::iterator it_volume;
264     for(MG_VOLUME* volume=forme_multi_volume_matrice->get_premier_mg_volume(it_volume);volume!=NULL;volume=forme_multi_volume_matrice->get_suivant_mg_volume(it_volume))
265     {
266     OT_GEOMETRIE::get_map_mg_ele_topo_sous_jacent(volume,tpl_map_sommet,tpl_map_arete,tpl_map_face);
267     }
268     }
269     }
270     TPL_MAP_ENTITE<MG_FACE*>::ITERATEUR it_face;
271     for(MG_FACE* face=tpl_map_face.get_premier(it_face);face!=NULL;face=tpl_map_face.get_suivant(it_face))
272     {
273     if(OT_GEOMETRIE::get_aire(face)<m_aire_min) return FAIL;
274     tpl_map_ele_topo_matrice.ajouter(face);
275     tpl_grille_ele_topo.inserer(face);
276     }
277     tpl_map_face.vide();
278     TPL_MAP_ENTITE<MG_ARETE*>::ITERATEUR it_arete;
279     for(MG_ARETE* arete=tpl_map_arete.get_premier(it_arete);arete!=NULL;arete=tpl_map_arete.get_suivant(it_arete))
280     {
281     if(OT_GEOMETRIE::get_longueur(arete)<m_longueur_min) return FAIL;
282     if(!angle_entre_face_acceptable(arete)) return FAIL;
283     tpl_map_ele_topo_matrice.ajouter(arete);
284     tpl_grille_ele_topo.inserer(arete);
285     }
286     tpl_map_arete.vide();
287     TPL_MAP_ENTITE<MG_SOMMET*>::ITERATEUR it_sommet;
288     for(MG_SOMMET* sommet=tpl_map_sommet.get_premier(it_sommet);sommet!=NULL;sommet=tpl_map_sommet.get_suivant(it_sommet))
289     {
290     tpl_map_ele_topo_matrice.ajouter(sommet);
291     tpl_grille_ele_topo.inserer(sommet);
292     }
293     tpl_map_sommet.vide();
294    
295     OT_GEOMETRIE::get_map_mg_ele_topo_sous_jacent(forme_volume_inclusion->get_mg_volume(),tpl_map_sommet,tpl_map_arete,tpl_map_face);
296     std::stack<MG_ELEMENT_TOPOLOGIQUE*> stack_ele_topo;
297     for(MG_FACE* face=tpl_map_face.get_premier(it_face);face!=NULL;face=tpl_map_face.get_suivant(it_face))
298     {
299     if(OT_GEOMETRIE::get_aire(face)<m_aire_min) return FAIL;
300     tpl_grille_ele_topo.inserer(face);
301     stack_ele_topo.push(face);
302     }
303     for(MG_ARETE* arete=tpl_map_arete.get_premier(it_arete);arete!=NULL;arete=tpl_map_arete.get_suivant(it_arete))
304     {
305     if(OT_GEOMETRIE::get_longueur(arete)<m_longueur_min) return FAIL;
306     if(!angle_entre_face_acceptable(arete)) return FAIL;
307    
308     tpl_grille_ele_topo.inserer(arete);
309     stack_ele_topo.push(arete);
310    
311     }
312     for(MG_SOMMET* sommet=tpl_map_sommet.get_premier(it_sommet);sommet!=NULL;sommet=tpl_map_sommet.get_suivant(it_sommet))
313     {
314     tpl_grille_ele_topo.inserer(sommet);
315     stack_ele_topo.push(sommet);
316     }
317     while(!stack_ele_topo.empty())
318     {
319     MG_ELEMENT_TOPOLOGIQUE* ele_topo=stack_ele_topo.top(); stack_ele_topo.pop();
320     TPL_MAP_ENTITE<MG_ELEMENT_TOPOLOGIQUE*> tpl_map_topo_trouve;
321     BOITE_3D boite_recherche(ele_topo->get_boite_3D().get_xmin()-m_distance_min,
322     ele_topo->get_boite_3D().get_ymin()-m_distance_min,
323     ele_topo->get_boite_3D().get_zmin()-m_distance_min,
324     ele_topo->get_boite_3D().get_xmax()+m_distance_min,
325     ele_topo->get_boite_3D().get_ymax()+m_distance_min,
326     ele_topo->get_boite_3D().get_zmax()+m_distance_min);
327     tpl_grille_ele_topo.rechercher(boite_recherche,tpl_map_topo_trouve);
328     TPL_MAP_ENTITE<MG_ELEMENT_TOPOLOGIQUE*>::ITERATEUR it_trouve;
329     for(MG_ELEMENT_TOPOLOGIQUE* ele_top_trouve=tpl_map_topo_trouve.get_premier(it_trouve);ele_top_trouve!=NULL;ele_top_trouve=tpl_map_topo_trouve.get_suivant(it_trouve))
330     {
331     if(ele_topo==ele_top_trouve) continue;
332     if(OT_GEOMETRIE::get_lien_topologique(ele_topo,ele_top_trouve)==OT_GEOMETRIE::TYPE_LIEN_TOPOLOGIQUE::AUCUN)
333     {
334     double distance;
335     if(OT_GEOMETRIE::get_distance_min_mg_eletopo_mg_eletopo(ele_topo,ele_top_trouve,distance)==FAIL) return FAIL;
336     if(distance<m_distance_min) return FAIL;
337     }
338     }
339     tpl_grille_ele_topo.supprimer(ele_topo);
340     }
341     return OK;
342     }
343    
344    
345    
346     void MSTRUCT_GENERATEUR::echec_insertion_inclusion(std::vector< MG_CG_FORME* >& vector_forme, std::vector< MG_CG_OPERATEUR* >& vector_operateur)
347     {
348     std::vector<MG_CG_FORME*>::iterator it_forme;
349     for(it_forme=vector_forme.begin();it_forme!=vector_forme.end();it_forme++) m_mgcg_modele->supprimer_mgcg_formeid((*it_forme)->get_id());
350     std::vector<MG_CG_OPERATEUR*>::iterator it_operateur;
351     for(it_operateur=vector_operateur.begin();it_operateur!=vector_operateur.end();it_operateur++) m_mgcg_modele->supprimer_mgcg_operateurid((*it_operateur)->get_id());
352     }
353    
354     int MSTRUCT_GENERATEUR::faces_correspondantes(MG_VOLUME* volume1, MG_VOLUME* volume2,double decallage,int nb_pas,double precision)
355     {
356     TPL_MAP_ENTITE<MG_FACE*> tpl_map_face_volume1;
357     TPL_MAP_ENTITE<MG_FACE*> tpl_map_face_volume2;
358     OT_GEOMETRIE::get_map_mg_face_sous_jacent(volume1,tpl_map_face_volume1);
359     OT_GEOMETRIE::get_map_mg_face_sous_jacent(volume2,tpl_map_face_volume2);
360     TPL_MAP_ENTITE<MG_FACE*>::ITERATEUR it_face1;
361     TPL_MAP_ENTITE<MG_FACE*>::ITERATEUR it_face2;
362     while(tpl_map_face_volume1.get_nb()>0 && tpl_map_face_volume2.get_nb()>0)
363     {
364     bool trouve=false;
365     MG_FACE* face1=tpl_map_face_volume1.get_premier(it_face1);
366     MG_FACE* face2=tpl_map_face_volume2.get_premier(it_face2);
367     while(trouve==false && face2!=NULL)
368     {
369     double distance_moyenne;
370     double ecart_type_distance;
371     if(OT_GEOMETRIE::declage_mg_face_mg_face(face1,face2,distance_moyenne,ecart_type_distance,nb_pas,precision)!=FAIL)
372     {
373     if(OPERATEUR::egal(distance_moyenne,decallage,precision) && ecart_type_distance<precision)
374     {
375     trouve=true;
376     tpl_map_face_volume1.supprimer(face1);
377     tpl_map_face_volume2.supprimer(face2);
378     }
379     else if(OPERATEUR::egal(distance_moyenne,0.0,precision) && ecart_type_distance<precision)
380     {
381     trouve=true;
382     tpl_map_face_volume1.supprimer(face1);
383     tpl_map_face_volume2.supprimer(face2);
384     }
385     }
386     face2=tpl_map_face_volume2.get_suivant(it_face2);
387     }
388     if(trouve==false) return FAIL;
389     }
390     return OK;
391     }
392    
393    
394     bool MSTRUCT_GENERATEUR::insertion_inclusion(MG_CG_FORME* forme_inclusion,
395     double fraction_volumique_cible,
396     double eps_fra_vol,
397     double &fraction_volumique_actuelle,
398     std::vector<MG_CG_FORME*> &vector_inclusion_genere,
399     int &position_relative,
400     std::vector< double > *vector_epaisseur_couche,
401     bool porosite)
402     {
403    
404     MG_CG_FORME_VOLUME* boite_frac_vol=MG_CG::creer_MG_CG_FORME_VOLUME_BOITE(m_mgcg_modele,
405     m_boite3d_distribution.get_xmin(),m_boite3d_distribution.get_ymin(),m_boite3d_distribution.get_zmin(),
406     m_boite3d_distribution.get_xmax(),m_boite3d_distribution.get_ymax(),m_boite3d_distribution.get_zmax());
407     boite_frac_vol->construire();
408     double volume_particule;
409     MG_CG_FORME* matrice=m_matrice;
410     std::vector<MG_CG_FORME*> vector_forme_genere;
411     std::vector<MG_CG_OPERATEUR*> vector_operateur_genere;
412     vector_forme_genere.push_back(forme_inclusion);
413     vector_forme_genere.push_back(boite_frac_vol);
414     if(forme_inclusion->get_type_forme()!=MG_CG_FORME::TYPE_FORME::VOLUME)
415     {
416     echec_insertion_inclusion(vector_forme_genere,vector_operateur_genere);
417     return FAIL;
418     }
419     MG_CG_FORME_VOLUME* forme_volume_inclusion = (MG_CG_FORME_VOLUME*)forme_inclusion;
420     MG_CG_FORME_VOLUME* forme_volume_verif_grille=NULL;
421     bool multicouche=false;
422     if(vector_epaisseur_couche!=NULL) multicouche=true;
423     if(multicouche)
424     {
425     // Creation de l'enveloppe externe
426     double decallage_max=0;
427     std::vector<double>::iterator it_epaisseur;
428     for(it_epaisseur=vector_epaisseur_couche->begin();it_epaisseur!=vector_epaisseur_couche->end();it_epaisseur++) decallage_max+=*it_epaisseur;
429     MG_CG_OP_TRANSF_DECALLAGE* op_decallage_couche_externe=MG_CG::creer_MG_CG_OP_TRANSF_DECALLAGE(m_mgcg_modele,MAGIC::SEMANTIQUECSG::DECALLAGE,forme_volume_inclusion);
430     vector_operateur_genere.push_back(op_decallage_couche_externe);
431     op_decallage_couche_externe->change_decallage(decallage_max);
432     if(op_decallage_couche_externe->construire()==FAIL)
433     {
434     echec_insertion_inclusion(vector_forme_genere,vector_operateur_genere);
435     return FAIL;
436     }
437     vector_forme_genere.push_back(op_decallage_couche_externe->get_mgcg_forme_sortie());
438     if(op_decallage_couche_externe->get_mgcg_forme_sortie()->get_type_forme()!=MG_CG_FORME::TYPE_FORME::VOLUME)
439     {
440     echec_insertion_inclusion(vector_forme_genere,vector_operateur_genere);
441     return FAIL;
442     }
443     MG_CG_FORME_VOLUME* forme_volume_couche_externe=(MG_CG_FORME_VOLUME*)op_decallage_couche_externe->get_mgcg_forme_sortie();
444     forme_volume_verif_grille=forme_volume_couche_externe;
445     MG_VOLUME* volume_inclusion=forme_volume_inclusion->get_mg_volume();
446     MG_VOLUME* volume_inclusion_couche_externe=forme_volume_couche_externe->get_mg_volume();
447     // Verification du positionnement
448     if(positionnement_acceptable(forme_volume_couche_externe->get_mg_volume())==FAIL)
449     {
450     echec_insertion_inclusion(vector_forme_genere,vector_operateur_genere);
451     return FAIL;
452     }
453     // Verification de la fraction volumique
454     volume_particule=OT_GEOMETRIE::get_volume_intersection(volume_inclusion_couche_externe,boite_frac_vol->get_mg_volume(),m_precision);
455     m_mgcg_modele->supprimer_mgcg_formeid(boite_frac_vol->get_id());
456     // volume_particule=OT_GEOMETRIE::get_volume(volume_inclusion_couche_externe,m_precision);
457     if((volume_particule+fraction_volumique_actuelle)>fraction_volumique_cible+eps_fra_vol)
458     {
459     echec_insertion_inclusion(vector_forme_genere,vector_operateur_genere);
460     return FAIL;
461     }
462     // Verification des plans de coupe (coeur VS enveloppe)
463     bool plan_x0yz_inc;bool plan_x1yz_inc;bool plan_xy0z_inc;bool plan_xy1z_inc;bool plan_xyz0_inc;bool plan_xyz1_inc;
464     bool plan_x0yz_ext;bool plan_x1yz_ext;bool plan_xy0z_ext;bool plan_xy1z_ext;bool plan_xyz0_ext;bool plan_xyz1_ext;
465     int position_relative_inc=position_relative_inclusion_ves(volume_inclusion,plan_x0yz_inc,plan_x1yz_inc,plan_xy0z_inc,plan_xy1z_inc,plan_xyz0_inc,plan_xyz1_inc);
466     int position_relative_ext=position_relative_inclusion_ves(volume_inclusion_couche_externe,plan_x0yz_ext,plan_x1yz_ext,plan_xy0z_ext,plan_xy1z_ext,plan_xyz0_ext,plan_xyz1_ext);
467     if(position_relative_inc!=position_relative_ext || position_relative_ext==POSITION_RELATIVE::EXTERIEUR)
468     {
469     echec_insertion_inclusion(vector_forme_genere,vector_operateur_genere);
470     return FAIL;
471     }
472     if(!m_intersection_bords_ves && position_relative_ext==POSITION_RELATIVE::AU_BORD)
473     {
474     echec_insertion_inclusion(vector_forme_genere,vector_operateur_genere);
475     return FAIL;
476     }
477     if(position_relative_ext==POSITION_RELATIVE::AU_BORD)
478     {
479     bool ok_plan=false;
480     if(plan_x0yz_inc==plan_x0yz_ext)
481     if(plan_x1yz_inc==plan_x1yz_ext)
482     if(plan_xy0z_inc==plan_xy0z_ext)
483     if(plan_xy1z_inc==plan_xy1z_ext)
484     if(plan_xyz0_inc==plan_xyz0_ext)
485     if(plan_xyz1_inc==plan_xyz1_ext) ok_plan=true;
486     if(!ok_plan)
487     {
488     echec_insertion_inclusion(vector_forme_genere,vector_operateur_genere);
489     return FAIL;
490     }
491     }
492     // Creation de l'intersection couche externe / matrice
493     MG_CG_OP_BOOL_INTERSECTION* op_inter_ext=MG_CG::creer_MG_CG_OP_BOOL_INTERSECTION(m_mgcg_modele,MAGIC::SEMANTIQUECSG::INTERSECTION,forme_volume_couche_externe,matrice);
494     vector_operateur_genere.push_back(op_inter_ext);
495     if(op_inter_ext->construire()==FAIL)
496     {
497     echec_insertion_inclusion(vector_forme_genere,vector_operateur_genere);
498     return FAIL;
499     }
500     vector_forme_genere.push_back(op_inter_ext->get_mgcg_forme_sortie());
501     if(op_inter_ext->get_mgcg_forme_sortie()->get_type_forme()!=MG_CG_FORME::TYPE_FORME::VOLUME)
502     {
503     echec_insertion_inclusion(vector_forme_genere,vector_operateur_genere);
504     return FAIL;
505     }
506     MG_CG_FORME_VOLUME* couche_externe_inter = (MG_CG_FORME_VOLUME*)op_inter_ext->get_mgcg_forme_sortie();
507     // Fragmentation de l'inclusion / matrice
508     MG_CG_OP_INCLUSION_MATRICE* op_inc_matrice_ext=MG_CG::creer_MG_CG_OP_INCLUSION_MATRICE(m_mgcg_modele,MAGIC::SEMANTIQUECSG::FRAGMENT,couche_externe_inter,matrice);
509     vector_operateur_genere.push_back(op_inc_matrice_ext);
510     if(op_inc_matrice_ext->construire()==FAIL)
511     {
512     echec_insertion_inclusion(vector_forme_genere,vector_operateur_genere);
513     return FAIL;
514     }
515     vector_forme_genere.push_back(op_inc_matrice_ext->get_mgcg_forme_matrice_sortie());
516     vector_forme_genere.push_back(op_inc_matrice_ext->get_mgcg_forme_inclusion_sortie());
517     if(op_inc_matrice_ext->get_mgcg_forme_inclusion_sortie()->get_type_forme()!=MG_CG_FORME::TYPE_FORME::VOLUME)
518     {
519     echec_insertion_inclusion(vector_forme_genere,vector_operateur_genere);
520     return FAIL;
521     }
522     matrice=op_inc_matrice_ext->get_mgcg_forme_matrice_sortie();
523     MG_CG_FORME_VOLUME* couche_externe_matrice=(MG_CG_FORME_VOLUME*)op_inc_matrice_ext->get_mgcg_forme_inclusion_sortie();
524     if(configuration_acceptable(couche_externe_matrice,matrice)==FAIL)
525     {
526     echec_insertion_inclusion(vector_forme_genere,vector_operateur_genere);
527     return FAIL;
528     }
529     // OK POUR GENERER
530     MG_CG_FORME_VOLUME* matrice_couche_partielle;
531     if(porosite)
532     {
533     MG_CG_OP_BOOL_DIFFERENCE* op_diff_porosite=MG_CG::creer_MG_CG_OP_BOOL_DIFFERENCE(m_mgcg_modele,MAGIC::SEMANTIQUECSG::SOUSTRACTION,couche_externe_matrice,forme_volume_inclusion);
534     vector_operateur_genere.push_back(op_diff_porosite);
535     if(op_diff_porosite->construire()==FAIL)
536     {
537     echec_insertion_inclusion(vector_forme_genere,vector_operateur_genere);
538     return FAIL;
539     }
540     vector_forme_genere.push_back(op_diff_porosite->get_mgcg_forme_sortie());
541     if(op_diff_porosite->get_mgcg_forme_sortie()->get_type_forme()!=MG_CG_FORME::TYPE_FORME::VOLUME)
542     {
543     echec_insertion_inclusion(vector_forme_genere,vector_operateur_genere);
544     return FAIL;
545     }
546     matrice_couche_partielle=(MG_CG_FORME_VOLUME*)op_diff_porosite->get_mgcg_forme_sortie();
547     }
548     else
549     {
550     MG_CG_OP_BOOL_INTERSECTION* op_inter_couche=MG_CG::creer_MG_CG_OP_BOOL_INTERSECTION(m_mgcg_modele,MAGIC::SEMANTIQUECSG::INTERSECTION,forme_volume_inclusion,couche_externe_matrice);
551     vector_operateur_genere.push_back(op_inter_couche);
552     if(op_inter_couche->construire()==FAIL)
553     {
554     echec_insertion_inclusion(vector_forme_genere,vector_operateur_genere);
555     return FAIL;
556     }
557     vector_forme_genere.push_back(op_inter_couche->get_mgcg_forme_sortie());
558     if(op_inter_couche->get_mgcg_forme_sortie()->get_type_forme()!=MG_CG_FORME::TYPE_FORME::VOLUME)
559     {
560     echec_insertion_inclusion(vector_forme_genere,vector_operateur_genere);
561     return FAIL;
562     }
563     MG_CG_OP_INCLUSION_MATRICE* op_inc_matrice_coeur=MG_CG::creer_MG_CG_OP_INCLUSION_MATRICE(m_mgcg_modele,MAGIC::SEMANTIQUECSG::FRAGMENT,op_inter_couche->get_mgcg_forme_sortie(),couche_externe_matrice);
564     vector_operateur_genere.push_back(op_inc_matrice_coeur);
565     if(op_inc_matrice_coeur->construire()==FAIL)
566     {
567     echec_insertion_inclusion(vector_forme_genere,vector_operateur_genere);
568     return FAIL;
569     }
570     vector_forme_genere.push_back(op_inc_matrice_coeur->get_mgcg_forme_inclusion_sortie());
571     vector_forme_genere.push_back(op_inc_matrice_coeur->get_mgcg_forme_matrice_sortie());
572     if(op_inc_matrice_coeur->get_mgcg_forme_inclusion_sortie()->get_type_forme()!=MG_CG_FORME::TYPE_FORME::VOLUME ||
573     op_inc_matrice_coeur->get_mgcg_forme_matrice_sortie()->get_type_forme()!=MG_CG_FORME::TYPE_FORME::VOLUME)
574     {
575     echec_insertion_inclusion(vector_forme_genere,vector_operateur_genere);
576     return FAIL;
577     }
578    
579     MG_CG_FORME_VOLUME* inclusion_coeur=(MG_CG_FORME_VOLUME*)op_inc_matrice_coeur->get_mgcg_forme_inclusion_sortie();
580     matrice_couche_partielle=(MG_CG_FORME_VOLUME*)op_inc_matrice_coeur->get_mgcg_forme_matrice_sortie();
581     if(configuration_acceptable(inclusion_coeur,NULL)==FAIL)
582     {
583     echec_insertion_inclusion(vector_forme_genere,vector_operateur_genere);
584     return FAIL;
585     }
586     // Verification des faces correspondantes
587     if(position_relative_ext==POSITION_RELATIVE::AU_BORD)
588     if(faces_correspondantes(inclusion_coeur->get_mg_volume(),couche_externe_matrice->get_mg_volume(),decallage_max,m_nb_pas_echantillonnage,m_precision)==FAIL)
589     {
590     echec_insertion_inclusion(vector_forme_genere,vector_operateur_genere);
591     return FAIL;
592     }
593     if(position_relative_ext==POSITION_RELATIVE::INTERIEUR) MG_CG::creer_MG_CG_INFO_STRING(op_inc_matrice_coeur->get_mgcg_forme_inclusion_sortie(),(char*)"POSITION_RELATIVE",(char*)"INTERIEUR");
594     if(position_relative_ext==POSITION_RELATIVE::AU_BORD) MG_CG::creer_MG_CG_INFO_STRING(op_inc_matrice_coeur->get_mgcg_forme_inclusion_sortie(),(char*)"POSITION_RELATIVE",(char*)"AU_BORD");
595     vector_inclusion_genere.push_back(inclusion_coeur);
596     }
597     double decallage=0;
598     double epaisseur;
599     for(it_epaisseur=vector_epaisseur_couche->begin();it_epaisseur!=vector_epaisseur_couche->end()-1;it_epaisseur++)
600     {
601     epaisseur=*it_epaisseur;
602     decallage+=*it_epaisseur;
603     MG_CG_OP_TRANSF_DECALLAGE* op_decallage_couche=MG_CG::creer_MG_CG_OP_TRANSF_DECALLAGE(m_mgcg_modele,MAGIC::SEMANTIQUECSG::DECALLAGE,forme_volume_inclusion);
604     vector_operateur_genere.push_back(op_decallage_couche);
605     op_decallage_couche->change_decallage(decallage);
606     if(op_decallage_couche->construire()==FAIL)
607     {
608     echec_insertion_inclusion(vector_forme_genere,vector_operateur_genere);
609     return FAIL;
610     }
611     vector_forme_genere.push_back(op_decallage_couche->get_mgcg_forme_sortie());
612     if(op_decallage_couche->get_mgcg_forme_sortie()->get_type_forme()!=MG_CG_FORME::TYPE_FORME::VOLUME)
613     {
614     echec_insertion_inclusion(vector_forme_genere,vector_operateur_genere);
615     return FAIL;
616     }
617     MG_CG_OP_BOOL_INTERSECTION* op_inter_couche=MG_CG::creer_MG_CG_OP_BOOL_INTERSECTION(m_mgcg_modele,MAGIC::SEMANTIQUECSG::INTERSECTION,op_decallage_couche->get_mgcg_forme_sortie(),matrice_couche_partielle);
618     vector_operateur_genere.push_back(op_inter_couche);
619     if(op_inter_couche->construire()==FAIL)
620     {
621     echec_insertion_inclusion(vector_forme_genere,vector_operateur_genere);
622     return FAIL;
623     }
624     vector_forme_genere.push_back(op_inter_couche->get_mgcg_forme_sortie());
625     if(op_inter_couche->get_mgcg_forme_sortie()->get_type_forme()!=MG_CG_FORME::TYPE_FORME::VOLUME)
626     {
627     echec_insertion_inclusion(vector_forme_genere,vector_operateur_genere);
628     return FAIL;
629     }
630     MG_CG_OP_INCLUSION_MATRICE* op_inc_matrice_couche=MG_CG::creer_MG_CG_OP_INCLUSION_MATRICE(m_mgcg_modele,MAGIC::SEMANTIQUECSG::FRAGMENT,op_inter_couche->get_mgcg_forme_sortie(),matrice_couche_partielle);
631     vector_operateur_genere.push_back(op_inc_matrice_couche);
632     if(op_inc_matrice_couche->construire()==FAIL)
633     {
634     echec_insertion_inclusion(vector_forme_genere,vector_operateur_genere);
635     return FAIL;
636     }
637     vector_forme_genere.push_back(op_inc_matrice_couche->get_mgcg_forme_inclusion_sortie());
638     vector_forme_genere.push_back(op_inc_matrice_couche->get_mgcg_forme_matrice_sortie());
639     if(op_inc_matrice_couche->get_mgcg_forme_inclusion_sortie()->get_type_forme()!=MG_CG_FORME::TYPE_FORME::VOLUME ||
640     op_inc_matrice_couche->get_mgcg_forme_matrice_sortie()->get_type_forme()!=MG_CG_FORME::TYPE_FORME::VOLUME)
641     {
642     echec_insertion_inclusion(vector_forme_genere,vector_operateur_genere);
643     return FAIL;
644     }
645     MG_CG_FORME_VOLUME* inclusion_couche = (MG_CG_FORME_VOLUME*)op_inc_matrice_couche->get_mgcg_forme_inclusion_sortie();
646     vector_inclusion_genere.push_back(inclusion_couche);
647     if(position_relative_ext==POSITION_RELATIVE::INTERIEUR) MG_CG::creer_MG_CG_INFO_STRING(inclusion_couche,(char*)"POSITION_RELATIVE",(char*)"INTERIEUR");
648     if(position_relative_ext==POSITION_RELATIVE::AU_BORD) MG_CG::creer_MG_CG_INFO_STRING(inclusion_couche,(char*)"POSITION_RELATIVE",(char*)"AU_BORD");
649     MG_CG::creer_MG_CG_INFO_DOUBLE(inclusion_couche,(char*)"EPAISSEUR",epaisseur);
650     matrice_couche_partielle=(MG_CG_FORME_VOLUME*)op_inc_matrice_couche->get_mgcg_forme_matrice_sortie();
651     configurer_volume_mince(inclusion_couche->get_mg_volume(),epaisseur);
652     }
653     epaisseur=*it_epaisseur;
654     if(position_relative_ext==POSITION_RELATIVE::INTERIEUR) MG_CG::creer_MG_CG_INFO_STRING(matrice_couche_partielle,(char*)"POSITION_RELATIVE",(char*)"INTERIEUR");
655     if(position_relative_ext==POSITION_RELATIVE::AU_BORD) MG_CG::creer_MG_CG_INFO_STRING(matrice_couche_partielle,(char*)"POSITION_RELATIVE",(char*)"AU_BORD");
656     MG_CG::creer_MG_CG_INFO_DOUBLE(matrice_couche_partielle,(char*)"EPAISSEUR",epaisseur);
657     vector_inclusion_genere.push_back(matrice_couche_partielle);
658     configurer_volume_mince(matrice_couche_partielle->get_mg_volume(),epaisseur);
659     position_relative=position_relative_ext;
660     }
661     else
662     {
663     // Verification du positionnement
664     if(positionnement_acceptable(forme_volume_inclusion->get_mg_volume())==FAIL)
665     {
666     echec_insertion_inclusion(vector_forme_genere,vector_operateur_genere);
667     return FAIL;
668     }
669     // Verification de la fraction volumique
670     volume_particule=OT_GEOMETRIE::get_volume_intersection(forme_volume_inclusion->get_mg_volume(),boite_frac_vol->get_mg_volume(),m_precision);
671     m_mgcg_modele->supprimer_mgcg_formeid(boite_frac_vol->get_id());
672     if((volume_particule+fraction_volumique_actuelle)>fraction_volumique_cible+eps_fra_vol)
673     {
674     echec_insertion_inclusion(vector_forme_genere,vector_operateur_genere);
675     return FAIL;
676     }
677     // Verification de la position relative
678     bool plan_x0yz;bool plan_x1yz;bool plan_xy0z;bool plan_xy1z;bool plan_xyz0;bool plan_xyz1;
679     position_relative=position_relative_inclusion_ves(forme_volume_inclusion->get_mg_volume(),plan_x0yz,plan_x1yz,plan_xy0z,plan_xy1z,plan_xyz0,plan_xyz1);
680     if(position_relative==POSITION_RELATIVE::EXTERIEUR)
681     {
682     echec_insertion_inclusion(vector_forme_genere,vector_operateur_genere);
683     return FAIL;
684     }
685     if(!m_intersection_bords_ves && position_relative==POSITION_RELATIVE::AU_BORD)
686     {
687     echec_insertion_inclusion(vector_forme_genere,vector_operateur_genere);
688     return FAIL;
689     }
690     if(position_relative==POSITION_RELATIVE::AU_BORD)
691     {
692     MG_CG_OP_BOOL_INTERSECTION* op_inter=MG_CG::creer_MG_CG_OP_BOOL_INTERSECTION(m_mgcg_modele,MAGIC::SEMANTIQUECSG::INTERSECTION,forme_volume_inclusion,m_matrice);
693     vector_operateur_genere.push_back(op_inter);
694     if(op_inter->construire()==FAIL)
695     {
696     echec_insertion_inclusion(vector_forme_genere,vector_operateur_genere);
697     return FAIL;
698     }
699     vector_forme_genere.push_back(op_inter->get_mgcg_forme_sortie());
700     if(op_inter->get_mgcg_forme_sortie()->get_type_forme()!=MG_CG_FORME::TYPE_FORME::VOLUME)
701     {
702     echec_insertion_inclusion(vector_forme_genere,vector_operateur_genere);
703     return FAIL;
704     }
705     forme_volume_inclusion=(MG_CG_FORME_VOLUME*)op_inter->get_mgcg_forme_sortie();
706     }
707     MG_CG_OP_INCLUSION_MATRICE* op_inc_matrice=MG_CG::creer_MG_CG_OP_INCLUSION_MATRICE(m_mgcg_modele,MAGIC::SEMANTIQUECSG::FRAGMENT,forme_volume_inclusion,m_matrice);
708     vector_operateur_genere.push_back(op_inc_matrice);
709     if(op_inc_matrice->construire()==FAIL)
710     {
711     echec_insertion_inclusion(vector_forme_genere,vector_operateur_genere);
712     return FAIL;
713     }
714     matrice=op_inc_matrice->get_mgcg_forme_matrice_sortie();
715     vector_forme_genere.push_back(matrice);
716     vector_forme_genere.push_back(op_inc_matrice->get_mgcg_forme_inclusion_sortie());
717     if(op_inc_matrice->get_mgcg_forme_inclusion_sortie()->get_type_forme()!=MG_CG_FORME::TYPE_FORME::VOLUME)
718     {
719     echec_insertion_inclusion(vector_forme_genere,vector_operateur_genere);
720     return FAIL;
721     }
722     forme_volume_inclusion=(MG_CG_FORME_VOLUME*)op_inc_matrice->get_mgcg_forme_inclusion_sortie();
723     if(configuration_acceptable(forme_volume_inclusion,matrice)==FAIL)
724     {
725     echec_insertion_inclusion(vector_forme_genere,vector_operateur_genere);
726     return FAIL;
727     }
728     if(position_relative==POSITION_RELATIVE::INTERIEUR) MG_CG::creer_MG_CG_INFO_STRING(forme_volume_inclusion,(char*)"POSITION_RELATIVE",(char*)"INTERIEUR");
729     if(position_relative==POSITION_RELATIVE::AU_BORD) MG_CG::creer_MG_CG_INFO_STRING(forme_volume_inclusion,(char*)"POSITION_RELATIVE",(char*)"AU_BORD");
730     vector_inclusion_genere.push_back(forme_volume_inclusion);
731     forme_volume_verif_grille=forme_volume_inclusion;
732     }
733     fraction_volumique_actuelle+=volume_particule;
734     m_grille_volume_inclusion.inserer(forme_volume_verif_grille->get_mg_volume());
735     m_matrice=matrice;
736     return OK;
737     }
738    
739     bool MSTRUCT_GENERATEUR::insertion_inclusion(MG_CG_FORME* forme_inclusion,
740     double& fraction_volumique_actuelle,
741     std::vector< MG_CG_FORME* >& vector_inclusion_genere,
742     int& position_relative)
743     {
744     double volume_particule;
745     MG_CG_FORME* matrice=m_matrice;
746     std::vector<MG_CG_FORME*> vector_forme_genere;
747     std::vector<MG_CG_OPERATEUR*> vector_operateur_genere;
748     vector_forme_genere.push_back(forme_inclusion);
749     if(forme_inclusion->get_type_forme()!=MG_CG_FORME::TYPE_FORME::VOLUME)
750     {
751     echec_insertion_inclusion(vector_forme_genere,vector_operateur_genere);
752     return FAIL;
753     }
754     MG_CG_FORME_VOLUME* forme_volume_inclusion = (MG_CG_FORME_VOLUME*)forme_inclusion;
755     MG_CG_FORME_VOLUME* forme_volume_verif_grille=NULL;
756     // Verification du positionnement
757     if(positionnement_acceptable(forme_volume_inclusion->get_mg_volume())==FAIL)
758     {
759     echec_insertion_inclusion(vector_forme_genere,vector_operateur_genere);
760     return FAIL;
761     }
762     // Verification de la position relative
763     bool plan_x0yz;bool plan_x1yz;bool plan_xy0z;bool plan_xy1z;bool plan_xyz0;bool plan_xyz1;
764     position_relative=position_relative_inclusion_ves(forme_volume_inclusion->get_mg_volume(),plan_x0yz,plan_x1yz,plan_xy0z,plan_xy1z,plan_xyz0,plan_xyz1);
765     if(position_relative==POSITION_RELATIVE::EXTERIEUR)
766     {
767     echec_insertion_inclusion(vector_forme_genere,vector_operateur_genere);
768     return FAIL;
769     }
770     if(!m_intersection_bords_ves && position_relative==POSITION_RELATIVE::AU_BORD)
771     {
772     echec_insertion_inclusion(vector_forme_genere,vector_operateur_genere);
773     return FAIL;
774     }
775     if(position_relative==POSITION_RELATIVE::AU_BORD)
776     {
777     MG_CG_OP_BOOL_INTERSECTION* op_inter=MG_CG::creer_MG_CG_OP_BOOL_INTERSECTION(m_mgcg_modele,MAGIC::SEMANTIQUECSG::INTERSECTION,forme_volume_inclusion,m_matrice);
778     vector_operateur_genere.push_back(op_inter);
779     if(op_inter->construire()==FAIL)
780     {
781     echec_insertion_inclusion(vector_forme_genere,vector_operateur_genere);
782     return FAIL;
783     }
784     vector_forme_genere.push_back(op_inter->get_mgcg_forme_sortie());
785     if(op_inter->get_mgcg_forme_sortie()->get_type_forme()!=MG_CG_FORME::TYPE_FORME::VOLUME)
786     {
787     echec_insertion_inclusion(vector_forme_genere,vector_operateur_genere);
788     return FAIL;
789     }
790     forme_volume_inclusion=(MG_CG_FORME_VOLUME*)op_inter->get_mgcg_forme_sortie();
791     }
792     MG_CG_OP_INCLUSION_MATRICE* op_inc_matrice=MG_CG::creer_MG_CG_OP_INCLUSION_MATRICE(m_mgcg_modele,MAGIC::SEMANTIQUECSG::FRAGMENT,forme_volume_inclusion,m_matrice);
793     vector_operateur_genere.push_back(op_inc_matrice);
794     if(op_inc_matrice->construire()==FAIL)
795     {
796     echec_insertion_inclusion(vector_forme_genere,vector_operateur_genere);
797     return FAIL;
798     }
799     matrice=op_inc_matrice->get_mgcg_forme_matrice_sortie();
800     vector_forme_genere.push_back(matrice);
801     vector_forme_genere.push_back(op_inc_matrice->get_mgcg_forme_inclusion_sortie());
802     if(op_inc_matrice->get_mgcg_forme_inclusion_sortie()->get_type_forme()!=MG_CG_FORME::TYPE_FORME::VOLUME)
803     {
804     echec_insertion_inclusion(vector_forme_genere,vector_operateur_genere);
805     return FAIL;
806     }
807     forme_volume_inclusion=(MG_CG_FORME_VOLUME*)op_inc_matrice->get_mgcg_forme_inclusion_sortie();
808     if(position_relative==POSITION_RELATIVE::INTERIEUR)
809     {
810     if(distance_acceptable_bord_ves(forme_volume_inclusion->get_mg_volume())==FAIL)
811     {
812     echec_insertion_inclusion(vector_forme_genere,vector_operateur_genere);
813     return FAIL;
814     }
815     }
816     if(configuration_acceptable(forme_volume_inclusion,matrice)==FAIL)
817     {
818     echec_insertion_inclusion(vector_forme_genere,vector_operateur_genere);
819     return FAIL;
820     }
821     if(position_relative==POSITION_RELATIVE::INTERIEUR) MG_CG::creer_MG_CG_INFO_STRING(forme_volume_inclusion,(char*)"POSITION_RELATIVE",(char*)"INTERIEUR");
822     if(position_relative==POSITION_RELATIVE::AU_BORD) MG_CG::creer_MG_CG_INFO_STRING(forme_volume_inclusion,(char*)"POSITION_RELATIVE",(char*)"AU_BORD");
823     vector_inclusion_genere.push_back(forme_volume_inclusion);
824     forme_volume_verif_grille=forme_volume_inclusion;
825     volume_particule=OT_GEOMETRIE::get_volume(forme_volume_inclusion->get_mg_volume());
826     fraction_volumique_actuelle+=volume_particule;
827     m_grille_volume_inclusion.inserer(forme_volume_verif_grille->get_mg_volume());
828     m_matrice=matrice;
829     return OK;
830     }
831    
832    
833     bool MSTRUCT_GENERATEUR::positionnement_acceptable(MG_VOLUME* inclusion)
834     {
835     TPL_MAP_ENTITE<MG_VOLUME*> map_volume_recherche;
836     BOITE_3D boite_recherche = inclusion->get_boite_3D();
837     boite_recherche.change_grosseur(1.1);
838     m_grille_volume_inclusion.rechercher(boite_recherche,map_volume_recherche);
839     TPL_MAP_ENTITE<MG_VOLUME*>::ITERATEUR it_volume;
840     MG_VOLUME* volume_recherche=map_volume_recherche.get_premier(it_volume);
841     double distance;
842     while(volume_recherche!=NULL)
843     {
844     if(OT_GEOMETRIE::get_distance_min_mg_volume_mg_volume(inclusion,volume_recherche,distance)==FAIL) return FAIL;
845     else
846     {
847     if(distance<=m_distance_min) return FAIL;
848     else volume_recherche=map_volume_recherche.get_suivant(it_volume);
849     }
850     }
851     return OK;
852     }
853    
854     bool MSTRUCT_GENERATEUR::lien_topologique_entre(MG_ELEMENT_TOPOLOGIQUE* mg_ele_topo1, MG_ELEMENT_TOPOLOGIQUE* mg_ele_topo2)
855     {
856     TPL_MAP_ENTITE<MG_SOMMET*> map_sommet_volume1;
857     OT_GEOMETRIE::get_map_mg_sommet_sous_jacent(mg_ele_topo1,map_sommet_volume1);
858     TPL_MAP_ENTITE<MG_SOMMET*> map_sommet_volume2;
859     OT_GEOMETRIE::get_map_mg_sommet_sous_jacent(mg_ele_topo2,map_sommet_volume2);
860     TPL_MAP_ENTITE<MG_SOMMET*>::ITERATEUR it_sommet;
861     for(MG_SOMMET* sommet=map_sommet_volume1.get_premier(it_sommet);sommet!=NULL;sommet=map_sommet_volume1.get_suivant(it_sommet))
862     {
863     if(map_sommet_volume2.existe(sommet)) return true;
864     }
865     return false;
866     }
867    
868     void MSTRUCT_GENERATEUR::cree_vector_groupe_forme(std::string nom, std::vector< MG_CG_GROUPE_FORME* >& vector_groupe_forme, std::vector< double >* vector_epaisseur_couches)
869     {
870     if(vector_epaisseur_couches==NULL)
871     {
872     MG_CG_GROUPE_FORME* groupe_forme;
873     groupe_forme=MG_CG::creer_MG_CG_GROUPE_FORME(m_mgcg_modele,nom);
874     m_vector_mgcg_groupe_forme_inclusion.push_back(groupe_forme);
875     vector_groupe_forme.push_back(groupe_forme);
876     }
877     else
878     {
879     MG_CG_GROUPE_FORME* groupe_forme;
880     groupe_forme=MG_CG::creer_MG_CG_GROUPE_FORME(m_mgcg_modele,nom);
881     m_vector_mgcg_groupe_forme_inclusion.push_back(groupe_forme);
882     vector_groupe_forme.push_back(groupe_forme);
883     for(int i=0;i<vector_epaisseur_couches->size();i++)
884     {
885     char nom_couche[500];
886     sprintf(nom_couche,"%s_%i",nom.c_str(),i);
887     groupe_forme=MG_CG::creer_MG_CG_GROUPE_FORME(m_mgcg_modele,nom_couche);
888     m_vector_mgcg_groupe_forme_inclusion.push_back(groupe_forme);
889     vector_groupe_forme.push_back(groupe_forme);
890     }
891     }
892     }
893    
894     int MSTRUCT_GENERATEUR::configurer_volume_mince(MG_CG_FORME* forme,double epaisseur)
895     {
896     if(forme->get_type_forme()==MG_CG_FORME::TYPE_FORME::VOLUME)
897     {
898     MG_VOLUME* volume = ((MG_CG_FORME_VOLUME*)forme)->get_mg_volume();
899     configurer_volume_mince(volume,epaisseur);
900     }
901     if(forme->get_type_forme()==MG_CG_FORME::TYPE_FORME::MULTI_VOLUME)
902     {
903     MG_CG_FORME_MULTI_VOLUME* forme_multi_volume = (MG_CG_FORME_MULTI_VOLUME*)forme;
904     std::map<long,MG_VOLUME*>::iterator it_volume;
905     for(MG_VOLUME* volume=forme_multi_volume->get_premier_mg_volume(it_volume);volume!=NULL;volume=forme_multi_volume->get_suivant_mg_volume(it_volume))
906     {
907     configurer_volume_mince(volume,epaisseur);
908     }
909     }
910     }
911    
912     int MSTRUCT_GENERATEUR::configurer_volume_mince(MG_VOLUME* volume, double epaisseur)
913     {
914     TPL_MAP_ENTITE<MG_FACE*> tpl_map_face;
915     OT_GEOMETRIE::get_map_mg_face_sous_jacent(volume,tpl_map_face);
916     if(tpl_map_face.get_nb()<2) return FAIL;
917     TPL_MAP_ENTITE<MG_FACE*>::ITERATEUR it_face1;
918     TPL_MAP_ENTITE<MG_FACE*>::ITERATEUR it_face2;
919     std::cout << "NB_FACE = " << tpl_map_face.get_nb() << std::endl;
920     while(tpl_map_face.get_nb()>=2)
921     {
922     bool trouve=false;
923     MG_FACE* face1=tpl_map_face.get_premier(it_face1);
924     it_face2=it_face1;
925     MG_FACE* face2=tpl_map_face.get_suivant(it_face2);
926     while(trouve==false && face1!=NULL && face2!=NULL)
927     {
928     double distance_moyenne;
929     double ecart_type_distance;
930     OT_GEOMETRIE::declage_mg_face_mg_face(face1,face2,distance_moyenne,ecart_type_distance,32);
931     std::cout << " Distance = " << distance_moyenne << " Std-dev = " << ecart_type_distance << " Epaisseur = " << epaisseur << std::endl;
932     if(OPERATEUR::egal(epaisseur,distance_moyenne,m_precision) && ecart_type_distance<m_precision)
933     {
934     std::cout << " Correspondance : " << face1->get_id() << " , " << face2->get_id() << std::endl;
935     trouve=true;
936     bool face1_base=false;
937     bool face2_base=false;
938     if(face1_base==true)
939     {
940     // std::cout << "Face1 base" << std::endl;
941     volume->ajouter_face_correspondante(face1,face2);
942     }
943     else if(face2_base==true)
944     {
945     // std::cout << "Face2 base" << std::endl;
946     volume->ajouter_face_correspondante(face2,face1);
947     }
948     else
949     {
950     double aire_face1 = OT_GEOMETRIE::get_aire(face1,m_precision);
951     double aire_face2 = OT_GEOMETRIE::get_aire(face2,m_precision);
952     if(aire_face1<aire_face2)
953     {
954     // std::cout << "Aire Face1 base" << std::endl;
955     volume->ajouter_face_correspondante(face1,face2);
956     }
957     else
958     {
959     // std::cout << "Aire Face2 base" << std::endl;
960     volume->ajouter_face_correspondante(face2,face1);
961     }
962     }
963     }
964     else face2=tpl_map_face.get_suivant(it_face2);
965     }
966     if(trouve)
967     {
968     tpl_map_face.supprimer(face1);
969     tpl_map_face.supprimer(face2);
970     }
971     else tpl_map_face.supprimer(face1);
972     }
973     }
974    
975    
976     void MSTRUCT_GENERATEUR::active_affichage(fonction_affiche* fonc)
977     {
978 couturad 929 fonc_affiche = fonc;
979     affichageactif = 1;
980     }
981    
982 couturad 951 void MSTRUCT_GENERATEUR::affiche(char* message)
983 couturad 929 {
984     if(affichageactif==1) fonc_affiche(message);
985     }