mg_face_element.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_face_element.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:54 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
 
 
#include "gestionversion.h"
#include "mg_face_element.h"
 
 
 
MG_FACE_ELEMENT::MG_FACE_ELEMENT(unsigned long num):MG_FACE("",num,NULL,0)
{
}
 
MG_FACE_ELEMENT::MG_FACE_ELEMENT():MG_FACE("",NULL,0)
{
}
 
MG_FACE_ELEMENT::MG_FACE_ELEMENT(MG_FACE_ELEMENT& mdd):MG_FACE(mdd)
{
}
 
MG_FACE_ELEMENT::~MG_FACE_ELEMENT()
{
}
 
bool MG_FACE_ELEMENT::est_une_topo_element(void)
{
return true;
}
 
 
void MG_FACE_ELEMENT::ajouter_element(MG_ELEMENT_MAILLAGE *ele)
{
ele->change_lien_topologie(this);
}
 
BOITE_3D MG_FACE_ELEMENT::get_boite_3D(void)
{
double xmax=-1e308,ymax=-1e308,zmax=-1e308;
double xmin=1e308,ymin=1e308,zmin=1e308;
TPL_SET<MG_ELEMENT_MAILLAGE*>::ITERATEUR it;   
for (MG_ELEMENT_MAILLAGE* ele=get_lien_maillage()->get_premier(it);ele!=NULL;ele=get_lien_maillage()->get_suivant(it))   
    {
    if (ele->get_type_entite()==MAGIC::TYPE_ENTITE::IDMG_TRIANGLE)
        {
        MG_TRIANGLE* tri=(MG_TRIANGLE*)ele;
        
        double x=tri->get_noeud1()->get_x();
        double y=tri->get_noeud1()->get_y();
        double z=tri->get_noeud1()->get_z();
        xmax=std::max(xmax,x);
        ymax=std::max(ymax,y);
        zmax=std::max(zmax,z);
        xmin=std::min(xmin,x);
        ymin=std::min(ymin,y);
        zmin=std::min(zmin,z);
        x=tri->get_noeud2()->get_x();
        y=tri->get_noeud2()->get_y();
        z=tri->get_noeud2()->get_z();
        xmax=std::max(xmax,x);
        ymax=std::max(ymax,y);
        zmax=std::max(zmax,z);
        xmin=std::min(xmin,x);
        ymin=std::min(ymin,y);
        zmin=std::min(zmin,z);
        x=tri->get_noeud3()->get_x();
        y=tri->get_noeud3()->get_y();
        z=tri->get_noeud3()->get_z();
        xmax=std::max(xmax,x);
        ymax=std::max(ymax,y);
        zmax=std::max(zmax,z);
        xmin=std::min(xmin,x);
        ymin=std::min(ymin,y);
        zmin=std::min(zmin,z);
        }
    if (ele->get_type_entite()==MAGIC::TYPE_ENTITE::IDMG_QUADRANGLE)
        {
        MG_QUADRANGLE* qua=(MG_QUADRANGLE*)ele;
        
        double x=qua->get_noeud1()->get_x();
        double y=qua->get_noeud1()->get_y();
        double z=qua->get_noeud1()->get_z();
        xmax=std::max(xmax,x);
        ymax=std::max(ymax,y);
        zmax=std::max(zmax,z);
        xmin=std::min(xmin,x);
        ymin=std::min(ymin,y);
        zmin=std::min(zmin,z);
        x=qua->get_noeud2()->get_x();
        y=qua->get_noeud2()->get_y();
        z=qua->get_noeud2()->get_z();
        xmax=std::max(xmax,x);
        ymax=std::max(ymax,y);
        zmax=std::max(zmax,z);
        xmin=std::min(xmin,x);
        ymin=std::min(ymin,y);
        zmin=std::min(zmin,z);
        x=qua->get_noeud3()->get_x();
        y=qua->get_noeud3()->get_y();
        z=qua->get_noeud3()->get_z();
        xmax=std::max(xmax,x);
        ymax=std::max(ymax,y);
        zmax=std::max(zmax,z);
        xmin=std::min(xmin,x);
        ymin=std::min(ymin,y);
        zmin=std::min(zmin,z);
        x=qua->get_noeud4()->get_x();
        y=qua->get_noeud4()->get_y();
        z=qua->get_noeud4()->get_z();
        xmax=std::max(xmax,x);
        ymax=std::max(ymax,y);
        zmax=std::max(zmax,z);
        xmin=std::min(xmin,x);
        ymin=std::min(ymin,y);
        zmin=std::min(zmin,z);
        }
        
    }
BOITE_3D boite(xmin,ymin,zmin,xmax,ymax,zmax);
return boite;    
}
 
 
MG_ELEMENT_MAILLAGE* MG_FACE_ELEMENT::get_element(int i)
{
return get_lien_maillage()->get(i);  
}
 
int MG_FACE_ELEMENT::get_nb_element(void)
{
return get_lien_maillage()->get_nb();  
}
 
void  MG_FACE_ELEMENT::enregistrer(std::ostream& o,double version)
{
    int nb=get_nb_ccf();
    int nbele=get_nb_element();
    o <<"%" << get_id() << "=FACE_ELEMENT(" << nbele << ",("; 
    /*for (int i=0;i<nbele;i++)
      if (i!=nbele-1) o << "$" << get_element(i)->get_id() << ","; 
    else o << "$" << get_element(nbele-1)->get_id() ;*/
      o << "),";
    if (version>2.35)
        {
        enregistrer_contraintes(o,version);
        o << ",";
        }
    enregistrer_ccf(o,version);
    o <<   ");" << std::endl;
}