ct_point.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/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       ct_point.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:53 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
 
 
#include "gestionversion.h"
 
#pragma hdrstop
 
#include "ct_point.h"
#include "ot_boite_3d.h"
 
template <int N> unsigned long CT_POINT<N>::idmax=0;   //definir ds CT_POINT la var idmax
 
template <int N> CT_POINT<N>::CT_POINT(double x,double y,double z,DOUBLEN<N> en,BOITE_3D boitetmp, double courbure):boite(boitetmp)
{
    point[0]=x;
    point[1]=y;
    point[2]=z;
    ecart=en;
    id=idmax;
    idmax++;
    c=courbure;
}
 
 
template <int N> CT_POINT<N>::~CT_POINT()
{
}
 
template <int N> void CT_POINT<N>::evaluer(double *xyz)
{
    xyz[0]=point[0];
    xyz[1]=point[1];
    xyz[2]=point[2];
}
 
template <int N> unsigned long CT_POINT<N>::get_id()
{
    return id;
}
 
template <int N> DOUBLEN<N> &CT_POINT<N>::get_valeur(void)
{
    return ecart;
}
 
template <int N> double CT_POINT<N>::get_fonction_influence(double r)
{
 
    if (!c) return (1-r);
    double f,b=0.1;              // la fonction est soit une droite (c=0) ou une courbe de courbure positive(f>0) ou n�gative (f<0) nulle en r=1
    if (c>0)                     // cette fonction d�termine le poids d'un point d'�chantillonnage en fonction de sa distance
    {                           // les deux fonctions c>0 et c<0 engendrent des �volutions oppos�es l'une � l'autre, et tendent toutes les deux vers 1-r quand c tend vers 0.
        f=(c+b)*b/(c*(c*r+b))-b/c;
        return f;
    }
    if (c<0)
    {
        f=(b-c)*b/((-c)*((-c)*(1-r)+b))-b/(-c);
        return r*(1-r)/f;
    }
 
 
    return 0;
}
 
template <int N> double CT_POINT<N>::get_fonction_influence(double r, double crbr)
{
 
    if (!crbr) return (1-r);
    double f,b=0.1;              // la fonction est soit une droite (c=0) ou une courbe de courbure positive(f>0) ou n�gative (f<0) nulle en r=1
    if (crbr>0)                     // cette fonction d�termine le poids d'un point d'�chantillonnage en fonction de sa distance
    {                           // les deux fonctions c>0 et c<0 engendrent des �volutions oppos�es l'une � l'autre, et tendent toutes les deux vers 1-r quand c tend vers 0.
        f=(crbr+b)*b/(crbr*(crbr*r+b))-b/crbr;
        return f;
    }
    if (crbr<0)
    {
        f=(b-crbr)*b/((-crbr)*((-crbr)*(1-r)+b))-b/(-crbr);
        return r*(1-r)/f;
    }
 
 
    return 0;
}
 
 
 
template <int N> BOITE_3D CT_POINT<N>::get_boite_3D(void)
{
    return boite;
}
 
 
 
#ifdef BORLANDCPP
template class __export CT_POINT<1>;
template class __export CT_POINT<4>;
#else
template class CT_POINT<1>;
template class CT_POINT<4>;
#endif