fct_taille_arete.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_arete.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:53 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
 
 
#include "gestionversion.h"
#include "fct_taille_arete.h"
#include <math.h>
 
 
 
 
FCT_TAILLE_ARETE::FCT_TAILLE_ARETE(double eps,double dist,MG_ARETE* mgarete):FCT_TAILLE_METRIQUE(),epsilon(eps),distance_maximale(dist),arete(mgarete)
{
 
    t1=arete->get_tmin();
    t2=arete->get_tmax();
    dens1=distance_maximale;
    dens2=distance_maximale;
    for (int i=0;i<arete->get_cosommet1()->get_sommet()->get_nb_mg_cosommet();i++)
    {
        MG_ARETE* mgarete=arete->get_cosommet1()->get_sommet()->get_mg_cosommet(i)->get_arete();
        double longueur=mgarete->get_longueur(mgarete->get_tmin(),mgarete->get_tmax());
        dens1=std::min(dens1,longueur);
    }
    for (int j=0;j<arete->get_cosommet2()->get_sommet()->get_nb_mg_cosommet();j++)
    {
        MG_ARETE* mgarete=arete->get_cosommet2()->get_sommet()->get_mg_cosommet(j)->get_arete();
        double longueur=mgarete->get_longueur(mgarete->get_tmin(),mgarete->get_tmax());
        dens2=std::min(dens2,longueur);
    }
}
 
 
FCT_TAILLE_ARETE::~FCT_TAILLE_ARETE()
{
}
 
void FCT_TAILLE_ARETE::evaluer(double *param,double *resultat)
{
 
 
 
    double t=*param;
    double tm=t1*dens1+t2*dens2;
    tm=tm/(dens1+dens2);
    double facteur2=(t-t1)*(t-t2)/(tm-t1)/(tm-t2)*distance_maximale+(t-tm)*(t-t2)/(t1-tm)/(t1-t2)*dens1+(t-t1)*(t-tm)/(t2-t1)/(t2-tm)*dens2;
    facteur2=1/facteur2/facteur2;
 
    double M=arete->get_M(*param);
    int nb_face=arete->get_nb_mg_coarete();
    double xyz[3];
    arete->evaluer(t,xyz);
    for (int j=0;j<nb_face;j++)
    {
        MG_FACE* face=arete->get_mg_coarete(j)->get_boucle()->get_mg_face();
        double uv[2];
        face->inverser(uv,xyz);
        double cmin,cmax;
        face->get_courbure(uv,cmax,cmin);
        cmax=fabs(cmax);
        if (cmax>M) M=cmax;
    }
    double facteur1=M/8./epsilon;
    double valeur=std::max(facteur1,facteur2);
    resultat[0]=valeur;
    resultat[1]=0.;
    resultat[2]=0.;
    resultat[3]=0.;
    resultat[4]=valeur;
    resultat[5]=0.;
    resultat[6]=0.;
    resultat[7]=0.;
    resultat[8]=valeur;
}
 
 
void  FCT_TAILLE_ARETE::deriver(double *param,double *resultat,int num_param)
{
    double t1=param[0];
    double t2=param[0]+1e-6;
    double coo1[9],coo2[9];
    evaluer(&t1,coo1);
    evaluer(&t2,coo2);
    resultat[0]=(coo2[0]-coo1[0])/1e-6;
    resultat[1]=(coo2[1]-coo1[1])/1e-6;
    resultat[2]=(coo2[2]-coo1[2])/1e-6;
    resultat[3]=(coo2[3]-coo1[3])/1e-6;
    resultat[4]=(coo2[4]-coo1[4])/1e-6;
    resultat[5]=(coo2[5]-coo1[5])/1e-6;
    resultat[6]=(coo2[6]-coo1[6])/1e-6;
    resultat[7]=(coo2[7]-coo1[7])/1e-6;
    resultat[8]=(coo2[8]-coo1[8])/1e-6;
}
 
void FCT_TAILLE_ARETE::evaluer_decompose(double *metrique_depart,double *metrique_decompose)
{
    metrique_decompose[0]=metrique_depart[0];
    arete->deriver(metrique_depart[1],metrique_decompose+1);
}
 
int FCT_TAILLE_ARETE::valide_parametre(double *param)
{
    return 1;
}
 
 
double FCT_TAILLE_ARETE::get_valeur_maximale(int num)
{
    return distance_maximale;
}