fem_hexa20.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_hexa20.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:55 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
 
 
#include "gestionversion.h"
#include "fem_hexa20.h"
#include "fem_maillage.h"
#include "fem_noeud.h"
#include "mg_element_maillage.h"
#include "ot_tenseur.h"
#include "ot_quadrature_gauss.h"
 
 
FEM_HEXA20::FEM_HEXA20(unsigned long num,class MG_ELEMENT_MAILLAGE* mai,class FEM_NOEUD** tabnoeud):FEM_ELEMENT3(num,mai),FEM_TEMPLATE_ELEMENT<20>(tabnoeud)
{
    if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==3) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
    tab[0]->get_lien_element3()->ajouter(this);
    tab[1]->get_lien_element3()->ajouter(this);
    tab[2]->get_lien_element3()->ajouter(this);
    tab[3]->get_lien_element3()->ajouter(this);
    tab[4]->get_lien_element3()->ajouter(this);
    tab[5]->get_lien_element3()->ajouter(this);
    tab[6]->get_lien_element3()->ajouter(this);
    tab[7]->get_lien_element3()->ajouter(this);
    tab[8]->get_lien_element3()->ajouter(this);
    tab[9]->get_lien_element3()->ajouter(this);
    tab[10]->get_lien_element3()->ajouter(this);
    tab[11]->get_lien_element3()->ajouter(this);
    tab[12]->get_lien_element3()->ajouter(this);
    tab[13]->get_lien_element3()->ajouter(this);
    tab[14]->get_lien_element3()->ajouter(this);
    tab[15]->get_lien_element3()->ajouter(this);
    tab[16]->get_lien_element3()->ajouter(this);
    tab[17]->get_lien_element3()->ajouter(this);
    tab[18]->get_lien_element3()->ajouter(this);
    tab[19]->get_lien_element3()->ajouter(this);
    get_fem_noeudpetitid()->get_lien_petit_element3()->ajouter(this);
}
 
FEM_HEXA20::FEM_HEXA20(class MG_ELEMENT_MAILLAGE* mai,FEM_NOEUD** tabnoeud):FEM_ELEMENT3(mai),FEM_TEMPLATE_ELEMENT<20>(tabnoeud)
{
    if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==3) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
    tab[0]->get_lien_element3()->ajouter(this);
    tab[1]->get_lien_element3()->ajouter(this);
    tab[2]->get_lien_element3()->ajouter(this);
    tab[3]->get_lien_element3()->ajouter(this);
    tab[4]->get_lien_element3()->ajouter(this);
    tab[5]->get_lien_element3()->ajouter(this);
    tab[6]->get_lien_element3()->ajouter(this);
    tab[7]->get_lien_element3()->ajouter(this);
    tab[8]->get_lien_element3()->ajouter(this);
    tab[9]->get_lien_element3()->ajouter(this);
    tab[10]->get_lien_element3()->ajouter(this);
    tab[11]->get_lien_element3()->ajouter(this);
    tab[12]->get_lien_element3()->ajouter(this);
    tab[13]->get_lien_element3()->ajouter(this);
    tab[14]->get_lien_element3()->ajouter(this);
    tab[15]->get_lien_element3()->ajouter(this);
    tab[16]->get_lien_element3()->ajouter(this);
    tab[17]->get_lien_element3()->ajouter(this);
    tab[18]->get_lien_element3()->ajouter(this);
    tab[19]->get_lien_element3()->ajouter(this);
    get_fem_noeudpetitid()->get_lien_petit_element3()->ajouter(this);
}
 
FEM_HEXA20::FEM_HEXA20(unsigned long num,class MG_ELEMENT_TOPOLOGIQUE* topo,class FEM_NOEUD** tabnoeud):FEM_ELEMENT3(num,topo),FEM_TEMPLATE_ELEMENT<20>(tabnoeud)
{
    if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==3) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
    tab[0]->get_lien_element3()->ajouter(this);
    tab[1]->get_lien_element3()->ajouter(this);
    tab[2]->get_lien_element3()->ajouter(this);
    tab[3]->get_lien_element3()->ajouter(this);
    tab[4]->get_lien_element3()->ajouter(this);
    tab[5]->get_lien_element3()->ajouter(this);
    tab[6]->get_lien_element3()->ajouter(this);
    tab[7]->get_lien_element3()->ajouter(this);
    tab[8]->get_lien_element3()->ajouter(this);
    tab[9]->get_lien_element3()->ajouter(this);
    tab[10]->get_lien_element3()->ajouter(this);
    tab[11]->get_lien_element3()->ajouter(this);
    tab[12]->get_lien_element3()->ajouter(this);
    tab[13]->get_lien_element3()->ajouter(this);
    tab[14]->get_lien_element3()->ajouter(this);
    tab[15]->get_lien_element3()->ajouter(this);
    tab[16]->get_lien_element3()->ajouter(this);
    tab[17]->get_lien_element3()->ajouter(this);
    tab[18]->get_lien_element3()->ajouter(this);
    tab[19]->get_lien_element3()->ajouter(this);
    get_fem_noeudpetitid()->get_lien_petit_element3()->ajouter(this);
}
 
FEM_HEXA20::FEM_HEXA20(class MG_ELEMENT_TOPOLOGIQUE* topo,FEM_NOEUD** tabnoeud):FEM_ELEMENT3(topo),FEM_TEMPLATE_ELEMENT<20>(tabnoeud)
{
    if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==3) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
    tab[0]->get_lien_element3()->ajouter(this);
    tab[1]->get_lien_element3()->ajouter(this);
    tab[2]->get_lien_element3()->ajouter(this);
    tab[3]->get_lien_element3()->ajouter(this);
    tab[4]->get_lien_element3()->ajouter(this);
    tab[5]->get_lien_element3()->ajouter(this);
    tab[6]->get_lien_element3()->ajouter(this);
    tab[7]->get_lien_element3()->ajouter(this);
    tab[8]->get_lien_element3()->ajouter(this);
    tab[9]->get_lien_element3()->ajouter(this);
    tab[10]->get_lien_element3()->ajouter(this);
    tab[11]->get_lien_element3()->ajouter(this);
    tab[12]->get_lien_element3()->ajouter(this);
    tab[13]->get_lien_element3()->ajouter(this);
    tab[14]->get_lien_element3()->ajouter(this);
    tab[15]->get_lien_element3()->ajouter(this);
    tab[16]->get_lien_element3()->ajouter(this);
    tab[17]->get_lien_element3()->ajouter(this);
    tab[18]->get_lien_element3()->ajouter(this);
    tab[19]->get_lien_element3()->ajouter(this);
    get_fem_noeudpetitid()->get_lien_petit_element3()->ajouter(this);
}
 
FEM_HEXA20::FEM_HEXA20(unsigned long num,class MG_ELEMENT_TOPOLOGIQUE* topo,class MG_ELEMENT_MAILLAGE* mai,class FEM_NOEUD** tabnoeud):FEM_ELEMENT3(num,topo,mai),FEM_TEMPLATE_ELEMENT<20>(tabnoeud)
{
    if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==3) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
    tab[0]->get_lien_element3()->ajouter(this);
    tab[1]->get_lien_element3()->ajouter(this);
    tab[2]->get_lien_element3()->ajouter(this);
    tab[3]->get_lien_element3()->ajouter(this);
    tab[4]->get_lien_element3()->ajouter(this);
    tab[5]->get_lien_element3()->ajouter(this);
    tab[6]->get_lien_element3()->ajouter(this);
    tab[7]->get_lien_element3()->ajouter(this);
    tab[8]->get_lien_element3()->ajouter(this);
    tab[9]->get_lien_element3()->ajouter(this);
    tab[10]->get_lien_element3()->ajouter(this);
    tab[11]->get_lien_element3()->ajouter(this);
    tab[12]->get_lien_element3()->ajouter(this);
    tab[13]->get_lien_element3()->ajouter(this);
    tab[14]->get_lien_element3()->ajouter(this);
    tab[15]->get_lien_element3()->ajouter(this);
    tab[16]->get_lien_element3()->ajouter(this);
    tab[17]->get_lien_element3()->ajouter(this);
    tab[18]->get_lien_element3()->ajouter(this);
    tab[19]->get_lien_element3()->ajouter(this);
    get_fem_noeudpetitid()->get_lien_petit_element3()->ajouter(this);
}
 
FEM_HEXA20::FEM_HEXA20(class MG_ELEMENT_TOPOLOGIQUE* topo,class MG_ELEMENT_MAILLAGE* mai,FEM_NOEUD** tabnoeud):FEM_ELEMENT3(topo,mai),FEM_TEMPLATE_ELEMENT<20>(tabnoeud)
{
    if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==3) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
    tab[0]->get_lien_element3()->ajouter(this);
    tab[1]->get_lien_element3()->ajouter(this);
    tab[2]->get_lien_element3()->ajouter(this);
    tab[3]->get_lien_element3()->ajouter(this);
    tab[4]->get_lien_element3()->ajouter(this);
    tab[5]->get_lien_element3()->ajouter(this);
    tab[6]->get_lien_element3()->ajouter(this);
    tab[7]->get_lien_element3()->ajouter(this);
    tab[8]->get_lien_element3()->ajouter(this);
    tab[9]->get_lien_element3()->ajouter(this);
    tab[10]->get_lien_element3()->ajouter(this);
    tab[11]->get_lien_element3()->ajouter(this);
    tab[12]->get_lien_element3()->ajouter(this);
    tab[13]->get_lien_element3()->ajouter(this);
    tab[14]->get_lien_element3()->ajouter(this);
    tab[15]->get_lien_element3()->ajouter(this);
    tab[16]->get_lien_element3()->ajouter(this);
    tab[17]->get_lien_element3()->ajouter(this);
    tab[18]->get_lien_element3()->ajouter(this);
    tab[19]->get_lien_element3()->ajouter(this);
    get_fem_noeudpetitid()->get_lien_petit_element3()->ajouter(this);
}
 
 
 
 
 
 
 
 
 
FEM_HEXA20::FEM_HEXA20(FEM_HEXA20& mdd):FEM_ELEMENT3(mdd),FEM_TEMPLATE_ELEMENT<20>(mdd)
{
    if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==3) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
    tab[0]->get_lien_element3()->ajouter(this);
    tab[1]->get_lien_element3()->ajouter(this);
    tab[2]->get_lien_element3()->ajouter(this);
    tab[3]->get_lien_element3()->ajouter(this);
    tab[4]->get_lien_element3()->ajouter(this);
    tab[5]->get_lien_element3()->ajouter(this);
    tab[6]->get_lien_element3()->ajouter(this);
    tab[7]->get_lien_element3()->ajouter(this);
    tab[8]->get_lien_element3()->ajouter(this);
    tab[9]->get_lien_element3()->ajouter(this);
    tab[10]->get_lien_element3()->ajouter(this);
    tab[11]->get_lien_element3()->ajouter(this);
    tab[12]->get_lien_element3()->ajouter(this);
    tab[13]->get_lien_element3()->ajouter(this);
    tab[14]->get_lien_element3()->ajouter(this);
    tab[15]->get_lien_element3()->ajouter(this);
    tab[16]->get_lien_element3()->ajouter(this);
    tab[17]->get_lien_element3()->ajouter(this);
    tab[18]->get_lien_element3()->ajouter(this);
    tab[19]->get_lien_element3()->ajouter(this);
    get_fem_noeudpetitid()->get_lien_petit_element3()->ajouter(this);
}
FEM_HEXA20::~FEM_HEXA20()
{
    if (liaison_topologique!=NULL) 
      if (liaison_topologique->get_dimension()==3) liaison_topologique->get_lien_fem_maillage()->supprimer(this);
    tab[0]->get_lien_element3()->supprimer(this);
    tab[1]->get_lien_element3()->supprimer(this);
    tab[2]->get_lien_element3()->supprimer(this);
    tab[3]->get_lien_element3()->supprimer(this);
    tab[4]->get_lien_element3()->supprimer(this);
    tab[5]->get_lien_element3()->supprimer(this);
    tab[6]->get_lien_element3()->supprimer(this);
    tab[7]->get_lien_element3()->supprimer(this);
    tab[8]->get_lien_element3()->supprimer(this);
    tab[9]->get_lien_element3()->supprimer(this);
    tab[10]->get_lien_element3()->supprimer(this);
    tab[11]->get_lien_element3()->supprimer(this);
    tab[12]->get_lien_element3()->supprimer(this);
    tab[13]->get_lien_element3()->supprimer(this);
    tab[14]->get_lien_element3()->supprimer(this);
    tab[15]->get_lien_element3()->supprimer(this);
    tab[16]->get_lien_element3()->supprimer(this);
    tab[17]->get_lien_element3()->supprimer(this);
    tab[18]->get_lien_element3()->supprimer(this);
    tab[19]->get_lien_element3()->supprimer(this);
    get_fem_noeudpetitid()->get_lien_petit_element3()->supprimer(this);
}
 
 
FEM_ELEMENT_MAILLAGE* FEM_HEXA20::dupliquer(FEM_MAILLAGE *femmai,long decalage)
{
    FEM_NOEUD* tabnoeud[20];
    tabnoeud[0]=femmai->get_fem_noeudid(tab[0]->get_id()+decalage);
    tabnoeud[1]=femmai->get_fem_noeudid(tab[1]->get_id()+decalage);
    tabnoeud[2]=femmai->get_fem_noeudid(tab[2]->get_id()+decalage);
    tabnoeud[3]=femmai->get_fem_noeudid(tab[3]->get_id()+decalage);
    tabnoeud[4]=femmai->get_fem_noeudid(tab[4]->get_id()+decalage);
    tabnoeud[5]=femmai->get_fem_noeudid(tab[5]->get_id()+decalage);
    tabnoeud[6]=femmai->get_fem_noeudid(tab[6]->get_id()+decalage);
    tabnoeud[7]=femmai->get_fem_noeudid(tab[7]->get_id()+decalage);
    tabnoeud[8]=femmai->get_fem_noeudid(tab[7]->get_id()+decalage);
    tabnoeud[9]=femmai->get_fem_noeudid(tab[7]->get_id()+decalage);
    tabnoeud[10]=femmai->get_fem_noeudid(tab[7]->get_id()+decalage);
    tabnoeud[11]=femmai->get_fem_noeudid(tab[7]->get_id()+decalage);
    tabnoeud[12]=femmai->get_fem_noeudid(tab[7]->get_id()+decalage);
    tabnoeud[13]=femmai->get_fem_noeudid(tab[7]->get_id()+decalage);
    tabnoeud[14]=femmai->get_fem_noeudid(tab[7]->get_id()+decalage);
    tabnoeud[15]=femmai->get_fem_noeudid(tab[7]->get_id()+decalage);
    tabnoeud[16]=femmai->get_fem_noeudid(tab[7]->get_id()+decalage);
    tabnoeud[17]=femmai->get_fem_noeudid(tab[7]->get_id()+decalage);
    tabnoeud[18]=femmai->get_fem_noeudid(tab[7]->get_id()+decalage);
    tabnoeud[19]=femmai->get_fem_noeudid(tab[7]->get_id()+decalage);
    FEM_HEXA20* quad=new FEM_HEXA20(get_id()+decalage,maillage,tabnoeud);
    femmai->ajouter_fem_element3(quad);
    return quad;
}
 
 
void FEM_HEXA20::get_voisin_noeud(class FEM_NOEUD* no,TPL_LISTE_ENTITE<FEM_NOEUD*> &voisin)
{
voisin.vide();
if (no==tab[0])
      {
      voisin.ajouter(tab[1]);
      voisin.ajouter(tab[7]);
      voisin.ajouter(tab[8]);
      }
if (no==tab[1])
      {
      voisin.ajouter(tab[0]);
      voisin.ajouter(tab[2]);
      }
if (no==tab[2])
      {
      voisin.ajouter(tab[1]);
      voisin.ajouter(tab[3]);
      voisin.ajouter(tab[9]);
      }
if (no==tab[3])
      {
      voisin.ajouter(tab[2]);
      voisin.ajouter(tab[4]);
      }
if (no==tab[4])
      {
      voisin.ajouter(tab[3]);
      voisin.ajouter(tab[5]);
      voisin.ajouter(tab[10]);
      }
if (no==tab[5])
      {
      voisin.ajouter(tab[4]);
      voisin.ajouter(tab[6]);
      }
if (no==tab[6])
      {
      voisin.ajouter(tab[5]);
      voisin.ajouter(tab[7]);
      voisin.ajouter(tab[11]);
      }
if (no==tab[7])
      {
      voisin.ajouter(tab[0]);
      voisin.ajouter(tab[6]);
      }
if (no==tab[8])
      {
      voisin.ajouter(tab[0]);
      voisin.ajouter(tab[12]);
      }
if (no==tab[9])
      {
      voisin.ajouter(tab[2]);
      voisin.ajouter(tab[14]);
      }
if (no==tab[10])
      {
      voisin.ajouter(tab[4]);
      voisin.ajouter(tab[16]);
      }
if (no==tab[11])
      {
      voisin.ajouter(tab[6]);
      voisin.ajouter(tab[18]);
      }
 
if (no==tab[12])
      {
      voisin.ajouter(tab[8]);
      voisin.ajouter(tab[13]);
      voisin.ajouter(tab[19]);
      }
if (no==tab[13])
      {
      voisin.ajouter(tab[12]);
      voisin.ajouter(tab[14]);
      }
if (no==tab[14])
      {
      voisin.ajouter(tab[9]);
      voisin.ajouter(tab[13]);
      voisin.ajouter(tab[15]);
      }
if (no==tab[15])
      {
      voisin.ajouter(tab[14]);
      voisin.ajouter(tab[16]);
      }
if (no==tab[16])
      {
      voisin.ajouter(tab[10]);
      voisin.ajouter(tab[15]);
      voisin.ajouter(tab[17]);
      }
if (no==tab[17])
      {
      voisin.ajouter(tab[16]);
      voisin.ajouter(tab[18]);
      }
if (no==tab[18])
      {
      voisin.ajouter(tab[11]);
      voisin.ajouter(tab[17]);
      voisin.ajouter(tab[19]);
      }
if (no==tab[19])
      {
      voisin.ajouter(tab[18]);
      voisin.ajouter(tab[12]);
      }      
}
 
 
int FEM_HEXA20::get_type_entite(void)
{
    return   MAGIC::TYPE_ENTITE::IDFEM_HEXA20;
}
 
int FEM_HEXA20::get_dimension(void)
{
    return 3;
}
 
 
 
int FEM_HEXA20::get_nb_fem_noeud(void)
{
    return FEM_TEMPLATE_ELEMENT<20>::get_nb_fem_noeud();
}
 
FEM_NOEUD* FEM_HEXA20::get_fem_noeud(int num)
{
    return FEM_TEMPLATE_ELEMENT<20>::get_fem_noeud(num);
}
 
void FEM_HEXA20::change_noeud(int num,FEM_NOEUD* noeud)
{
    FEM_TEMPLATE_ELEMENT<20>::change_noeud(num,noeud);
}
 
BOITE_3D& FEM_HEXA20::get_boite_3D(void)
{
    return FEM_TEMPLATE_ELEMENT<20>::get_boite_3D();
}
 
 
void FEM_HEXA20::enregistrer(std::ostream& o,double version)
{
    if (maillage!=NULL)
      if (get_lien_topologie()!=NULL) o << "%" << get_id() << "=FEM_HEXA20($"<< get_lien_topologie()->get_id() << ",$" << maillage->get_id() << ",$" << tab[0]->get_id() << ",$" << tab[1]->get_id() << ",$" << tab[2]->get_id()<<  ",$" << tab[3]->get_id() << ",$" << tab[4]->get_id() << ",$" << tab[5]->get_id() << ",$" << tab[6]->get_id()<<  ",$" << tab[7]->get_id()<< ",$"  << tab[8]->get_id() << ",$" << tab[9]->get_id() << ",$" << tab[10]->get_id()<<  ",$" << tab[11]->get_id() << ",$" << tab[12]->get_id() << ",$" << tab[13]->get_id() << ",$" << tab[14]->get_id() <<  ",$" << tab[15]->get_id() <<  ",$" << tab[16]->get_id()<<  ",$" << tab[17]->get_id()<<  ",$" << tab[18]->get_id()<<  ",$" << tab[19]->get_id()<<  ");" << std::endl;
      else o << "%" << get_id() << "=FEM_HEXA20(NULL,$" << maillage->get_id() << ",$" << tab[0]->get_id() << ",$" << tab[1]->get_id() << ",$" << tab[2]->get_id()<<  ",$" << tab[3]->get_id() << ",$" << tab[4]->get_id() << ",$" << tab[5]->get_id() << ",$" << tab[6]->get_id()<<  ",$" << tab[7]->get_id()<< ",$"  << tab[8]->get_id() << ",$" << tab[9]->get_id() << ",$" << tab[10]->get_id()<<  ",$" << tab[11]->get_id() << ",$" << tab[12]->get_id() << ",$" << tab[13]->get_id() << ",$" << tab[14]->get_id() <<  ",$" << tab[15]->get_id() <<  ",$" << tab[16]->get_id()<<  ",$" << tab[17]->get_id()<<  ",$" << tab[18]->get_id()<<  ",$" << tab[19]->get_id()<<");" << std::endl;
    else
      if (get_lien_topologie()!=NULL) o << "%" << get_id() << "=FEM_HEXA20($"<< get_lien_topologie()->get_id() << ",NULL,$" << tab[0]->get_id() << ",$" << tab[1]->get_id() << ",$" << tab[2]->get_id()<<  ",$" << tab[3]->get_id() << ",$" << tab[4]->get_id() << ",$" << tab[5]->get_id() << ",$" << tab[6]->get_id()<<  ",$" << tab[7]->get_id()<< ",$"  << tab[8]->get_id() << ",$" << tab[9]->get_id() << ",$" << tab[10]->get_id()<<  ",$" << tab[11]->get_id() << ",$" << tab[12]->get_id() << ",$" << tab[13]->get_id() << ",$" << tab[14]->get_id() <<  ",$" << tab[15]->get_id() <<  ",$" << tab[16]->get_id()<<  ",$" << tab[17]->get_id()<<  ",$" << tab[18]->get_id()<<  ",$" << tab[19]->get_id()<<  ");" << std::endl;
      else o << "%" << get_id() << "=FEM_HEXA20(NULL,NULL,$" << tab[0]->get_id() << ",$" << tab[1]->get_id() << ",$" << tab[2]->get_id()<<  ",$" << tab[3]->get_id() << ",$" << tab[4]->get_id() << ",$" << tab[5]->get_id() << ",$" << tab[6]->get_id()<<  ",$" << tab[7]->get_id()<< ",$"  << tab[8]->get_id() << ",$" << tab[9]->get_id() << ",$" << tab[10]->get_id()<<  ",$" << tab[11]->get_id() << ",$" << tab[12]->get_id() << ",$" << tab[13]->get_id() << ",$" << tab[14]->get_id() <<  ",$" << tab[15]->get_id() <<  ",$" << tab[16]->get_id()<<  ",$" << tab[17]->get_id()<<  ",$" << tab[18]->get_id()<<  ",$" << tab[19]->get_id()<<");" << std::endl;
      
}
 
 
 
int FEM_HEXA20::get_nb_pt_gauss(int degre)
{
return OT_POINTS_GAUSS::get_nb_point_hexa_prod(degre);
}
 
void FEM_HEXA20::get_pt_gauss(int degre,int num,double &w,double *uvw)
{
    return OT_POINTS_GAUSS::get_pt_gauss_hex_prod(degre,num,w,uvw);
}
int FEM_HEXA20::get_nb_fonction_interpolation(void)
{
    return 20;
}
int FEM_HEXA20::get_degremax_fonction_interpolation(void)
{
     return 4;  
}
 
int FEM_HEXA20::get_degre_gauss(int num)
{
return OT_POINTS_GAUSS::get_degre_gauss_hexa_prod(num);  
}
 
 
double FEM_HEXA20::get_fonction_interpolation(int num,double *uv)
{
    double val;
    switch (num)
    {
    case 1:
        {double a=-1;
        double b=-1;
        double c=-1;
        val=0.125*(1.+uv[0]*a)*(1.+uv[1]*b)*(1.+uv[2]*c)*(-2.+uv[0]*a+uv[1]*b+uv[2]*c);
        break;}
    case 3:
        {double a=1;
        double b=-1;
        double c=-1;
        val=0.125*(1.+uv[0]*a)*(1.+uv[1]*b)*(1.+uv[2]*c)*(-2.+uv[0]*a+uv[1]*b+uv[2]*c);
        break;}
    case 5:
        {double a=1;
        double b=1;
        double c=-1;
        val=0.125*(1.+uv[0]*a)*(1.+uv[1]*b)*(1.+uv[2]*c)*(-2.+uv[0]*a+uv[1]*b+uv[2]*c);
        break;}
    case 7:
       {double a=-1;
        double b=1;
        double c=-1;
        val=0.125*(1.+uv[0]*a)*(1.+uv[1]*b)*(1.+uv[2]*c)*(-2.+uv[0]*a+uv[1]*b+uv[2]*c);
        break;}
    case 13:
        {double a=-1;
        double b=-1;
        double c=1;
        val=0.125*(1.+uv[0]*a)*(1.+uv[1]*b)*(1.+uv[2]*c)*(-2.+uv[0]*a+uv[1]*b+uv[2]*c);
        break;}
    case 15:
        {double a=1;
        double b=-1;
        double c=1;
        val=0.125*(1.+uv[0]*a)*(1.+uv[1]*b)*(1.+uv[2]*c)*(-2.+uv[0]*a+uv[1]*b+uv[2]*c);
        break;}
    case 17:
        {double a=1;
        double b=1;
        double c=1;
        val=0.125*(1.+uv[0]*a)*(1.+uv[1]*b)*(1.+uv[2]*c)*(-2.+uv[0]*a+uv[1]*b+uv[2]*c);
        break;}
    case 19:
        {double a=-1;
        double b=1;
        double c=1;
        val=0.125*(1.+uv[0]*a)*(1.+uv[1]*b)*(1.+uv[2]*c)*(-2.+uv[0]*a+uv[1]*b+uv[2]*c);
        break;}
    case 2:
        {double a=0.;
        double b=-1;
        double c=-1;
        val=0.25*(1.-uv[0]*uv[0])*(1.+uv[1]*b)*(1.+uv[2]*c);
        break;}
    case 6:
        {double a=0.;
        double b=1;
        double c=-1;
        val=0.25*(1.-uv[0]*uv[0])*(1.+uv[1]*b)*(1.+uv[2]*c);
        break;}
    case 14:
        {double a=0.;
        double b=-1;
        double c=1;
        val=0.25*(1.-uv[0]*uv[0])*(1.+uv[1]*b)*(1.+uv[2]*c);
        break;}
    case 18:
        {double a=0.;
        double b=1;
        double c=1;
        val=0.25*(1.-uv[0]*uv[0])*(1.+uv[1]*b)*(1.+uv[2]*c);
        break;}
    case 4:
        {double a=1.;
        double b=0.;
        double c=-1;
        val=0.25*(1.+uv[0]*a)*(1.-uv[1]*uv[1])*(1.+uv[2]*c);
        break;}
    case 8:
        {double a=-1.;
        double b=0.;
        double c=-1;
        val=0.25*(1.+uv[0]*a)*(1.-uv[1]*uv[1])*(1.+uv[2]*c);
        break;}
    case 16:
        {double a=1.;
        double b=0.;
        double c=1;
        val=0.25*(1.+uv[0]*a)*(1.-uv[1]*uv[1])*(1.+uv[2]*c);
        break;}
    case 20:
        {double a=-1.;
        double b=0.;
        double c=1;
        val=0.25*(1.+uv[0]*a)*(1.-uv[1]*uv[1])*(1.+uv[2]*c);
        break;}
    case 9:
        {double a=-1.;
        double b=-1.;
        double c=0.;
        val=0.25*(1.+uv[0]*a)*(1.+uv[1]*b)*(1.-uv[2]*uv[2]);
        break;}
    case 10:
        {double a=1.;
        double b=-1.;
        double c=0.;
        val=0.25*(1.+uv[0]*a)*(1.+uv[1]*b)*(1.-uv[2]*uv[2]);
        break;}
    case 11:
        {double a=1.;
        double b=1.;
        double c=0.;
        val=0.25*(1.+uv[0]*a)*(1.+uv[1]*b)*(1.-uv[2]*uv[2]);
        break;}
    case 12:
        {double a=-1.;
        double b=1.;
        double c=0.;
        val=0.25*(1.+uv[0]*a)*(1.+uv[1]*b)*(1.-uv[2]*uv[2]);
        break;}
        
 
    }
    return val;
}
void FEM_HEXA20::reinit_boite_3D(void)
{
FEM_TEMPLATE_ELEMENT<20>::reinit_boite_3D();  
}
 
double FEM_HEXA20::get_fonction_derive_interpolation(int num,int num_variable,double *uv)
{
    double val;
    switch (num)
    {
    case 1:
        {
        double a=-1;
        double b=-1;
        double c=-1;
        switch (num_variable)
        {
        case 1:
            val=0.125*a*(1+uv[1]*b)*(1+uv[2]*c)*(-1+2*uv[0]*a+uv[1]*b+uv[2]*c);
            break;
        case 2:
            val=0.125*b*(1+uv[0]*a)*(1+uv[2]*c)*(-1+uv[0]*a+2*uv[1]*b+uv[2]*c);
            break;
        case 3:
            val=0.125*c*(1+uv[0]*a)*(1+uv[1]*b)*(-1+uv[0]*a+uv[1]*b+2*uv[2]*c);
            break;
        } break;}
    case 3:
        {
        double a=1;
        double b=-1;
        double c=-1;
        switch (num_variable)
        {
        case 1:
            val=0.125*a*(1+uv[1]*b)*(1+uv[2]*c)*(-1+2*uv[0]*a+uv[1]*b+uv[2]*c);
            break;
        case 2:
            val=0.125*b*(1+uv[0]*a)*(1+uv[2]*c)*(-1+uv[0]*a+2*uv[1]*b+uv[2]*c);
            break;
        case 3:
            val=0.125*c*(1+uv[0]*a)*(1+uv[1]*b)*(-1+uv[0]*a+uv[1]*b+2*uv[2]*c);
            break;
        } break;}
    case 5:
        {
        double a=1;
        double b=1;
        double c=-1;
        switch (num_variable)
        {
        case 1:
            val=0.125*a*(1+uv[1]*b)*(1+uv[2]*c)*(-1+2*uv[0]*a+uv[1]*b+uv[2]*c);
            break;
        case 2:
            val=0.125*b*(1+uv[0]*a)*(1+uv[2]*c)*(-1+uv[0]*a+2*uv[1]*b+uv[2]*c);
            break;
        case 3:
            val=0.125*c*(1+uv[0]*a)*(1+uv[1]*b)*(-1+uv[0]*a+uv[1]*b+2*uv[2]*c);
            break;
        } break;}
    case 7:
        {
        double a=-1;
        double b=1;
        double c=-1;
        switch (num_variable)
        {
        case 1:
            val=0.125*a*(1+uv[1]*b)*(1+uv[2]*c)*(-1+2*uv[0]*a+uv[1]*b+uv[2]*c);
            break;
        case 2:
            val=0.125*b*(1+uv[0]*a)*(1+uv[2]*c)*(-1+uv[0]*a+2*uv[1]*b+uv[2]*c);
            break;
        case 3:
            val=0.125*c*(1+uv[0]*a)*(1+uv[1]*b)*(-1+uv[0]*a+uv[1]*b+2*uv[2]*c);
            break;
        } break;}
    case 13:
        {
        double a=-1;
        double b=-1;
        double c=1;
        switch (num_variable)
        {
        case 1:
            val=0.125*a*(1+uv[1]*b)*(1+uv[2]*c)*(-1+2*uv[0]*a+uv[1]*b+uv[2]*c);
            break;
        case 2:
            val=0.125*b*(1+uv[0]*a)*(1+uv[2]*c)*(-1+uv[0]*a+2*uv[1]*b+uv[2]*c);
            break;
        case 3:
            val=0.125*c*(1+uv[0]*a)*(1+uv[1]*b)*(-1+uv[0]*a+uv[1]*b+2*uv[2]*c);
            break;
        } break;}
    case 15:
        {
        double a=1;
        double b=-1;
        double c=1;
        switch (num_variable)
        {
        case 1:
            val=0.125*a*(1+uv[1]*b)*(1+uv[2]*c)*(-1+2*uv[0]*a+uv[1]*b+uv[2]*c);
            break;
        case 2:
            val=0.125*b*(1+uv[0]*a)*(1+uv[2]*c)*(-1+uv[0]*a+2*uv[1]*b+uv[2]*c);
            break;
        case 3:
            val=0.125*c*(1+uv[0]*a)*(1+uv[1]*b)*(-1+uv[0]*a+uv[1]*b+2*uv[2]*c);
            break;
        } break;}
    case 17:
        {
        double a=1;
        double b=1;
        double c=1;
        switch (num_variable)
        {
        case 1:
            val=0.125*a*(1+uv[1]*b)*(1+uv[2]*c)*(-1+2*uv[0]*a+uv[1]*b+uv[2]*c);
            break;
        case 2:
            val=0.125*b*(1+uv[0]*a)*(1+uv[2]*c)*(-1+uv[0]*a+2*uv[1]*b+uv[2]*c);
            break;
        case 3:
            val=0.125*c*(1+uv[0]*a)*(1+uv[1]*b)*(-1+uv[0]*a+uv[1]*b+2*uv[2]*c);
            break;
        } break;}
    case 19:
        {
        double a=-1;
        double b=1;
        double c=1;
        switch (num_variable)
        {
        case 1:
            val=0.125*a*(1+uv[1]*b)*(1+uv[2]*c)*(-1+2*uv[0]*a+uv[1]*b+uv[2]*c);
            break;
        case 2:
            val=0.125*b*(1+uv[0]*a)*(1+uv[2]*c)*(-1+uv[0]*a+2*uv[1]*b+uv[2]*c);
            break;
        case 3:
            val=0.125*c*(1+uv[0]*a)*(1+uv[1]*b)*(-1+uv[0]*a+uv[1]*b+2*uv[2]*c);
            break;
        } break;}
    case 2:
        {
        double a=0;
        double b=-1;
        double c=-1;
        switch (num_variable)
        {
        case 1:
            val=-0.5*uv[0]*(1+uv[1]*b)*(1+uv[2]*c);
            break;
        case 2:
            val=0.25*b*(1-uv[0]*uv[0])*(1+uv[2]*c);
            break;
        case 3:
            val=0.25*c*(1-uv[0]*uv[0])*(1+uv[1]*b);
            break;
        } break;}
    case 18:
        {
        double a=0;
        double b=1;
        double c=1;
        switch (num_variable)
        {
        case 1:
            val=-0.5*uv[0]*(1+uv[1]*b)*(1+uv[2]*c);
            break;
        case 2:
            val=0.25*b*(1-uv[0]*uv[0])*(1+uv[2]*c);
            break;
        case 3:
            val=0.25*c*(1-uv[0]*uv[0])*(1+uv[1]*b);
            break;
        } break;}
    case 14:
        {
        double a=0;
        double b=-1;
        double c=1;
        switch (num_variable)
        {
        case 1:
            val=-0.5*uv[0]*(1+uv[1]*b)*(1+uv[2]*c);
            break;
        case 2:
            val=0.25*b*(1-uv[0]*uv[0])*(1+uv[2]*c);
            break;
        case 3:
            val=0.25*c*(1-uv[0]*uv[0])*(1+uv[1]*b);
            break;
        } break;}
    case 6:
        {
        double a=0;
        double b=1;
        double c=-1;
        switch (num_variable)
        {
        case 1:
            val=-0.5*uv[0]*(1+uv[1]*b)*(1+uv[2]*c);
            break;
        case 2:
            val=0.25*b*(1-uv[0]*uv[0])*(1+uv[2]*c);
            break;
        case 3:
            val=0.25*c*(1-uv[0]*uv[0])*(1+uv[1]*b);
            break;
        } break;}
    case 4:
        {
        double a=1;
        double b=0;
        double c=-1;
        switch (num_variable)
        {
        case 1:
            val=0.25*a*(1-uv[1]*uv[1])*(1+uv[2]*c);
            break;
        case 2:
            val=-0.5*uv[1]*(1+uv[0]*a)*(1+uv[2]*c);
            break;
        case 3:
            val=0.25*c*(1+uv[0]*a)*(1-uv[1]*uv[1]);
            break;
        } break;}
    case 8:
        {
        double a=-1;
        double b=0;
        double c=-1;
        switch (num_variable)
        {
        case 1:
            val=0.25*a*(1-uv[1]*uv[1])*(1+uv[2]*c);
            break;
        case 2:
            val=-0.5*uv[1]*(1+uv[0]*a)*(1+uv[2]*c);
            break;
        case 3:
            val=0.25*c*(1+uv[0]*a)*(1-uv[1]*uv[1]);
            break;
        } break;}
    case 16:
        {
        double a=1;
        double b=0;
        double c=1;
        switch (num_variable)
        {
        case 1:
            val=0.25*a*(1-uv[1]*uv[1])*(1+uv[2]*c);
            break;
        case 2:
            val=-0.5*uv[1]*(1+uv[0]*a)*(1+uv[2]*c);
            break;
        case 3:
            val=0.25*c*(1+uv[0]*a)*(1-uv[1]*uv[1]);
            break;
        } break;}
    case 20:
        {
        double a=-1;
        double b=0;
        double c=1;
        switch (num_variable)
        {
        case 1:
            val=0.25*a*(1-uv[1]*uv[1])*(1+uv[2]*c);
            break;
        case 2:
            val=-0.5*uv[1]*(1+uv[0]*a)*(1+uv[2]*c);
            break;
        case 3:
            val=0.25*c*(1+uv[0]*a)*(1-uv[1]*uv[1]);
            break;
        } break;}
    case 9:
        {
        double a=-1;
        double b=-1;
        double c=0.;
        switch (num_variable)
        {
        case 1:
            val=0.25*a*(1+uv[1]*b)*(1-uv[2]*uv[2]);
            break;
        case 2:
            val=0.25*b*(1+uv[0]*a)*(1-uv[2]*uv[2]);
            break;
        case 3:
            val=-0.5*uv[2]*(1+uv[0]*a)*(1+uv[1]*b);
            break;
        } break;}
    case 10:
        {
        double a=1;
        double b=-1;
        double c=0.;
        switch (num_variable)
        {
        case 1:
            val=0.25*a*(1+uv[1]*b)*(1-uv[2]*uv[2]);
            break;
        case 2:
            val=0.25*b*(1+uv[0]*a)*(1-uv[2]*uv[2]);
            break;
        case 3:
            val=-0.5*uv[2]*(1+uv[0]*a)*(1+uv[1]*b);
            break;
        } break;}
    case 11:
        {
        double a=1;
        double b=1;
        double c=0.;
        switch (num_variable)
        {
        case 1:
            val=0.25*a*(1+uv[1]*b)*(1-uv[2]*uv[2]);
            break;
        case 2:
            val=0.25*b*(1+uv[0]*a)*(1-uv[2]*uv[2]);
            break;
        case 3:
            val=-0.5*uv[2]*(1+uv[0]*a)*(1+uv[1]*b);
            break;
        } break;}
    case 12:
        {
        double a=-1;
        double b=1;
        double c=0.;
        switch (num_variable)
        {
        case 1:
            val=0.25*a*(1+uv[1]*b)*(1-uv[2]*uv[2]);
            break;
        case 2:
            val=0.25*b*(1+uv[0]*a)*(1-uv[2]*uv[2]);
            break;
        case 3:
            val=-0.5*uv[2]*(1+uv[0]*a)*(1+uv[1]*b);
            break;
        } break;}
 
    }
    return val;
}
 
 
 
bool FEM_HEXA20::valide_parametre_EF(double* uvw)
{
if (uvw[0]>=-1.-1e-10)
   if (uvw[1]>=-1.-1e-10)
     if (uvw[2]>=-1.-1e-10)
        if (uvw[0]<=1.+1e-10)  
            if (uvw[1]<=1.+1e-10)  
               if (uvw[2]<=1.+1e-10)  
                     return true;
return false;
}
 
int FEM_HEXA20::verifie_validite_decoupage_xfem(double *vol)
{
return 0;  
}