mg_segment.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/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       mg_segment.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:54 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
 
 
#include "gestionversion.h"
#include <math.h>
#include "mg_segment.h"
#include "mg_sommet.h"
#include "ot_decalage_parametre.h"
#include "mg_maillage.h"
 
 
MG_SEGMENT::MG_SEGMENT(unsigned long num,MG_ELEMENT_TOPOLOGIQUE* topo,MG_NOEUD *mgnoeud1,MG_NOEUD *mgnoeud2,int origine,double longue):MG_ELEMENT_MAILLAGE(num,topo,origine),OT_REFERENCE(),noeud1(mgnoeud1),noeud2(mgnoeud2),longueur(longue),dimension(0),distance_metrique_valide(false)
{
    if (longueur==0.) calcule_longueur();
    if (noeud1!=NULL) noeud1->get_lien_segment()->ajouter(this);
    if (noeud2!=NULL) noeud2->get_lien_segment()->ajouter(this);
    if ((noeud1!=NULL) && (noeud2!=NULL))
        if (noeud1->get_id()<noeud2->get_id()) noeud1->get_lien_petit_segment()->ajouter(this);
        else noeud2->get_lien_petit_segment()->ajouter(this);
    noeud1->incrementer();
    noeud2->incrementer();
    if (liaison_topologique==NULL) return;
    if (liaison_topologique->get_dimension()==1) liaison_topologique->get_lien_maillage()->ajouter(this);
}
 
MG_SEGMENT::MG_SEGMENT(MG_ELEMENT_TOPOLOGIQUE* topo,MG_NOEUD *mgnoeud1,MG_NOEUD *mgnoeud2,int origine,double longue):MG_ELEMENT_MAILLAGE(topo,origine),OT_REFERENCE(),noeud1(mgnoeud1),noeud2(mgnoeud2),longueur(longue),dimension(0),distance_metrique_valide(false)
{
    if (longueur==0.) calcule_longueur();
    if (noeud1!=NULL) noeud1->get_lien_segment()->ajouter(this);
    if (noeud2!=NULL) noeud2->get_lien_segment()->ajouter(this);
    if ((noeud1!=NULL) && (noeud2!=NULL))
        if (noeud1->get_id()<noeud2->get_id()) noeud1->get_lien_petit_segment()->ajouter(this);
        else noeud2->get_lien_petit_segment()->ajouter(this);
    noeud1->incrementer();
    noeud2->incrementer();
    if (liaison_topologique==NULL) return;
    if (liaison_topologique->get_dimension()==1) liaison_topologique->get_lien_maillage()->ajouter(this);
}
 
MG_SEGMENT::MG_SEGMENT(MG_SEGMENT& mdd):MG_ELEMENT_MAILLAGE(mdd),OT_REFERENCE(),noeud1(mdd.noeud1),noeud2(mdd.noeud2),longueur(mdd.longueur),dimension(mdd.dimension),distance_metrique_valide(mdd.distance_metrique_valide),distance_metrique(mdd.distance_metrique)
{
    if (noeud1!=NULL) noeud1->get_lien_segment()->ajouter(this);
    if (noeud2!=NULL) noeud2->get_lien_segment()->ajouter(this);
    if ((noeud1!=NULL) && (noeud2!=NULL))
        if (noeud1->get_id()<noeud2->get_id()) noeud1->get_lien_petit_segment()->ajouter(this);
        else noeud2->get_lien_petit_segment()->ajouter(this);
    noeud1->incrementer();
    noeud2->incrementer();
    if (liaison_topologique==NULL) return;
    if (liaison_topologique->get_dimension()==1) liaison_topologique->get_lien_maillage()->ajouter(this);
}
 
 
MG_SEGMENT* MG_SEGMENT::dupliquer(MG_MAILLAGE *mgmai,long decalage)
{
    MG_NOEUD* nd1=mgmai->get_mg_noeudid(noeud1->get_id()+decalage);
    MG_NOEUD* nd2=mgmai->get_mg_noeudid(noeud2->get_id()+decalage);
    if (nd1==NULL)
        noeud1->get_id();
    if (nd2==NULL)
        noeud2->get_id();
    MG_SEGMENT* seg=new MG_SEGMENT(this->get_id()+decalage,liaison_topologique,nd1,nd2,MAGIC::ORIGINE::DUPLIQUER,longueur);
    mgmai->ajouter_mg_segment(seg);
    return seg;
}
 
 
MG_SEGMENT::~MG_SEGMENT()
{
    if (noeud1!=NULL) noeud1->get_lien_segment()->supprimer(this);
    if (noeud2!=NULL) noeud2->get_lien_segment()->supprimer(this);
    if ((noeud1!=NULL) && (noeud2!=NULL))
        if (noeud1->get_id()<noeud2->get_id()) noeud1->get_lien_petit_segment()->supprimer(this);
        else noeud2->get_lien_petit_segment()->supprimer(this);
    noeud1->decrementer();
    noeud2->decrementer();
    if (liaison_topologique!=NULL) 
        if (liaison_topologique->get_dimension()==1) liaison_topologique->get_lien_maillage()->supprimer(this);
}
 
 
int MG_SEGMENT::get_type_entite(void)
{
    return MAGIC::TYPE_ENTITE::IDMG_SEGMENT;
}
 
 
MG_NOEUD* MG_SEGMENT::get_noeud1(void)
{
    return noeud1;
}
 
MG_NOEUD* MG_SEGMENT::get_noeud2(void)
{
    return noeud2;
}
 
void MG_SEGMENT::inverse_sens(void)
{
MG_NOEUD* noeud=noeud1;
noeud1=noeud2;
noeud2=noeud1;
}
 
double MG_SEGMENT::get_longueur(void)
{
    return longueur;
}
 
void MG_SEGMENT::change_noeud1(MG_NOEUD* mgnoeud,double longue)
{
    if (noeud1!=NULL) noeud1->get_lien_segment()->supprimer(this);
    if (noeud1!=NULL)
        if (noeud1->get_id()<noeud2->get_id()) noeud1->get_lien_petit_segment()->supprimer(this);
    noeud1=mgnoeud;
    longueur=longue;
    if (longueur==0.) calcule_longueur();
    if (noeud1!=NULL) noeud1->get_lien_segment()->ajouter(this);
    if (noeud1!=NULL)
        if (noeud1->get_id()<noeud2->get_id()) noeud1->get_lien_petit_segment()->ajouter(this);
}
 
 
void MG_SEGMENT::change_noeud2(MG_NOEUD* mgnoeud,double longue)
{
    if (noeud2!=NULL) noeud2->get_lien_segment()->supprimer(this);
    if (noeud2!=NULL)
        if (noeud2->get_id()<noeud1->get_id()) noeud2->get_lien_petit_segment()->supprimer(this);
    noeud2=mgnoeud;
    longueur=longue;
    if (longueur==0.) calcule_longueur();
    if (noeud2!=NULL) noeud2->get_lien_segment()->ajouter(this);
    if (noeud2!=NULL)
        if (noeud2->get_id()<noeud1->get_id()) noeud2->get_lien_petit_segment()->ajouter(this);
}
 
 
void MG_SEGMENT::enregistrer(std::ostream& o,double version)
{
    if ((liaison_topologique==NULL) && (MG_MAILLAGE::DIMENSIONMAILLAGESANSTOPO!=1)) return;
    if (liaison_topologique==NULL) o << "%" << get_id() << "=SEGMENT(NULL,$" << noeud1->get_id() << ",$" << noeud2->get_id() << "," << origine_creation << ");" << std::endl;
    else if (liaison_topologique->get_dimension()==1) o << "%" << get_id() << "=SEGMENT($"<< get_lien_topologie()->get_id() << ",$" << noeud1->get_id() << ",$" << noeud2->get_id() << "," << origine_creation << ");" << std::endl;
}
 
int MG_SEGMENT::get_dimension(void)
{
    return 1;
}
 
 
void MG_SEGMENT::calcule_longueur(void)
{
    double *coo1=noeud1->get_coord();
    double *coo2=noeud2->get_coord();
    longueur=sqrt( (coo1[0]-coo2[0])*(coo1[0]-coo2[0])+(coo1[1]-coo2[1])*(coo1[1]-coo2[1])+(coo1[2]-coo2[2])*(coo1[2]-coo2[2]));
}
 
 
 
 
BOITE_2D MG_SEGMENT::get_boite_2D(double periodeu,double periodev)
{
    OT_DECALAGE_PARAMETRE decalage(periodeu,periodev);
    double coouv1[2];
    double coouvd1[2];
    double coouv2[2];
    double coouvd2[2];
 
    coouv1[0]=noeud1->get_u();
    coouv1[1]=noeud1->get_v();
    coouv2[0]=noeud2->get_u();
    coouv2[1]=noeud2->get_v();
    double du=decalage.calcul_decalage_parametre_u(coouv1[0]);
    double dv=decalage.calcul_decalage_parametre_v(coouv1[1]);
    coouvd1[0]=decalage.decalage_parametre_u(coouv1[0],du);
    coouvd1[1]=decalage.decalage_parametre_v(coouv1[1],dv);
    coouvd2[0]=decalage.decalage_parametre_u(coouv2[0],du);
    coouvd2[1]=decalage.decalage_parametre_v(coouv2[1],dv);
    double umin,vmin,umax,vmax;
    umin=std::min(coouvd1[0],coouvd2[0]);
    umax=std::max(coouvd1[0],coouvd2[0]);
    vmin=std::min(coouvd1[1],coouvd2[1]);
    vmax=std::max(coouvd1[1],coouvd2[1]);
    umin=umin-du;
    umax=umax-du;
    vmin=vmin-dv;
    vmax=vmax-dv;
    BOITE_2D boite(umin,vmin,umax,vmax);
    return boite;
}
 
BOITE_3D MG_SEGMENT::get_boite_3D(void)
{
    double coo1[3];
    double coo2[3];
 
    coo1[0]=noeud1->get_x();
    coo1[1]=noeud1->get_y();
    coo1[2]=noeud1->get_z();
    coo2[0]=noeud2->get_x();
    coo2[1]=noeud2->get_y();
    coo2[2]=noeud2->get_z();
    double xmin,ymin,zmin,xmax,ymax,zmax;
    xmin=std::min(coo1[0],coo2[0]);
    xmax=std::max(coo1[0],coo2[0]);
    ymin=std::min(coo1[1],coo2[1]);
    ymax=std::max(coo1[1],coo2[1]);
    zmin=std::min(coo1[2],coo2[2]);
    zmax=std::max(coo1[2],coo2[2]);
    BOITE_3D boite(xmin,ymin,zmin,xmax,ymax,zmax);
    return boite;
}
 
void MG_SEGMENT::extrapoler_solution_noeud(void)
{
    for (int j=0;j<MAX_TYPE_SOLUTION;j++)  
      {
        noeud1->change_solution(solution[j],j);
        noeud2->change_solution(solution[j],j);
      }
}
 
TPL_LISTE_ENTITE<class MG_TRIANGLE*>* MG_SEGMENT::get_lien_triangle(void)
{
    return &triangle;
}
TPL_LISTE_ENTITE<class MG_QUADRANGLE*>* MG_SEGMENT::get_lien_quadrangle(void)
{
    return &quadrangle;
}
void MG_SEGMENT::change_dimension_topo_null(int num)
{
dimension=num;
}
 
 
bool MG_SEGMENT::get_distance_metrique(double *l)
{
*l=distance_metrique;
return distance_metrique_valide;  
}
 
void MG_SEGMENT::change_distance_metrique(double val)
{
distance_metrique=val;
distance_metrique_valide=true;
}
 
int MG_SEGMENT::get_dimension_topo_null(void)
{
return dimension;
}