fem_graph_noeud.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/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       fem_graph_noeud.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:53 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
 
 
#pragma hdrstop
#include "gestionversion.h"
#include "fem_graph_noeud.h"
#include "fem_noeud.h"
#include "fem_element3.h"
#include "fem_maillage.h"
 
 
 
FEM_GRAPH_NOEUD::FEM_GRAPH_NOEUD(FEM_NOEUD* rac,FEM_MAILLAGE* mai):racine(rac),maillage(mai)
{
    construire();
}
 
FEM_GRAPH_NOEUD::~FEM_GRAPH_NOEUD()
{
}
 
 
void FEM_GRAPH_NOEUD::construire(void)
{
    int nbnoeudtotal=maillage->get_nb_fem_noeud();
    int *tab=new int[nbnoeudtotal+1];
    for (int i=0;i<nbnoeudtotal+1;i++) tab[i]=0;
    tab[racine->get_numero()]=1;
    TPL_MAP_ENTITE<FEM_NOEUD*> liste;
    graph.insert(graph.end(),liste);
    int niveau=0;
    graph[0].ajouter(racine);
    while (graph[niveau].get_nb()!=0)
    {
        graph.insert(graph.end(),liste);
        int nb=graph[niveau].get_nb();
        for (int i=0;i<nb;i++)
        {
            FEM_NOEUD* noeud=graph[niveau].get(i);
            TPL_MAP_ENTITE<FEM_NOEUD*> liste_voisin;
            get_voisin_noeud(noeud,liste_voisin);
            int nbvoisin=liste_voisin.get_nb();
            for (int j=0;j<nbvoisin;j++)
            {
                FEM_NOEUD* no=liste_voisin.get(j);
                if (tab[no->get_numero()]!=1)
                {
                    tab[no->get_numero()]=1;
                    graph[niveau+1].ajouter(no);
                }
            }
        }
        niveau++;
    }
    delete [] tab;
}
 
int FEM_GRAPH_NOEUD::get_excentricite(void)
{
    return graph.size()-2;
}
 
TPL_MAP_ENTITE<FEM_NOEUD*>& FEM_GRAPH_NOEUD::get_dernier_niveau(void)
{
    int num=graph.size()-2;
    return graph[num];
}
 
TPL_MAP_ENTITE<FEM_NOEUD*>& FEM_GRAPH_NOEUD::get_niveau(int num)
{
    return graph[num];
}
 
void FEM_GRAPH_NOEUD::get_voisin_noeud(FEM_NOEUD* no,TPL_MAP_ENTITE<FEM_NOEUD*> &liste)
{
    int dim;
    int nb_ele=no->get_lien_element3()->get_nb();
    if (nb_ele!=0) dim=3;
    else 
        {
          nb_ele=no->get_lien_element2()->get_nb();
          if (nb_ele!=0) dim=2;
          else 
            {
            nb_ele=no->get_lien_element1()->get_nb();
            dim=1;  
            }
        }
    for (int i=0;i<nb_ele;i++)
    {
        TPL_LISTE_ENTITE<FEM_NOEUD*> voisin;
        if (dim==3)
            {
            FEM_ELEMENT3* ele=no->get_lien_element3()->get(i);
            ele->get_voisin_noeud(no,voisin);  
            }
        if (dim==2)
            {
            FEM_ELEMENT2* ele=no->get_lien_element2()->get(i);
            ele->get_voisin_noeud(no,voisin);  
            }
        if (dim==1)
            {
            FEM_ELEMENT1* ele=no->get_lien_element1()->get(i);
            ele->get_voisin_noeud(no,voisin);  
            }
        for (int i=0;i<voisin.get_nb();i++)
          liste.ajouter(voisin.get(i));
      
        }
    
}
 
#pragma package(smart_init)