CAD4FE_mg_utils.h

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/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       CAD4FE_mg_utils.h
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:56 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
 
#ifndef CAD4FE_mg_utilsH
#define CAD4FE_mg_utilsH
 
#include "CAD4FE_Common_platform.h"
 
#include "mg_maillage.h"
#include "mg_geometrie.h"
#include "tpl_grille.h"
#include "mg_volume.h"
 
 
#define MG_TRIANGLE_GET_NOEUDS(triang,vec_nodes) \
vec_nodes[0] = triang->get_noeud1();  \
vec_nodes[1] = triang->get_noeud2(); \
vec_nodes[2] = triang->get_noeud3();
 
#define MG_TRIANGLE_GET_SEGMENTS(triang,vec_segs) \
vec_segs[0] = triang->get_segment1();  \
vec_segs[1] = triang->get_segment2(); \
vec_segs[2] = triang->get_segment3();
                                                   
#define MG_SEGMENT_GET_NOEUDS(seg,vec_nodes)\
vec_nodes[0] = seg->get_noeud1();\
vec_nodes[1] = seg->get_noeud2();
 
#define MG_SEGMENT_GET_2TRIANGLES(seg,vec_triangles)\
vec_triangles[0]=
 
#define MG_SOMMET_GET_XYZ(s,xyz) xyz[0]=s->get_point()->evaluer(xyz)
 
class  MG_UTILS {
public:
              
 
static void MG_NOEUD_GET_XYZ(MG_NOEUD* node, double * xyz){
 xyz[0]=node->get_x(); xyz[1]=node->get_y(); xyz[2]=node->get_z();
}
static void MG_NOEUD_GET_XYZ(MG_NOEUD* node, float * xyz){
 xyz[0]=node->get_x(); xyz[1]=node->get_y(); xyz[2]=node->get_z();
}
 
static void GetFEVCount(MG_VOLUME * __volume, int * __faceCount,int * __edgeCount,int * __vertexCount);
 
 
        static void MG_MAILLAGE_GetBoundingBox (MG_MAILLAGE * __mesh, double __bbox [6] );
 
        
        static double MG_MAILLAGE_GetVolume (MG_MAILLAGE * __mesh);
        static double MG_MAILLAGE_GetSurface (MG_MAILLAGE * __mesh);
 
 
 
template <class A>
void
TPL_GRILLE_Construct (double __bbox[6], int __nb_steps[3], TPL_GRILLE<A> & __grid)
{
        // Compute the center of the bounding box
        double center[3];
        for (int i=0; i<3; i++)
                center[i] = .5*(__bbox[i] + __bbox[3+i]);
 
        // Compute the size of each cell of the grid
        double step[3];
        for (int i=0; i<3; i++)
                step[i] =  (__bbox[3+i] - __bbox[i])/__nb_steps[i];
 
        // Enlarge the bounding box of the grid
        // by 10 % of cell's size
        double bbox[6];
        for (int i=0; i<6; i++) bbox[i] = __bbox[i];
        for (int i=0; i<3; i++)
                bbox[i+3] += .1*step[i];
 
        // Construct the grid with the enlarged bounding box
        __grid.initialiser(bbox[0], bbox[1], bbox[2], bbox[3], bbox[4], bbox[5], __nb_steps[0], __nb_steps[1], __nb_steps[2]);
}
 
 
};
 
#endif