ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/REPOS_ERICCA/magic/lib/mailleur/src/mailleur3d.h
Revision: 78
Committed: Thu Apr 10 00:05:47 2008 UTC (17 years, 1 month ago) by francois
Content type: text/plain
Original Path: magic/lib/mailleur/mailleur/src/mailleur3d.h
File size: 9792 byte(s)
Log Message: 
Ajout du mailleur octal pour XFEM et ajout de fusion de face pour obtenir les cofaces avec vectorisation ameliore

File Contents

# User Rev Content
1 5 //------------------------------------------------------------
2     //------------------------------------------------------------
3     // MAGiC
4     // Jean Christophe Cuilli�e et Vincent FRANCOIS
5     // D�artement de G�ie M�anique - UQTR
6     //------------------------------------------------------------
7     // Le projet MAGIC est un projet de recherche du d�artement
8     // de g�ie m�anique de l'Universit�du Qu�ec �
9     // Trois Rivi�es
10     // Les librairies ne peuvent �re utilis�s sans l'accord
11     // des auteurs (contact : francois@uqtr.ca)
12     //------------------------------------------------------------
13     //------------------------------------------------------------
14     //
15     // mailleur3d.h
16     //
17     //------------------------------------------------------------
18     //------------------------------------------------------------
19     // COPYRIGHT 2000
20     // Version du 02/03/2006 �11H23
21     //------------------------------------------------------------
22     //------------------------------------------------------------
23    
24     #ifndef mailleur3dH
25     #define mailleur3dH
26    
27    
28    
29     #ifdef WINDOWS_VERSION
30     #ifdef BUILT_DLL_MAILLEUR
31     #define DLLPORTMAIL __declspec(dllexport)
32     #else
33     #define DLLPORTMAIL __declspec(dllimport)
34     #endif
35     #else
36     #define DLLPORTMAIL
37     #endif
38    
39    
40    
41     #ifdef BORLANDCPP
42     #include <map.h>
43     #else
44     #include <map>
45     #define CLK_TCK CLOCKS_PER_SEC
46     #endif
47    
48     #include <vector>
49     #include "mailleur.h"
50     #include "mg_front_3D.h"
51     #include "mg_maillage.h"
52     #include "mg_geometrie.h"
53     #include "mg_volume.h"
54     #include "fct_taille_volume.h"
55 francois 78 #include "fct_generateur_3D.h"
56 5 #include "tpl_octree.h"
57    
58    
59     const unsigned int NONFORCE=159;
60     const unsigned int ATTENTE=160;
61 francois 35 const unsigned int PRIORITAIRE=160;
62 5
63     const unsigned int GENERATION_NOEUD=505;
64     const unsigned int NOEUD_EXISTANT=506;
65     const unsigned int FACE_EXISTANTE=506;
66    
67     const unsigned int EXISTE=700;
68     const unsigned int MULTIEXISTE=701;
69     const unsigned int CREATION=702;
70    
71 francois 78 const unsigned int FRONTAL=1000;
72     const unsigned int OCTAL=1001;
73 5
74    
75 francois 78
76 5 class DLLPORTMAIL MAILLEUR3D:public MAILLEUR
77     {
78     public:
79     typedef std::multimap<double,MG_FRONT_3D*,std::less<double> > FRONT;
80 francois 35 typedef std::multimap<double,class M3D_TETRA*,std::less<double> > ORDRE_TETRA;
81     typedef std::map<unsigned long,ORDRE_TETRA::iterator,std::less<unsigned long> > ORDRE_TETRA_PARID;
82 5
83    
84    
85     MAILLEUR3D(MG_MAILLAGE* mgmai,MG_GEOMETRIE *mggeo,MG_VOLUME* mgvol=NULL,FCT_TAILLE* fct_taille=NULL);
86 francois 78 MAILLEUR3D(FCT_GENERATEUR_3D<4> *carte,MG_MAILLAGE* mgmai);
87 5 virtual ~MAILLEUR3D();
88    
89     void maille(void);
90 francois 54 void maille(MG_VOLUME* mgvol,TPL_LISTE_ENTITE<MG_TRIANGLE*> *lsttri=NULL,TPL_LISTE_ENTITE<MG_TETRA*> *lsttet=NULL);
91 5 double cpu[100]; //debug
92     long nbfois[100]; //debug
93    
94     void active_log(char* nomfile);
95     void desactive_log(void);
96 francois 35 void analyse_maillage_obtenu(double &vol);
97 francois 61 void change_niveau_optimisation(int num);
98     int get_niveau_optimisation(void);
99 francois 78 void change_coef(double val);
100     void verifie_conformite_octale(MG_MAILLAGE* mai,MG_MAILLAGE* maisur);
101 5
102     private:
103 francois 78 void maille(FCT_GENERATEUR_3D<4> *carte);
104 5
105 francois 54 void insere_contrainte_tetra(MG_VOLUME* mgvol,TPL_LISTE_ENTITE<MG_TETRA*> *lsttet);
106     void insere_contrainte_triangle(MG_VOLUME* mgvol,TPL_LISTE_ENTITE<MG_TRIANGLE*> *lsttri);
107    
108    
109    
110 5 void initialise_frontiere(MG_VOLUME* mgvol);
111     void cree_octree(void);
112     void initialise_front(MG_VOLUME* mgvol);
113     void progresse_front(MG_VOLUME* mgvol);
114 francois 35 int genere_noeud(MG_VOLUME* mgvol,MG_FRONT_3D* ft,TPL_MAP_ENTITE<MG_NOEUD*> &liste_noeud,std::vector<class CAS_FRONT*> &liste_cas);
115 5 void detruit_element_inactif(void);
116    
117     int traite_front(class MG_VOLUME* mgvol,int type,MG_FRONT_3D* ft,std::vector<class CAS_FRONT*>& liste_cas,std::vector<class CAS_FRONT*>& liste_classe_cas,class CAS_FRONT **cas);
118     int traite_front_adjacent3(class MG_VOLUME* mgvol,MG_FRONT_3D* ft,std::vector<CAS_FRONT*> &liste_cas,std::vector<CAS_FRONT*> &liste_classe_cas,CAS_FRONT **cas);
119     int traite_front_adjacent2(class MG_VOLUME* mgvol,MG_FRONT_3D* ft,std::vector<CAS_FRONT*> &liste_cas,std::vector<CAS_FRONT*> &liste_classe_cas,CAS_FRONT **cas);
120     int traite_front_adjacent1(class MG_VOLUME* mgvol,MG_FRONT_3D* ft,std::vector<CAS_FRONT*> &liste_cas,std::vector<CAS_FRONT*> &liste_classe_cas,CAS_FRONT **cas);
121     int traite_front_lie(class MG_VOLUME* mgvol,MG_FRONT_3D* ft,std::vector<CAS_FRONT*> &liste_cas,std::vector<CAS_FRONT*> &liste_classe_cas,CAS_FRONT **cas);
122     int traite_front_general(class MG_VOLUME* mgvol,MG_FRONT_3D* ft,std::vector<CAS_FRONT*> &liste_cas,std::vector<CAS_FRONT*> &liste_classe_cas,CAS_FRONT **cas);
123    
124     int traite_destruction(MG_FRONT_3D* ft,TPL_MAP_ENTITE<class MG_TRIANGLE*>& liste_intersection);
125     void detruit_tetra(MG_TRIANGLE* tri);
126    
127 francois 54 class MG_FRONT_3D* mise_a_jour_front(MG_TETRA* tet,int cote,int numfront=ATTENTE);
128 5 void mise_a_jour_voisin_front(MG_FRONT_3D* ft);
129    
130     int test_front(void);
131 francois 61 void genere_log(int fincouche);
132     void ini_log(void);
133     void fin_log(void);
134 5
135    
136     double angle_front(MG_FRONT_3D* ft1,MG_FRONT_3D* ft2);
137     double angle_front(MG_FRONT_3D* ft1,MG_SEGMENT* seg);
138     double evaluer_distance_noeud_triangle(double x,double y,double z,MG_TRIANGLE* tri);
139     void insere_cas_front(int type,MG_NOEUD* noeud,int numero_cote,double angle,std::vector<CAS_FRONT*>& liste_cas);
140     void efface_cas_front(std::vector<CAS_FRONT*>& liste_cas,std::vector<CAS_FRONT*> liste_classe_cas[8]);
141    
142    
143     int inter_tetra_triangle(MG_NOEUD* noeud1,MG_NOEUD* noeud2,MG_NOEUD* noeud3,MG_NOEUD* noeud4,MG_NOEUD* noeud5,MG_NOEUD* noeud6,MG_NOEUD* noeud7);
144     int inter_triangle_triangle(MG_NOEUD* noeud1,MG_NOEUD* noeud2,MG_NOEUD* noeud3,MG_NOEUD* noeud4,MG_NOEUD* noeud5,MG_NOEUD* noeud6);
145     int inter_segment_triangle(MG_NOEUD* noeud1,MG_NOEUD* noeud2,MG_NOEUD* noeud3,MG_NOEUD* noeud4,MG_NOEUD* noeud5);
146     int inter_segment_segment(MG_NOEUD* noeud1,MG_NOEUD* noeud2,MG_NOEUD* noeud3,MG_NOEUD* noeud4);
147     int examine_solution(double sol1,double sol2,int type);
148     int noeud_est_triangle(MG_NOEUD* noeud1,MG_NOEUD* noeud2,MG_NOEUD* noeud3,double x,double y,double z);
149    
150    
151     double calcule_volume(MG_NOEUD* noeud1,MG_NOEUD* noeud2,MG_NOEUD* noeud3,MG_NOEUD* noeud4);
152     double calcule_longueur_caracteristique(MG_NOEUD* noeud1,MG_NOEUD* noeud2,MG_NOEUD* noeud3=NULL,MG_NOEUD* noeud4=NULL);
153    
154    
155     class MG_TETRA* insere_tetra(class MG_VOLUME* mgvol,class MG_FRONT_3D* ft,MG_NOEUD* noeud4,int type,TPL_MAP_ENTITE<class MG_TRIANGLE*>& liste_intersection);
156 francois 54 class MG_TRIANGLE* insere_triangle(class MG_VOLUME* mgvol,MG_NOEUD* noeud1,MG_NOEUD* noeud2,MG_NOEUD* noeud3,int origine);
157     class MG_SEGMENT* insere_segment(class MG_VOLUME* mgvol,MG_NOEUD* noeud1,MG_NOEUD* noeud2,int origine);
158     class MG_TETRA* cree_tetra(class MG_VOLUME* mgvol,MG_NOEUD* noeud1,MG_NOEUD* noeud2,MG_NOEUD* noeud3,MG_NOEUD* noeud4,MG_TRIANGLE* triangle1,MG_TRIANGLE* triangle2,MG_TRIANGLE* triangle3,MG_TRIANGLE* triangle4,int origine);
159     class MG_TRIANGLE* cree_triangle(class MG_VOLUME* mgvol,MG_NOEUD* noeud1,MG_NOEUD* noeud2,MG_NOEUD* noeud3,MG_SEGMENT* segment1,MG_SEGMENT* segment2,MG_SEGMENT* segment3,int origine);
160     class MG_SEGMENT* cree_segment(class MG_VOLUME* mgvol,MG_NOEUD* noeud1,MG_NOEUD* noeud2,int origine);
161 5
162 francois 78 void insere_tetra(MG_MAILLAGE* mai,MG_NOEUD* noeud1,MG_NOEUD* noeud2,MG_NOEUD* noeud3,MG_NOEUD* noeud4);
163    
164    
165    
166 5 MG_FRONT_3D* ajouter_front_courant(int numero_front,MG_NOEUD* noeud1,MG_NOEUD* noeud2,MG_NOEUD* noeud3,MG_TRIANGLE* triangle);
167     MG_FRONT_3D* ajouter_front(FRONT& front,MG_NOEUD* noeud1,MG_NOEUD* noeud2,MG_NOEUD* noeud3,MG_TRIANGLE* triangle);
168     void ajouter_front(FRONT& front,MG_FRONT_3D *ft);
169     MG_FRONT_3D* get_front(FRONT& front,unsigned int num);
170     unsigned int get_nb_front(FRONT& front);
171     unsigned int get_nb_front_courant(void);
172     void supprimer_front_en_avancant(MG_FRONT_3D* ft);
173     void supprimer_front_en_reculant(MG_FRONT_3D* ft);
174     void echange_de_front(FRONT& front_original,FRONT& front_destination,MG_FRONT_3D* ft);
175     void echange_de_front(MG_FRONT_3D* ft,int num);
176     void swap_front(FRONT& front_original,FRONT& front_destination);
177     int premier_front_courant(void);
178    
179 francois 35 double calcul_distance_metrique(MG_SEGMENT* seg,int pas=32);
180     double calcul_distance_metrique(double *xyz1,double* xyz2,int pas=32);
181     void ajuste_distance_metrique(double *xyz1,double *xyz2,double longueur_desiree,int pas=32);
182     double calcul_volume_tetra_metrique(MG_TETRA* tet);
183 5
184    
185 francois 54
186 francois 61 void optimise(MG_VOLUME* mgvol,int& nbmauvais);
187 francois 35 void o3d_data(void);
188     void o3d_data2(void);
189     int bouge_point(MG_VOLUME* mgvol,class M3D_NOEUD* noeud,double& crit,double& x,double& y, double& z);
190     void remaille_coquille(MG_NOEUD* noeud1,MG_NOEUD* noeud2, double& crit, class COQUILLE& coque);
191     void ajouter_ordre_tetra(M3D_TETRA* tet,int num=1);
192     void supprimer_ordre_tetra(M3D_TETRA* tet);
193 5
194 francois 35
195    
196 5 MG_MAILLAGE* mg_maillage;
197     MG_GEOMETRIE* mg_geometrie;
198     MG_VOLUME* mg_volume;
199     FCT_TAILLE *metrique;
200 francois 78 int type_mailleur;
201 5
202     FRONT front_courant[11];
203     FRONT front_attente;
204    
205     int activelog;
206     std::string nomlog;
207    
208     TPL_OCTREE<MG_TRIANGLE*,MG_NOEUD*> *octree_de_triangle;
209     TPL_OCTREE<MG_TRIANGLE*,MG_NOEUD*> *octree_de_frontiere;
210     TPL_OCTREE<MG_FRONT_3D*,MG_NOEUD*> *octree_de_front;
211     TPL_MAP_ENTITE<MG_TRIANGLE*> triangle_frontiere;
212     TPL_MAP_ENTITE<MG_NOEUD*> noeud_frontiere;
213     TPL_MAP_ENTITE<MG_NOEUD*> noeud_arete_frontiere;
214    
215 francois 35 ORDRE_TETRA lst_tetra[2];
216     ORDRE_TETRA_PARID lst_tetraid[2];
217     int tab_solution[7][1430][8];
218     int tab_face[7][120][3];
219 francois 61 int niveau_optimisation;
220     FILE* inlog;
221 francois 78 double coef;
222 5
223     };
224    
225    
226    
227    
228 francois 35 class COQUILLE
229     {
230     public:
231     COQUILLE():volume(0) {};
232     ~COQUILLE() {};
233 5
234 francois 35 int taille;
235     MG_TETRA* tet[10];
236     MG_NOEUD* new_tetra[64];
237     double volume;
238     };
239 5
240    
241    
242    
243    
244    
245    
246    
247    
248     //---------------------------------------------------------------------------
249     #endif