mailleur/src/mailleur_adaptatif.cpp

//------------------------------------------------------------
//------------------------------------------------------------
//  MAGiC
//  Jean Christophe Cuilli�re et Vincent FRANCOIS
//  D�partement de G�nie M�canique - UQTR
//------------------------------------------------------------
//  Le projet MAGIC est un projet de recherche du d�partement
//  de g�nie m�canique de l'Universit� du Qu�bec �
//  Trois Rivi�res
//  Les librairies ne peuvent �tre utilis�es sans l'accord
//  des auteurs (contact : francois@uqtr.ca)
//------------------------------------------------------------
//------------------------------------------------------------
//
//       mailleur0d.cpp
//
//------------------------------------------------------------
//------------------------------------------------------------
//  COPYRIGHT 2000
//  Version du 02/03/2006 � 11H23
//------------------------------------------------------------
//------------------------------------------------------------


#include "gestionversion.h"
#include "mailleur_adaptatif.h"
#include "mg_solution.h"
#include "mg_maillage.h"
#include "mg_gestionnaire.h"
#include "fct_taille.h"
#include "mailleur0d.h"
#include "mailleur1d.h"
#include "mailleur2d.h"
#include "mailleur3d.h"
#include <math.h>

#define DETRUIT 2000
#define NONDETRUIT 3000


MAILLEUR_ADAPTATIF::MAILLEUR_ADAPTATIF(MG_MAILLAGE* mai,MG_MAILLAGE* nvmai,FCT_TAILLE* carte,double coef):MAILLEUR(false),mg_maillage(mai),metrique(carte),limite(coef),nv_maillage(nvmai),activelog(0)
{
}

MAILLEUR_ADAPTATIF::~MAILLEUR_ADAPTATIF()
{
}

void MAILLEUR_ADAPTATIF::maille(MG_GROUPE_TOPOLOGIQUE* mggt)
{
affiche((char*)"  Adaptation de maillage");
char mess[255];
sprintf(mess,"     Destruction des mailles ne respectant pas la carte de taille a %.2f %%",limite*100);
affiche(mess);
LISTE_MG_NOEUD::iterator itnoeud;
for (MG_NOEUD* noeud=mg_maillage->get_premier_noeud(itnoeud);noeud;noeud=mg_maillage->get_suivant_noeud(itnoeud))
        noeud->change_nouveau_numero(NONDETRUIT);
LISTE_MG_SEGMENT::iterator itsegment;
for (MG_SEGMENT* seg=mg_maillage->get_premier_segment(itsegment);seg;seg=mg_maillage->get_suivant_segment(itsegment))
        seg->change_nouveau_numero(NONDETRUIT);
LISTE_MG_TRIANGLE::iterator ittriangle;
for (MG_TRIANGLE* tri=mg_maillage->get_premier_triangle(ittriangle);tri;tri=mg_maillage->get_suivant_triangle(ittriangle))
        tri->change_nouveau_numero(NONDETRUIT);
LISTE_MG_TETRA::iterator ittetra;
for (MG_TETRA* tet=mg_maillage->get_premier_tetra(ittetra);tet;tet=mg_maillage->get_suivant_tetra(ittetra))
        tet->change_nouveau_numero(NONDETRUIT);


MG_SOLUTION* sol=mg_maillage->calculer_carte_taille_reelle((char*)"temp.soltmp");
sol->active_solution(0);
for (MG_NOEUD* noeud=mg_maillage->get_premier_noeud(itnoeud);noeud;noeud=mg_maillage->get_suivant_noeud(itnoeud))
        {
        double *xyz=noeud->get_coord();
        double tenseur[9];
        metrique->evaluer(xyz,tenseur);
        double val=1./sqrt(tenseur[0]);
        double err=fabs((val-noeud->get_solution())/noeud->get_solution());
        if (err>limite)
                {
                int nb_tetra=noeud->get_lien_tetra()->get_nb();
                if (nb_tetra!=0)
                        {
                        for (int j=0;j<nb_tetra;j++)
                                noeud->get_lien_tetra()->get(j)->change_nouveau_numero(DETRUIT);
                        }
                int nb_tri=noeud->get_lien_triangle()->get_nb();
                if (nb_tri!=0)
                        {
                        for (int j=0;j<nb_tri;j++)
                                if (noeud->get_lien_triangle()->get(j)->get_lien_topologie()->get_dimension()==2)
                                        noeud->get_lien_triangle()->get(j)->change_nouveau_numero(DETRUIT);
                        }
                int nb_seg=noeud->get_lien_segment()->get_nb();
                if (nb_seg!=0)
                        {
                        for (int j=0;j<nb_seg;j++)
                                if (noeud->get_lien_segment()->get(j)->get_lien_topologie()->get_dimension()==1)
                                        noeud->get_lien_segment()->get(j)->change_nouveau_numero(DETRUIT);
                        }

                }
        }
delete sol;
for (MG_NOEUD* noeud=mg_maillage->get_premier_noeud(itnoeud);noeud;noeud=mg_maillage->get_suivant_noeud(itnoeud))
        {
        if (noeud->get_lien_topologie()->get_dimension()==1)
                {
                int nb_seg_sur_arete=0;
                int nb_seg=noeud->get_lien_segment()->get_nb();
                for (int j=0;j<nb_seg;j++)
                        if ((noeud->get_lien_segment()->get(j)->get_lien_topologie()->get_dimension()==1) &&
                             (noeud->get_lien_segment()->get(j)->get_nouveau_numero()==NONDETRUIT)) 
                                         nb_seg_sur_arete++;

                if (nb_seg_sur_arete==0)
                        {
                        int nb_tri=noeud->get_lien_triangle()->get_nb();
                        for (int j=0;j<nb_tri;j++)
                        if (noeud->get_lien_triangle()->get(j)->get_lien_topologie()->get_dimension()==2)
                                noeud->get_lien_triangle()->get(j)->change_nouveau_numero(DETRUIT);
                        int nb_tetra=noeud->get_lien_tetra()->get_nb();
                        for (int j=0;j<nb_tetra;j++)
                                noeud->get_lien_tetra()->get(j)->change_nouveau_numero(DETRUIT);
                        }
                }
        if (noeud->get_lien_topologie()->get_dimension()==2)
                {
                int nb_tri_sur_face=0;
                int nb_tri=noeud->get_lien_triangle()->get_nb();
                for (int j=0;j<nb_tri;j++)
                        if ((noeud->get_lien_triangle()->get(j)->get_lien_topologie()->get_dimension()==2) &&
                            (noeud->get_lien_triangle()->get(j)->get_nouveau_numero()==NONDETRUIT) )
                                        nb_tri_sur_face++;

                if (nb_tri_sur_face==0)
                        {
                        int nb_tetra=noeud->get_lien_tetra()->get_nb();
                        for (int j=0;j<nb_tetra;j++)
                                noeud->get_lien_tetra()->get(j)->change_nouveau_numero(DETRUIT);
                        }
                }
}
int nb_sommet=mg_maillage->get_mg_geometrie()->get_nb_mg_sommet();
for (int i=0;i<nb_sommet;i++)
        {
        MG_SOMMET* som=mg_maillage->get_mg_geometrie()->get_mg_sommet(i);
        int nb_noeud=som->get_lien_maillage()->get_nb();
        for (int j=0;j<nb_noeud;j++)
                {
                MG_NOEUD* noeud=(MG_NOEUD*)som->get_lien_maillage()->get(j);
                if (mg_maillage->get_mg_noeudid(noeud->get_id())==NULL) continue;
                if (noeud->get_nouveau_numero()==DETRUIT) continue;
                MG_NOEUD* nv_noeud=nv_maillage->ajouter_mg_noeud(som,noeud->get_x(),noeud->get_y(),noeud->get_z(),MAGIC::ORIGINE::IMPOSE);
                noeud->change_element_reference(nv_noeud);
                nv_noeud->change_element_reference(noeud);
                }

        }
affiche((char*)"     Adaptation maillage 0D");
MAILLEUR0D m0d(nv_maillage,nv_maillage->get_mg_geometrie());
m0d.adapte();
affiche((char*)"     Adaptation maillage 1D");
int nb_arete=mg_maillage->get_mg_geometrie()->get_nb_mg_arete();
for (int i=0;i<nb_arete;i++)
        {
        MG_ARETE* arete=mg_maillage->get_mg_geometrie()->get_mg_arete(i);
        int nb_seg=arete->get_lien_maillage()->get_nb();
        for (int j=0;j<nb_seg;j++)
                {
                MG_SEGMENT* seg=(MG_SEGMENT*)arete->get_lien_maillage()->get(j);
                if (mg_maillage->get_mg_segmentid(seg->get_id())==NULL) continue;
                if (seg->get_nouveau_numero()==DETRUIT) continue;
                MG_NOEUD* noeud1=(MG_NOEUD*)seg->get_noeud1()->get_element_reference();
                MG_NOEUD* noeud2=(MG_NOEUD*)seg->get_noeud2()->get_element_reference();
                if (noeud1==NULL)
                        {
                        noeud1=nv_maillage->ajouter_mg_noeud(arete,seg->get_noeud1()->get_x(),seg->get_noeud1()->get_y(),seg->get_noeud1()->get_z(),MAGIC::ORIGINE::IMPOSE);
                        noeud1->change_element_reference(seg->get_noeud1());
                        seg->get_noeud1()->change_element_reference(noeud1);
                        }
                if (noeud2==NULL)
                        {
                        noeud2=nv_maillage->ajouter_mg_noeud(arete,seg->get_noeud2()->get_x(),seg->get_noeud2()->get_y(),seg->get_noeud2()->get_z(),MAGIC::ORIGINE::IMPOSE);
                        noeud2->change_element_reference(seg->get_noeud2());
                        seg->get_noeud2()->change_element_reference(noeud2);
                        }
                nv_maillage->ajouter_mg_segment(arete,noeud1,noeud2,MAGIC::ORIGINE::IMPOSE);
                }
        MAILLEUR1D m1d(nv_maillage,nv_maillage->get_mg_geometrie(),metrique,arete);
        m1d.change_pas_integration(pas);
        m1d.adapte();
        }
affiche((char*)"     Adaptation maillage 2D");
int nb_face=mg_maillage->get_mg_geometrie()->get_nb_mg_face();
for (int i=0;i<nb_face;i++)
        {
        MG_FACE* face=mg_maillage->get_mg_geometrie()->get_mg_face(i);
        int nb_tri=face->get_lien_maillage()->get_nb();
        TPL_LISTE_ENTITE<MG_TRIANGLE*> lsttri;
        for (int j=0;j<nb_tri;j++)
                {
                MG_TRIANGLE* tri=(MG_TRIANGLE*)face->get_lien_maillage()->get(j);
                if (mg_maillage->get_mg_triangleid(tri->get_id())==NULL) continue;
                if (tri->get_nouveau_numero()==DETRUIT) continue;
                MG_NOEUD* noeud1=tri->get_noeud1();
                MG_NOEUD* noeud2=tri->get_noeud2();
                MG_NOEUD* noeud3=tri->get_noeud3();
                double *xyz1,*xyz2,*xyz3;
                xyz1=noeud1->get_coord();
                xyz2=noeud2->get_coord();
                xyz3=noeud3->get_coord();
                MG_NOEUD* noeudtmp1=new MG_NOEUD(NULL,xyz1[0],xyz1[1],xyz1[2],MAGIC::ORIGINE::IMPOSE);
                MG_NOEUD* noeudtmp2=new MG_NOEUD(NULL,xyz2[0],xyz2[1],xyz2[2],MAGIC::ORIGINE::IMPOSE);
                MG_NOEUD* noeudtmp3=new MG_NOEUD(NULL,xyz3[0],xyz3[1],xyz3[2],MAGIC::ORIGINE::IMPOSE);
                MG_TRIANGLE* tritmp=new MG_TRIANGLE(NULL,noeudtmp1,noeudtmp2,noeudtmp3,NULL,NULL,NULL,MAGIC::ORIGINE::IMPOSE);
                lsttri.ajouter(tritmp);
                }
        MAILLEUR2D m2d(nv_maillage,nv_maillage->get_mg_geometrie(),metrique,face);
        m2d.change_pas_integration(pas);
        m2d.maille(face,NULL,&lsttri);
        nb_tri=lsttri.get_nb();
        for (int j=0;j<nb_tri;j++)
                {
                MG_TRIANGLE* tri=lsttri.get(0);
                MG_NOEUD* noeud1=tri->get_noeud1();
                MG_NOEUD* noeud2=tri->get_noeud2();
                MG_NOEUD* noeud3=tri->get_noeud3();
                lsttri.supprimer(tri);
                delete tri;
                delete noeud1;
                delete noeud2;
                delete noeud3;
                }

        }
affiche((char*)"     Adaptation maillage 3D");
int nb_vol=mg_maillage->get_mg_geometrie()->get_nb_mg_volume();
for (int i=0;i<nb_vol;i++)
        {
        MG_VOLUME* vol=mg_maillage->get_mg_geometrie()->get_mg_volume(i);
        int nb_tet=vol->get_lien_maillage()->get_nb();
        TPL_LISTE_ENTITE<MG_TETRA*> lsttet;
        for (int j=0;j<nb_tet;j++)
                {
                MG_TETRA* tet=(MG_TETRA*)vol->get_lien_maillage()->get(j);
                if (mg_maillage->get_mg_tetraid(tet->get_id())==NULL) continue;
                if (tet->get_nouveau_numero()==DETRUIT) continue;
                MG_NOEUD* noeud1=tet->get_noeud1();
                MG_NOEUD* noeud2=tet->get_noeud2();
                MG_NOEUD* noeud3=tet->get_noeud3();
                MG_NOEUD* noeud4=tet->get_noeud4();
                double *xyz1,*xyz2,*xyz3,*xyz4;
                xyz1=noeud1->get_coord();
                xyz2=noeud2->get_coord();
                xyz3=noeud3->get_coord();
                xyz4=noeud4->get_coord();
                MG_NOEUD* noeudtmp1=new MG_NOEUD(NULL,xyz1[0],xyz1[1],xyz1[2],MAGIC::ORIGINE::IMPOSE);
                MG_NOEUD* noeudtmp2=new MG_NOEUD(NULL,xyz2[0],xyz2[1],xyz2[2],MAGIC::ORIGINE::IMPOSE);
                MG_NOEUD* noeudtmp3=new MG_NOEUD(NULL,xyz3[0],xyz3[1],xyz3[2],MAGIC::ORIGINE::IMPOSE);
                MG_NOEUD* noeudtmp4=new MG_NOEUD(NULL,xyz4[0],xyz4[1],xyz4[2],MAGIC::ORIGINE::IMPOSE);
                MG_TETRA* tetra=new MG_TETRA(NULL,noeudtmp1,noeudtmp2,noeudtmp3,noeudtmp4,NULL,NULL,NULL,NULL,MAGIC::ORIGINE::IMPOSE);
                lsttet.ajouter(tetra);
                }
        MAILLEUR3D m3d(nv_maillage,nv_maillage->get_mg_geometrie(),metrique,vol);
        m3d.active_affichage(affiche2); 
        m3d.change_pas_integration(pas);
        if (activelog) m3d.active_log(nomlog); 
        m3d.maille(vol,NULL,&lsttet);
        nb_tet=lsttet.get_nb();
        for (int j=0;j<nb_tet;j++)
                {
                MG_TETRA* tet=lsttet.get(0);
                MG_NOEUD* noeud1=tet->get_noeud1();
                MG_NOEUD* noeud2=tet->get_noeud2();
                MG_NOEUD* noeud3=tet->get_noeud3();
                MG_NOEUD* noeud4=tet->get_noeud4();
                lsttet.supprimer(tet);
                delete tet;
                delete noeud1;
                delete noeud2;
                delete noeud3;
                delete noeud4;
                }

        }

affiche((char*)"  Fin de l'adaptation de maillage");
}


void MAILLEUR_ADAPTATIF::active_log(char* nomfile)
{
activelog=1;
nomlog=nomfile;
}

void MAILLEUR_ADAPTATIF::desactive_log(void)
{
activelog=0;
}


Revision:	791
Committed:	Fri Mar 18 21:23:00 2016 UTC (9 years, 1 month ago) by francois
Original Path:	*magic/lib/mailleur_auto/src/mailleur_adaptatif.cpp*
File size:	14181 byte(s)
Log Message:	Transferts de la definition des origines dans le namespace MAGiC
#	User	Rev	Content
1	francois	35	//------------------------------------------------------------
2			//------------------------------------------------------------
3			// MAGiC
4	francois	210	// Jean Christophe Cuilli�re et Vincent FRANCOIS
5			// D�partement de G�nie M�canique - UQTR
6	francois	35	//------------------------------------------------------------
7	francois	210	// Le projet MAGIC est un projet de recherche du d�partement
8			// de g�nie m�canique de l'Universit� du Qu�bec �
9			// Trois Rivi�res
10			// Les librairies ne peuvent �tre utilis�es sans l'accord
11	francois	35	// des auteurs (contact : francois@uqtr.ca)
12			//------------------------------------------------------------
13			//------------------------------------------------------------
14			//
15			// mailleur0d.cpp
16			//
17			//------------------------------------------------------------
18			//------------------------------------------------------------
19			// COPYRIGHT 2000
20	francois	210	// Version du 02/03/2006 � 11H23
21	francois	35	//------------------------------------------------------------
22			//------------------------------------------------------------
23
24
25			#include "gestionversion.h"
26			#include "mailleur_adaptatif.h"
27			#include "mg_solution.h"
28			#include "mg_maillage.h"
29	francois	61	#include "mg_gestionnaire.h"
30	francois	35	#include "fct_taille.h"
31			#include "mailleur0d.h"
32			#include "mailleur1d.h"
33			#include "mailleur2d.h"
34	francois	54	#include "mailleur3d.h"
35	francois	35	#include <math.h>
36
37	francois	54	#define DETRUIT 2000
38			#define NONDETRUIT 3000
39	francois	35
40
41	francois	494	MAILLEUR_ADAPTATIF::MAILLEUR_ADAPTATIF(MG_MAILLAGE* mai,MG_MAILLAGE* nvmai,FCT_TAILLE* carte,double coef):MAILLEUR(false),mg_maillage(mai),metrique(carte),limite(coef),nv_maillage(nvmai),activelog(0)
42	francois	35	{
43			}
44
45			MAILLEUR_ADAPTATIF::~MAILLEUR_ADAPTATIF()
46			{
47			}
48
49	francois	210	void MAILLEUR_ADAPTATIF::maille(MG_GROUPE_TOPOLOGIQUE* mggt)
50	francois	35	{
51	francois	479	affiche((char*)" Adaptation de maillage");
52	francois	425	char mess[255];
53			sprintf(mess," Destruction des mailles ne respectant pas la carte de taille a %.2f %%",limite*100);
54			affiche(mess);
55	francois	54	LISTE_MG_NOEUD::iterator itnoeud;
56			for (MG_NOEUD* noeud=mg_maillage->get_premier_noeud(itnoeud);noeud;noeud=mg_maillage->get_suivant_noeud(itnoeud))
57			noeud->change_nouveau_numero(NONDETRUIT);
58			LISTE_MG_SEGMENT::iterator itsegment;
59			for (MG_SEGMENT* seg=mg_maillage->get_premier_segment(itsegment);seg;seg=mg_maillage->get_suivant_segment(itsegment))
60			seg->change_nouveau_numero(NONDETRUIT);
61			LISTE_MG_TRIANGLE::iterator ittriangle;
62			for (MG_TRIANGLE* tri=mg_maillage->get_premier_triangle(ittriangle);tri;tri=mg_maillage->get_suivant_triangle(ittriangle))
63			tri->change_nouveau_numero(NONDETRUIT);
64			LISTE_MG_TETRA::iterator ittetra;
65			for (MG_TETRA* tet=mg_maillage->get_premier_tetra(ittetra);tet;tet=mg_maillage->get_suivant_tetra(ittetra))
66			tet->change_nouveau_numero(NONDETRUIT);
67
68
69	francois	272	MG_SOLUTION* sol=mg_maillage->calculer_carte_taille_reelle((char*)"temp.soltmp");
70	francois	35	sol->active_solution(0);
71	francois	54	for (MG_NOEUD* noeud=mg_maillage->get_premier_noeud(itnoeud);noeud;noeud=mg_maillage->get_suivant_noeud(itnoeud))
72	francois	35	{
73			double *xyz=noeud->get_coord();
74			double tenseur[9];
75			metrique->evaluer(xyz,tenseur);
76			double val=1./sqrt(tenseur[0]);
77			double err=fabs((val-noeud->get_solution())/noeud->get_solution());
78			if (err>limite)
79			{
80			int nb_tetra=noeud->get_lien_tetra()->get_nb();
81			if (nb_tetra!=0)
82			{
83			for (int j=0;j<nb_tetra;j++)
84	francois	54	noeud->get_lien_tetra()->get(j)->change_nouveau_numero(DETRUIT);
85	francois	35	}
86			int nb_tri=noeud->get_lien_triangle()->get_nb();
87			if (nb_tri!=0)
88			{
89			for (int j=0;j<nb_tri;j++)
90			if (noeud->get_lien_triangle()->get(j)->get_lien_topologie()->get_dimension()==2)
91	francois	54	noeud->get_lien_triangle()->get(j)->change_nouveau_numero(DETRUIT);
92	francois	35	}
93			int nb_seg=noeud->get_lien_segment()->get_nb();
94			if (nb_seg!=0)
95			{
96			for (int j=0;j<nb_seg;j++)
97			if (noeud->get_lien_segment()->get(j)->get_lien_topologie()->get_dimension()==1)
98	francois	54	noeud->get_lien_segment()->get(j)->change_nouveau_numero(DETRUIT);
99	francois	35	}
100
101			}
102			}
103			delete sol;
104	francois	54	for (MG_NOEUD* noeud=mg_maillage->get_premier_noeud(itnoeud);noeud;noeud=mg_maillage->get_suivant_noeud(itnoeud))
105	francois	35	{
106			if (noeud->get_lien_topologie()->get_dimension()==1)
107			{
108			int nb_seg_sur_arete=0;
109			int nb_seg=noeud->get_lien_segment()->get_nb();
110			for (int j=0;j<nb_seg;j++)
111	francois	54	if ((noeud->get_lien_segment()->get(j)->get_lien_topologie()->get_dimension()==1) &&
112			(noeud->get_lien_segment()->get(j)->get_nouveau_numero()==NONDETRUIT))
113			nb_seg_sur_arete++;
114	francois	35
115			if (nb_seg_sur_arete==0)
116			{
117			int nb_tri=noeud->get_lien_triangle()->get_nb();
118			for (int j=0;j<nb_tri;j++)
119			if (noeud->get_lien_triangle()->get(j)->get_lien_topologie()->get_dimension()==2)
120	francois	54	noeud->get_lien_triangle()->get(j)->change_nouveau_numero(DETRUIT);
121	francois	35	int nb_tetra=noeud->get_lien_tetra()->get_nb();
122			for (int j=0;j<nb_tetra;j++)
123	francois	54	noeud->get_lien_tetra()->get(j)->change_nouveau_numero(DETRUIT);
124	francois	35	}
125			}
126			if (noeud->get_lien_topologie()->get_dimension()==2)
127			{
128			int nb_tri_sur_face=0;
129			int nb_tri=noeud->get_lien_triangle()->get_nb();
130			for (int j=0;j<nb_tri;j++)
131	francois	54	if ((noeud->get_lien_triangle()->get(j)->get_lien_topologie()->get_dimension()==2) &&
132			(noeud->get_lien_triangle()->get(j)->get_nouveau_numero()==NONDETRUIT) )
133	francois	35	nb_tri_sur_face++;
134
135			if (nb_tri_sur_face==0)
136			{
137			int nb_tetra=noeud->get_lien_tetra()->get_nb();
138			for (int j=0;j<nb_tetra;j++)
139	francois	54	noeud->get_lien_tetra()->get(j)->change_nouveau_numero(DETRUIT);
140	francois	35	}
141			}
142			}
143	francois	54	int nb_sommet=mg_maillage->get_mg_geometrie()->get_nb_mg_sommet();
144			for (int i=0;i<nb_sommet;i++)
145			{
146			MG_SOMMET* som=mg_maillage->get_mg_geometrie()->get_mg_sommet(i);
147			int nb_noeud=som->get_lien_maillage()->get_nb();
148			for (int j=0;j<nb_noeud;j++)
149			{
150			MG_NOEUD* noeud=(MG_NOEUD*)som->get_lien_maillage()->get(j);
151			if (mg_maillage->get_mg_noeudid(noeud->get_id())==NULL) continue;
152			if (noeud->get_nouveau_numero()==DETRUIT) continue;
153	francois	791	MG_NOEUD* nv_noeud=nv_maillage->ajouter_mg_noeud(som,noeud->get_x(),noeud->get_y(),noeud->get_z(),MAGIC::ORIGINE::IMPOSE);
154	francois	54	noeud->change_element_reference(nv_noeud);
155			nv_noeud->change_element_reference(noeud);
156			}
157
158			}
159	francois	479	affiche((char*)" Adaptation maillage 0D");
160	francois	54	MAILLEUR0D m0d(nv_maillage,nv_maillage->get_mg_geometrie());
161	francois	35	m0d.adapte();
162	francois	479	affiche((char*)" Adaptation maillage 1D");
163	francois	35	int nb_arete=mg_maillage->get_mg_geometrie()->get_nb_mg_arete();
164			for (int i=0;i<nb_arete;i++)
165			{
166			MG_ARETE* arete=mg_maillage->get_mg_geometrie()->get_mg_arete(i);
167	francois	54	int nb_seg=arete->get_lien_maillage()->get_nb();
168			for (int j=0;j<nb_seg;j++)
169			{
170			MG_SEGMENT* seg=(MG_SEGMENT*)arete->get_lien_maillage()->get(j);
171			if (mg_maillage->get_mg_segmentid(seg->get_id())==NULL) continue;
172			if (seg->get_nouveau_numero()==DETRUIT) continue;
173			MG_NOEUD* noeud1=(MG_NOEUD*)seg->get_noeud1()->get_element_reference();
174			MG_NOEUD* noeud2=(MG_NOEUD*)seg->get_noeud2()->get_element_reference();
175			if (noeud1==NULL)
176			{
177	francois	791	noeud1=nv_maillage->ajouter_mg_noeud(arete,seg->get_noeud1()->get_x(),seg->get_noeud1()->get_y(),seg->get_noeud1()->get_z(),MAGIC::ORIGINE::IMPOSE);
178	francois	54	noeud1->change_element_reference(seg->get_noeud1());
179			seg->get_noeud1()->change_element_reference(noeud1);
180			}
181			if (noeud2==NULL)
182			{
183	francois	791	noeud2=nv_maillage->ajouter_mg_noeud(arete,seg->get_noeud2()->get_x(),seg->get_noeud2()->get_y(),seg->get_noeud2()->get_z(),MAGIC::ORIGINE::IMPOSE);
184	francois	54	noeud2->change_element_reference(seg->get_noeud2());
185			seg->get_noeud2()->change_element_reference(noeud2);
186			}
187	francois	791	nv_maillage->ajouter_mg_segment(arete,noeud1,noeud2,MAGIC::ORIGINE::IMPOSE);
188	francois	54	}
189	francois	446	MAILLEUR1D m1d(nv_maillage,nv_maillage->get_mg_geometrie(),metrique,arete);
190	francois	532	m1d.change_pas_integration(pas);
191	francois	35	m1d.adapte();
192			}
193	francois	479	affiche((char*)" Adaptation maillage 2D");
194	francois	35	int nb_face=mg_maillage->get_mg_geometrie()->get_nb_mg_face();
195			for (int i=0;i<nb_face;i++)
196			{
197			MG_FACE* face=mg_maillage->get_mg_geometrie()->get_mg_face(i);
198			int nb_tri=face->get_lien_maillage()->get_nb();
199			TPL_LISTE_ENTITE<MG_TRIANGLE*> lsttri;
200			for (int j=0;j<nb_tri;j++)
201			{
202			MG_TRIANGLE* tri=(MG_TRIANGLE*)face->get_lien_maillage()->get(j);
203	francois	54	if (mg_maillage->get_mg_triangleid(tri->get_id())==NULL) continue;
204			if (tri->get_nouveau_numero()==DETRUIT) continue;
205	francois	35	MG_NOEUD* noeud1=tri->get_noeud1();
206			MG_NOEUD* noeud2=tri->get_noeud2();
207			MG_NOEUD* noeud3=tri->get_noeud3();
208			double xyz1,xyz2,*xyz3;
209			xyz1=noeud1->get_coord();
210			xyz2=noeud2->get_coord();
211			xyz3=noeud3->get_coord();
212	francois	791	MG_NOEUD* noeudtmp1=new MG_NOEUD(NULL,xyz1[0],xyz1[1],xyz1[2],MAGIC::ORIGINE::IMPOSE);
213			MG_NOEUD* noeudtmp2=new MG_NOEUD(NULL,xyz2[0],xyz2[1],xyz2[2],MAGIC::ORIGINE::IMPOSE);
214			MG_NOEUD* noeudtmp3=new MG_NOEUD(NULL,xyz3[0],xyz3[1],xyz3[2],MAGIC::ORIGINE::IMPOSE);
215			MG_TRIANGLE* tritmp=new MG_TRIANGLE(NULL,noeudtmp1,noeudtmp2,noeudtmp3,NULL,NULL,NULL,MAGIC::ORIGINE::IMPOSE);
216	francois	35	lsttri.ajouter(tritmp);
217			}
218	francois	446	MAILLEUR2D m2d(nv_maillage,nv_maillage->get_mg_geometrie(),metrique,face);
219	francois	532	m2d.change_pas_integration(pas);
220	francois	35	m2d.maille(face,NULL,&lsttri);
221	francois	54	nb_tri=lsttri.get_nb();
222	francois	35	for (int j=0;j<nb_tri;j++)
223			{
224			MG_TRIANGLE* tri=lsttri.get(0);
225			MG_NOEUD* noeud1=tri->get_noeud1();
226			MG_NOEUD* noeud2=tri->get_noeud2();
227			MG_NOEUD* noeud3=tri->get_noeud3();
228			lsttri.supprimer(tri);
229			delete tri;
230			delete noeud1;
231			delete noeud2;
232			delete noeud3;
233			}
234
235			}
236	francois	479	affiche((char*)" Adaptation maillage 3D");
237	francois	54	int nb_vol=mg_maillage->get_mg_geometrie()->get_nb_mg_volume();
238			for (int i=0;i<nb_vol;i++)
239			{
240			MG_VOLUME* vol=mg_maillage->get_mg_geometrie()->get_mg_volume(i);
241			int nb_tet=vol->get_lien_maillage()->get_nb();
242			TPL_LISTE_ENTITE<MG_TETRA*> lsttet;
243			for (int j=0;j<nb_tet;j++)
244			{
245			MG_TETRA* tet=(MG_TETRA*)vol->get_lien_maillage()->get(j);
246			if (mg_maillage->get_mg_tetraid(tet->get_id())==NULL) continue;
247			if (tet->get_nouveau_numero()==DETRUIT) continue;
248			MG_NOEUD* noeud1=tet->get_noeud1();
249			MG_NOEUD* noeud2=tet->get_noeud2();
250			MG_NOEUD* noeud3=tet->get_noeud3();
251			MG_NOEUD* noeud4=tet->get_noeud4();
252			double xyz1,xyz2,xyz3,xyz4;
253			xyz1=noeud1->get_coord();
254			xyz2=noeud2->get_coord();
255			xyz3=noeud3->get_coord();
256			xyz4=noeud4->get_coord();
257	francois	791	MG_NOEUD* noeudtmp1=new MG_NOEUD(NULL,xyz1[0],xyz1[1],xyz1[2],MAGIC::ORIGINE::IMPOSE);
258			MG_NOEUD* noeudtmp2=new MG_NOEUD(NULL,xyz2[0],xyz2[1],xyz2[2],MAGIC::ORIGINE::IMPOSE);
259			MG_NOEUD* noeudtmp3=new MG_NOEUD(NULL,xyz3[0],xyz3[1],xyz3[2],MAGIC::ORIGINE::IMPOSE);
260			MG_NOEUD* noeudtmp4=new MG_NOEUD(NULL,xyz4[0],xyz4[1],xyz4[2],MAGIC::ORIGINE::IMPOSE);
261			MG_TETRA* tetra=new MG_TETRA(NULL,noeudtmp1,noeudtmp2,noeudtmp3,noeudtmp4,NULL,NULL,NULL,NULL,MAGIC::ORIGINE::IMPOSE);
262	francois	54	lsttet.ajouter(tetra);
263			}
264	francois	446	MAILLEUR3D m3d(nv_maillage,nv_maillage->get_mg_geometrie(),metrique,vol);
265	francois	425	m3d.active_affichage(affiche2);
266	francois	532	m3d.change_pas_integration(pas);
267	francois	61	if (activelog) m3d.active_log(nomlog);
268	francois	54	m3d.maille(vol,NULL,&lsttet);
269			nb_tet=lsttet.get_nb();
270			for (int j=0;j<nb_tet;j++)
271			{
272			MG_TETRA* tet=lsttet.get(0);
273			MG_NOEUD* noeud1=tet->get_noeud1();
274			MG_NOEUD* noeud2=tet->get_noeud2();
275			MG_NOEUD* noeud3=tet->get_noeud3();
276			MG_NOEUD* noeud4=tet->get_noeud4();
277			lsttet.supprimer(tet);
278			delete tet;
279			delete noeud1;
280			delete noeud2;
281			delete noeud3;
282			delete noeud4;
283			}
284	francois	35
285	francois	54	}
286
287	francois	479	affiche((char*)" Fin de l'adaptation de maillage");
288	francois	35	}
289
290
291
292	francois	61	void MAILLEUR_ADAPTATIF::active_log(char* nomfile)
293			{
294			activelog=1;
295			nomlog=nomfile;
296			}
297
298			void MAILLEUR_ADAPTATIF::desactive_log(void)
299			{
300			activelog=0;
301			}
302
303