fct_generateur_constante.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_generateur_constante.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:52 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
 
 
#include "gestionversion.h"
 
#pragma hdrstop
 
#include "fct_generateur_constante.h"
#include "mg_gestionnaire.h"
#include "ct_face.h"
#include <math.h>
 
FCT_GENERATEUR_CONSTANTE::FCT_GENERATEUR_CONSTANTE(MG_GESTIONNAIRE& mggest,double dn):FCT_GENERATEUR_3D<4>(),gest(mggest),dng(dn),pas_echantillon(20),pas_estimation(20)
{
}
 
FCT_GENERATEUR_CONSTANTE::~FCT_GENERATEUR_CONSTANTE()
{
}
 
 
 
 
void FCT_GENERATEUR_CONSTANTE::construit(int pas_echan,int pas_estim,MG_GROUPE_TOPOLOGIQUE* mggt)
{
    TPL_MAP_ENTITE<MG_ELEMENT_TOPOLOGIQUE*> lst;
    if (mggt!=NULL)
    {
        std::map<MG_ELEMENT_TOPOLOGIQUE*,MG_ELEMENT_TOPOLOGIQUE*>::iterator it;
        for(MG_ELEMENT_TOPOLOGIQUE* ele=mggt->get_premier(it);ele!=NULL;ele=mggt->get_suivant(it))
        {
          lst.ajouter(ele);
          ele->get_topologie_sousjacente(&lst);
        }
    }
    MG_ELEMENT_TOPOLOGIQUE::change_pas_echantillon(pas_echan);
    pas_estimation=pas_estim;
    int nb_geo=gest.get_nb_mg_geometrie();
    BOITE_3D boite;
    bool affect=false;
    for (int i=0;i<nb_geo;i++)
        {
        MG_GEOMETRIE* geo=gest.get_mg_geometrie(i);
        int nbvol=geo->get_nb_mg_volume();
        for (int j=0;j<nbvol;j++)
            {
            MG_VOLUME* vol=geo->get_mg_volume(j);
            if (mggt!=NULL) 
              if (!mggt->est_dans_laliste(vol)) continue;
            BOITE_3D tmp=vol->get_boite_3D();
            if (affect==false) boite=tmp;
            else boite=boite+tmp;
            affect=true;
            }
        int nbcoque=geo->get_nb_mg_coque();
        for (int j=0;j<nbcoque;j++)
            {
            MG_COQUE* coq=geo->get_mg_coque(j);
            if (mggt!=NULL) 
              if (!mggt->est_dans_laliste(coq)) continue;
            BOITE_3D tmp=coq->get_boite_3D();
            if (affect==false) boite=tmp;
            else boite=boite+tmp;
            affect=true;
            }
        }
    MG_ELEMENT_TOPOLOGIQUE::restaure_pas_echantillon();
    boite.change_grosseur(1.25);
    grille.initialiser(boite.get_xmin(),boite.get_ymin(),boite.get_zmin(),boite.get_xmax(),boite.get_ymax(),boite.get_zmax(),pas_estimation,pas_estimation,pas_estimation);
    double valmax=-1e308;
    for (int i=0;i<pas_estimation;i++)
        for (int j=0;j<pas_estimation;j++)
            for (int k=0;k<pas_estimation;k++)
            {
                TPL_CELLULE_GRILLE_INFO<CT_POINT<1> *,DOUBLEN<4> > *cell=grille.get_cellule(i,j,k);
                double tab[4];
                tab[0]=dng;
                tab[1]=0.;
                tab[2]=0.;
                tab[3]=0.;
                DOUBLEN<4> ecart0(tab);
                if (tab[0]>valmax)
                {
                    valeurmax=ecart0;
                    valmax=tab[0];
                }
                cell->change_info(0,ecart0);
                cell->change_info(1,ecart0);
                cell->change_info(2,ecart0);
                cell->change_info(3,ecart0);
                cell->change_info(4,ecart0);
                cell->change_info(5,ecart0);
                cell->change_info(6,ecart0);
                cell->change_info(7,ecart0);
            }
    legende[0]="F(x,y,z)";
    legende[1]="dF/dx";
    legende[2]="dF/dy";
    legende[3]="dF/dz";
}
 
void FCT_GENERATEUR_CONSTANTE::construit(double xmin,double ymin,double zmin,double xmax,double ymax,double zmax,int pas_echan,int pas_estim)
{
    pas_echantillon=pas_echan;
    pas_estimation=pas_estim;
    BOITE_3D boite(xmin,ymin,zmin,xmax,ymax,zmax);
    double xc=boite.get_xcentre();
    double yc=boite.get_ycentre();
    double zc=boite.get_zcentre();
    double rayonx=1.125*(boite.get_xmax()-xc);
    double rayony=1.125*(boite.get_ymax()-yc);
    double rayonz=1.125*(boite.get_zmax()-zc);
    boite.reinit(xc-rayonx,yc-rayony,zc-rayonz,xc+rayonx,yc+rayony,zc+rayonz);
    grille.initialiser(boite.get_xmin(),boite.get_ymin(),boite.get_zmin(),boite.get_xmax(),boite.get_ymax(),boite.get_zmax(),pas_estimation,pas_estimation,pas_estimation);
    for (int i=0;i<pas_estimation;i++)
        for (int j=0;j<pas_estimation;j++)
            for (int k=0;k<pas_estimation;k++)
            {
                TPL_CELLULE_GRILLE_INFO<CT_POINT<1> *,DOUBLEN<4> > *cell=grille.get_cellule(i,j,k);
                double tab[4];
                tab[0]=dng;
                tab[1]=0.;
                tab[2]=0.;
                tab[3]=0.;
                DOUBLEN<4> ecart0(tab);
                cell->change_info(0,ecart0);
                cell->change_info(1,ecart0);
                cell->change_info(2,ecart0);
                cell->change_info(3,ecart0);
                cell->change_info(4,ecart0);
                cell->change_info(5,ecart0);
                cell->change_info(6,ecart0);
                cell->change_info(7,ecart0);
            }
    legende[0]="F(x,y,z)";
    legende[1]="dF/dx";
    legende[2]="dF/dy";
    legende[3]="dF/dz";
}