fct_taille_fem_solution_generateur_fichier.cpp

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/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       fct_taille_fem_solution_generateur_fichier.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:52 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
#include "gestionversion.h"
#include "fct_taille_fem_solution_generateur_fichier.h"
#include <string.h>
#include <math.h>
#include "fct_taille_fem_solution_generateur_echantillon_frontiere.h"
 
 
FCT_TAILLE_FEM_SOLUTION_GENERATEUR_FICHIER::FCT_TAILLE_FEM_SOLUTION_GENERATEUR_FICHIER(MG_GESTIONNAIRE *gesttmp,MG_GEOMETRIE* geotmp,double dg,int fechantillon,int nxtmp,
 int nytmp,int nztmp,char *nomsolution,char* nomfich,double ind1,double ind2):
 FCT_TAILLE_FEM_SOLUTION_ECHANTILLONNAGE_ESTIMATION(gesttmp, geotmp, dg,fechantillon,nxtmp,
nytmp, nztmp, nomsolution, ind1, ind2),nomficher(nomfich)
{
 
}
FCT_TAILLE_FEM_SOLUTION_GENERATEUR_FICHIER::FCT_TAILLE_FEM_SOLUTION_GENERATEUR_FICHIER(FCT_TAILLE_FEM_SOLUTION_GENERATEUR_FICHIER& mdd):
FCT_TAILLE_FEM_SOLUTION_ECHANTILLONNAGE_ESTIMATION(mdd)
{
  
}    
 
FCT_TAILLE_FEM_SOLUTION_GENERATEUR_FICHIER::~FCT_TAILLE_FEM_SOLUTION_GENERATEUR_FICHIER()
{
TPL_MAP_ENTITE<CT_POINT_ECHANTILLONNAGE* >::ITERATEUR it0; 
for(CT_POINT_ECHANTILLONNAGE* PT=lstpoint.get_premier(it0);PT!=NULL;PT=lstpoint.get_suivant(it0))
 delete PT;
}
  
 
 
void FCT_TAILLE_FEM_SOLUTION_GENERATEUR_FICHIER::echantillonnage(void )
{
  double dng;
  double xmin,ymin,zmin,xmax,ymax,zmax;
   FILE *in =fopen(nomficher,"rt");
    if (in!=NULL)
    { 
      char message[500];
        char* res=fgets(message,255,in);
        sscanf(message,"%lf %lf %lf %lf %lf %lf ",&xmin,&ymin,&zmin,&xmax,&ymax,&zmax);
        res=fgets(message,255,in);
        sscanf(message,"%lf ",&dng);
        ecart_nodal=dng;
        while (!feof(in))
        {
          double x,y,z,dn;
          res=fgets(message,255,in);
      double lambda=zonezoom;
      int n=sscanf(message,"%lf %lf %lf %lf %lf",&x,&y,&z,&dn,&lambda); 
          if (n>=4)
          {
          double dist=dng*(1.-dn)*lambda;
          if (!(OPERATEUR::egal(dist,0.0,1e-8)))
          {
            BOITE_3D boite(x-dist,y-dist,z-dist,x+dist,y+dist,z+dist);
            double ecart=dn*dng;
            CT_POINT_ECHANTILLONNAGE* pt=new CT_POINT_ECHANTILLONNAGE(x,y,z,ecart,boite);
            lstpoint.ajouter(pt);
 
 
              }
          }
          
        }
              
   fclose(in);
        TPL_MAP_ENTITE<CT_POINT_ECHANTILLONNAGE* >::ITERATEUR it;
        BOITE_3D boi(xmin,ymin,zmin,xmax,ymax,zmax); 
        boi.change_grosseur(facteur_grossissement*1.01);
        xmin=boi.get_xmin();
        ymin=boi.get_ymin();
        zmin=boi.get_zmin();
        xmax=boi.get_xmax();
        ymax=boi.get_ymax();
        zmax=boi.get_zmax();
        TPL_MAP_ENTITE<CT_POINT_ECHANTILLONNAGE* > lstpoint2;
        for(CT_POINT_ECHANTILLONNAGE* PT=lstpoint.get_premier(it);PT!=NULL;PT=lstpoint.get_suivant(it))
              {
                 double xyz[3];
                 PT->evaluer(xyz);
                 double ecart=PT->get_valeur();
                 BOITE_3D boite(xyz[0],xyz[1],xyz[2],xyz[0],xyz[1],xyz[2]); 
                 CT_POINT_ECHANTILLONNAGE* pt=new CT_POINT_ECHANTILLONNAGE(xyz[0],xyz[1],xyz[2],ecart,boite); 
                 lstpoint2.ajouter(pt);
              }
        octre.initialiser(&lstpoint2,1,xmin,ymin,zmin,xmax, ymax, zmax); 
        for(CT_POINT_ECHANTILLONNAGE* PT=lstpoint2.get_premier(it);PT!=NULL;PT=lstpoint2.get_suivant(it))
          delete PT;
        int nb_point=lstpoint.get_nb();
        for(CT_POINT_ECHANTILLONNAGE* PT=lstpoint.get_premier(it);PT!=NULL;PT=lstpoint.get_suivant(it))
        {
          octre.inserer(PT);
        }
    }
 
   }