pir_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/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       pir_courbe.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:53 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
 
 
#include "gestionversion.h"
 
#ifdef BREP_PIRATE
 
#include "pir_courbe.h"
 
 
 
 
 
PIR_COURBE::PIR_COURBE(unsigned long num,PG_COURBE *crb):MG_COURBE(num),courbe(crb)
{
    pos=new PG_POSITION1D_COURBE((*courbe));
    t_min=pos->get_tmin();
    t_max=pos->get_tmax();
}
 
PIR_COURBE::PIR_COURBE(PG_COURBE *crb):MG_COURBE(),courbe(crb)
{
    pos=new PG_POSITION1D_COURBE((*courbe));
    t_min=pos->get_tmin();
    t_max=pos->get_tmax();
}
 
PIR_COURBE::PIR_COURBE(PIR_COURBE& mdd):MG_COURBE(mdd),courbe(mdd.courbe)
{
    pos=new PG_POSITION1D_COURBE((*courbe));
    t_min=pos->get_tmin();
    t_max=pos->get_tmax();
}
 
 
 
PIR_COURBE::~PIR_COURBE()
{
    delete pos;
}
 
void  PIR_COURBE::evaluer(double t,double *xyz)
{
    courbe->evaluer(&t,xyz);
}
 
void  PIR_COURBE::deriver(double t,double *xyz)
{
    courbe->deriver(&t,xyz);
}
 
void  PIR_COURBE::deriver_seconde(double t,double *ddxyz,double *dxyz,double *xyz)
{
    double pt[3]={0,0,0};
    double dt[3]={0,0,0};
    double dtt[3]={0,0,0};
    double* d2C[]={pt, dt, dtt};
    courbe->deriver_seconde(&t,d2C);
    if (xyz!=NULL)
    {
        xyz[0]=d2C[0][0];
        xyz[1]=d2C[0][1];
        xyz[2]=d2C[0][2];
    }
    if (dxyz!=NULL)
    {
        dxyz[0]=d2C[1][0];
        dxyz[1]=d2C[1][1];
        dxyz[2]=d2C[1][2];
    }
    ddxyz[0]=d2C[2][0];
    ddxyz[1]=d2C[2][1];
    ddxyz[2]=d2C[2][2];
}
 
void  PIR_COURBE::inverser(double& t,double *xyz,double precision)
{
    if (precision>0.0001*pos->get_longueur()) pos->localiser(xyz,&t,0.0001*pos->get_longueur());
    else pos->localiser(xyz,&t,precision);
}
 
 
 
 
double  PIR_COURBE::get_longueur(double t1,double t2,double precis)
{
    int niveau;
    if ( (t1==t_min) && (t2==t_max) ) return pos->get_longueur();
    return pos->calculer_longueur_segment(t1,t2,pos->evaluer_integrant_longueur(t1),pos->evaluer_integrant_longueur(t2),precis,niveau);
}
 
int PIR_COURBE::est_periodique(void)
{
    return courbe->estfermee();
}
 
double PIR_COURBE::get_periode(void)
{
    return courbe->get_periode();
}
 
void PIR_COURBE::enregistrer(std::ostream& o,double version)
{
    o << "%" << get_id() << "=COURBE_PIRATE("<< courbe->get_id()<< ");" << std::endl;
}
 
#endif