mg_fast_marching.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/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       mg_fast_marching.h
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:56 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
#ifndef _MGFASTMARCHING_
#define _MGFASTMARCHING_
#include <vector>
#include <mg_definition.h>
#include "fem_maillage.h"
#include "tpl_octree.h"
#include "magic_affiche.h"
 
class FEM_SOLUTION;
class FEM_MAILLAGE;
class FEM_NOEUD;
class FEM_ELEMENT3;
class FEM_ELEMENT2;
class MG_GESTIONNAIRE;
 
 
namespace MAGIC
{
 
const int VERSION_CONSERVE_GRADT=10;
const int VERSION_CONSERVE_DISTANCE=11;
 
class MG_FAST_MARCHING : public MAGIC_AFFICHE
{
public :
    MG_FAST_MARCHING(FEM_SOLUTION* solu,int numchamp);
    MG_FAST_MARCHING(MG_FAST_MARCHING &mdd);
    virtual ~MG_FAST_MARCHING();
    
    //virtual void init(int versiongrad,int versiondistance);
    virtual void init(void);
    virtual void init_noeud(std::vector<FEM_NOEUD*> &lst);
    virtual void init_noeud(char* entite);
    virtual void init_coord(double x,double y,double z)=0;
    virtual void propage(double vitesse=1.0)=0;
    virtual void finalise(void);
    virtual double get_valeur(double x,double y,double z)=0;
protected:
    //virtual void propage2(double vitesse);
    //virtual bool  traite_particulier2(FEM_ELEMENT2 *ele,FEM_NOEUD* no,double v,double *res);
    //virtual void propage3(double vitesse);
    FEM_SOLUTION *sol;
    //virtual int resoudgradTtriangle(FEM_ELEMENT2* ele,double v,FEM_NOEUD** noupdate,double *res);
    //virtual int resoudgradTtriangle_sethian(FEM_ELEMENT2* ele,double v,FEM_NOEUD** noupdate,double *res);
    //virtual int resoudgradTtriangle_perso(FEM_ELEMENT2* ele,double v,FEM_NOEUD** noupdate,double *res);
    //bool noeud_est_dans_espace_element(FEM_ELEMENT2* ele,FEM_NOEUD* noele,FEM_NOEUD* no);
    
    virtual int sign ( double x );
     virtual double fdiv ( double a, double b );
    virtual  double dt2d ( double xi,double xb,double yb,double zb,double xd,double yd,double zd,double tb,double v );
    virtual int Solve_dt2d ( double ax,double ay,double az,double bx,double by,double bz,double dx,double dy,double dz,double at,double bt,double v,double &dt,double &alpha ); 
    virtual int val2d ( double ax,double ay,double az,double bx,double by,double bz,double dx,double dy,double dz,double at,double bt,double v,double& dt,double &alpha );
    virtual double dt3d ( double xi,double eta,double xb,double yb,double zb,double xc,double yc,double zc,double xd,double yd,double zd,double tb,double tc,double v );
    virtual int Solve_dt3d ( double ax,double ay,double az,double bx,double by,double bz,double cx,double cy,double cz,double dx,double dy,double dz,double at,double bt,double ct,double v,double &dt,double &alpha,double &beta );
    virtual int val3d ( double ax,double ay,double az,double bx,double by,double bz,double cx,double cy,double cz,double dx,double dy,double dz,double at,double bt,double ct,double v,double& dt,double &alpha,double &beta );
 
   
   
   
   
    FEM_MAILLAGE *fem;
    int num;
    int nbnoeudinit;
    FEM_NOEUD* last;
    
    
    //int versiongradTtriangle;
    //int respectgrad;
    TPL_OCTREE<FEM_ELEMENT2*,FEM_NOEUD*> oc;
};
 
 
};
 
 
#endif