sat_plane.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_plane.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:54:00 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
 
 
#include "sat_plane.h"
#include "constantegeo.h"
#include "ot_mathematique.h"
 
 
SAT_PLANE::SAT_PLANE(unsigned long num):SAT_SURFACE(num)
{
}
 
SAT_PLANE::SAT_PLANE():SAT_SURFACE()
{
}
 
SAT_PLANE::~SAT_PLANE()
{
}
 
 
 
void SAT_PLANE::calcule_parametre(void)
{
    OT_VECTEUR_3D vec_normal(normal);
    vec_normal.norme();
    if (!(OPERATEUR::egal(vec_normal.get_x(),0.,0.0001)))
    {
        dir1[0]=vec_normal.get_y();
        dir1[1]=(-vec_normal.get_x());
        dir1[2]=0.;
    }
    else if (!(OPERATEUR::egal(vec_normal.get_y(),0.,0.0001)))
    {
        dir1[0]=0.;
        dir1[1]=(-vec_normal.get_z());
        dir1[2]=vec_normal.get_y();
    }
    else
    {
        dir1[0]=vec_normal.get_z();
        dir1[1]=0.;
        dir1[2]=0.;
    }
    OT_VECTEUR_3D vec_dir1(dir1);
    vec_dir1.norme();
    OT_VECTEUR_3D vec_dir2=vec_normal&vec_dir1;
    dir2[0]=vec_dir2.get_x();
    dir2[1]=vec_dir2.get_y();
    dir2[2]=vec_dir2.get_z();
}
 
 
void  SAT_PLANE::evaluer(double *uv,double *xyz)
{
    xyz[0]=root[0]+uv[0]*dir1[0]+uv[1]*dir2[0];
    xyz[1]=root[1]+uv[0]*dir1[1]+uv[1]*dir2[1];
    xyz[2]=root[2]+uv[0]*dir1[2]+uv[1]*dir2[2];
}
 
void  SAT_PLANE::deriver(double *uv,double *xyzdu, double *xyzdv)
{
    xyzdu[0]=dir1[0];
    xyzdu[1]=dir1[1];
    xyzdu[2]=dir1[2];
    xyzdv[0]=dir2[0];
    xyzdv[1]=dir2[1];
    xyzdv[2]=dir2[2];
}
 
void  SAT_PLANE::deriver_seconde(double *uv,double* xyzduu,double* xyzduv,double* xyzdvv,double *xyz, double *xyzdu , double *xyzdv )
{
    xyzduu[0]=0.;
    xyzduu[1]=0.;
    xyzduu[2]=0.;
    xyzduv[0]=0.;
    xyzduv[1]=0.;
    xyzduv[2]=0.;
    xyzdvv[0]=0.;
    xyzdvv[1]=0.;
    xyzdvv[2]=0.;
    evaluer(uv,xyz);
    deriver(uv,xyzdu,xyzdv);
}
 
void  SAT_PLANE::inverser(double *uv,double *xyz,double precision)
{
    double coord[3];
    int n1,n2;
    coord[0]=xyz[0]-root[0];
    coord[1]=xyz[1]-root[1];
    coord[2]=xyz[2]-root[2];;
    double det=dir1[0]*dir2[1]-dir1[1]*dir2[0];
    if (!(OPERATEUR::egal(det,0.0,0.0001)))
    {
        n1=0;
        n2=1;
    }
    else
    {
        det=dir1[0]*dir2[2]-dir1[2]*dir2[0];
        if (!(OPERATEUR::egal(det,0.0,0.0001)))
        {
            n1=0;
            n2=2;
        }
        else
        {
            n1=1;
            n2=2;
            det=dir1[1]*dir2[2]-dir1[2]*dir2[1];
        }
    }
    uv[0]=(coord[n1]*dir2[n2]-coord[n2]*dir2[n1])/det;
    uv[1]=(dir1[n1]*coord[n2]-dir1[n2]*coord[n1])/det;
 
}
 
int SAT_PLANE::est_periodique_u(void)
{
    return 0;
}
 
int SAT_PLANE::est_periodique_v(void)
{
    return 0;
}
 
double SAT_PLANE::get_periode_u(void)
{
    return 0.;
}
 
double SAT_PLANE::get_periode_v(void)
{
    return 0.;
}
 
double  SAT_PLANE::get_umin(void)
{
    return -1e300;
}
 
double  SAT_PLANE::get_umax(void)
{
    return 1e300;
}
 
double  SAT_PLANE::get_vmin(void)
{
    return -1e300;
}
 
double  SAT_PLANE::get_vmax(void)
{
    return 1e300;
}
 
 
int SAT_PLANE::get_type_geometrique(TPL_LISTE_ENTITE<double> &param)
{
    param.ajouter(root[0]);
    param.ajouter(root[1]);
    param.ajouter(root[2]);
    param.ajouter(normal[0]);
    param.ajouter(normal[1]);
    param.ajouter(normal[2]);
    return GEOMETRIE::CONST::Co_PLAN;
}
 
 
 
void SAT_PLANE::get_param_NURBS(int& indx_premier_ptctr,TPL_LISTE_ENTITE<double> &param)
{
 
    double uv[2];
    double xyz[3];
 
    param.ajouter(2);
 
    param.ajouter(2);
    param.ajouter(2);
 
 
    param.ajouter(2);
    param.ajouter(2);
 
 
    param.ajouter(0);
    param.ajouter(0);
    param.ajouter(1);
    param.ajouter(1);
 
 
    param.ajouter(0);
    param.ajouter(0);
    param.ajouter(1);
    param.ajouter(1);
 
 
 
    double inf_val=10e6;
 
 
    xyz[0]=root[0]+inf_val*dir1[0]+inf_val*dir2[0];
    xyz[1]=root[1]+inf_val*dir1[1]+inf_val*dir2[1];
    xyz[2]=root[2]+inf_val*dir1[2]+inf_val*dir2[2];
 
    param.ajouter(xyz[0]);
    param.ajouter(xyz[1]);
    param.ajouter(xyz[2]);
    param.ajouter(1.);
 
 
 
    xyz[0]=root[0]-inf_val*dir1[0]+inf_val*dir2[0];
    xyz[1]=root[1]-inf_val*dir1[1]+inf_val*dir2[1];
    xyz[2]=root[2]-inf_val*dir1[2]+inf_val*dir2[2];
 
    param.ajouter(xyz[0]);
    param.ajouter(xyz[1]);
    param.ajouter(xyz[2]);
    param.ajouter(1.);
 
 
 
    xyz[0]=root[0]+inf_val*dir1[0]-inf_val*dir2[0];
    xyz[1]=root[1]+inf_val*dir1[1]-inf_val*dir2[1];
    xyz[2]=root[2]+inf_val*dir1[2]-inf_val*dir2[2];
 
    param.ajouter(xyz[0]);
    param.ajouter(xyz[1]);
    param.ajouter(xyz[2]);
    param.ajouter(1.);
 
 
    xyz[0]=root[0]-inf_val*dir1[0]-inf_val*dir2[0];
    xyz[1]=root[1]-inf_val*dir1[1]-inf_val*dir2[1];
    xyz[2]=root[2]-inf_val*dir1[2]-inf_val*dir2[2];
 
    param.ajouter(xyz[0]);
    param.ajouter(xyz[1]);
    param.ajouter(xyz[2]);
    param.ajouter(1.);
 
    indx_premier_ptctr=13;
 
}