fem_template_element.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/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       fem_template_element.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:55 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
 
 
#include "gestionversion.h"
#include <math.h>
#include "fem_template_element.h"
#include "fem_noeud.h"
#include "mg_element_maillage.h"
 
template <int N> FEM_TEMPLATE_ELEMENT<N>::FEM_TEMPLATE_ELEMENT()
{
}
 
 
template <int N> FEM_TEMPLATE_ELEMENT<N>::FEM_TEMPLATE_ELEMENT(FEM_NOEUD** tabnoeud)
{
    double xmax=-1e308,ymax=-1e308,zmax=-1e308;
    double xmin=1e308,ymin=1e308,zmin=1e308;
    for (int i=0;i<N;i++)
    {
        tab[i]=tabnoeud[i];
        FEM_NOEUD* noeud=tab[i];
        double x=noeud->get_x();
        if (x<xmin) xmin=x;
        if (x>xmax) xmax=x;
        double y=noeud->get_y();
        if (y<ymin) ymin=y;
        if (y>ymax) ymax=y;
        double z=noeud->get_z();
        if (z<zmin) zmin=z;
        if (z>zmax) zmax=z;
    }
    boite.reinit(xmin,ymin,zmin,xmax,ymax,zmax);
}
 
 
template <int N> void FEM_TEMPLATE_ELEMENT<N>::reinit_boite_3D(void)
{
    double xmax=-1e308,ymax=-1e308,zmax=-1e308;
    double xmin=1e308,ymin=1e308,zmin=1e308;
    for (int i=0;i<N;i++)
    {
        FEM_NOEUD* noeud=tab[i];
        double x=noeud->get_x();
        if (x<xmin) xmin=x;
        if (x>xmax) xmax=x;
        double y=noeud->get_y();
        if (y<ymin) ymin=y;
        if (y>ymax) ymax=y;
        double z=noeud->get_z();
        if (z<zmin) zmin=z;
        if (z>zmax) zmax=z;
    }
    boite.reinit(xmin,ymin,zmin,xmax,ymax,zmax);
}
 
 
 
template <int N> FEM_TEMPLATE_ELEMENT<N>::FEM_TEMPLATE_ELEMENT(FEM_TEMPLATE_ELEMENT& mdd)
{
    for (int i=0;i<N;i++)
        tab[i]=mdd.tab[i];
    boite=mdd.boite;
}
 
 
 
template <int N> FEM_TEMPLATE_ELEMENT<N>::~FEM_TEMPLATE_ELEMENT()
{
}
 
 
 
template <int N> int FEM_TEMPLATE_ELEMENT<N>::get_nb_fem_noeud(void)
{
    return N;
}
 
 
template <int N> FEM_NOEUD* FEM_TEMPLATE_ELEMENT<N>::get_fem_noeud(int num)
{
    return tab[num];
}
 
 
template <int N> void FEM_TEMPLATE_ELEMENT<N>::change_noeud(int num,FEM_NOEUD* noeud)
{
    tab[num]=noeud;
}
 
template <int N> FEM_NOEUD* FEM_TEMPLATE_ELEMENT<N>::get_fem_noeudpetitid(void)
{
    long id=tab[0]->get_id();
    FEM_NOEUD* noeud=tab[0];
    for (int i=0;i<N;i++)
    {
        long id2=tab[i]->get_id();
        if (id2<id)
        {
            id=id2;
            noeud=tab[i];
        }
    }
    return noeud;
}
 
template <int N> BOITE_3D& FEM_TEMPLATE_ELEMENT<N>::get_boite_3D(void)
{
    return boite;
}
 
#ifdef BORLANDCPP
template class __export FEM_TEMPLATE_ELEMENT<1>;
template class __export FEM_TEMPLATE_ELEMENT<2>;
template class __export FEM_TEMPLATE_ELEMENT<3>;
template class __export FEM_TEMPLATE_ELEMENT<6>;
template class __export FEM_TEMPLATE_ELEMENT<4>;
template class __export FEM_TEMPLATE_ELEMENT<10>;
template class __export FEM_TEMPLATE_ELEMENT<8>;
template class __export FEM_TEMPLATE_ELEMENT<20>;
template class __export FEM_TEMPLATE_ELEMENT<15>;
 
#else
#ifdef GCC
template class FEM_TEMPLATE_ELEMENT<1>;
template class FEM_TEMPLATE_ELEMENT<2>;
template class FEM_TEMPLATE_ELEMENT<3>;
template class FEM_TEMPLATE_ELEMENT<6>;
template class FEM_TEMPLATE_ELEMENT<4>;
template class FEM_TEMPLATE_ELEMENT<10>;
template class FEM_TEMPLATE_ELEMENT<8>;
template class FEM_TEMPLATE_ELEMENT<20>;
template class FEM_TEMPLATE_ELEMENT<15>;
#else
#endif
#endif