vct_arete.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_arete.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:55 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
#include "gestionversion.h"
 
 
#pragma hdrstop
#include "vct_arete.h"
#include "vct_courbe.h"
#include "mg_courbe.h"
#include "mg_arete.h"
#include <iomanip>
#include<math.h>
#include "ot_mathematique.h"
 
#pragma package(smart_init)
 
 
 
VCT_ARETE::VCT_ARETE(MG_ARETE* arete):VCT_ELEMENT_TOPOLOGIQUE(arete)
{
    double2 ZERO=0.;
    int indx_premier_ptctr;
    OT_VECTEUR_4DD V1,V2,V;
 
    TPL_LISTE_ENTITE<double> nurbs_params;
    arete->get_param_NURBS(indx_premier_ptctr,nurbs_params);
    int nb_points=1/4.*(nurbs_params.get_nb()-indx_premier_ptctr);
 
    OT_VECTEUR_4DD VCT_NUL(0.,0.,0.,0.);
 
    for (int s=0;s<nb_points-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;
/*
        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.est_nul_3d()) continue;
        V.vecteur_norme_3d();
        if (V==VCT_NUL)
            lst_vecteurs.insert(lst_vecteurs.end(),VCT_NUL);
 
        //if (norm!=ZERO) {
          //  V=1./norm*V;
          else  lst_vecteurs.insert(lst_vecteurs.end(),V);
        
       // }
 
    }
 
    construire_forme_tensorielle();
}
 
VCT_ARETE::VCT_ARETE(VCT_ARETE& mdd):VCT_ELEMENT_TOPOLOGIQUE(mdd.elem_topo)
{
}
 
 
VCT_ARETE::~ VCT_ARETE()
{
 
}