vct_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/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       vct_courbe.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:54 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
#include "gestionversion.h"
 
 
 
#pragma hdrstop
#include<math.h>
#include "vct_courbe.h"
#include "mg_courbe.h"
#include "mg_arete.h"
#include <iomanip>
 
#pragma package(smart_init)
 
 
 
VCT_COURBE::VCT_COURBE(MG_COURBE* elemgeo):VCT_ELEMENT_GEOMETRIQUE(elemgeo)
{
    int indx_premier_ptctr;
    OT_VECTEUR_4DD V1,V2,V;
    TPL_LISTE_ENTITE<double> nurbs_params;
    elem_geo->get_param_NURBS(indx_premier_ptctr,nurbs_params);
    int nb_points=1/4.*(nurbs_params.get_nb()-indx_premier_ptctr);
    double2 ZERO=0.;
    OT_VECTEUR_4DD VCT_NUL(0.,0.,0.,0.);
 
 
    int nb_u=nurbs_params.get(3);
 
 
    for (int s=0;s<nb_u-1;s++)
    {
        V1[0]=nurbs_params.get(indx_premier_ptctr+  4 * s);
        V1[1]=nurbs_params.get(indx_premier_ptctr+  4 * s + 1);
        V1[2]=nurbs_params.get(indx_premier_ptctr+  4 * s + 2);
        V1[3]=nurbs_params.get(indx_premier_ptctr+  4 * s + 3);
        if (s==0)
            lst_points.insert(lst_points.end(),V1);
 
        V2[0]=nurbs_params.get(indx_premier_ptctr+  4 * (s + 1));
        V2[1]=nurbs_params.get(indx_premier_ptctr+  4 * (s + 1) + 1);
        V2[2]=nurbs_params.get(indx_premier_ptctr+  4 * (s + 1) + 2);
        V2[3]=nurbs_params.get(indx_premier_ptctr+  4 * (s + 1) + 3);
 
        lst_points.insert(lst_points.end(),V2);
 
        V=V2-V1;
        if (V.est_nul_3d()) continue;
        V.vecteur_norme_3d();
/*
        double2 norm_au_carre=(V[0]*V[0])+(V[1]*V[1])+(V[2]*V[2])+(V[3]*V[3]);
        double2 norm=norm_au_carre^0.5;
*/
        if (V==VCT_NUL)
            lst_vecteurs.insert(lst_vecteurs.end(),VCT_NUL);
 
        //if (norm!=ZERO) //{
          else  lst_vecteurs.insert(lst_vecteurs.end(),V);
        //}
 
 
    }
    construire_forme_tensorielle();
}
 
 
VCT_COURBE::VCT_COURBE(VCT_COURBE& mdd):VCT_ELEMENT_GEOMETRIQUE(mdd.elem_geo)
{
}
 
 
VCT_COURBE::~VCT_COURBE()
{
}