mg_courbe.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/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       mg_courbe.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:54 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
 
 
#include "gestionversion.h"
#include <math.h>
#include "mg_courbe.h"
#include "vct_courbe.h"
 
 
 
MG_COURBE::MG_COURBE(unsigned long num):MG_ELEMENT_GEOMETRIQUE(num),vect(NULL)
{
}
 
 
MG_COURBE::MG_COURBE():MG_ELEMENT_GEOMETRIQUE(),vect(NULL)
{
}
 
 
 
MG_COURBE::MG_COURBE(MG_COURBE& mdd):MG_ELEMENT_GEOMETRIQUE(),vect(NULL)
{
}
 
 
MG_COURBE::~MG_COURBE()
{
    if (vect!=NULL) delete vect;
}
 
int MG_COURBE::get_type(void)
{
  return TYPE_ELEMENT_GEOMETRIQUE::COURBE;
}
 
 
double  MG_COURBE::get_tmin()
{
    return t_min;
}
double  MG_COURBE::get_tmax()
{
    return t_max;
}
 
double MG_COURBE::get_M(double t)
{
    double ddxyz[3];
    double dxyz[3];
    double xyz[3];
 
    deriver_seconde(t,ddxyz,dxyz,xyz);
    double facteur1=(dxyz[0]*dxyz[0]+dxyz[1]*dxyz[1]+dxyz[2]*dxyz[2]);
    double facteur2=(ddxyz[0]*dxyz[0]+ddxyz[1]*dxyz[1]+ddxyz[2]*dxyz[2])/facteur1/facteur1;
    double m[3];
    m[0]=ddxyz[0]/facteur1-dxyz[0]*facteur2;
    m[1]=ddxyz[1]/facteur1-dxyz[1]*facteur2;
    m[2]=ddxyz[2]/facteur1-dxyz[2]*facteur2;
    return sqrt(m[0]*m[0]+m[1]*m[1]+m[2]*m[2]);
 
}
 
VCT& MG_COURBE::get_vectorisation(void)
{
    if (vect==NULL) vect=new VCT_COURBE(this);
    return *vect;
 
}