fct_taille_3d.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_3d.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:52 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
 
 
#include "gestionversion.h"
 
 
#pragma hdrstop
 
#include "fct_taille_3d.h"
 
 
#pragma package(smart_init)
 
 
FCT_TAILLE_3D::FCT_TAILLE_3D(FCT_GENERATEUR_3D<4> *carte):FCT_TAILLE_METRIQUE(),fctgen(carte)
{
}
 
FCT_TAILLE_3D::~FCT_TAILLE_3D()
{
}
 
 
void FCT_TAILLE_3D::evaluer(double *xyz,double *resultat)
{
    double res;
    fctgen->evaluer(xyz,&res);
    resultat[0]=1./res/res;
    resultat[1]=0.;
    resultat[2]=0.;
    resultat[3]=0.;
    resultat[4]=1./res/res;
    resultat[5]=0.;
    resultat[6]=0.;
    resultat[7]=0.;
    resultat[8]=1./res/res;
}
 
void FCT_TAILLE_3D::deriver(double *xyz,double *resultat,int num_param)
{
    double res1,res2;
    fctgen->evaluer(xyz,&res1);
    fctgen->deriver(xyz,&res2,num_param);
    resultat[0]=-2.*res2/res1/res1/res1;
    resultat[1]=0.;
    resultat[2]=0.;
    resultat[3]=0.;
    resultat[4]=-2.*res2/res1/res1/res1;
    resultat[5]=0.;
    resultat[6]=0.;
    resultat[7]=0.;
    resultat[8]=-2.*res2/res1/res1/res1;
}
 
void evaluer_decompose(double *metrique_depart,double *metrique_decompose)
{
}