geometrie/src/fem_solution.cpp

//------------------------------------------------------------
//------------------------------------------------------------
//  MAGiC
//  Jean Christophe Cuilli�re et Vincent FRANCOIS
//  D�partement de G�nie M�canique - UQTR
//------------------------------------------------------------
//  Le projet MAGIC est un projet de recherche du d�partement
//  de g�nie m�canique de l'Universit� du Qu�bec �
//  Trois Rivi�res
//  Les librairies ne peuvent �tre utilis�es sans l'accord
//  des auteurs (contact : francois@uqtr.ca)
//------------------------------------------------------------
//------------------------------------------------------------
//
//       fem_solution.cpp
//
//------------------------------------------------------------
//------------------------------------------------------------
//  COPYRIGHT 2000
//  Version du 02/03/2006 � 11H22
//------------------------------------------------------------
//------------------------------------------------------------


#include "gestionversion.h"
#include <math.h>
#include <iostream>
#include <stdlib.h>
#include <string.h>

#include "fem_solution.h"
#include "fem_maillage.h"
#include "mg_geometrie.h"
#include "ot_mathematique.h"


#include <stdio.h> // pour la fonction remove() qui remplace DeleteFile


FEM_SOLUTION::FEM_SOLUTION(FEM_MAILLAGE* mai,int nb,char* chemin,int code,std::string nomsol,int entite,int type_solution):MG_IDENTIFICATEUR(),nb_champs(nb),femmai(mai),nomsolution(nomsol),typeentite(entite),dim_solution(type_solution)
{
inisolution(chemin,code);    
}

FEM_SOLUTION::FEM_SOLUTION(unsigned long num,FEM_MAILLAGE* mai,int nb,char* chemin,int code,std::string nomsol,int entite,int type_solution):MG_IDENTIFICATEUR(num),nb_champs(nb),femmai(mai),nomsolution(nomsol),typeentite(entite),dim_solution(type_solution)
{
inisolution(chemin,code);
 
}

void FEM_SOLUTION::inisolution(char* chemin,int code)
{
  nom_fichier=new char[strlen(chemin)+2];
  legende=new std::string[nb_champs];
  for (int i=0;i<nb_champs;i++)
  {
      char mess[15];
      sprintf(mess,"Champs_%d",i);
      legende[i]=mess ;
  }
  strcpy(nom_fichier,chemin);
  if (code!=0)
    {
        in=fopen(nom_fichier,"wb");
        unsigned long nbele=0;
        if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_NOEUD) nbele=femmai->get_nb_fem_noeud();
        if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT1) nbele=femmai->get_nb_fem_element1();
        if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT2) nbele=femmai->get_nb_fem_element2();
        if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT3) nbele=femmai->get_nb_fem_element3();
        if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT1_NOEUD) nbele=femmai->get_nb_fem_element1();
        if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT2_NOEUD) nbele=femmai->get_nb_fem_element2();
        if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT3_NOEUD) nbele=femmai->get_nb_fem_element3();
        if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT1_GAUSS) nbele=femmai->get_nb_fem_element1();
        if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT2_GAUSS) nbele=femmai->get_nb_fem_element2();
        if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT3_GAUSS) nbele=femmai->get_nb_fem_element3();
        fwrite(&nbele,sizeof(unsigned long),1,in);
        for (unsigned long i=0;i<nbele+1;i++)
          {
          unsigned long val=0;
          fwrite(&val,sizeof(unsigned long),1,in);
          }
        unsigned long *tab=new unsigned long[nbele+1];  
        unsigned long i=0;
        unsigned long cumul_iter=0;
        if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_NOEUD)
          {
          LISTE_FEM_NOEUD::iterator ite;
          for (FEM_NOEUD* ele=femmai->get_premier_noeud(ite);ele;ele=femmai->get_suivant_noeud(ite))
            {
            tab[i]=cumul_iter;
            unsigned long iter=1*nb_champs*dim_solution;
            for (unsigned long j=0;j<iter;j++)
                {
                double val=0;
                fwrite(&val,sizeof(double),1,in);
                }
            i++;
            cumul_iter=cumul_iter+iter;
            }
                    
          }
        if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT1)
          {
          LISTE_FEM_ELEMENT1::iterator ite;
          for (FEM_ELEMENT1* ele=femmai->get_premier_element1(ite);ele;ele=femmai->get_suivant_element1(ite))
            {
            tab[i]=cumul_iter;
            unsigned long iter=1*nb_champs*dim_solution;
            for (unsigned long j=0;j<iter;j++)
                {
                double val=0;
                fwrite(&val,sizeof(double),1,in);
                }
            i++;
            cumul_iter=cumul_iter+iter;
            }
                    
          }
        if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT2)
          {
          LISTE_FEM_ELEMENT2::iterator ite;
          for (FEM_ELEMENT2* ele=femmai->get_premier_element2(ite);ele;ele=femmai->get_suivant_element2(ite))
            {
            tab[i]=cumul_iter;
            unsigned long iter=1*nb_champs*dim_solution;
            for (unsigned long j=0;j<iter;j++)
                {
                double val=0;
                fwrite(&val,sizeof(double),1,in);
                }
            i++;
            cumul_iter=cumul_iter+iter;
            }
                    
          }
        if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT3)
          {
          LISTE_FEM_ELEMENT3::iterator ite;
          for (FEM_ELEMENT3* ele=femmai->get_premier_element3(ite);ele;ele=femmai->get_suivant_element3(ite))
            {
            tab[i]=cumul_iter;
            unsigned long iter=1*nb_champs*dim_solution;
            for (unsigned long j=0;j<iter;j++)
                {
                double val=0;
                fwrite(&val,sizeof(double),1,in);
                }
            i++;
            cumul_iter=cumul_iter+iter;
            }
                    
          }
        if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT3_NOEUD)
          {
          LISTE_FEM_ELEMENT3::iterator ite;
          for (FEM_ELEMENT3* ele=femmai->get_premier_element3(ite);ele;ele=femmai->get_suivant_element3(ite))
            {
            tab[i]=cumul_iter;
            unsigned long iter=ele->get_nb_fem_noeud()*nb_champs*dim_solution;
            for (unsigned long j=0;j<iter;j++)
                {
                double val=0;
                fwrite(&val,sizeof(double),1,in);
                }
            i++;
            cumul_iter=cumul_iter+iter;
            }
                    
          }
       if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT1_NOEUD)
          {
          LISTE_FEM_ELEMENT1::iterator ite;
          for (FEM_ELEMENT1* ele=femmai->get_premier_element1(ite);ele;ele=femmai->get_suivant_element1(ite))
            {
            tab[i]=cumul_iter;
            unsigned long iter=ele->get_nb_fem_noeud()*nb_champs*dim_solution;
            for (unsigned long j=0;j<iter;j++)
                {
                double val=0;
                fwrite(&val,sizeof(double),1,in);
                }
            i++;
            cumul_iter=cumul_iter+iter;
            }
                    
          }
        if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT2_NOEUD)
          {
          LISTE_FEM_ELEMENT2::iterator ite;
          for (FEM_ELEMENT2* ele=femmai->get_premier_element2(ite);ele;ele=femmai->get_suivant_element2(ite))
            {
            tab[i]=cumul_iter;
            unsigned long iter=ele->get_nb_fem_noeud()*nb_champs*dim_solution;
            for (unsigned long j=0;j<iter;j++)
                {
                double val=0;
                fwrite(&val,sizeof(double),1,in);
                }
            i++;
            cumul_iter=cumul_iter+iter;
            }
                    
          }
       if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT3_GAUSS)
          {
          LISTE_FEM_ELEMENT3::iterator ite;
          for (FEM_ELEMENT3* ele=femmai->get_premier_element3(ite);ele;ele=femmai->get_suivant_element3(ite))
            {
            tab[i]=cumul_iter;
            unsigned long iter=1*nb_champs*dim_solution;
            /*for (unsigned long j=0;j<iter;j++)
                {
                double val=0;
                fwrite(&val,sizeof(double),1,in);
                }*/
            i++;
            cumul_iter=cumul_iter+iter;
            }
                    
          }
       if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT2_GAUSS)
          {
          LISTE_FEM_ELEMENT2::iterator ite;
          for (FEM_ELEMENT2* ele=femmai->get_premier_element2(ite);ele;ele=femmai->get_suivant_element2(ite))
            {
            tab[i]=cumul_iter;
            unsigned long iter=1*nb_champs*dim_solution;
            /*for (unsigned long j=0;j<iter;j++)
                {
                double val=0;
                fwrite(&val,sizeof(double),1,in);
                }*/
            i++;
            cumul_iter=cumul_iter+iter;
            }
                    
          }
       if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT1_GAUSS)
          {
          LISTE_FEM_ELEMENT1::iterator ite;
          for (FEM_ELEMENT1* ele=femmai->get_premier_element1(ite);ele;ele=femmai->get_suivant_element1(ite))
            {
            tab[i]=cumul_iter;
            unsigned long iter=1*nb_champs*dim_solution;
            /*for (unsigned long j=0;j<iter;j++)
                {
                double val=0;
                fwrite(&val,sizeof(double),1,in);
                }*/
            i++;
            cumul_iter=cumul_iter+iter;
            }
                    
          }
        tab[nbele]=cumul_iter;
        fseek(in,sizeof(unsigned long),SEEK_SET);
        for (unsigned long i=0;i<nbele+1;i++)
          {
          unsigned long val=0;
          fwrite(&(tab[i]),sizeof(unsigned long),1,in);
          }
        delete [] tab;
        fclose(in);
    }
    in=fopen(nom_fichier,"r+b");
    hash.clear();
    fseek(in,0,SEEK_SET);
    unsigned long nb;
    fread(&nb,sizeof(unsigned long),1,in);
    hash.insert(hash.end(),nb);
    for (unsigned long i=0;i<nb+1;i++)
      {
      unsigned long val;
      fread(&val,sizeof(unsigned long),1,in);
      hash.insert(hash.end(),val);
      }
}


FEM_SOLUTION::~FEM_SOLUTION()
{
    delete [] nom_fichier;
    delete [] legende;
    if (in!=NULL) fclose(in);
}

void FEM_SOLUTION::efface(void)
{
    fclose(in);
    remove(nom_fichier);
// DeleteFile(nom_fichier);
    in=NULL;
}

int FEM_SOLUTION::get_entite_solution(void)
{
    return typeentite;
}

double FEM_SOLUTION::lire(int i,int j,int coord,int num_no)
{
    unsigned long nb=hash[0];
    unsigned long pos=hash[i+1];
    unsigned long pospro=hash[i+2];
    int nbval=(pospro-pos)/nb_champs/dim_solution;
    long posval=(nb+2)*sizeof(unsigned long)+(pos+j*dim_solution*nbval+num_no*dim_solution+coord)*sizeof(double);
    fseek(in,posval,SEEK_SET);
    double val;
    fread(&val,sizeof(double),1,in);
    return (val);
}

void FEM_SOLUTION::ecrire(double val,int i,int j,int coord,int num_no)
{
    unsigned long nb=hash[0];
    unsigned long pos=hash[i+1];
    unsigned long pospro=hash[i+2];
    int nbval=(pospro-pos)/nb_champs/dim_solution;
    long posval=(nb+2)*sizeof(unsigned long)+(pos+j*dim_solution*nbval+num_no*dim_solution+coord)*sizeof(double);
    fseek(in,posval,SEEK_SET);
    fwrite(&val,sizeof(double),1,in);
    fflush(in);
}

void FEM_SOLUTION::change_nb_gauss(int i,int num)
{
if ((typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT1_GAUSS)||(typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT2_GAUSS)||(typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT3_GAUSS))
  {
  unsigned long nb=hash[0];
  unsigned long pos=hash[i+1];
  unsigned long oldpospro=hash[i+2];
  unsigned long pospro=num*nb_champs*dim_solution+pos;
  hash[i+2]=pospro;
  unsigned long decalage=pospro-oldpospro;
  for (unsigned long j=i+3;j<nb+1;j++)
    hash[j]=hash[j]+decalage;
  fseek(in,0,SEEK_SET);
  for (unsigned long j=0;j<nb+2;j++)
      fwrite(&(hash[j]),sizeof(unsigned long),1,in);
  }
  
}

int FEM_SOLUTION::get_nb_gauss(int i)
{
if ((typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT1_GAUSS)||(typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT2_GAUSS)||(typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT3_GAUSS))
  {
  unsigned long nb=hash[0];
  unsigned long pos=hash[i+1];
  unsigned long pospro=hash[i+2];
  int nbval=(pospro-pos)/nb_champs/dim_solution;
  return nbval;  
  }
return 0;  
}


void FEM_SOLUTION::change_legende(int num,std::string val)
{
    legende[num]=val;
}

std::string FEM_SOLUTION::get_legende(int num)
{
    return legende[num];
}

int FEM_SOLUTION::get_nb_champ(void)
{
    return nb_champs;
}

FEM_MAILLAGE* FEM_SOLUTION::get_maillage(void)
{
    return femmai;
}

double FEM_SOLUTION::get_legende_min(int num)
{
    return solmin[num];
}

double FEM_SOLUTION::get_legende_max(int num)
{
    return solmax[num];
}
int FEM_SOLUTION::get_dim_solution(void)
{
    return dim_solution;
}
void FEM_SOLUTION::active_solution(int num)
{
    for (int j=0;j<9;j++)
      {
      solmax[j]=-1e308;
      solmin[j]=1e308;
      }
    int i=0;
    if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_NOEUD)
    {
        LISTE_FEM_NOEUD::iterator it;
        for (FEM_NOEUD* no=femmai->get_premier_noeud(it);no!=NULL;no=femmai->get_suivant_noeud(it))
        {
            for (int k=0;k<dim_solution;k++)
              {
              double val=lire(i,num,k);
              no->change_solution(val,k);
              if (val<solmin[k]) solmin[k]=val;
              if (val>solmax[k]) solmax[k]=val;
              }
            i++;
        }
    }
    if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT1)
    {
        LISTE_FEM_ELEMENT1::iterator it;
        for (FEM_ELEMENT1* seg=femmai->get_premier_element1(it);seg!=NULL;seg=femmai->get_suivant_element1(it))
        {
            for (int k=0;k<dim_solution;k++)
              {
                double val=lire(i,num,k);
                seg->change_solution(val,k);
                if (val<solmin[k]) solmin[k]=val;
                if (val>solmax[k]) solmax[k]=val;
              }
            i++;
        }
    }
    if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT2)
    {
        LISTE_FEM_ELEMENT2::iterator it;
        for (FEM_ELEMENT2* tri=femmai->get_premier_element2(it);tri!=NULL;tri=femmai->get_suivant_element2(it))
        {
            for (int k=0;k<dim_solution;k++)
              {
              double val=lire(i,num);
              tri->change_solution(val);
              if (val<solmin[k]) solmin[k]=val;
              if (val>solmax[k]) solmax[k]=val;
              }
            i++;
        }
    }
    if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT3)
    {
        LISTE_FEM_ELEMENT3::iterator it;
        for (FEM_ELEMENT3* tet=femmai->get_premier_element3(it);tet!=NULL;tet=femmai->get_suivant_element3(it))
        {
             for (int k=0;k<dim_solution;k++)
              {
              double val=lire(i,num);
              tet->change_solution(val);
              if (val<solmin[k]) solmin[k]=val;
              if (val>solmax[k]) solmax[k]=val;
              }
            i++;
        }
    }
    for (int k=0;k<dim_solution;k++)
      if (OPERATEUR::egal(fabs(solmin[k]-solmax[k]),0.,1e-10)==1)
        solmax[k]=solmin[k]+1e-10;
}


double FEM_SOLUTION::get_moyenne_volumique_champs(MG_VOLUME* mgvol,int numchamps,int coord)
{
if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT3_GAUSS)  
    {
    LISTE_FEM_ELEMENT3::iterator it;
    double unite=femmai->get_mg_geometrie()->get_valeur_unite();
    int i=0;
    double val=0.;
    double vol=0.;
    for (FEM_ELEMENT3* ele=femmai->get_premier_element3(it);ele!=NULL;ele=femmai->get_suivant_element3(it))
      {
      if (mgvol!=NULL)
        if (ele->get_lien_topologie()!=mgvol) {i++;continue;}
      int nbgauss=get_nb_gauss(i);
      double valele=0.0;
      for (int r=0;r<nbgauss;r++)
        {
        int degre=ele->get_degre_gauss(nbgauss);
        double w;
        double uvw[3];
        ele->get_pt_gauss(degre,r,w,uvw);
        double jac[9];
        int l,c;
        double det=ele->get_jacobien(jac,uvw,l,c,unite);
        double valeur=lire(i,numchamps,coord,r)*fabs(det);
        valele=valele+w*valeur;
        vol=vol+w*fabs(det);
        }
      val=val+valele;
      i++;
      }  
    return val/vol;
    }
return 0;
}


double FEM_SOLUTION::get_moyenne_surfacique_champs(MG_FACE* fac,int numchamps,int coord)
{
if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT2_GAUSS)  
    {
    LISTE_FEM_ELEMENT2::iterator it;
    double unite=femmai->get_mg_geometrie()->get_valeur_unite();
    int i=0;
    double val=0.;
    double vol=0.;
    for (FEM_ELEMENT2* ele=femmai->get_premier_element2(it);ele!=NULL;ele=femmai->get_suivant_element2(it))
      {
      if (fac!=NULL)
        if (ele->get_lien_topologie()!=fac) {i++;continue;}
      int nbgauss=get_nb_gauss(i);
      double valele=0.0;
      for (int r=0;r<nbgauss;r++)
        {
        int degre=ele->get_degre_gauss(nbgauss);
        double w;
        double uvw[3];
        ele->get_pt_gauss(degre,r,w,uvw);
        double jac[9];
        int l,c;
        double det=ele->get_jacobien(jac,uvw,l,c,unite);
        double valeur=lire(i,numchamps,coord,r)*fabs(det);
        valele=valele+w*valeur;
        vol=vol+w*fabs(det);
        }
      val=val+valele;
      i++;
      }  
    return val/vol;
    }
return 0;
}


double FEM_SOLUTION::get_moyenne_lineique_champs(MG_ARETE* are,int numchamps,int coord)
{
if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT1_GAUSS)  
    {
    LISTE_FEM_ELEMENT1::iterator it;
    double unite=femmai->get_mg_geometrie()->get_valeur_unite();
    int i=0;
    double val=0.;
    double vol=0.;
    for (FEM_ELEMENT1* ele=femmai->get_premier_element1(it);ele!=NULL;ele=femmai->get_suivant_element1(it))
      {
      int nbgauss=get_nb_gauss(i);
      double valele=0.0;
      for (int r=0;r<nbgauss;r++)
        {
        if (are!=NULL)
          if (ele->get_lien_topologie()!=are) {i++;continue;}
        int degre=ele->get_degre_gauss(nbgauss);
        double w;
        double uvw[3];
        ele->get_pt_gauss(degre,r,w,uvw);
        double jac[9];
        int l,c;
        double det=ele->get_jacobien(jac,uvw,l,c,unite);
        double valeur=lire(i,numchamps,coord,r)*fabs(det);
        valele=valele+w*valeur;
        vol=vol+w*fabs(det);
        }
      val=val+valele;
      i++;
      }  
    return val/vol;
    }
return 0;
}
double FEM_SOLUTION::get_moyenne_volumique_elementaire_champs(int numchamps,FEM_ELEMENT3* ele,int num,int coord)
{
if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT3_GAUSS)  
    {
    double unite=femmai->get_mg_geometrie()->get_valeur_unite();
    double val=0.;
    double vol=0.;
    int nbgauss=get_nb_gauss(num);
    for (int r=0;r<nbgauss;r++)
        {
        int degre=ele->get_degre_gauss(nbgauss);
        double w;
        double uvw[3];
        ele->get_pt_gauss(degre,r,w,uvw);
        double jac[9];
        int l,c;
        double det=ele->get_jacobien(jac,uvw,l,c,unite);
        double valeur=lire(num,numchamps,coord,r)*fabs(det);
        val=val+w*valeur;
        vol=vol+w*fabs(det);
        }
    return val/vol;
    }
return 0;
}

double FEM_SOLUTION::get_moyenne_surfacique_elementaire_champs(int numchamps,FEM_ELEMENT2* ele,int num,int coord)
{
if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT2_GAUSS)  
    {
    double unite=femmai->get_mg_geometrie()->get_valeur_unite();
    double val=0.;
    double vol=0.;
    int nbgauss=get_nb_gauss(num);
    for (int r=0;r<nbgauss;r++)
        {
        int degre=ele->get_degre_gauss(nbgauss);
        double w;
        double uvw[3];
        ele->get_pt_gauss(degre,r,w,uvw);
        double jac[9];
        int l,c;
        double det=ele->get_jacobien(jac,uvw,l,c,unite);
        double valeur=lire(num,numchamps,coord,r)*fabs(det);
        val=val+w*valeur;
        vol=vol+w*fabs(det);
        }
    return val/vol;
    }
return 0;
}

double FEM_SOLUTION::get_moyenne_lineique_elementaire_champs(int numchamps,FEM_ELEMENT1* ele,int num,int coord)
{
if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT1_GAUSS)  
    {
    double unite=femmai->get_mg_geometrie()->get_valeur_unite();
    double val=0.;
    double vol=0.;
    int nbgauss=get_nb_gauss(num);
    for (int r=0;r<nbgauss;r++)
        {
        int degre=ele->get_degre_gauss(nbgauss);
        double w;
        double uvw[3];
        ele->get_pt_gauss(degre,r,w,uvw);
        double jac[9];
        int l,c;
        double det=ele->get_jacobien(jac,uvw,l,c,unite);
        double valeur=lire(num,numchamps,coord,r)*fabs(det);
        val=val+w*valeur;
        vol=vol+w*fabs(det);
        }
    return val/vol;
    }
return 0;
}


std::string FEM_SOLUTION::get_nom(void)
{
    return nomsolution;
}

std::string FEM_SOLUTION::get_nom_fichier(void)
{
    return nom_fichier;
}
void FEM_SOLUTION::enregistrer(std::ostream& o)
{
    o << "%" << get_id() << "=FEM_SOLUTION(" << typeentite << "," << dim_solution   << "," << nomsolution << ",$" << femmai->get_id() << "," << nb_champs << "," << nom_fichier << ",(";
    for (int i=0;i<nb_champs;i++)
    {
        o << legende[i];
        if (i!=nb_champs-1) o<<",";
        else o<<")";
    }


    o << ");" << std::endl;
}


std::ostream& operator << (std::ostream& o,FEM_SOLUTION& sol)
{
    sol.enregistrer(o);
    return o;
}


Revision:	762
Committed:	Mon Nov 23 15:31:53 2015 UTC (9 years, 5 months ago) by francois
File size:	19212 byte(s)
Log Message:	calcul de moyenne volumique integrale dans fem_solution a partir des solutions aux points de gauss+ bug orthographe dans les conditions aux limites
#	User	Rev	Content
1	francois	283	//------------------------------------------------------------
2			//------------------------------------------------------------
3			// MAGiC
4			// Jean Christophe Cuilli�re et Vincent FRANCOIS
5			// D�partement de G�nie M�canique - UQTR
6			//------------------------------------------------------------
7			// Le projet MAGIC est un projet de recherche du d�partement
8			// de g�nie m�canique de l'Universit� du Qu�bec �
9			// Trois Rivi�res
10			// Les librairies ne peuvent �tre utilis�es sans l'accord
11			// des auteurs (contact : francois@uqtr.ca)
12			//------------------------------------------------------------
13			//------------------------------------------------------------
14			//
15			// fem_solution.cpp
16			//
17			//------------------------------------------------------------
18			//------------------------------------------------------------
19			// COPYRIGHT 2000
20			// Version du 02/03/2006 � 11H22
21			//------------------------------------------------------------
22			//------------------------------------------------------------
23
24
25			#include "gestionversion.h"
26			#include <math.h>
27			#include <iostream>
28			#include <stdlib.h>
29			#include <string.h>
30
31			#include "fem_solution.h"
32			#include "fem_maillage.h"
33	francois	762	#include "mg_geometrie.h"
34	francois	283	#include "ot_mathematique.h"
35
36
37			#include <stdio.h> // pour la fonction remove() qui remplace DeleteFile
38
39
40
41
42
43
44	francois	375	FEM_SOLUTION::FEM_SOLUTION(FEM_MAILLAGE* mai,int nb,char* chemin,int code,std::string nomsol,int entite,int type_solution):MG_IDENTIFICATEUR(),nb_champs(nb),femmai(mai),nomsolution(nomsol),typeentite(entite),dim_solution(type_solution)
45	francois	283	{
46	francois	375	inisolution(chemin,code);
47	francois	283	}
48
49	francois	375	FEM_SOLUTION::FEM_SOLUTION(unsigned long num,FEM_MAILLAGE* mai,int nb,char* chemin,int code,std::string nomsol,int entite,int type_solution):MG_IDENTIFICATEUR(num),nb_champs(nb),femmai(mai),nomsolution(nomsol),typeentite(entite),dim_solution(type_solution)
50	francois	283	{
51	francois	375	inisolution(chemin,code);
52
53			}
54
55			void FEM_SOLUTION::inisolution(char* chemin,int code)
56			{
57			nom_fichier=new char[strlen(chemin)+2];
58			legende=new std::string[nb_champs];
59			for (int i=0;i<nb_champs;i++)
60			{
61			char mess[15];
62			sprintf(mess,"Champs_%d",i);
63			legende[i]=mess ;
64			}
65			strcpy(nom_fichier,chemin);
66			if (code!=0)
67	francois	283	{
68			in=fopen(nom_fichier,"wb");
69	francois	375	unsigned long nbele=0;
70			if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_NOEUD) nbele=femmai->get_nb_fem_noeud();
71			if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT1) nbele=femmai->get_nb_fem_element1();
72			if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT2) nbele=femmai->get_nb_fem_element2();
73			if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT3) nbele=femmai->get_nb_fem_element3();
74			if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT1_NOEUD) nbele=femmai->get_nb_fem_element1();
75			if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT2_NOEUD) nbele=femmai->get_nb_fem_element2();
76			if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT3_NOEUD) nbele=femmai->get_nb_fem_element3();
77	francois	757	if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT1_GAUSS) nbele=femmai->get_nb_fem_element1();
78			if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT2_GAUSS) nbele=femmai->get_nb_fem_element2();
79			if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT3_GAUSS) nbele=femmai->get_nb_fem_element3();
80	francois	375	fwrite(&nbele,sizeof(unsigned long),1,in);
81			for (unsigned long i=0;i<nbele+1;i++)
82			{
83			unsigned long val=0;
84			fwrite(&val,sizeof(unsigned long),1,in);
85			}
86			unsigned long *tab=new unsigned long[nbele+1];
87			unsigned long i=0;
88			unsigned long cumul_iter=0;
89			if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_NOEUD)
90			{
91			LISTE_FEM_NOEUD::iterator ite;
92			for (FEM_NOEUD* ele=femmai->get_premier_noeud(ite);ele;ele=femmai->get_suivant_noeud(ite))
93			{
94			tab[i]=cumul_iter;
95			unsigned long iter=1nb_champsdim_solution;
96			for (unsigned long j=0;j<iter;j++)
97			{
98			double val=0;
99			fwrite(&val,sizeof(double),1,in);
100			}
101			i++;
102			cumul_iter=cumul_iter+iter;
103			}
104
105			}
106			if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT1)
107			{
108			LISTE_FEM_ELEMENT1::iterator ite;
109			for (FEM_ELEMENT1* ele=femmai->get_premier_element1(ite);ele;ele=femmai->get_suivant_element1(ite))
110			{
111			tab[i]=cumul_iter;
112			unsigned long iter=1nb_champsdim_solution;
113			for (unsigned long j=0;j<iter;j++)
114			{
115			double val=0;
116			fwrite(&val,sizeof(double),1,in);
117			}
118			i++;
119			cumul_iter=cumul_iter+iter;
120			}
121
122			}
123			if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT2)
124			{
125			LISTE_FEM_ELEMENT2::iterator ite;
126			for (FEM_ELEMENT2* ele=femmai->get_premier_element2(ite);ele;ele=femmai->get_suivant_element2(ite))
127			{
128			tab[i]=cumul_iter;
129			unsigned long iter=1nb_champsdim_solution;
130			for (unsigned long j=0;j<iter;j++)
131			{
132			double val=0;
133			fwrite(&val,sizeof(double),1,in);
134			}
135			i++;
136			cumul_iter=cumul_iter+iter;
137			}
138
139			}
140			if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT3)
141			{
142			LISTE_FEM_ELEMENT3::iterator ite;
143			for (FEM_ELEMENT3* ele=femmai->get_premier_element3(ite);ele;ele=femmai->get_suivant_element3(ite))
144			{
145			tab[i]=cumul_iter;
146			unsigned long iter=1nb_champsdim_solution;
147			for (unsigned long j=0;j<iter;j++)
148			{
149			double val=0;
150			fwrite(&val,sizeof(double),1,in);
151			}
152			i++;
153			cumul_iter=cumul_iter+iter;
154			}
155
156			}
157			if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT3_NOEUD)
158			{
159			LISTE_FEM_ELEMENT3::iterator ite;
160			for (FEM_ELEMENT3* ele=femmai->get_premier_element3(ite);ele;ele=femmai->get_suivant_element3(ite))
161			{
162			tab[i]=cumul_iter;
163			unsigned long iter=ele->get_nb_fem_noeud()nb_champsdim_solution;
164			for (unsigned long j=0;j<iter;j++)
165			{
166			double val=0;
167			fwrite(&val,sizeof(double),1,in);
168			}
169			i++;
170			cumul_iter=cumul_iter+iter;
171			}
172
173			}
174			if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT1_NOEUD)
175			{
176			LISTE_FEM_ELEMENT1::iterator ite;
177			for (FEM_ELEMENT1* ele=femmai->get_premier_element1(ite);ele;ele=femmai->get_suivant_element1(ite))
178			{
179			tab[i]=cumul_iter;
180			unsigned long iter=ele->get_nb_fem_noeud()nb_champsdim_solution;
181			for (unsigned long j=0;j<iter;j++)
182			{
183			double val=0;
184			fwrite(&val,sizeof(double),1,in);
185			}
186			i++;
187			cumul_iter=cumul_iter+iter;
188			}
189
190			}
191			if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT2_NOEUD)
192			{
193			LISTE_FEM_ELEMENT2::iterator ite;
194			for (FEM_ELEMENT2* ele=femmai->get_premier_element2(ite);ele;ele=femmai->get_suivant_element2(ite))
195			{
196			tab[i]=cumul_iter;
197			unsigned long iter=ele->get_nb_fem_noeud()nb_champsdim_solution;
198			for (unsigned long j=0;j<iter;j++)
199			{
200			double val=0;
201			fwrite(&val,sizeof(double),1,in);
202			}
203			i++;
204			cumul_iter=cumul_iter+iter;
205			}
206
207			}
208	francois	757	if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT3_GAUSS)
209			{
210			LISTE_FEM_ELEMENT3::iterator ite;
211			for (FEM_ELEMENT3* ele=femmai->get_premier_element3(ite);ele;ele=femmai->get_suivant_element3(ite))
212			{
213			tab[i]=cumul_iter;
214			unsigned long iter=1nb_champsdim_solution;
215			/*for (unsigned long j=0;j<iter;j++)
216			{
217			double val=0;
218			fwrite(&val,sizeof(double),1,in);
219			}*/
220			i++;
221			cumul_iter=cumul_iter+iter;
222			}
223
224			}
225			if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT2_GAUSS)
226			{
227			LISTE_FEM_ELEMENT2::iterator ite;
228			for (FEM_ELEMENT2* ele=femmai->get_premier_element2(ite);ele;ele=femmai->get_suivant_element2(ite))
229			{
230			tab[i]=cumul_iter;
231			unsigned long iter=1nb_champsdim_solution;
232			/*for (unsigned long j=0;j<iter;j++)
233			{
234			double val=0;
235			fwrite(&val,sizeof(double),1,in);
236			}*/
237			i++;
238			cumul_iter=cumul_iter+iter;
239			}
240
241			}
242			if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT1_GAUSS)
243			{
244			LISTE_FEM_ELEMENT1::iterator ite;
245			for (FEM_ELEMENT1* ele=femmai->get_premier_element1(ite);ele;ele=femmai->get_suivant_element1(ite))
246			{
247			tab[i]=cumul_iter;
248			unsigned long iter=1nb_champsdim_solution;
249			/*for (unsigned long j=0;j<iter;j++)
250			{
251			double val=0;
252			fwrite(&val,sizeof(double),1,in);
253			}*/
254			i++;
255			cumul_iter=cumul_iter+iter;
256			}
257
258			}
259	francois	375	tab[nbele]=cumul_iter;
260			fseek(in,sizeof(unsigned long),SEEK_SET);
261			for (unsigned long i=0;i<nbele+1;i++)
262			{
263			unsigned long val=0;
264			fwrite(&(tab[i]),sizeof(unsigned long),1,in);
265			}
266			delete [] tab;
267	francois	283	fclose(in);
268			}
269			in=fopen(nom_fichier,"r+b");
270	francois	375	hash.clear();
271			fseek(in,0,SEEK_SET);
272			unsigned long nb;
273			fread(&nb,sizeof(unsigned long),1,in);
274			hash.insert(hash.end(),nb);
275			for (unsigned long i=0;i<nb+1;i++)
276			{
277			unsigned long val;
278			fread(&val,sizeof(unsigned long),1,in);
279			hash.insert(hash.end(),val);
280			}
281	francois	283	}
282
283	francois	375
284
285	francois	283	FEM_SOLUTION::~FEM_SOLUTION()
286			{
287			delete [] nom_fichier;
288			delete [] legende;
289			if (in!=NULL) fclose(in);
290			}
291
292			void FEM_SOLUTION::efface(void)
293			{
294			fclose(in);
295			remove(nom_fichier);
296			// DeleteFile(nom_fichier);
297			in=NULL;
298			}
299
300	francois	635	int FEM_SOLUTION::get_entite_solution(void)
301	francois	283	{
302			return typeentite;
303			}
304
305	francois	375	double FEM_SOLUTION::lire(int i,int j,int coord,int num_no)
306	francois	283	{
307	francois	375	unsigned long nb=hash[0];
308			unsigned long pos=hash[i+1];
309			unsigned long pospro=hash[i+2];
310			int nbval=(pospro-pos)/nb_champs/dim_solution;
311			long posval=(nb+2)sizeof(unsigned long)+(pos+jdim_solutionnbval+num_nodim_solution+coord)*sizeof(double);
312			fseek(in,posval,SEEK_SET);
313	francois	283	double val;
314	francois	375	fread(&val,sizeof(double),1,in);
315	francois	283	return (val);
316			}
317
318	francois	375	void FEM_SOLUTION::ecrire(double val,int i,int j,int coord,int num_no)
319	francois	283	{
320	francois	375	unsigned long nb=hash[0];
321			unsigned long pos=hash[i+1];
322			unsigned long pospro=hash[i+2];
323			int nbval=(pospro-pos)/nb_champs/dim_solution;
324			long posval=(nb+2)sizeof(unsigned long)+(pos+jdim_solutionnbval+num_nodim_solution+coord)*sizeof(double);
325			fseek(in,posval,SEEK_SET);
326	francois	283	fwrite(&val,sizeof(double),1,in);
327			fflush(in);
328			}
329
330	francois	757	void FEM_SOLUTION::change_nb_gauss(int i,int num)
331			{
332			if ((typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT1_GAUSS)\|\|(typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT2_GAUSS)\|\|(typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT3_GAUSS))
333			{
334			unsigned long nb=hash[0];
335			unsigned long pos=hash[i+1];
336			unsigned long oldpospro=hash[i+2];
337			unsigned long pospro=numnb_champsdim_solution+pos;
338			hash[i+2]=pospro;
339			unsigned long decalage=pospro-oldpospro;
340			for (unsigned long j=i+3;j<nb+1;j++)
341			hash[j]=hash[j]+decalage;
342			fseek(in,0,SEEK_SET);
343			for (unsigned long j=0;j<nb+2;j++)
344			fwrite(&(hash[j]),sizeof(unsigned long),1,in);
345			}
346
347			}
348	francois	283
349	francois	757	int FEM_SOLUTION::get_nb_gauss(int i)
350			{
351			if ((typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT1_GAUSS)\|\|(typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT2_GAUSS)\|\|(typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT3_GAUSS))
352			{
353			unsigned long nb=hash[0];
354			unsigned long pos=hash[i+1];
355			unsigned long pospro=hash[i+2];
356			int nbval=(pospro-pos)/nb_champs/dim_solution;
357			return nbval;
358			}
359			return 0;
360			}
361
362
363	francois	283	void FEM_SOLUTION::change_legende(int num,std::string val)
364			{
365			legende[num]=val;
366			}
367
368			std::string FEM_SOLUTION::get_legende(int num)
369			{
370			return legende[num];
371			}
372
373			int FEM_SOLUTION::get_nb_champ(void)
374			{
375			return nb_champs;
376			}
377
378			FEM_MAILLAGE* FEM_SOLUTION::get_maillage(void)
379			{
380			return femmai;
381			}
382
383	francois	375	double FEM_SOLUTION::get_legende_min(int num)
384	francois	283	{
385	francois	375	return solmin[num];
386	francois	283	}
387
388	francois	375	double FEM_SOLUTION::get_legende_max(int num)
389	francois	283	{
390	francois	375	return solmax[num];
391	francois	283	}
392	francois	375	int FEM_SOLUTION::get_dim_solution(void)
393			{
394			return dim_solution;
395			}
396	francois	283	void FEM_SOLUTION::active_solution(int num)
397			{
398	francois	375	for (int j=0;j<9;j++)
399			{
400			solmax[j]=-1e308;
401			solmin[j]=1e308;
402			}
403	francois	283	int i=0;
404	francois	375	if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_NOEUD)
405	francois	283	{
406			LISTE_FEM_NOEUD::iterator it;
407			for (FEM_NOEUD* no=femmai->get_premier_noeud(it);no!=NULL;no=femmai->get_suivant_noeud(it))
408			{
409	francois	375	for (int k=0;k<dim_solution;k++)
410			{
411			double val=lire(i,num,k);
412			no->change_solution(val,k);
413			if (val<solmin[k]) solmin[k]=val;
414			if (val>solmax[k]) solmax[k]=val;
415			}
416	francois	283	i++;
417			}
418			}
419	francois	375	if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT1)
420	francois	283	{
421	francois	309	LISTE_FEM_ELEMENT1::iterator it;
422			for (FEM_ELEMENT1* seg=femmai->get_premier_element1(it);seg!=NULL;seg=femmai->get_suivant_element1(it))
423	francois	283	{
424	francois	375	for (int k=0;k<dim_solution;k++)
425			{
426			double val=lire(i,num,k);
427			seg->change_solution(val,k);
428			if (val<solmin[k]) solmin[k]=val;
429			if (val>solmax[k]) solmax[k]=val;
430			}
431	francois	283	i++;
432			}
433			}
434	francois	375	if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT2)
435	francois	283	{
436	francois	309	LISTE_FEM_ELEMENT2::iterator it;
437			for (FEM_ELEMENT2* tri=femmai->get_premier_element2(it);tri!=NULL;tri=femmai->get_suivant_element2(it))
438	francois	283	{
439	francois	375	for (int k=0;k<dim_solution;k++)
440			{
441			double val=lire(i,num);
442			tri->change_solution(val);
443			if (val<solmin[k]) solmin[k]=val;
444			if (val>solmax[k]) solmax[k]=val;
445			}
446	francois	283	i++;
447			}
448			}
449	francois	375	if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT3)
450	francois	283	{
451	francois	309	LISTE_FEM_ELEMENT3::iterator it;
452			for (FEM_ELEMENT3* tet=femmai->get_premier_element3(it);tet!=NULL;tet=femmai->get_suivant_element3(it))
453	francois	283	{
454	francois	375	for (int k=0;k<dim_solution;k++)
455			{
456			double val=lire(i,num);
457			tet->change_solution(val);
458			if (val<solmin[k]) solmin[k]=val;
459			if (val>solmax[k]) solmax[k]=val;
460			}
461	francois	283	i++;
462			}
463			}
464	francois	375	for (int k=0;k<dim_solution;k++)
465			if (OPERATEUR::egal(fabs(solmin[k]-solmax[k]),0.,1e-10)==1)
466			solmax[k]=solmin[k]+1e-10;
467	francois	283	}
468
469	francois	762
470			double FEM_SOLUTION::get_moyenne_volumique_champs(MG_VOLUME* mgvol,int numchamps,int coord)
471			{
472			if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT3_GAUSS)
473			{
474			LISTE_FEM_ELEMENT3::iterator it;
475			double unite=femmai->get_mg_geometrie()->get_valeur_unite();
476			int i=0;
477			double val=0.;
478			double vol=0.;
479			for (FEM_ELEMENT3* ele=femmai->get_premier_element3(it);ele!=NULL;ele=femmai->get_suivant_element3(it))
480			{
481			if (mgvol!=NULL)
482			if (ele->get_lien_topologie()!=mgvol) {i++;continue;}
483			int nbgauss=get_nb_gauss(i);
484			double valele=0.0;
485			for (int r=0;r<nbgauss;r++)
486			{
487			int degre=ele->get_degre_gauss(nbgauss);
488			double w;
489			double uvw[3];
490			ele->get_pt_gauss(degre,r,w,uvw);
491			double jac[9];
492			int l,c;
493			double det=ele->get_jacobien(jac,uvw,l,c,unite);
494			double valeur=lire(i,numchamps,coord,r)*fabs(det);
495			valele=valele+w*valeur;
496			vol=vol+w*fabs(det);
497			}
498			val=val+valele;
499			i++;
500			}
501			return val/vol;
502			}
503			return 0;
504			}
505
506
507			double FEM_SOLUTION::get_moyenne_surfacique_champs(MG_FACE* fac,int numchamps,int coord)
508			{
509			if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT2_GAUSS)
510			{
511			LISTE_FEM_ELEMENT2::iterator it;
512			double unite=femmai->get_mg_geometrie()->get_valeur_unite();
513			int i=0;
514			double val=0.;
515			double vol=0.;
516			for (FEM_ELEMENT2* ele=femmai->get_premier_element2(it);ele!=NULL;ele=femmai->get_suivant_element2(it))
517			{
518			if (fac!=NULL)
519			if (ele->get_lien_topologie()!=fac) {i++;continue;}
520			int nbgauss=get_nb_gauss(i);
521			double valele=0.0;
522			for (int r=0;r<nbgauss;r++)
523			{
524			int degre=ele->get_degre_gauss(nbgauss);
525			double w;
526			double uvw[3];
527			ele->get_pt_gauss(degre,r,w,uvw);
528			double jac[9];
529			int l,c;
530			double det=ele->get_jacobien(jac,uvw,l,c,unite);
531			double valeur=lire(i,numchamps,coord,r)*fabs(det);
532			valele=valele+w*valeur;
533			vol=vol+w*fabs(det);
534			}
535			val=val+valele;
536			i++;
537			}
538			return val/vol;
539			}
540			return 0;
541			}
542
543
544			double FEM_SOLUTION::get_moyenne_lineique_champs(MG_ARETE* are,int numchamps,int coord)
545			{
546			if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT1_GAUSS)
547			{
548			LISTE_FEM_ELEMENT1::iterator it;
549			double unite=femmai->get_mg_geometrie()->get_valeur_unite();
550			int i=0;
551			double val=0.;
552			double vol=0.;
553			for (FEM_ELEMENT1* ele=femmai->get_premier_element1(it);ele!=NULL;ele=femmai->get_suivant_element1(it))
554			{
555			int nbgauss=get_nb_gauss(i);
556			double valele=0.0;
557			for (int r=0;r<nbgauss;r++)
558			{
559			if (are!=NULL)
560			if (ele->get_lien_topologie()!=are) {i++;continue;}
561			int degre=ele->get_degre_gauss(nbgauss);
562			double w;
563			double uvw[3];
564			ele->get_pt_gauss(degre,r,w,uvw);
565			double jac[9];
566			int l,c;
567			double det=ele->get_jacobien(jac,uvw,l,c,unite);
568			double valeur=lire(i,numchamps,coord,r)*fabs(det);
569			valele=valele+w*valeur;
570			vol=vol+w*fabs(det);
571			}
572			val=val+valele;
573			i++;
574			}
575			return val/vol;
576			}
577			return 0;
578			}
579			double FEM_SOLUTION::get_moyenne_volumique_elementaire_champs(int numchamps,FEM_ELEMENT3* ele,int num,int coord)
580			{
581			if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT3_GAUSS)
582			{
583			double unite=femmai->get_mg_geometrie()->get_valeur_unite();
584			double val=0.;
585			double vol=0.;
586			int nbgauss=get_nb_gauss(num);
587			for (int r=0;r<nbgauss;r++)
588			{
589			int degre=ele->get_degre_gauss(nbgauss);
590			double w;
591			double uvw[3];
592			ele->get_pt_gauss(degre,r,w,uvw);
593			double jac[9];
594			int l,c;
595			double det=ele->get_jacobien(jac,uvw,l,c,unite);
596			double valeur=lire(num,numchamps,coord,r)*fabs(det);
597			val=val+w*valeur;
598			vol=vol+w*fabs(det);
599			}
600			return val/vol;
601			}
602			return 0;
603			}
604
605			double FEM_SOLUTION::get_moyenne_surfacique_elementaire_champs(int numchamps,FEM_ELEMENT2* ele,int num,int coord)
606			{
607			if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT2_GAUSS)
608			{
609			double unite=femmai->get_mg_geometrie()->get_valeur_unite();
610			double val=0.;
611			double vol=0.;
612			int nbgauss=get_nb_gauss(num);
613			for (int r=0;r<nbgauss;r++)
614			{
615			int degre=ele->get_degre_gauss(nbgauss);
616			double w;
617			double uvw[3];
618			ele->get_pt_gauss(degre,r,w,uvw);
619			double jac[9];
620			int l,c;
621			double det=ele->get_jacobien(jac,uvw,l,c,unite);
622			double valeur=lire(num,numchamps,coord,r)*fabs(det);
623			val=val+w*valeur;
624			vol=vol+w*fabs(det);
625			}
626			return val/vol;
627			}
628			return 0;
629			}
630
631			double FEM_SOLUTION::get_moyenne_lineique_elementaire_champs(int numchamps,FEM_ELEMENT1* ele,int num,int coord)
632			{
633			if (typeentite==MAGIC::ENTITE_SOLUTION::ENTITE_ELEMENT1_GAUSS)
634			{
635			double unite=femmai->get_mg_geometrie()->get_valeur_unite();
636			double val=0.;
637			double vol=0.;
638			int nbgauss=get_nb_gauss(num);
639			for (int r=0;r<nbgauss;r++)
640			{
641			int degre=ele->get_degre_gauss(nbgauss);
642			double w;
643			double uvw[3];
644			ele->get_pt_gauss(degre,r,w,uvw);
645			double jac[9];
646			int l,c;
647			double det=ele->get_jacobien(jac,uvw,l,c,unite);
648			double valeur=lire(num,numchamps,coord,r)*fabs(det);
649			val=val+w*valeur;
650			vol=vol+w*fabs(det);
651			}
652			return val/vol;
653			}
654			return 0;
655			}
656
657
658
659	francois	283	std::string FEM_SOLUTION::get_nom(void)
660			{
661			return nomsolution;
662			}
663
664	francois	326	std::string FEM_SOLUTION::get_nom_fichier(void)
665			{
666			return nom_fichier;
667			}
668	francois	283	void FEM_SOLUTION::enregistrer(std::ostream& o)
669			{
670	francois	375	o << "%" << get_id() << "=FEM_SOLUTION(" << typeentite << "," << dim_solution << "," << nomsolution << ",$" << femmai->get_id() << "," << nb_champs << "," << nom_fichier << ",(";
671	francois	283	for (int i=0;i<nb_champs;i++)
672			{
673			o << legende[i];
674			if (i!=nb_champs-1) o<<",";
675			else o<<")";
676			}
677
678
679			o << ");" << std::endl;
680			}
681
682
683
684
685			std::ostream& operator << (std::ostream& o,FEM_SOLUTION& sol)
686			{
687			sol.enregistrer(o);
688			return o;
689			}
690
691