dly_triangle.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/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       dly_triangle.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:55 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
 
 
#include "gestionversion.h"
#include <stdlib.h>
 
#include "dly_noeud.h"
#include "dly_triangle.h"
#include "dly_tetra.h"
#include "ot_mathematique.h"
#include <math.h>
 
 
 
DLY_TRIANGLE::DLY_TRIANGLE(class DLY_NOEUD *no1,class DLY_NOEUD *no2,class DLY_NOEUD *no3):noeud1(no1),noeud2(no2),noeud3(no3),frontiere(-1)
{
}
DLY_TRIANGLE::DLY_TRIANGLE(DLY_TRIANGLE& mdd):noeud1(mdd.noeud1),noeud2(mdd.noeud2),noeud3(mdd.noeud3),frontiere(mdd.frontiere)
{
}
DLY_TRIANGLE::~DLY_TRIANGLE()
{
}
 
DLY_NOEUD* DLY_TRIANGLE::get_noeud1(void)
{
return noeud1;
}
DLY_NOEUD* DLY_TRIANGLE::get_noeud2(void)
{
return noeud2;
}
DLY_NOEUD* DLY_TRIANGLE::get_noeud3(void)
{
return noeud3;
}
 
void DLY_TRIANGLE::insere_lien_tetra(DLY_TETRA* tet)
{
lsttet.push_back(tet);
}
 
 
int DLY_TRIANGLE::tet_est_present(DLY_TETRA* tet)
{
for (std::vector<DLY_TETRA*>::iterator it=lsttet.begin();it!=lsttet.end();it++)
    if ((*it)==tet) return 1;return 0;
}
 
int DLY_TRIANGLE::get_nb_tetra(void)
{
return lsttet.size();
}
 
DLY_TETRA* DLY_TRIANGLE::get_tetra(int i)
{
return lsttet[i];
}
 
int DLY_TRIANGLE::est_coplanaire(double x,double y, double z)
{
double xyz1[3],xyz2[3],xyz3[3];
noeud1->get_coord(xyz1);
noeud2->get_coord(xyz2);
noeud3->get_coord(xyz3);
double xyz[3]={x,y,z};
OT_VECTEUR_3D vec1(xyz1,xyz2);
OT_VECTEUR_3D vec2(xyz1,xyz3);
OT_VECTEUR_3D n=vec1&vec2;
OT_VECTEUR_3D vec(xyz1,xyz);
n.norme();
vec.norme();
double ps=n*vec;
int resultat=0;
if (fabs(ps)<0.0000000001) resultat=1;
return resultat;
}
 
int DLY_TRIANGLE::est_frontiere(void )
{
if (frontiere==-1)
  {
  int nbtetra=0;
  for (int j=0;j<lsttet.size();j++)
  if ((lsttet[j])->get_feuille()==1) nbtetra++;
  if (nbtetra==1) frontiere=1; else frontiere=0;
  }
return frontiere;
}