mg_pyramide.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/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       mg_pyramide.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:53 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
 
 
#include "gestionversion.h"
#include "mg_pyramide.h"
#include "mg_maillage.h"
 
MG_PYRAMIDE::MG_PYRAMIDE(unsigned long num,MG_ELEMENT_TOPOLOGIQUE* topo,class MG_NOEUD *mgnoeud1,class MG_NOEUD *mgnoeud2,class MG_NOEUD *mgnoeud3,class MG_NOEUD *mgnoeud4,class MG_NOEUD *mgnoeud5,class MG_TRIANGLE* mgtriangle1,MG_TRIANGLE* mgtriangle2,MG_TRIANGLE* mgtriangle3,MG_TRIANGLE* mgtriangle4,class MG_QUADRANGLE* mgquadrangle1,int origine):MG_ELEMENT_MAILLAGE(num,topo,origine),noeud1(mgnoeud1),noeud2(mgnoeud2),noeud3(mgnoeud3),noeud4(mgnoeud4),noeud5(mgnoeud5),triangle1(mgtriangle1),triangle2(mgtriangle2),triangle3(mgtriangle3),triangle4(mgtriangle4),quadrangle1(mgquadrangle1)
{
    init_pyramide();
}
 
MG_PYRAMIDE::MG_PYRAMIDE(MG_ELEMENT_TOPOLOGIQUE* topo,class MG_NOEUD *mgnoeud1,class MG_NOEUD *mgnoeud2,class MG_NOEUD *mgnoeud3,class MG_NOEUD *mgnoeud4,class MG_NOEUD *mgnoeud5,class MG_TRIANGLE* mgtriangle1,MG_TRIANGLE* mgtriangle2,MG_TRIANGLE* mgtriangle3,MG_TRIANGLE* mgtriangle4,class MG_QUADRANGLE* mgquadrangle1,int origine):MG_ELEMENT_MAILLAGE(topo,origine),noeud1(mgnoeud1),noeud2(mgnoeud2),noeud3(mgnoeud3),noeud4(mgnoeud4),noeud5(mgnoeud5),triangle1(mgtriangle1),triangle2(mgtriangle2),triangle3(mgtriangle3),triangle4(mgtriangle4),quadrangle1(mgquadrangle1)
{
    init_pyramide();
}
 
 
MG_PYRAMIDE::MG_PYRAMIDE(MG_PYRAMIDE& mdd):MG_ELEMENT_MAILLAGE(mdd),noeud1(mdd.noeud1),noeud2(mdd.noeud2),noeud3(mdd.noeud3),noeud4(mdd.noeud4),noeud5(mdd.noeud5),triangle1(mdd.triangle1),triangle2(mdd.triangle2),triangle3(mdd.triangle3),triangle4(mdd.triangle4),quadrangle1(mdd.quadrangle1)
{
    init_pyramide();
}
 
 
MG_PYRAMIDE* MG_PYRAMIDE::dupliquer(MG_MAILLAGE *mgmai,long decalage)
{
    MG_NOEUD* nd1=mgmai->get_mg_noeudid(noeud1->get_id()+decalage);
    MG_NOEUD* nd2=mgmai->get_mg_noeudid(noeud2->get_id()+decalage);
    MG_NOEUD* nd3=mgmai->get_mg_noeudid(noeud3->get_id()+decalage);
    MG_NOEUD* nd4=mgmai->get_mg_noeudid(noeud4->get_id()+decalage);
    MG_NOEUD* nd5=mgmai->get_mg_noeudid(noeud5->get_id()+decalage);
    MG_TRIANGLE* tri1=mgmai->get_mg_triangleid(triangle1->get_id()+decalage);
    MG_TRIANGLE* tri2=mgmai->get_mg_triangleid(triangle2->get_id()+decalage);
    MG_TRIANGLE* tri3=mgmai->get_mg_triangleid(triangle3->get_id()+decalage);
    MG_TRIANGLE* tri4=mgmai->get_mg_triangleid(triangle4->get_id()+decalage);
    MG_QUADRANGLE* quad1=mgmai->get_mg_quadrangleid(quadrangle1->get_id()+decalage);
    MG_PYRAMIDE* pyr=new MG_PYRAMIDE(this->get_id()+decalage,liaison_topologique,nd1,nd2,nd3,nd4,nd5,tri1,tri2,tri3,tri4,quad1,MAGIC::ORIGINE::DUPLIQUER);
    mgmai->ajouter_mg_pyramide(pyr);
    return pyr;
}
 
MG_PYRAMIDE::~MG_PYRAMIDE()
{
    if (noeud1!=NULL) noeud1->get_lien_pyramide()->supprimer(this);
    if (noeud2!=NULL) noeud2->get_lien_pyramide()->supprimer(this);
    if (noeud3!=NULL) noeud3->get_lien_pyramide()->supprimer(this);
    if (noeud4!=NULL) noeud4->get_lien_pyramide()->supprimer(this);
    if (noeud5!=NULL) noeud5->get_lien_pyramide()->supprimer(this);
    unsigned long noeudminid=std::min(noeud1->get_id(),noeud2->get_id());
    noeudminid=std::min(noeudminid,noeud3->get_id());
    noeudminid=std::min(noeudminid,noeud4->get_id());
    noeudminid=std::min(noeudminid,noeud5->get_id());
    if (noeud1!=NULL)
        if (noeud1->get_id()==noeudminid) noeud1->get_lien_petit_pyramide()->supprimer(this);
    if (noeud2!=NULL)
        if (noeud2->get_id()==noeudminid) noeud2->get_lien_petit_pyramide()->supprimer(this);
    if (noeud3!=NULL)
        if (noeud3->get_id()==noeudminid) noeud3->get_lien_petit_pyramide()->supprimer(this);
    if (noeud4!=NULL)
        if (noeud4->get_id()==noeudminid) noeud4->get_lien_petit_pyramide()->supprimer(this);
    if (noeud5!=NULL)
        if (noeud5->get_id()==noeudminid) noeud5->get_lien_petit_pyramide()->supprimer(this);
    if (triangle1!=NULL) {
        triangle1->decrementer();
        triangle1->get_lien_pyramide()->supprimer(this);
    }
    if (triangle2!=NULL) {
        triangle2->decrementer();
        triangle2->get_lien_pyramide()->supprimer(this);
    }
    if (triangle3!=NULL) {
        triangle3->decrementer();
        triangle3->get_lien_pyramide()->supprimer(this);
    }
    if (triangle4!=NULL) {
        triangle4->decrementer();
        triangle4->get_lien_pyramide()->supprimer(this);
    }
    if (quadrangle1!=NULL) {
        quadrangle1->decrementer();
        quadrangle1->get_lien_pyramide()->supprimer(this);
    }
    if (liaison_topologique!=NULL)
      if (liaison_topologique->get_dimension()==3) liaison_topologique->get_lien_maillage()->supprimer(this);
}
 
void MG_PYRAMIDE::init_pyramide(void)
{
    if (noeud1!=NULL) noeud1->get_lien_pyramide()->ajouter(this);
    if (noeud2!=NULL) noeud2->get_lien_pyramide()->ajouter(this);
    if (noeud3!=NULL) noeud3->get_lien_pyramide()->ajouter(this);
    if (noeud4!=NULL) noeud4->get_lien_pyramide()->ajouter(this);
    if (noeud5!=NULL) noeud5->get_lien_pyramide()->ajouter(this);
    unsigned long noeudminid=std::min(noeud1->get_id(),noeud2->get_id());
    noeudminid=std::min(noeudminid,noeud3->get_id());
    noeudminid=std::min(noeudminid,noeud4->get_id());
    noeudminid=std::min(noeudminid,noeud5->get_id());
    if (noeud1!=NULL)
        if (noeud1->get_id()==noeudminid) noeud1->get_lien_petit_pyramide()->ajouter(this);
    if (noeud2!=NULL)
        if (noeud2->get_id()==noeudminid) noeud2->get_lien_petit_pyramide()->ajouter(this);
    if (noeud3!=NULL)
        if (noeud3->get_id()==noeudminid) noeud3->get_lien_petit_pyramide()->ajouter(this);
    if (noeud4!=NULL)
        if (noeud4->get_id()==noeudminid) noeud4->get_lien_petit_pyramide()->ajouter(this);
    if (noeud5!=NULL)
        if (noeud5->get_id()==noeudminid) noeud5->get_lien_petit_pyramide()->ajouter(this);
    if (triangle1!=NULL) {
        triangle1->incrementer();
        triangle1->get_lien_pyramide()->ajouter(this);
    }
    if (triangle2!=NULL) {
        triangle2->incrementer();
        triangle2->get_lien_pyramide()->ajouter(this);
    }
    if (triangle3!=NULL) {
        triangle3->incrementer();
        triangle3->get_lien_pyramide()->ajouter(this);
    }
    if (triangle4!=NULL) {
        triangle4->incrementer();
        triangle4->get_lien_pyramide()->ajouter(this);
    }
    if (quadrangle1!=NULL) {
        quadrangle1->incrementer();
        quadrangle1->get_lien_pyramide()->ajouter(this);
    }
    
    if (liaison_topologique==NULL) return;
    if (liaison_topologique->get_dimension()==3) liaison_topologique->get_lien_maillage()->ajouter(this);
}
 
 
int MG_PYRAMIDE::get_type_entite(void)
{
    return MAGIC::TYPE_ENTITE::IDMG_PYRAMIDE;
}
 
 
 
 
 
MG_NOEUD* MG_PYRAMIDE::get_noeud1(void)
{
    return noeud1;
}
 
MG_NOEUD* MG_PYRAMIDE::get_noeud2(void)
{
    return noeud2;
}
 
MG_NOEUD* MG_PYRAMIDE::get_noeud3(void)
{
    return noeud3;
}
 
MG_NOEUD* MG_PYRAMIDE::get_noeud4(void)
{
    return noeud4;
}
MG_NOEUD* MG_PYRAMIDE::get_noeud5(void)
{
    return noeud5;
}
 
 
MG_QUADRANGLE* MG_PYRAMIDE::get_quadrangle1(void)
{
    return quadrangle1;
}
 
MG_TRIANGLE* MG_PYRAMIDE::get_triangle1(void)
{
    return triangle1;
}
 
MG_TRIANGLE* MG_PYRAMIDE::get_triangle2(void)
{
    return triangle2;
}
MG_TRIANGLE* MG_PYRAMIDE::get_triangle3(void)
{
    return triangle3;
}
 
MG_TRIANGLE* MG_PYRAMIDE::get_triangle4(void)
{
    return triangle4;
}
 
void MG_PYRAMIDE::inverse_sens(void)
{
MG_NOEUD* noeud=noeud4;
noeud4=noeud2;
noeud2=noeud;
MG_TRIANGLE* triangle=triangle1;
triangle1=triangle4;
triangle4=triangle;
triangle=triangle2;
triangle2=triangle3;
triangle3=triangle;
}
 
 
int MG_PYRAMIDE::get_dimension(void)
{
    return 3;
}
 
 
BOITE_3D MG_PYRAMIDE::get_boite_3D(void)
{
    double coo1[3];
    double coo2[3];
    double coo3[3];
    double coo4[3];
    double coo5[3];
    coo1[0]=noeud1->get_x();
    coo1[1]=noeud1->get_y();
    coo1[2]=noeud1->get_z();
    coo2[0]=noeud2->get_x();
    coo2[1]=noeud2->get_y();
    coo2[2]=noeud2->get_z();
    coo3[0]=noeud3->get_x();
    coo3[1]=noeud3->get_y();
    coo3[2]=noeud3->get_z();
    coo4[0]=noeud4->get_x();
    coo4[1]=noeud4->get_y();
    coo4[2]=noeud4->get_z();
    coo5[0]=noeud5->get_x();
    coo5[1]=noeud5->get_y();
    coo5[2]=noeud5->get_z();
    double xmin,ymin,zmin,xmax,ymax,zmax;
    xmin=std::min(coo1[0],coo2[0]);
    xmin=std::min(xmin,coo3[0]);
    xmin=std::min(xmin,coo4[0]);
    xmin=std::min(xmin,coo5[0]);
    xmax=std::max(coo1[0],coo2[0]);
    xmax=std::max(xmax,coo3[0]);
    xmax=std::max(xmax,coo4[0]);
    xmax=std::max(xmax,coo5[0]);
    ymin=std::min(coo1[1],coo2[1]);
    ymin=std::min(ymin,coo3[1]);
    ymin=std::min(ymin,coo4[1]);
    ymin=std::min(ymin,coo5[1]);
    ymax=std::max(coo1[1],coo2[1]);
    ymax=std::max(ymax,coo3[1]);
    ymax=std::max(ymax,coo4[1]);
    ymax=std::max(ymax,coo5[1]);
    zmin=std::min(coo1[2],coo2[2]);
    zmin=std::min(zmin,coo3[2]);
    zmin=std::min(zmin,coo4[2]);
    zmin=std::min(zmin,coo5[2]);
    zmax=std::max(coo1[2],coo2[2]);
    zmax=std::max(zmax,coo3[2]);
    zmax=std::max(zmax,coo4[2]);
    zmax=std::max(zmax,coo5[2]);
    BOITE_3D boite(xmin,ymin,zmin,xmax,ymax,zmax);
    return boite;
}
void MG_PYRAMIDE::extrapoler_solution_noeud(void)
{
    for (int j=0;j<MAX_TYPE_SOLUTION;j++)  
      {
        noeud1->change_solution(solution[j],j);
        noeud2->change_solution(solution[j],j);
        noeud3->change_solution(solution[j],j);
        noeud4->change_solution(solution[j],j);
        noeud5->change_solution(solution[j],j);
          }
}
 
void MG_PYRAMIDE::enregistrer(std::ostream& o,double version)
{
    if ((liaison_topologique==NULL) && (MG_MAILLAGE::DIMENSIONMAILLAGESANSTOPO!=3)) return;
    if (liaison_topologique==NULL) o << "%" << get_id() << "=PYRAMIDE(NULL,$" << noeud1->get_id() << ",$" << noeud2->get_id() << ",$" << noeud3->get_id()  << ",$" << noeud4->get_id()  << ",$" << noeud5->get_id()  << "," << origine_creation << ");" << std::endl;
    else if (liaison_topologique->get_dimension()==3)  o << "%" << get_id() << "=PYRAMIDE($"<< get_lien_topologie()->get_id() << ",$" << noeud1->get_id() << ",$" << noeud2->get_id() << ",$" << noeud3->get_id()  << ",$" << noeud4->get_id()  << ",$" <<noeud5->get_id() <<  "," << origine_creation << ");" << std::endl;
}