pir_import.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/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       pir_import.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:53 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
 
 
#include "gestionversion.h"
#ifdef BREP_PIRATE
 
 
#include "affiche.h"
#include "pir_import.h"
#include "magic.h"
#include "message.h"
#include "lc_point.h"
#include "mg_identificateur.h"
#include "mg_point.h"
#include "mg_courbe.h"
#include "mg_surface.h"
#include "mg_sommet.h"
#include "mg_cosommet.h"
#include "mg_arete.h"
#include "mg_coarete.h"
#include "mg_boucle.h"
#include "mg_coface.h"
#include "mg_face.h"
#include "mg_coquille.h"
#include "mg_volume.h"
#include "mg_geometrie.h"
#include "mg_gestionnaire.h"
#include "pgpoint.h"
#include "pgcourbe.h"
#include "pgface.h"
#include "pgboucle.h"
#include "pgsurfac.h"
#include "pggeo.h"
#include "pgarete.h"
#include "pgvolume.h"
#include "pgfacass.h"
#include "pgareass.h"
#include "pgcoquil.h"
#include "pgsomass.h"
#include "pgsommet.h"
#include "pgpos0d.h"
#include "pgpos2d.h"
#include "pifich.h"
#include "pir_point.h"
#include "pir_courbe.h"
#include "pir_surface.h"
#include "pmmai.h"
#include "ptlagr.h"
#include "pmzone.h"
 
 
 
 
 
 
 
PIR_IMPORT::PIR_IMPORT()
{
}
 
 
MG_GEOMETRIE* PIR_IMPORT::importer(MG_GESTIONNAIRE& gest,PI_FICHIER_IMPORT& fichier,char *geo,int *volume)
{
    if (fichier.analyser_syntaxe()!=0) return NULL;
    if (geo!=NULL)
    {
        P__IDENTIFICATEUR idgeo(geo);
        PG_GEOMETRIE *pggeo=fichier.trouver_geometrie(idgeo);
        if (pggeo==NULL)
        {
            afficheur << ERRGEOMETRIE << enderr;
            return NULL;
        }
        MG_GEOMETRIE *mggeo=new MG_GEOMETRIE("PIRATE",fichier.get_id(),pggeo->get_id());
        gest.ajouter_mg_geometrie(mggeo);
        if (volume!=NULL)
        {
            PG_VOLUME *pgvol=pggeo->trouver_volume(*volume);
            importer_volume(mggeo,pggeo,pgvol);
        }
        else importer_volume(mggeo,pggeo);
        return mggeo;
    }
    return NULL;
}
 
 
 
 
PIR_IMPORT::~PIR_IMPORT()
{
}
 
 
 
void PIR_IMPORT::importer_volume(MG_GEOMETRIE* mggeo,PG_GEOMETRIE* pggeo,PG_VOLUME* pgvoldesire)
{
    for (int nbvol=0;nbvol<pggeo->get_nbvolumes();nbvol++)
    {
        PG_VOLUME* pgvol=pggeo->get_volume(nbvol);
        if ((pgvoldesire==NULL)||((pgvoldesire!=NULL)&&(pgvol==pgvoldesire)))
        {
            MG_VOLUME* mgvol=new MG_VOLUME;
            mggeo->ajouter_mg_volume(mgvol);
            for (int ii=0;ii<pgvol->get_nbcoquilles();ii++)
            {
                const PG_COQUILLE* pgcoq=pgvol->get_coquille(ii);
                MG_COQUILLE* mgcoq=new MG_COQUILLE(mgvol);
                mggeo->ajouter_mg_coquille(mgcoq,mgvol);
                for (int jj=0;jj<pgcoq->get_nbfacesassociees();jj++)
                {
                    const PG_FACE_ASSOCIEE* pgfaceass=pgcoq->get_faceassociee(jj);
                    const PG_FACE& pgface=pgfaceass->get_face();
                    const PG_SURFACE& pgsurface=pgface.get_position()->get_surface();
                    MG_SURFACE* mgsrf=trouver((PG_SURFACE*)&pgsurface);
                    if (mgsrf==NULL)
                    {
                        PIR_SURFACE* pirsrf=new PIR_SURFACE((PG_SURFACE*)&pgsurface);
                        mgsrf=pirsrf;
                        mggeo->ajouter_mg_surface(mgsrf);
                        relation((PG_SURFACE*)&pgsurface,mgsrf);
                    }
                    MG_FACE* mgface=trouver((PG_FACE*)&pgface);
                    if (mgface==NULL)
                    {
                        mgface=new MG_FACE(mgsrf,MEME_SENS);
                        mggeo->ajouter_mg_face(mgface);
                        relation((PG_FACE*)&pgface,mgface);
                    }
                    int sens;
                    if (pgfaceass->get_orientation()==1) sens=1;
                    else sens=-1;
                    MG_COFACE *mgcoface=mggeo->ajouter_mg_coface(mgface,mgcoq,sens);
                    mgcoq->ajouter_mg_coface(mgcoface);
                    for (int i=0;i<pgface.get_nbboucles();i++)
                    {
                        const PG_BOUCLE* pgboucle=pgface.get_boucle(i);
                        MG_BOUCLE* mgboucle=new MG_BOUCLE(mgface);
                        mggeo->ajouter_mg_boucle(mgboucle,mgface);
                        for (int j=0;j<pgboucle->get_nbaretesassociees();j++)
                        {
                            const PG_ARETE_ASSOCIEE* pgareteass=pgboucle->get_areteassociee(j);
                            const PG_ARETE& pgarete=pgareteass->get_arete();
                            const PG_POSITION1D *pos1d=pgarete.get_position();
                            const PG_COURBE& pcrb=pos1d->get_courbe();
                            MG_COURBE* mgcrb=trouver((PG_COURBE*)&pcrb);
                            if (mgcrb==NULL)
                            {
                                PIR_COURBE *pircrb=new PIR_COURBE((PG_COURBE* )&pcrb);
                                mgcrb=pircrb;
                                mggeo->ajouter_mg_courbe(mgcrb);
                                relation((PG_COURBE*)&pcrb,mgcrb);
                            }
                            const PG_SOMMET& psom1=pgarete.get_sommetassocie1()->get_sommet();
                            const PG_SOMMET& psom2=pgarete.get_sommetassocie2()->get_sommet();
                            double coo[3];
                            MG_SOMMET* mgsom=trouver((PG_SOMMET*)&psom1);
                            if (mgsom==NULL)
                            {
                                psom1.get_position()->get_coord(coo);
                                LC_POINT* lcpt=new LC_POINT(coo);
                                MG_POINT* mgpt=lcpt;
                                mggeo->ajouter_mg_point(mgpt);
                                mgsom=new MG_SOMMET(mgpt);
                                mggeo->ajouter_mg_sommet(mgsom);
                                relation((PG_SOMMET*)&psom1,mgsom);
                            }
                            MG_SOMMET* mgsom2=trouver((PG_SOMMET*)&psom2);
                            if (mgsom2==NULL)
                            {
                                psom2.get_position()->get_coord(coo);
                                LC_POINT* lcpt=new LC_POINT(coo);
                                MG_POINT* mgpt2=lcpt;
                                mggeo->ajouter_mg_point(mgpt2);
                                mgsom2=new MG_SOMMET(mgpt2);
                                mggeo->ajouter_mg_sommet(mgsom2);
                                relation((PG_SOMMET*)&psom2,mgsom2);
                            }
                            MG_ARETE* mgarete=trouver((PG_ARETE*)&pgarete);
                            if (mgarete==NULL)
                            {
                                mgarete=new MG_ARETE(mgcrb,MEME_SENS);
                                mggeo->ajouter_mg_arete(mgarete);
                                relation((PG_ARETE*)&pgarete,mgarete);
                                MG_COSOMMET* mgcosom1=mggeo->ajouter_mg_cosommet(mgarete,mgsom);
                                MG_COSOMMET* mgcosom2=mggeo->ajouter_mg_cosommet(mgarete,mgsom2);
                                mgarete->changer_cosommet1(mgcosom1);
                                mgarete->changer_cosommet2(mgcosom2);
                            }
                            int sens;
                            if (pgareteass->get_orientation_arete()==1) sens=1;
                            else sens=-1;
                            MG_COARETE* mgcoarete=mggeo->ajouter_mg_coarete(mgarete,mgboucle,sens);
                            mgboucle->ajouter_mg_coarete(mgcoarete);
                        }
                    }
                }
            }
        }
    }
 
 
}
 
 
void PIR_IMPORT::transferer_maillage_arete(MG_GEOMETRIE& mggeo,PM_MAILLAGE& mai,int *arete)
{
    MG_ARETE* mgarete=NULL;
    if (arete!=NULL) mgarete=mggeo.get_mg_areteid(*arete);
    for (int i=0;i<mggeo.get_nb_mg_arete();i++)
    {
        MG_ARETE* mgarete_courante=mggeo.get_mg_arete(i);
        if ((mgarete==NULL) || (mgarete_courante==mgarete))
        {
            PG_ARETE* pgarete=trouver(mgarete_courante);
            P__IDENTIFICATEUR nomZone("dicretisationDeLArete");
            nomZone += pgarete->get_id();
            P_TMPL_CHAMP_MEMOIRE_REF<int> cnc;
            PM_ZONE* pmzone=mai.ajouter_zone( nomZone,*PT_INTERPOLANT_LAGRANGE::getLagrLigne02(), &cnc, pgarete->get_id() );
            TPL_MAP_ENTITE<MG_ELEMENT_MAILLAGE*> * lien_maillage = mgarete_courante->get_lien_maillage();
            TPL_MAP_ENTITE<MG_ELEMENT_MAILLAGE*>::ITERATEUR it;
            MG_ELEMENT_MAILLAGE* element;
            for (element = lien_maillage->get_premier(it); element; element = lien_maillage->get_suivant(it) )
            {
                MG_SEGMENT* mgsegment=(MG_SEGMENT*)element;
                MG_NOEUD* noeud1=mgsegment->get_noeud1();
                MG_NOEUD* noeud2=mgsegment->get_noeud2();
                int cncDeElm[2]={noeud1->get_nouveau_numero(),noeud2->get_nouveau_numero()};
                pmzone->ajouter_cncelement( cncDeElm );
            }
 
        }
    }
}
 
void PIR_IMPORT::transferer_maillage_face(MG_GEOMETRIE& mggeo,PM_MAILLAGE& mai,int *face)
{
    MG_FACE* mgface=NULL;
    if (face!=NULL) mgface=mggeo.get_mg_faceid(*face);
    for (int i=0;i<mggeo.get_nb_mg_face();i++)
    {
        MG_FACE* mgface_courante=mggeo.get_mg_face(i);
        if ((mgface==NULL) || (mgface_courante==mgface))
        {
            PG_FACE* pgface=trouver(mgface_courante);
            P__IDENTIFICATEUR nomZone("dicretisationDeLaFace");
            nomZone += pgface->get_id();
            P_TMPL_CHAMP_MEMOIRE_REF<int> cnc;
            PM_ZONE* pmzone=mai.ajouter_zone( nomZone,*PT_INTERPOLANT_LAGRANGE::getLagrTrian03(), &cnc, pgface->get_id() );
            TPL_MAP_ENTITE<MG_ELEMENT_MAILLAGE*> * lien_maillage = mface_courante->get_lien_maillage();
            TPL_MAP_ENTITE<MG_ELEMENT_MAILLAGE*>::ITERATEUR it;
            MG_ELEMENT_MAILLAGE* element;
            for (element = lien_maillage->get_premier(it); element; element = lien_maillage->get_suivant(it) )
            {
                MG_TRIANGLE* mgtri=(MG_TRIANGLE*)element;
                MG_NOEUD* noeud1=mgtri->get_noeud1();
                MG_NOEUD* noeud2=mgtri->get_noeud2();
                MG_NOEUD* noeud3=mgtri->get_noeud3();
                int cncDeElm[3]={noeud1->get_nouveau_numero(),noeud2->get_nouveau_numero(),noeud3->get_nouveau_numero()};
                pmzone->ajouter_cncelement( cncDeElm );
            }
 
        }
    }
}
 
 
void PIR_IMPORT::transferer_noeud(MG_MAILLAGE& mgmai,PM_MAILLAGE& pmmai)
{
    P_TMPL_CHAMP_MEMOIRE_REF<double> coord("", P__CHAMP_REF_DATA(3));
    P_TMPL_CHAMP_MEMOIRE_REF<int> codes("", P__CHAMP_REF_DATA(1));
    for (int i=0;i<mgmai.get_nb_mg_noeud();i++)
    {
        MG_NOEUD* mgnoeud=mgmai.get_mg_noeud(i);
        double coordDuNoeud[3];
        mgnoeud->get_coord(coordDuNoeud);
        coord.ajouter_groupes(coordDuNoeud);
        MG_ELEMENT_TOPOLOGIQUE* mg_topo=mgnoeud->get_lien_topologie();
        int dim=mg_topo->get_dimension();
        int code;
        switch (dim)
        {
        case 0:
            code=trouver((MG_SOMMET*)mg_topo)->get_id();
            break;
        case 1:
            code=trouver((MG_ARETE*)mg_topo)->get_id();
            break;
        case 2:
            code=trouver((MG_FACE*)mg_topo)->get_id();
            break;
            // case 3:code=trouver((MG_VOLUME*)mg_topo)->get_id();break;   pour plus tard;
        };
        codes.ajouter_groupes(&code);
        mgnoeud->change_nouveau_numero(i+1);
        pmmai.changer_champ_coordonnees(&coord);
        pmmai.changer_champ_codesnoeuds(&codes);
    }
 
}
 
 
#endif