mailleur3d_octal.cpp

Aller à la documentation de ce fichier.

//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//  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/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       mailleur3d_octal.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:55 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
 
 
#include "gestionversion.h"
#include "mailleur3d.h"
 
 
int MAILLEUR3D::maille(FCT_TAILLE *metri)
{
    metrique=metri;
    BOITE_3D boite=metri->get_boite_3D();
    double xmin=boite.get_xmin();
    double ymin=boite.get_ymin();
    double zmin=boite.get_zmin();
    double xmax=boite.get_xmax();
    double ymax=boite.get_ymax();
    double zmax=boite.get_zmax();
    TPL_OCTREE_FCT<MG_NOEUD*,FCT_TAILLE > octree;
    octree.change_coefficent_multiplicateur(coef);
    octree.initialiser(*metri,xmin,ymin,zmin,xmax,ymax,zmax);
    octree.equilibre();
    int nb_feuille=octree.get_nb_feuille();
    std::map<int,MG_NOEUD*,std::less<unsigned long> > lst;
    for (int i=0;i<nb_feuille;i++)
    {
        for (int j=0;j<8;j++)
        {
            unsigned long keycode=octree.get_feuille_keycode(i,j);
            if (lst[keycode]==NULL)
            {
                double x,y,z;
                octree.get_xyzsommet_feuille(i,j,x,y,z);
                MG_NOEUD* noeud=new MG_NOEUD(NULL,x,y,z,MAGIC::ORIGINE::MAILLEUR_AUTO);
                mg_maillage->ajouter_mg_noeud(noeud);
                lst[keycode]=noeud;
            }
        }
    }
    for (int i=0;i<nb_feuille;i++)
    {
        TPL_CELLULE_OCTREE<MG_NOEUD*,MG_NOEUD*> *feuille=octree.get_feuille(i);
        double xyz[3];
        feuille->boite.get_centre(xyz);
        MG_NOEUD* noeud=new MG_NOEUD(NULL,xyz[0],xyz[1],xyz[2],MAGIC::ORIGINE::MAILLEUR_AUTO);
        mg_maillage->ajouter_mg_noeud(noeud);
        int face[6]={0,0,0,0,0,0};
        for (int j=0;j<6;j++)
        {
            //OUTIL_OCTREE oo;
            TPL_CELLULE_OCTREE<MG_NOEUD*,MG_NOEUD*> *cell=octree.get_cellule_voisin(feuille,j);
            if (cell!=NULL)
                if (cell->feuille==0)
                    face[j]=1;
        }
        for (int j=0;j<6;j++)
        {
            OUTIL_OCTREE oo;
            int arete[5]={0,0,0,0};
            if (face[j]==1)
            {
                arete[0]=1;
                arete[1]=1;
                arete[2]=1;
                arete[3]=1;
            }
            for (int k=0;k<4;k++)
            {
                if (arete[k]==0)
                {
                    if (face[oo.faceadjacente(j,k)]==1)
                        arete[k]=1;
                    else
                    {
                        int numarete=oo.areteface(j,k);
                        TPL_CELLULE_OCTREE<MG_NOEUD*,MG_NOEUD*> *cell=octree.get_cellule_adjacent(feuille,numarete);
                        if (cell!=NULL)
                            if (cell->feuille==0) arete[k]=1;
                    }
                }
            }
            arete[4]=arete[0]+arete[1]+arete[2]+arete[3];
            if (arete[4]==0)
            {
                unsigned long key1=octree.get_feuille_keycode(i,oo.sommetface(j,0));
                unsigned long key2=octree.get_feuille_keycode(i,oo.sommetface(j,1));
                unsigned long key3=octree.get_feuille_keycode(i,oo.sommetface(j,2));
                unsigned long key4=octree.get_feuille_keycode(i,oo.sommetface(j,3));
                if ( (j==2) || (j==4) )
                {
                    if ( (feuille->get_sontype()==0) ||(feuille->get_sontype()==1) ||(feuille->get_sontype()==6) ||(feuille->get_sontype()==7))
                    {
                        insere_tetra(mg_maillage,lst[key1],lst[key3],lst[key2],noeud);
                        insere_tetra(mg_maillage,lst[key1],lst[key4],lst[key3],noeud);
                    }
                    else
                    {
                        insere_tetra(mg_maillage,lst[key1],lst[key2],lst[key4],noeud);
                        insere_tetra(mg_maillage,lst[key2],lst[key3],lst[key4],noeud);
                    }
                }
                if ( (j==1) || (j==3) )
                {
                    if ( (feuille->get_sontype()==0) ||(feuille->get_sontype()==2) ||(feuille->get_sontype()==5) ||(feuille->get_sontype()==7))
                    {
                        insere_tetra(mg_maillage,lst[key1],lst[key3],lst[key2],noeud);
                        insere_tetra(mg_maillage,lst[key1],lst[key4],lst[key3],noeud);
                    }
                    else
                    {
                        insere_tetra(mg_maillage,lst[key1],lst[key2],lst[key4],noeud);
                        insere_tetra(mg_maillage,lst[key2],lst[key3],lst[key4],noeud);
                    }
                }
                if ( (j==0) || (j==5) )
                {
                    if ( (feuille->get_sontype()==0) ||(feuille->get_sontype()==4) ||(feuille->get_sontype()==7) ||(feuille->get_sontype()==3))
                    {
                        insere_tetra(mg_maillage,lst[key1],lst[key3],lst[key2],noeud);
                        insere_tetra(mg_maillage,lst[key1],lst[key4],lst[key3],noeud);
                    }
                    else
                    {
                        insere_tetra(mg_maillage,lst[key1],lst[key2],lst[key4],noeud);
                        insere_tetra(mg_maillage,lst[key2],lst[key3],lst[key4],noeud);
                    }
                }
            }
            else
            {
                OUTIL_OCTREE oo;
                double xcentre=0.,ycentre=0.,zcentre=0.;
                double x,y,z;
                octree.get_xyzsommet_feuille(i,oo.sommetface(j,0),x,y,z);
                xcentre=xcentre+x;
                ycentre=ycentre+y;
                zcentre=zcentre+z;
                octree.get_xyzsommet_feuille(i,oo.sommetface(j,1),x,y,z);
                xcentre=xcentre+x;
                ycentre=ycentre+y;
                zcentre=zcentre+z;
                octree.get_xyzsommet_feuille(i,oo.sommetface(j,2),x,y,z);
                xcentre=xcentre+x;
                ycentre=ycentre+y;
                zcentre=zcentre+z;
                octree.get_xyzsommet_feuille(i,oo.sommetface(j,3),x,y,z);
                xcentre=xcentre+x;
                ycentre=ycentre+y;
                zcentre=zcentre+z;
                xcentre=0.25*xcentre;
                ycentre=0.25*ycentre;
                zcentre=0.25*zcentre;
                unsigned long key=octree.get_keycode(xcentre,ycentre,zcentre);
                if (lst[key]==NULL)
                {
                    MG_NOEUD* noeud=new MG_NOEUD(NULL,xcentre,ycentre,zcentre,MAGIC::ORIGINE::MAILLEUR_AUTO);
                    mg_maillage->ajouter_mg_noeud(noeud);
                    lst[key]=noeud;
                }
                MG_NOEUD* centroideface=lst[key];
                for (int k=0;k<4;k++)
                {
                    unsigned long key1=octree.get_feuille_keycode(i,oo.sommetarete(oo.areteface(j,k),0));
                    unsigned long key2=octree.get_feuille_keycode(i,oo.sommetarete(oo.areteface(j,k),1));
                    if (arete[k]==0)
                    {
                        insere_tetra(mg_maillage,lst[key1],lst[key2],centroideface,noeud);
                    }
                    else
                    {
                        double x1,y1,z1,x2,y2,z2;
                        octree.get_xyzsommet_feuille(i,oo.sommetarete(oo.areteface(j,k),0),x1,y1,z1);
                        octree.get_xyzsommet_feuille(i,oo.sommetarete(oo.areteface(j,k),1),x2,y2,z2);
                        unsigned long keymilieu=octree.get_keycode(0.5*(x1+x2),0.5*(y1+y2),0.5*(z1+z2));
                        insere_tetra(mg_maillage,lst[key1],lst[keymilieu],centroideface,noeud);
                        insere_tetra(mg_maillage,lst[keymilieu],lst[key2],centroideface,noeud);
 
                    }
                }
            }
 
        }
 
 
    }
return OK;
}
 
 
 
void MAILLEUR3D::insere_tetra(MG_MAILLAGE* mai,MG_NOEUD* noeud1,MG_NOEUD* noeud2,MG_NOEUD* noeud3,MG_NOEUD* noeud4)
{
    double q=OPERATEUR::qualite_tetra(noeud1->get_coord(),noeud2->get_coord(),noeud3->get_coord(),noeud4->get_coord());
    if (q>0.) mai->ajouter_mg_tetra(NULL,noeud1,noeud2,noeud3,noeud4,MAGIC::ORIGINE::MAILLEUR_AUTO);
    else  mai->ajouter_mg_tetra(NULL,noeud2,noeud1,noeud3,noeud4,MAGIC::ORIGINE::MAILLEUR_AUTO);
}
 
 
 
void MAILLEUR3D::change_coef(double val)
{
    coef=val;
}
 
void MAILLEUR3D::verifie_conformite_octale(MG_MAILLAGE* mai,MG_MAILLAGE* maisur)
{
    int nbtri=mai->get_nb_mg_triangle();
    for (int i=0;i<nbtri;i++)
    {
        MG_TRIANGLE* tri=mai->get_mg_triangle(i);
        if (tri->get_lien_tetra()->get_nb()!=1) continue;
        MG_NOEUD* noeud[3];
        noeud[0]=tri->get_noeud1();
        noeud[1]=tri->get_noeud2();
        noeud[2]=tri->get_noeud3();
        MG_NOEUD* n1=new MG_NOEUD(NULL,noeud[0]->get_x(),noeud[0]->get_y(),noeud[0]->get_z(),MAGIC::ORIGINE::MAILLEUR_AUTO);
        MG_NOEUD* n2=new MG_NOEUD(NULL,noeud[1]->get_x(),noeud[1]->get_y(),noeud[1]->get_z(),MAGIC::ORIGINE::MAILLEUR_AUTO);
        MG_NOEUD* n3=new MG_NOEUD(NULL,noeud[2]->get_x(),noeud[2]->get_y(),noeud[2]->get_z(),MAGIC::ORIGINE::MAILLEUR_AUTO);
        maisur->ajouter_mg_noeud(n1);
        maisur->ajouter_mg_noeud(n2);
        maisur->ajouter_mg_noeud(n3);
        maisur->ajouter_mg_triangle(NULL,n1,n2,n3,MAGIC::ORIGINE::MAILLEUR_AUTO);
    }
}