d1/d6e/toibrep_8h_source.html

//####//------------------------------------------------------------

//####//------------------------------------------------------------

//####//  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/

//####//------------------------------------------------------------

//####//------------------------------------------------------------

//####//

//####//       toibrep.h

//####//

//####//------------------------------------------------------------

//####//------------------------------------------------------------

//####//    COPYRIGHT 2000-2024

//####//  jeu 13 jun 2024 11:58:56 EDT

//####//------------------------------------------------------------

//####//------------------------------------------------------------

#ifndef _TOIBREP_

#define _TOIBREP_


#include <string>

#include "tpl_octree.h"

#include "mg_geometrie_outils.h"

#ifdef ENABLE_IBREP

  #include "IBrep.h"

#endif

#include "fem_maillage_outils.h"


class XFEM_TRIANGLE3;


class TOIBREP

{

public:

enum element_classe {INCONNU=1,EXTERIEUR=3,INTERIEUR=2,FRONTIERE=4,XINCONNU=11,XEXTERIEUR=13,XINTERIEUR=12,XFRONTIERE=14};

TOIBREP(class MG_GESTIONNAIRE *g,class MG_GEOMETRIE *ge,class FEM_MAILLAGE* femm,int nbpas,class OT_CPU* compt);

TOIBREP(class MG_GESTIONNAIRE *g,class MG_GEOMETRIE *ge,class FEM_MAILLAGE* femm,int nbpas);

TOIBREP(class MG_GESTIONNAIRE *g,class MG_GEOMETRIE *ge,class FEM_MAILLAGE* femm);

TOIBREP(class MG_GESTIONNAIRE *g,class MG_GEOMETRIE *ge,class FEM_MAILLAGE* femm,class OT_CPU* compt);

~TOIBREP();

#ifdef ENABLE_IBREP

IBrep importer(std::string nomfichier,std::string nomfichieribrep,class MG_GROUPE_TOPOLOGIQUE* mggt=NULL);

#endif

void active_affichage(void (*fonc)(char*));


#ifdef ENABLE_IBREP

IBrep exporter_IBrep(string chemin,class FEM_SOLUTION* solution,class FEM_SOLUTION* solution_ele,MG_GROUPE_TOPOLOGIQUE* mggt=NULL);

#endif

void importer_et_decouper(MG_GROUPE_TOPOLOGIQUE* mggt);


private:

typedef std::multimap<double,FEM_ELEMENT3*,std::less<double> > LISTE_FM_TRI;

typedef std::map<unsigned long,LISTE_FM_TRI::iterator,std::less<unsigned long> > LISTE_FM_TRI_ID;

typedef std::vector<FEM_ELEMENT3*> LISTE_FM;


FEM_MAILLAGE_OUTILS outilfem;

OT_CPU* compteur;

void recherche_arete_tangeante(TPL_MAP_ENTITE<class MG_ARETE*> &lst,TPL_MAP_ENTITE<class MG_ELEMENT_TOPOLOGIQUE*> &lsttopo);

void remplir_trou_tangeant(std::vector<class TOIBREP_POINT*> *lst,MG_ARETE* are,FEM_SOLUTION* solution_ele);

void remplir_trou(std::vector<class TOIBREP_POINT*> *lst,MG_FACE* face,int orientation,TPL_MAP_ENTITE<MG_ELEMENT_TOPOLOGIQUE*> *lstentitetopoface);

void decoupe_tetra_noeud(FEM_ELEMENT3* tet,double* lst,int *nbtet);

int decoupe_tetra(FEM_ELEMENT3* tet,MG_FACE* face);

void traite_arete_tangeante(class MG_ARETE* are,int numsol,class FEM_SOLUTION* solution,FEM_SOLUTION* solution_ele);

void ajouter_liste(LISTE_FM_TRI& lst,LISTE_FM_TRI_ID &lstid,FEM_ELEMENT3* tet,double val);

void supprimer_liste(LISTE_FM_TRI& lst,LISTE_FM_TRI_ID &lstid,FEM_ELEMENT3* tet);

void ajouter_liste(LISTE_FM& lst,FEM_ELEMENT3* tet);

void supprimer_liste(LISTE_FM& lst,FEM_ELEMENT3* tet);

void decoupe_maillage(int num,FEM_SOLUTION* solution,FEM_SOLUTION* solution_ele,int avecsuppression);

void cree_tag_maillage(void);

double get_volume(FEM_ELEMENT3* tet);

void calcul_valeur_sous_element(FEM_ELEMENT3* ele,FEM_ELEMENT3 *xele);

double calcul_distance(class FEM_NOEUD* noeud,class MG_ARETE* are,class TOIBREP_POINT *pt,double &tii,double *xyz,double precision=1e-6);

double calcul_distance(FEM_NOEUD* noeud,class MG_FACE* are,TOIBREP_POINT* pt,double *xyz,double *uv,double precision=1e-6);

double calcul_distance_level_ortho(FEM_NOEUD* no,MG_ARETE* are,TOIBREP_POINT* pt,double precision=1e-6);

double calcul_distance_level_ortho(FEM_NOEUD* no,MG_ARETE* are,TOIBREP_POINT* pt,double &t,double *xyz,double precision=1e-6);

double calculdist(double *n,double x,double y,double z,class FEM_NOEUD* noeud);

void levelsetn(TPL_MAP_ENTITE<class MG_ELEMENT_TOPOLOGIQUE*> *lst,MG_FACE *face,TPL_MAP_ENTITE<MG_ELEMENT_TOPOLOGIQUE*> *lstentitetopo);

void echantillonne_aretes(std::vector<TOIBREP_POINT*> &lst,MG_FACE* face);

void echantillonne_sommets(std::vector<TOIBREP_POINT*> &lst,MG_FACE* face);

void calcullevelsetpremierepasse(MG_FACE* face,int sens,std::vector<TOIBREP_POINT*> *lst,int n1,int n2);

void calcullevelsetdeuxiemepasse(std::vector<TOIBREP_POINT*> *lst);

void etendrelevelset(FEM_SOLUTION* sol,int numsol);

double resoudgradT(FEM_ELEMENT3* tet,int *signe);

int signe_tetra(FEM_ELEMENT3* tet);

TPL_OCTREE<class FEM_ELEMENT3*,class FEM_NOEUD*> octree_tetra;

TPL_OCTREE<class FEM_NOEUD*,class FEM_NOEUD*> octree_noeud;

FEM_NOEUD* inserer_noeud(MG_ELEMENT_TOPOLOGIQUE* topo,double x,double y, double z);

XFEM_TRIANGLE3* inserer_xtriangle(FEM_ELEMENT_MAILLAGE* tet,MG_FACE* face,FEM_NOEUD** tabnoeud);

MG_GEOMETRIE_OUTILS ot;

MG_GEOMETRIE* geo;

FEM_MAILLAGE* mai;

MG_GESTIONNAIRE* gest;

int NPAS;

double longueur_caracteristique;

int affichageactif;

void (*affiche)(char *mess);

void recherche_interieur_face(FEM_ELEMENT3* tet,MG_FACE* face,TPL_MAP_ENTITE<MG_ELEMENT_TOPOLOGIQUE*> *lstentitetopo);

void decoupe_element2_par_element1et0(FEM_ELEMENT3* ele,TPL_MAP_ENTITE<MG_ELEMENT_TOPOLOGIQUE*> *lstentitetopo);

int decoupe_element2(FEM_ELEMENT3* ele,FEM_NOEUD* no);

int decoupe_element(XFEM_ELEMENT2 *ele,FEM_NOEUD* no);

void decoupe_element(int status,FEM_NOEUD* no,FEM_ELEMENT3* ele);

int decoupe_element(int status,FEM_ELEMENT3 *ele,FEM_NOEUD* no,FEM_NOEUD* no1,FEM_NOEUD* no2,FEM_NOEUD* no3,FEM_NOEUD* no4);

void decoupe_element(XFEM_ELEMENT1* seg,class inter_ele_arete *ele_inter);

void decoupe_xtri(FEM_ELEMENT3* tet,int nb_xtri);

 int intersection_arete_triangle(MG_ARETE* are,FEM_NOEUD* no1,FEM_NOEUD *no2,FEM_NOEUD *no3,double *uv,double *xyz);

int decoupe_arete_tetra(MG_ARETE* are,double *xyz,double *xyzcible,double &t,class inter_ele_arete& ele_inter);

int decoupe_segment_xtetra(FEM_ELEMENT3* ele,XFEM_ELEMENT1* seg,double *xyz,double *xyzcible,double &t,inter_ele_arete& xele_inter);

int intersection_segment_triangle(FEM_NOEUD* nos1,FEM_NOEUD* nos2,FEM_NOEUD* no1,FEM_NOEUD *no2,FEM_NOEUD *no3,double *uv,double *xyz);

int test_du_point_milieu(FEM_NOEUD* no1,FEM_NOEUD* no2,FEM_ELEMENT3* tet);

void oriente_tri(XFEM_ELEMENT2* tri,MG_FACE* face);

int inter_segment_segment_plan(XFEM_ELEMENT2 *xfem2,FEM_NOEUD *no1,FEM_NOEUD *no2,FEM_NOEUD *nos1,FEM_NOEUD *nos2,double *xyz);

void testtopo(void);

};


#endif