m3d_tetra.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/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       m3d_tetra.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:55 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
 
 
#include "gestionversion.h"
 
 
#include <math.h>
#include "m3d_tetra.h"
#include "ot_mathematique.h"
#include "mg_definition.h"
 
 
 
 
 
 
M3D_TETRA::M3D_TETRA(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_TRIANGLE* mgtriangle1,class MG_TRIANGLE* mgtriangle2,class MG_TRIANGLE* mgtriangle3,class MG_TRIANGLE* mgtriangle4,int origine):MG_TETRA(num,topo,mgnoeud1,mgnoeud2,mgnoeud3,mgnoeud4,mgtriangle1,mgtriangle2,mgtriangle3,mgtriangle4,origine),etat(MAGIC::MAILLEURFRONTALETAT::ACTIF)//,optimisation(0)
{
calcule_volume();
}
 
M3D_TETRA::M3D_TETRA(MG_ELEMENT_TOPOLOGIQUE* topo,class MG_NOEUD *mgnoeud1,class MG_NOEUD *mgnoeud2,class MG_NOEUD *mgnoeud3,class MG_NOEUD *mgnoeud4,class MG_TRIANGLE* mgtriangle1,class MG_TRIANGLE* mgtriangle2,class MG_TRIANGLE* mgtriangle3,class MG_TRIANGLE* mgtriangle4,int origine):MG_TETRA(topo,mgnoeud1,mgnoeud2,mgnoeud3,mgnoeud4,mgtriangle1,mgtriangle2,mgtriangle3,mgtriangle4,origine),etat(MAGIC::MAILLEURFRONTALETAT::ACTIF)//,optimisation(0)
{
calcule_volume();
}
 
M3D_TETRA::M3D_TETRA(M3D_TETRA& mdd):MG_TETRA(mdd),etat(mdd.etat),volume(mdd.volume)//,optimisation(mdd.optimisation)
{
}
 
M3D_TETRA::~M3D_TETRA()
{
}
 
void M3D_TETRA::etat_detruit(void)
{
etat=MAGIC::MAILLEURFRONTALETAT::DETRUIT;
}
 
void M3D_TETRA::etat_actif(void)
{
etat=MAGIC::MAILLEURFRONTALETAT::ACTIF;
}
 
int M3D_TETRA::get_etat(void)
{
 
return etat;
}
 
int M3D_TETRA::get_type_entite(void)
{
return MAGIC::TYPE_ENTITE::IDM3D_TETRA;
}
 
 
double M3D_TETRA::get_volume(void)
{
return volume;
}
 
 
void M3D_TETRA::change_volume(double vol)
{
volume=vol;
}
 
 
void M3D_TETRA::calcule_volume(void)
{
OT_VECTEUR_3D ab(noeud2->get_x()-noeud1->get_x(),noeud2->get_y()-noeud1->get_y(),noeud2->get_z()-noeud1->get_z());
OT_VECTEUR_3D ac(noeud3->get_x()-noeud1->get_x(),noeud3->get_y()-noeud1->get_y(),noeud3->get_z()-noeud1->get_z());
OT_VECTEUR_3D da(noeud1->get_x()-noeud4->get_x(),noeud1->get_y()-noeud4->get_y(),noeud1->get_z()-noeud4->get_z());
 
OT_VECTEUR_3D pvec=ab&ac;
volume=pvec*da;
volume=fabs(volume)/3.;
 
}
/*
int M3D_TETRA::est_optimise(void)
{
return optimisation;
}
 
void M3D_TETRA::est_optimise(int num)
{
optimisation=num;
}
 */
double M3D_TETRA::get_qualite(void)
{
return qualite;
}
 
void M3D_TETRA::change_qualite(double qu)
{
qualite=qu;
}
 
void M3D_TETRA::enregistrer(std::ostream& o,double version)
{
if (etat!=MAGIC::MAILLEURFRONTALETAT::DETRUIT) MG_TETRA::enregistrer(o,version);  
}