sat_straight.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/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       sat_straight.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:54:00 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
 
 
#include "sat_straight.h"
#include "ot_mathematique.h"
#include "constantegeo.h"
 
#include <math.h>
 
SAT_STRAIGHT::SAT_STRAIGHT(unsigned long num):SAT_COURBE(num)
{
}
 
SAT_STRAIGHT::SAT_STRAIGHT():SAT_COURBE()
{
}
 
SAT_STRAIGHT::~SAT_STRAIGHT()
{
}
 
 
void  SAT_STRAIGHT::evaluer(double t,double *xyz)
{
    xyz[0]=root[0]+t*dir[0];
    xyz[1]=root[1]+t*dir[1];
    xyz[2]=root[2]+t*dir[2];
}
 
void  SAT_STRAIGHT::deriver(double t,double *xyz)
{
    xyz[0]=dir[0];
    xyz[1]=dir[1];
    xyz[2]=dir[2];
}
 
void  SAT_STRAIGHT::deriver_seconde(double t,double *ddxyz,double* dxyz,double* xyz)
{
    ddxyz[0]=0.;
    ddxyz[1]=0.;
    ddxyz[2]=0.;
    if (dxyz!=NULL) deriver(t,dxyz);
    if (xyz!=NULL) evaluer(t,xyz);
}
 
void  SAT_STRAIGHT::inverser(double& t,double *xyz,double precision)
{
    int num;
    if (!(OPERATEUR::egal(dir[0],0.0,0.0001))) num=0;
    else if (!(OPERATEUR::egal(dir[1],0.0,0.0001))) num=1;
    else num=2;
 
    t=(xyz[num]-root[num])/dir[num];
}
 
double  SAT_STRAIGHT::get_tmin()
{
    return -1e300;
}
 
double  SAT_STRAIGHT::get_tmax()
{
    return 1e300;
}
 
double SAT_STRAIGHT::get_longueur(double t1,double t2,double precis)
{
    return (sqrt(dir[0]*dir[0]+dir[1]*dir[1]+dir[2]*dir[2])*(t2-t1));
}
 
int SAT_STRAIGHT::est_periodique(void)
{
    return 0;
}
 
double SAT_STRAIGHT::get_periode(void)
{
    return 0.;
}
 
int SAT_STRAIGHT::get_type_geometrique(TPL_LISTE_ENTITE<double> &param)
{
    param.ajouter(root[0]);
    param.ajouter(root[1]);
    param.ajouter(root[2]);
    param.ajouter(dir[0]);
    param.ajouter(dir[1]);
    param.ajouter(dir[2]);
    return GEOMETRIE::CONST::Co_LINE;
}
 
 
 
void SAT_STRAIGHT::get_param_NURBS(int& indx_premier_ptctr,TPL_LISTE_ENTITE<double> &param)
{
    param.ajouter(1);
 
    param.ajouter(2);
    param.ajouter(0);
 
 
    param.ajouter(2);
    param.ajouter(0);
 
 
    param.ajouter(0);
    param.ajouter(0);
    param.ajouter(1);
    param.ajouter(1);
 
    double u_inf=10e6;
    double xyz[3];
 
 
    xyz[0]=root[0]+u_inf*dir[0];
    xyz[1]=root[1]+u_inf*dir[1];
    xyz[2]=root[2]+u_inf*dir[2];
 
    param.ajouter(xyz[0]);
    param.ajouter(xyz[1]);
    param.ajouter(xyz[2]);
    param.ajouter(1.);
 
    xyz[0]=root[0]-u_inf*dir[0];
    xyz[1]=root[1]-u_inf*dir[1];
    xyz[2]=root[2]-u_inf*dir[2];
 
    param.ajouter(xyz[0]);
    param.ajouter(xyz[1]);
    param.ajouter(xyz[2]);
    param.ajouter(1.);
 
    indx_premier_ptctr=9;
}