fem_solution.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/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       fem_solution.h
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:55 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
#ifndef __FEMSOLUTION_
#define __FEMSOLUTION_
 
 
 
 
#include <map>
#include <set>
#include <vector>
 
 
#include "mg_identificateur.h"
#include "mg_definition.h"
 
 
 std::ostream& operator << (std::ostream& o,class MG_SOLUTION& sol);
 
 
class FEM_SOLUTION : public MG_IDENTIFICATEUR
{
public:
    FEM_SOLUTION(class FEM_MAILLAGE *mai,int nb,char* chemin,int code,std::string nomsol,int entite=MAGIC::ENTITE_SOLUTION::ENTITE_NOEUD,int type_solution=MAGIC::TYPE_SOLUTION::SCALAIRE);
    FEM_SOLUTION(unsigned long num,FEM_MAILLAGE *mai,int nb,char* chemin,int code,std::string nomsol,int entite=MAGIC::ENTITE_SOLUTION::ENTITE_NOEUD,int type_solution=MAGIC::TYPE_SOLUTION::SCALAIRE);
    virtual ~FEM_SOLUTION();
    double lire(int i,int j,int coord=0,int num_no=0);
    void ecrire(double val,int i,int j,int coord=0,int num_no=0);
    void change_nb_gauss(int i,int num);
    int get_nb_gauss(int i);
    void active_solution(int num);
    void change_legende(int num,std::string val);
    std::string get_legende(int num);
    int get_nb_champ(void);
    FEM_MAILLAGE* get_maillage(void);
    double get_legende_min(int num=0);
    double get_legende_max(int num=0);
    void efface(void);
    std::string get_nom(void);
    std::string get_nom_fichier(void);
    
    virtual void enregistrer(std::ostream& o,double version);
    virtual void get_fichier_dependant(std::vector<std::string> &liste_fichier);
 
 
    friend std::ostream& operator << (std::ostream& o,FEM_SOLUTION& maillage);
    int get_entite_solution(void);
    int get_dim_solution(void);
    
    double get_moyenne_volumique_champs(class MG_VOLUME* vol,int numchamps,int coord=0); 
    double get_moyenne_volumique_champs(class BOITE_3D &boite_analyse,class MG_VOLUME* vol,int numchamps,int coord=0);// pour les solutions aux point de gauss
    double get_moyenne_surfacique_champs(class MG_FACE* face,int numchamps,int coord=0);
    double get_moyenne_lineique_champs(class MG_ARETE* are,int numchamps,int coord);
    double get_moyenne_volumique_elementaire_champs(int numchamps,class FEM_ELEMENT3* ele,int num,int coord=0);// pour les solutions aux point de gauss
    double get_moyenne_volumique_elementaire_champs(int numchamps,class FEM_ELEMENT3* ele,int num,double& volume,int coord=0);
    double get_moyenne_surfacique_elementaire_champs(int numchamps,class FEM_ELEMENT2* ele,int num,int coord);
    double get_moyenne_lineique_elementaire_champs(int numchamps,class FEM_ELEMENT1* ele,int num,int coord);
    double get_volume_element(class FEM_ELEMENT3* ele,int num);
private:
    void inisolution(char* chemin,int code);
    char *nom_fichier;
    int nb_champs;
    FILE* in;
    FEM_MAILLAGE* femmai;
    double solmin[MAX_TYPE_SOLUTION];
    double solmax[MAX_TYPE_SOLUTION];
    std::string* legende;
    std::string nomsolution;
    int typeentite;
    std::vector<unsigned long> hash;
    int dim_solution;
    bool change_gauss;
};
 
 
 
#endif