mailleur/src/mailleur1d.cpp

//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//  MAGiC
//####//  Jean Christophe Cuilliere et Vincent FRANCOIS
//####//  Departement de Genie Mecanique - UQTR
//####//------------------------------------------------------------
//####//  MAGIC est un projet de recherche de l equipe ERICCA
//####//  du departement de genie mecanique de l Universite du Quebec a Trois Rivieres
//####//  http://www.uqtr.ca/ericca
//####//  http://www.uqtr.ca/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       mailleur1d.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:55 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------


#include "gestionversion.h"
#include <math.h>
#include <vector>
#include "mailleur1d.h"
#include "tpl_set.h"
#include <fstream>
struct integrale
{
    double ti;
    double li;
};


MAILLEUR1D::MAILLEUR1D(MG_MAILLAGE* mgmai,MG_GEOMETRIE *mggeo,FCT_TAILLE* fct_taille,MG_ARETE* mgarete):MAILLEUR(false),mg_maillage(mgmai),mg_geometrie(mggeo),mg_arete(mgarete),metrique(fct_taille)
{
}


MAILLEUR1D::~MAILLEUR1D()
{
}


int MAILLEUR1D::maille(MG_GROUPE_TOPOLOGIQUE* mggt)
{
    if (mg_arete!=NULL) 
    {
      if(maille(mg_arete)==FAIL) return FAIL;        
    }
    else
    {
        TPL_MAP_ENTITE<MG_ELEMENT_TOPOLOGIQUE*> lst;
        if (mggt!=NULL)
        {
          std::map<MG_ELEMENT_TOPOLOGIQUE*,MG_ELEMENT_TOPOLOGIQUE*>::iterator it;
          for(MG_ELEMENT_TOPOLOGIQUE* ele=mggt->get_premier(it);ele!=NULL;ele=mggt->get_suivant(it))
          {
            lst.ajouter(ele);
            ele->get_topologie_sousjacente(&lst);
          }
        }
        int nb_arete=mg_geometrie->get_nb_mg_arete();
        for (int i=0;i<nb_arete;i++)
        {
            MG_ARETE* mgarete=mg_geometrie->get_mg_arete(i);
            if (mggt!=NULL)
                if (lst.existe(mgarete)==0) continue;
            if(maille(mgarete)==FAIL) return FAIL;
        }
    }
    return OK;
}


int MAILLEUR1D::maille(MG_ARETE* mgarete,double t1,MG_NOEUD* noeud_depart,double t2,MG_NOEUD* noeud_arrivee)
{
    /* calcul de la longueur de la arete dans la metrique */
    if (fabs(t1)<1e-14)
        if (fabs(t2)<1e-14)
        {
            t1=mgarete->get_tmin();
            t2=mgarete->get_tmax();
        }

    std::vector<integrale> tab_integrale;
    double ti,tii,t;
    double longueur_metrique=0.;
    tii=t1;
    double dt=1./pas;
    do
    {
        refresh();
        double coefficient_metrique_ti[9];
        double coefficient_metrique_tii[9];
        double coefficient_metrique_derive_ti[9];
        ti=tii;
        tii=ti+(t2-t1)*dt;
        if (tii>t2) tii=t2;

#define ON_A_ABANDONNE_ET_ON_PREND_DT_EGALE_T2_MOINS_T1_DIVISE_PAR_1000
#ifndef ON_A_ABANDONNE_ET_ON_PREND_DT_EGALE_T2_MOINS_T1_DIVISE_PAR_1000
        int bon_dt=0;
        do
        {
            tii=ti+(t2-t1)*dt;
            if (tii>t2) tii=t2;
            double ddxyz[3];
            double dxyz[3];
            double xyz[3];
            double dxyzii[3];
            double xyzii[3];
            mgarete->deriver_seconde(ti,ddxyz,dxyz,xyz);
            mgarete->evaluer(tii,xyzii);
            mgarete->deriver(tii,dxyzii);
            if (creation_metrique)
            {
                metrique->evaluer(&ti,coefficient_metrique_ti);
                metrique->evaluer(&tii,coefficient_metrique_tii);
                metrique->deriver(&ti,coefficient_metrique_derive_ti);
            }
            else
            {
                metrique->evaluer(xyz,coefficient_metrique_ti);
                metrique->evaluer(xyzii,coefficient_metrique_tii);
                metrique->deriver(xyz,coefficient_metrique_derive_ti);
            }
            double facteur1=dxyz[0]*dxyz[0]*coefficient_metrique_ti[0]+dxyz[0]*dxyz[1]*coefficient_metrique_ti[3]+dxyz[0]*dxyz[2]*coefficient_metrique_ti[6]+dxyz[1]*dxyz[0]*coefficient_metrique_ti[1]+dxyz[1]*dxyz[1]*coefficient_metrique_ti[4]+dxyz[1]*dxyz[2]*coefficient_metrique_ti[7]+dxyz[2]*dxyz[0]*coefficient_metrique_ti[2]+dxyz[2]*dxyz[1]*coefficient_metrique_ti[4]+dxyz[2]*dxyz[2]*coefficient_metrique_ti[8];
            double facteur2=ddxyz[0]*dxyz[0]*coefficient_metrique_ti[0]+dxyz[0]*ddxyz[0]*coefficient_metrique_ti[0]+dxyz[0]*dxyz[0]*coefficient_metrique_derive_ti[0]+
                            ddxyz[0]*dxyz[1]*coefficient_metrique_ti[3]+dxyz[0]*ddxyz[1]*coefficient_metrique_ti[3]+dxyz[0]*dxyz[1]*coefficient_metrique_derive_ti[3]+
                            ddxyz[0]*dxyz[2]*coefficient_metrique_ti[6]+dxyz[0]*ddxyz[2]*coefficient_metrique_ti[6]+dxyz[0]*dxyz[2]*coefficient_metrique_derive_ti[6]+
                            ddxyz[1]*dxyz[0]*coefficient_metrique_ti[1]+dxyz[1]*ddxyz[0]*coefficient_metrique_ti[1]+dxyz[1]*dxyz[0]*coefficient_metrique_derive_ti[1]+
                            ddxyz[1]*dxyz[1]*coefficient_metrique_ti[4]+dxyz[1]*ddxyz[1]*coefficient_metrique_ti[4]+dxyz[1]*dxyz[1]*coefficient_metrique_derive_ti[4]+
                            ddxyz[1]*dxyz[2]*coefficient_metrique_ti[7]+dxyz[1]*ddxyz[2]*coefficient_metrique_ti[7]+dxyz[1]*dxyz[2]*coefficient_metrique_derive_ti[7]+
                            ddxyz[2]*dxyz[0]*coefficient_metrique_ti[2]+dxyz[2]*ddxyz[0]*coefficient_metrique_ti[2]+dxyz[2]*dxyz[0]*coefficient_metrique_derive_ti[2]+
                            ddxyz[2]*dxyz[1]*coefficient_metrique_ti[4]+dxyz[2]*ddxyz[1]*coefficient_metrique_ti[4]+dxyz[2]*dxyz[1]*coefficient_metrique_derive_ti[4]+
                            ddxyz[2]*dxyz[2]*coefficient_metrique_ti[8]+dxyz[2]*ddxyz[2]*coefficient_metrique_ti[8]+dxyz[2]*dxyz[2]*coefficient_metrique_derive_ti[8];
            double facteur3=dxyzii[0]*dxyzii[0]*coefficient_metrique_tii[0]+dxyzii[0]*dxyzii[1]*coefficient_metrique_tii[3]+dxyzii[0]*dxyzii[2]*coefficient_metrique_tii[6]+dxyzii[1]*dxyzii[0]*coefficient_metrique_tii[1]+dxyzii[1]*dxyzii[1]*coefficient_metrique_tii[4]+dxyzii[1]*dxyzii[2]*coefficient_metrique_tii[7]+dxyzii[2]*dxyzii[0]*coefficient_metrique_tii[2]+dxyzii[2]*dxyzii[1]*coefficient_metrique_tii[4]+dxyzii[2]*dxyzii[2]*coefficient_metrique_tii[8];
            double residu=fabs(sqrt(facteur3)-sqrt(facteur1)-1./2./sqrt(facteur1)*facteur2*(tii-ti));
            //if (residu>0.01) dt=dt/2.; else {dt=dt*2;bon_dt=1;}
            //if (dt<1e-4) {dt=1e-4;bon_dt=1;}
            //if (dt>0.25) {dt=0.25;bon_dt=1;}
            bon_dt=1;

        }
        while (bon_dt!=1);
#endif // ON_A_ABANDONNE_ET_ON_PREND_DT_EGALE_T2_MOINS_T1_DIVISE_PAR_1000

        t=0.7886751345*ti+0.2113248654*tii;
        double coord[3];
        double coefficient_metrique[9];
        mgarete->deriver(t,coord);
            double xyz[3];
            mgarete->evaluer(t,xyz);
            metrique->evaluer(xyz,coefficient_metrique);
        double facteur1=coord[0]*coord[0]*coefficient_metrique[0]+coord[0]*coord[1]*coefficient_metrique[3]+coord[0]*coord[2]*coefficient_metrique[6]+coord[1]*coord[0]*coefficient_metrique[1]+coord[1]*coord[1]*coefficient_metrique[4]+coord[1]*coord[2]*coefficient_metrique[7]+coord[2]*coord[0]*coefficient_metrique[2]+coord[2]*coord[1]*coefficient_metrique[4]+coord[2]*coord[2]*coefficient_metrique[8];
        longueur_metrique=longueur_metrique+0.5*(tii-ti)*sqrt(facteur1);
        t=0.7886751345*tii+0.2113248654*ti;
        mgarete->deriver(t,coord);
       //     double xyz[3];
            mgarete->evaluer(t,xyz);
            metrique->evaluer(xyz,coefficient_metrique);
        facteur1=coord[0]*coord[0]*coefficient_metrique[0]+coord[0]*coord[1]*coefficient_metrique[3]+coord[0]*coord[2]*coefficient_metrique[6]+coord[1]*coord[0]*coefficient_metrique[1]+coord[1]*coord[1]*coefficient_metrique[4]+coord[1]*coord[2]*coefficient_metrique[7]+coord[2]*coord[0]*coefficient_metrique[2]+coord[2]*coord[1]*coefficient_metrique[4]+coord[2]*coord[2]*coefficient_metrique[8];
        longueur_metrique=longueur_metrique+0.5*(tii-ti)*sqrt(facteur1);
        integrale pas;
        pas.ti=ti;
        pas.li=longueur_metrique;
        tab_integrale.insert(tab_integrale.end(),pas);
    }
    while (tii<t2);
    int nombre_de_segment=(int)floor(longueur_metrique);
    if (longueur_metrique-floor(longueur_metrique)>0.5) nombre_de_segment++;
    if (nombre_de_segment<1)  nombre_de_segment=1;

    /* discretisation */

    double valeur_cible=longueur_metrique/nombre_de_segment;
    int bon_noeud;
    if (noeud_depart==NULL) bon_noeud=0;
    else bon_noeud=1;
    int numnoeud=0;
    while (bon_noeud==0)
    {
        noeud_depart=(MG_NOEUD*)mgarete->get_cosommet1()->get_sommet()->get_lien_maillage()->get(numnoeud);
        MG_NOEUD* noeudtemp=mg_maillage->get_mg_noeudid(noeud_depart->get_id());
        if (noeudtemp==NULL) numnoeud++;
        else bon_noeud=1;
    }
    if (noeud_arrivee==NULL) bon_noeud=0;
    else bon_noeud=1;
    numnoeud=0;
    while (bon_noeud==0)
    {
        noeud_arrivee=(MG_NOEUD*)mgarete->get_cosommet2()->get_sommet()->get_lien_maillage()->get(numnoeud);
        MG_NOEUD* noeudtemp=mg_maillage->get_mg_noeudid(noeud_arrivee->get_id());
        if (noeudtemp==NULL) numnoeud++;
        else bon_noeud=1;
    }


    MG_NOEUD* noeud_precedent=noeud_depart;

    int nb_segment_cree=0;
    int nb_pas_integrale=tab_integrale.size();
    for (int i=0;i<nb_pas_integrale;i++)
    {
        while ((tab_integrale[i].li>(nb_segment_cree+1)*valeur_cible) && (nb_segment_cree<nombre_de_segment-1))
        {
            refresh();
            double ti=tab_integrale[i].ti;
            double tii=t2;
            double li=0.;
            if (i!=0) li=tab_integrale[i-1].li;
            double lii=tab_integrale[i].li;
            if (i!=nb_pas_integrale-1) tii=tab_integrale[i+1].ti;
            t=ti+(tii-ti)/(lii-li)*((nb_segment_cree+1)*valeur_cible-li);
            double coo[3];
            mgarete->evaluer(t,coo);
            MG_NOEUD* nouveau_noeud=mg_maillage->ajouter_mg_noeud(mgarete,coo[0],coo[1],coo[2],MAGIC::ORIGINE::MAILLEUR_AUTO);
            mg_maillage->ajouter_mg_segment(mgarete,noeud_precedent,nouveau_noeud,MAGIC::ORIGINE::MAILLEUR_AUTO);
            noeud_precedent=nouveau_noeud;
            nb_segment_cree++;
            //ofstream o2("test_1D.mai",ios::out|ios::trunc);
            //o2.precision(16);
            //o2.setf(ios::showpoint);
            //mg_maillage->enregistrer_sous_mesh_1D(o2);
        }
    }
    mg_maillage->ajouter_mg_segment(mgarete,noeud_precedent,noeud_arrivee,MAGIC::ORIGINE::MAILLEUR_AUTO);


    return OK;
}


void MAILLEUR1D::adapte(void)
{
    MG_NOEUD* noeud_dep;
    MG_NOEUD* noeud_arr;
    int nb=mg_arete->get_cosommet1()->get_sommet()->get_lien_maillage()->get_nb();
    for (int i=0;i<nb;i++)
        if (mg_maillage->get_mg_noeudid(mg_arete->get_cosommet1()->get_sommet()->get_lien_maillage()->get(i)->get_id())!=NULL)
            noeud_dep=(MG_NOEUD*)mg_arete->get_cosommet1()->get_sommet()->get_lien_maillage()->get(i);
    nb=mg_arete->get_cosommet2()->get_sommet()->get_lien_maillage()->get_nb();
    for (int i=0;i<nb;i++)
        if (mg_maillage->get_mg_noeudid(mg_arete->get_cosommet2()->get_sommet()->get_lien_maillage()->get(i)->get_id())!=NULL)
            noeud_arr=(MG_NOEUD*)mg_arete->get_cosommet2()->get_sommet()->get_lien_maillage()->get(i);
    int nb_segment=mg_arete->get_lien_maillage()->get_nb();
    double tdeb=mg_arete->get_tmin();
    double tfin=mg_arete->get_tmax();
    std::map<double,MG_NOEUD*,std::less<double> > debut;
    std::map<double,MG_NOEUD*,std::less<double> > fin;
    std::pair<double,MG_NOEUD*> tmp(tfin,noeud_arr);
    debut.insert(tmp);
    std::pair<double,MG_NOEUD*> tmp2(tdeb,noeud_dep);
    fin.insert(tmp2);
    for (int i=0;i<nb_segment;i++)
    {
        MG_SEGMENT* seg=(MG_SEGMENT*)mg_arete->get_lien_maillage()->get(i);
        if (mg_maillage->get_mg_segmentid(seg->get_id())==NULL) continue;
        MG_NOEUD* noeud1=seg->get_noeud1();
        MG_NOEUD* noeud2=seg->get_noeud2();
        double t1,t2;
        if (noeud1==noeud_dep) t1=tdeb;
        else if (noeud1==noeud_arr) t1=tfin;
        else
        {
            double *xyz1=noeud1->get_coord();
            mg_arete->inverser(t1,xyz1);
        }
        if (noeud2==noeud_dep) t2=tdeb;
        else if (noeud2==noeud_arr) t2=tfin;
        else
        {
            double *xyz2=noeud2->get_coord();
            mg_arete->inverser(t2,xyz2);
        }
        if (mg_arete->get_courbe()->est_periodique())
        {
            if (t1<tdeb) t1=t1+mg_arete->get_courbe()->get_periode();
            if (t2<t1) t2=t2+mg_arete->get_courbe()->get_periode();
        }
        std::pair<double,MG_NOEUD*> tmp(t1,noeud1);
        std::pair<double,MG_NOEUD*> tmp2(t2,noeud2);
        debut.insert(tmp);
        fin.insert(tmp2);
    }
    std::map<double,MG_NOEUD*,std::less<double> > :: iterator itdebut;
    std::map<double,MG_NOEUD*,std::less<double> > :: iterator itfin;
    itdebut=debut.begin();
    itfin=fin.begin();
    do
    {
        if ((itdebut->second!=itfin->second) || (nb_segment==0) ) maille(mg_arete,(*itfin).first,(*itfin).second,(*itdebut).first,(*itdebut).second);
        itdebut++;
        itfin++;
    }
    while (itdebut!=debut.end());
}


Revision:	1158
Committed:	Thu Jun 13 22:18:49 2024 UTC (10 months, 4 weeks ago) by francois
Original Path:	*magic/lib/mailleur_auto/src/mailleur1d.cpp*
File size:	13487 byte(s)
Log Message:	compatibilité Ubuntu 22.04 Suppression des refeences à Windows Ajout d'une banière
#	User	Rev	Content
1	francois	1158	//####//------------------------------------------------------------
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			//####// mailleur1d.cpp
15			//####//
16			//####//------------------------------------------------------------
17			//####//------------------------------------------------------------
18			//####// COPYRIGHT 2000-2024
19			//####// jeu 13 jun 2024 11:58:55 EDT
20			//####//------------------------------------------------------------
21			//####//------------------------------------------------------------
22	francois	283
23
24			#include "gestionversion.h"
25			#include <math.h>
26			#include <vector>
27			#include "mailleur1d.h"
28			#include "tpl_set.h"
29			#include <fstream>
30			struct integrale
31			{
32			double ti;
33			double li;
34			};
35
36
37	francois	494	MAILLEUR1D::MAILLEUR1D(MG_MAILLAGE* mgmai,MG_GEOMETRIE mggeo,FCT_TAILLE fct_taille,MG_ARETE* mgarete):MAILLEUR(false),mg_maillage(mgmai),mg_geometrie(mggeo),mg_arete(mgarete),metrique(fct_taille)
38	francois	283	{
39			}
40
41
42
43			MAILLEUR1D::~MAILLEUR1D()
44			{
45			}
46
47
48
49
50	couturad	966	int MAILLEUR1D::maille(MG_GROUPE_TOPOLOGIQUE* mggt)
51	francois	283	{
52	couturad	966	if (mg_arete!=NULL)
53			{
54			if(maille(mg_arete)==FAIL) return FAIL;
55			}
56	francois	283	else
57			{
58			TPL_MAP_ENTITE<MG_ELEMENT_TOPOLOGIQUE*> lst;
59			if (mggt!=NULL)
60			{
61	couturad	906	std::map<MG_ELEMENT_TOPOLOGIQUE,MG_ELEMENT_TOPOLOGIQUE>::iterator it;
62			for(MG_ELEMENT_TOPOLOGIQUE* ele=mggt->get_premier(it);ele!=NULL;ele=mggt->get_suivant(it))
63			{
64			lst.ajouter(ele);
65			ele->get_topologie_sousjacente(&lst);
66			}
67	francois	283	}
68			int nb_arete=mg_geometrie->get_nb_mg_arete();
69			for (int i=0;i<nb_arete;i++)
70			{
71			MG_ARETE* mgarete=mg_geometrie->get_mg_arete(i);
72			if (mggt!=NULL)
73			if (lst.existe(mgarete)==0) continue;
74	couturad	967	if(maille(mgarete)==FAIL) return FAIL;
75	francois	283	}
76			}
77	couturad	966	return OK;
78	francois	283	}
79
80
81	couturad	966	int MAILLEUR1D::maille(MG_ARETE* mgarete,double t1,MG_NOEUD* noeud_depart,double t2,MG_NOEUD* noeud_arrivee)
82	francois	283	{
83			/* calcul de la longueur de la arete dans la metrique */
84			if (fabs(t1)<1e-14)
85			if (fabs(t2)<1e-14)
86			{
87			t1=mgarete->get_tmin();
88			t2=mgarete->get_tmax();
89			}
90
91			std::vector<integrale> tab_integrale;
92			double ti,tii,t;
93			double longueur_metrique=0.;
94			tii=t1;
95	francois	532	double dt=1./pas;
96	francois	283	do
97			{
98			refresh();
99			double coefficient_metrique_ti[9];
100			double coefficient_metrique_tii[9];
101			double coefficient_metrique_derive_ti[9];
102			ti=tii;
103			tii=ti+(t2-t1)*dt;
104			if (tii>t2) tii=t2;
105
106			#define ON_A_ABANDONNE_ET_ON_PREND_DT_EGALE_T2_MOINS_T1_DIVISE_PAR_1000
107			#ifndef ON_A_ABANDONNE_ET_ON_PREND_DT_EGALE_T2_MOINS_T1_DIVISE_PAR_1000
108			int bon_dt=0;
109			do
110			{
111			tii=ti+(t2-t1)*dt;
112			if (tii>t2) tii=t2;
113			double ddxyz[3];
114			double dxyz[3];
115			double xyz[3];
116			double dxyzii[3];
117			double xyzii[3];
118			mgarete->deriver_seconde(ti,ddxyz,dxyz,xyz);
119			mgarete->evaluer(tii,xyzii);
120			mgarete->deriver(tii,dxyzii);
121			if (creation_metrique)
122			{
123			metrique->evaluer(&ti,coefficient_metrique_ti);
124			metrique->evaluer(&tii,coefficient_metrique_tii);
125			metrique->deriver(&ti,coefficient_metrique_derive_ti);
126			}
127			else
128			{
129			metrique->evaluer(xyz,coefficient_metrique_ti);
130			metrique->evaluer(xyzii,coefficient_metrique_tii);
131			metrique->deriver(xyz,coefficient_metrique_derive_ti);
132			}
133			double facteur1=dxyz[0]dxyz[0]coefficient_metrique_ti[0]+dxyz[0]dxyz[1]coefficient_metrique_ti[3]+dxyz[0]dxyz[2]coefficient_metrique_ti[6]+dxyz[1]dxyz[0]coefficient_metrique_ti[1]+dxyz[1]dxyz[1]coefficient_metrique_ti[4]+dxyz[1]dxyz[2]coefficient_metrique_ti[7]+dxyz[2]dxyz[0]coefficient_metrique_ti[2]+dxyz[2]dxyz[1]coefficient_metrique_ti[4]+dxyz[2]dxyz[2]coefficient_metrique_ti[8];
134			double facteur2=ddxyz[0]dxyz[0]coefficient_metrique_ti[0]+dxyz[0]ddxyz[0]coefficient_metrique_ti[0]+dxyz[0]dxyz[0]coefficient_metrique_derive_ti[0]+
135			ddxyz[0]dxyz[1]coefficient_metrique_ti[3]+dxyz[0]ddxyz[1]coefficient_metrique_ti[3]+dxyz[0]dxyz[1]coefficient_metrique_derive_ti[3]+
136			ddxyz[0]dxyz[2]coefficient_metrique_ti[6]+dxyz[0]ddxyz[2]coefficient_metrique_ti[6]+dxyz[0]dxyz[2]coefficient_metrique_derive_ti[6]+
137			ddxyz[1]dxyz[0]coefficient_metrique_ti[1]+dxyz[1]ddxyz[0]coefficient_metrique_ti[1]+dxyz[1]dxyz[0]coefficient_metrique_derive_ti[1]+
138			ddxyz[1]dxyz[1]coefficient_metrique_ti[4]+dxyz[1]ddxyz[1]coefficient_metrique_ti[4]+dxyz[1]dxyz[1]coefficient_metrique_derive_ti[4]+
139			ddxyz[1]dxyz[2]coefficient_metrique_ti[7]+dxyz[1]ddxyz[2]coefficient_metrique_ti[7]+dxyz[1]dxyz[2]coefficient_metrique_derive_ti[7]+
140			ddxyz[2]dxyz[0]coefficient_metrique_ti[2]+dxyz[2]ddxyz[0]coefficient_metrique_ti[2]+dxyz[2]dxyz[0]coefficient_metrique_derive_ti[2]+
141			ddxyz[2]dxyz[1]coefficient_metrique_ti[4]+dxyz[2]ddxyz[1]coefficient_metrique_ti[4]+dxyz[2]dxyz[1]coefficient_metrique_derive_ti[4]+
142			ddxyz[2]dxyz[2]coefficient_metrique_ti[8]+dxyz[2]ddxyz[2]coefficient_metrique_ti[8]+dxyz[2]dxyz[2]coefficient_metrique_derive_ti[8];
143			double facteur3=dxyzii[0]dxyzii[0]coefficient_metrique_tii[0]+dxyzii[0]dxyzii[1]coefficient_metrique_tii[3]+dxyzii[0]dxyzii[2]coefficient_metrique_tii[6]+dxyzii[1]dxyzii[0]coefficient_metrique_tii[1]+dxyzii[1]dxyzii[1]coefficient_metrique_tii[4]+dxyzii[1]dxyzii[2]coefficient_metrique_tii[7]+dxyzii[2]dxyzii[0]coefficient_metrique_tii[2]+dxyzii[2]dxyzii[1]coefficient_metrique_tii[4]+dxyzii[2]dxyzii[2]coefficient_metrique_tii[8];
144			double residu=fabs(sqrt(facteur3)-sqrt(facteur1)-1./2./sqrt(facteur1)facteur2(tii-ti));
145			//if (residu>0.01) dt=dt/2.; else {dt=dt*2;bon_dt=1;}
146			//if (dt<1e-4) {dt=1e-4;bon_dt=1;}
147			//if (dt>0.25) {dt=0.25;bon_dt=1;}
148			bon_dt=1;
149
150			}
151			while (bon_dt!=1);
152			#endif // ON_A_ABANDONNE_ET_ON_PREND_DT_EGALE_T2_MOINS_T1_DIVISE_PAR_1000
153
154			t=0.7886751345ti+0.2113248654tii;
155			double coord[3];
156			double coefficient_metrique[9];
157			mgarete->deriver(t,coord);
158			double xyz[3];
159			mgarete->evaluer(t,xyz);
160			metrique->evaluer(xyz,coefficient_metrique);
161			double facteur1=coord[0]coord[0]coefficient_metrique[0]+coord[0]coord[1]coefficient_metrique[3]+coord[0]coord[2]coefficient_metrique[6]+coord[1]coord[0]coefficient_metrique[1]+coord[1]coord[1]coefficient_metrique[4]+coord[1]coord[2]coefficient_metrique[7]+coord[2]coord[0]coefficient_metrique[2]+coord[2]coord[1]coefficient_metrique[4]+coord[2]coord[2]coefficient_metrique[8];
162			longueur_metrique=longueur_metrique+0.5(tii-ti)sqrt(facteur1);
163			t=0.7886751345tii+0.2113248654ti;
164			mgarete->deriver(t,coord);
165	francois	446	// double xyz[3];
166	francois	283	mgarete->evaluer(t,xyz);
167			metrique->evaluer(xyz,coefficient_metrique);
168			facteur1=coord[0]coord[0]coefficient_metrique[0]+coord[0]coord[1]coefficient_metrique[3]+coord[0]coord[2]coefficient_metrique[6]+coord[1]coord[0]coefficient_metrique[1]+coord[1]coord[1]coefficient_metrique[4]+coord[1]coord[2]coefficient_metrique[7]+coord[2]coord[0]coefficient_metrique[2]+coord[2]coord[1]coefficient_metrique[4]+coord[2]coord[2]coefficient_metrique[8];
169			longueur_metrique=longueur_metrique+0.5(tii-ti)sqrt(facteur1);
170			integrale pas;
171			pas.ti=ti;
172			pas.li=longueur_metrique;
173			tab_integrale.insert(tab_integrale.end(),pas);
174			}
175			while (tii<t2);
176			int nombre_de_segment=(int)floor(longueur_metrique);
177			if (longueur_metrique-floor(longueur_metrique)>0.5) nombre_de_segment++;
178			if (nombre_de_segment<1) nombre_de_segment=1;
179
180			/* discretisation */
181
182			double valeur_cible=longueur_metrique/nombre_de_segment;
183			int bon_noeud;
184			if (noeud_depart==NULL) bon_noeud=0;
185			else bon_noeud=1;
186			int numnoeud=0;
187			while (bon_noeud==0)
188			{
189			noeud_depart=(MG_NOEUD*)mgarete->get_cosommet1()->get_sommet()->get_lien_maillage()->get(numnoeud);
190			MG_NOEUD* noeudtemp=mg_maillage->get_mg_noeudid(noeud_depart->get_id());
191			if (noeudtemp==NULL) numnoeud++;
192			else bon_noeud=1;
193			}
194			if (noeud_arrivee==NULL) bon_noeud=0;
195			else bon_noeud=1;
196			numnoeud=0;
197			while (bon_noeud==0)
198			{
199			noeud_arrivee=(MG_NOEUD*)mgarete->get_cosommet2()->get_sommet()->get_lien_maillage()->get(numnoeud);
200			MG_NOEUD* noeudtemp=mg_maillage->get_mg_noeudid(noeud_arrivee->get_id());
201			if (noeudtemp==NULL) numnoeud++;
202			else bon_noeud=1;
203			}
204
205
206
207			MG_NOEUD* noeud_precedent=noeud_depart;
208
209			int nb_segment_cree=0;
210			int nb_pas_integrale=tab_integrale.size();
211			for (int i=0;i<nb_pas_integrale;i++)
212			{
213			while ((tab_integrale[i].li>(nb_segment_cree+1)*valeur_cible) && (nb_segment_cree<nombre_de_segment-1))
214			{
215			refresh();
216			double ti=tab_integrale[i].ti;
217			double tii=t2;
218			double li=0.;
219			if (i!=0) li=tab_integrale[i-1].li;
220			double lii=tab_integrale[i].li;
221			if (i!=nb_pas_integrale-1) tii=tab_integrale[i+1].ti;
222			t=ti+(tii-ti)/(lii-li)((nb_segment_cree+1)valeur_cible-li);
223			double coo[3];
224			mgarete->evaluer(t,coo);
225	francois	791	MG_NOEUD* nouveau_noeud=mg_maillage->ajouter_mg_noeud(mgarete,coo[0],coo[1],coo[2],MAGIC::ORIGINE::MAILLEUR_AUTO);
226			mg_maillage->ajouter_mg_segment(mgarete,noeud_precedent,nouveau_noeud,MAGIC::ORIGINE::MAILLEUR_AUTO);
227	francois	283	noeud_precedent=nouveau_noeud;
228			nb_segment_cree++;
229			//ofstream o2("test_1D.mai",ios::out\|ios::trunc);
230			//o2.precision(16);
231			//o2.setf(ios::showpoint);
232			//mg_maillage->enregistrer_sous_mesh_1D(o2);
233			}
234			}
235	francois	791	mg_maillage->ajouter_mg_segment(mgarete,noeud_precedent,noeud_arrivee,MAGIC::ORIGINE::MAILLEUR_AUTO);
236	francois	283
237
238
239	couturad	966	return OK;
240	francois	283	}
241
242
243
244			void MAILLEUR1D::adapte(void)
245			{
246			MG_NOEUD* noeud_dep;
247			MG_NOEUD* noeud_arr;
248			int nb=mg_arete->get_cosommet1()->get_sommet()->get_lien_maillage()->get_nb();
249			for (int i=0;i<nb;i++)
250			if (mg_maillage->get_mg_noeudid(mg_arete->get_cosommet1()->get_sommet()->get_lien_maillage()->get(i)->get_id())!=NULL)
251			noeud_dep=(MG_NOEUD*)mg_arete->get_cosommet1()->get_sommet()->get_lien_maillage()->get(i);
252			nb=mg_arete->get_cosommet2()->get_sommet()->get_lien_maillage()->get_nb();
253			for (int i=0;i<nb;i++)
254			if (mg_maillage->get_mg_noeudid(mg_arete->get_cosommet2()->get_sommet()->get_lien_maillage()->get(i)->get_id())!=NULL)
255			noeud_arr=(MG_NOEUD*)mg_arete->get_cosommet2()->get_sommet()->get_lien_maillage()->get(i);
256			int nb_segment=mg_arete->get_lien_maillage()->get_nb();
257			double tdeb=mg_arete->get_tmin();
258			double tfin=mg_arete->get_tmax();
259			std::map<double,MG_NOEUD*,std::less<double> > debut;
260			std::map<double,MG_NOEUD*,std::less<double> > fin;
261			std::pair<double,MG_NOEUD*> tmp(tfin,noeud_arr);
262			debut.insert(tmp);
263			std::pair<double,MG_NOEUD*> tmp2(tdeb,noeud_dep);
264			fin.insert(tmp2);
265			for (int i=0;i<nb_segment;i++)
266			{
267			MG_SEGMENT* seg=(MG_SEGMENT*)mg_arete->get_lien_maillage()->get(i);
268			if (mg_maillage->get_mg_segmentid(seg->get_id())==NULL) continue;
269			MG_NOEUD* noeud1=seg->get_noeud1();
270			MG_NOEUD* noeud2=seg->get_noeud2();
271			double t1,t2;
272			if (noeud1==noeud_dep) t1=tdeb;
273			else if (noeud1==noeud_arr) t1=tfin;
274			else
275			{
276			double *xyz1=noeud1->get_coord();
277			mg_arete->inverser(t1,xyz1);
278			}
279			if (noeud2==noeud_dep) t2=tdeb;
280			else if (noeud2==noeud_arr) t2=tfin;
281			else
282			{
283			double *xyz2=noeud2->get_coord();
284			mg_arete->inverser(t2,xyz2);
285			}
286			if (mg_arete->get_courbe()->est_periodique())
287			{
288			if (t1<tdeb) t1=t1+mg_arete->get_courbe()->get_periode();
289			if (t2<t1) t2=t2+mg_arete->get_courbe()->get_periode();
290			}
291			std::pair<double,MG_NOEUD*> tmp(t1,noeud1);
292			std::pair<double,MG_NOEUD*> tmp2(t2,noeud2);
293			debut.insert(tmp);
294			fin.insert(tmp2);
295			}
296			std::map<double,MG_NOEUD*,std::less<double> > :: iterator itdebut;
297			std::map<double,MG_NOEUD*,std::less<double> > :: iterator itfin;
298			itdebut=debut.begin();
299			itfin=fin.begin();
300			do
301			{
302			if ((itdebut->second!=itfin->second) \|\| (nb_segment==0) ) maille(mg_arete,(itfin).first,(itfin).second,(itdebut).first,(itdebut).second);
303			itdebut++;
304			itfin++;
305			}
306			while (itdebut!=debut.end());
307			}
308
309
310
311