CAD4FE/src/CAD4FE_MCFile.cpp

//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//  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/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       CAD4FE_MCFile.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:56 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------


#pragma hdrstop 
                                  
#include "gestionversion.h"

#include "CAD4FE_MCFile.h"


#pragma package(smart_init)

#include "mg_file.h"
#include "CAD4FE_MCFile.h"
#include "pars_argument.h"
#include "parse.h"

#include "CAD4FE_MCEdge.h"
#include "CAD4FE_MCFace.h"
#include "CAD4FE_MCVertex.h"
#include "CAD4FE_PolySurface.h"
#include "CAD4FE_PolyCurve.h"
#include "CAD4FE_MCNode.h"
#include "CAD4FE_MCSegment.h"
#include "CAD4FE_MCTriangle.h"
#include "CAD4FE_m3d_MCTriangle.h"


#pragma package(smart_init)

using namespace CAD4FE;

MCFile::MCFile():MG_FILE()
{
}

MCFile::MCFile(char* chemin):MG_FILE(chemin)
{     
FILE *in;
char ligne[3000];
in=fopen(chemin,"rt");
if (in!=NULL)
{
PARS_ARGUMENT param[100];
PARSE parse;           

MG_GEOMETRIE *mggeo;
MG_MAILLAGE *mgmai;


do
        {
        int ierr;
        std::string ligne=parse.lire(in,';',&ierr);
        parse.decode(ligne.c_str(),"@;",param);
        const char *chaine=param[0].argument[0].c_str();
        if ((chaine[0]=='/') && (chaine[1]=='/'))
                {
                if ((chaine[2]=='*')&&(chaine[3]=='i'))
                        {
                        unsigned long id;
                        sscanf(chaine,"//*i:%lu;",&id);
                        entiteidmax=id;
                        }
                }
        else if (param[0].argument[0]!="FIN")
                {
                parse.decode(ligne.c_str(),"%@=@(@);",param);
                std::string entite=param[1].argument[0];
                long id=atol(param[0].argument[0].c_str());
                std::string data=param[2].argument[0] ;
                MG_MAILLAGE *mgmai;
                FEM_MAILLAGE *femmai;
                
                if (entite=="MAILLAGE")
                        {
                        parse.decode(data.c_str(),"@",param+2);
                        long idgeo=cid(param[2].argument[0]);
                        mggeo=get_mg_geometrieid(idgeo);
                        mgmai=get_mg_maillageid(id);
                        }         
                if (entite=="TETRAEDRE")
                        {
                        if (mgmai->get_mg_tetraid(id) != NULL) continue;
                        parse.decode(data.c_str(),"@,@,@,@,@",param+2);
                        long idtopo=cid(param[2].argument[0]);
                        long idn1=cid(param[3].argument[0]);
                        long idn2=cid(param[4].argument[0]);
                        long idn3=cid(param[5].argument[0]);
                        long idn4=cid(param[6].argument[0]);
                        MG_NOEUD* noeud1=mgmai->get_mg_noeudid(idn1);
                        MG_NOEUD* noeud2=mgmai->get_mg_noeudid(idn2);
                        MG_NOEUD* noeud3=mgmai->get_mg_noeudid(idn3);
                        MG_NOEUD* noeud4=mgmai->get_mg_noeudid(idn4);
                        MG_ELEMENT_TOPOLOGIQUE* topo;
                        if (idtopo>-1) topo=mggeo->get_mg_volumeid(idtopo);
                        else topo=NULL;
                        if (noeud1&&noeud2&&noeud3&&noeud4) mgmai->ajouter_mg_tetra(topo,noeud1,noeud2,noeud3,noeud4,id);
                        }    
                /*if (entite=="TRIANGLE")
                        {
                        if (mgmai->get_mg_triangleid(id) != NULL) continue;
                        parse.decode(data.c_str(),"@,@,@,@",param+2);
                        long idtopo=cid(param[2].argument[0]);
                        long idn1=cid(param[3].argument[0]);
                        long idn2=cid(param[4].argument[0]);
                        long idn3=cid(param[5].argument[0]);
                        MG_NOEUD* noeud1=mgmai->get_mg_noeudid(idn1);
                        MG_NOEUD* noeud2=mgmai->get_mg_noeudid(idn2);
                        MG_NOEUD* noeud3=mgmai->get_mg_noeudid(idn3);
                        MG_ELEMENT_TOPOLOGIQUE* topo;
                        if (idtopo>-1)
                                {
                                topo=mggeo->get_mg_faceid(idtopo);
                                if (topo==NULL) topo=mggeo->get_mg_volumeid(idtopo);
                                }
                        else topo=NULL;
                        if (noeud1&&noeud2&&noeud3)mgmai->ajouter_mg_triangle(topo,noeud1,noeud2,noeud3,id);
                        }
                if (entite=="SEGMENT")
                        {         
                        if (mgmai->get_mg_segmentid(id) != NULL) continue;
                        parse.decode(data.c_str(),"@,@,@",param+2);
                        long idtopo=cid(param[2].argument[0]);
                        long idn1=cid(param[3].argument[0]);
                        long idn2=cid(param[4].argument[0]);
                        MG_NOEUD* noeud1=mgmai->get_mg_noeudid(idn1);
                        MG_NOEUD* noeud2=mgmai->get_mg_noeudid(idn2);
                        MG_ELEMENT_TOPOLOGIQUE* topo;
                        if (idtopo>-1)
                                {
                                topo=mggeo->get_mg_areteid(idtopo);
                                if (topo==NULL) topo=mggeo->get_mg_faceid(idtopo);
                                if (topo==NULL) topo=mggeo->get_mg_volumeid(idtopo);
                                }
                        else topo=NULL;
                        MG_SEGMENT* seg=new MG_SEGMENT(id,topo,noeud1,noeud2);
                        if (noeud1&&noeud2) mgmai->ajouter_mg_segment(seg);
                        }             */
                if (entite=="CAD4FE_MCSEGMENT")
                        {
                        parse.decode(data.c_str(),"@,@,@",param+2);
                        long idtopo=cid(param[2].argument[0]);
                        long idn1=cid(param[3].argument[0]);
                        long idn2=cid(param[4].argument[0]);
                        MCNode* noeud1=(MCNode*)mgmai->get_mg_noeudid(idn1);
                        MCNode* noeud2=(MCNode*)mgmai->get_mg_noeudid(idn2);
                        MG_ELEMENT_TOPOLOGIQUE* topo;
                        if (idtopo>-1)
                                {
                                topo=mggeo->get_mg_areteid(idtopo);
                                if (topo==NULL) topo=mggeo->get_mg_faceid(idtopo);
                                if (topo==NULL) topo=mggeo->get_mg_volumeid(idtopo);
                                }
                        else topo=NULL;
                        MCSegment * seg = new MCSegment(id,topo,noeud1,noeud2);
                        mgmai->ajouter_mg_segment(seg);
                        }
                if (entite=="CAD4FE_MCNODE")
                        {
                        parse.decode(data.c_str(),"@,@,@,@,@",param+2);
                        long idRefTopo=cid(param[2].argument[0]);
                        long idMCTopo=cid(param[3].argument[0]);
                        double x=atof(param[4].argument[0].c_str());
                        double y=atof(param[5].argument[0].c_str());
                        double z=atof(param[6].argument[0].c_str());
                        MG_ELEMENT_TOPOLOGIQUE* refTopo=NULL;
                        if (refTopo==NULL) refTopo=mggeo->get_mg_sommetid(idRefTopo);
                        if (refTopo==NULL) refTopo=mggeo->get_mg_areteid(idRefTopo);
                        if (refTopo==NULL) refTopo=mggeo->get_mg_faceid(idRefTopo);
                        if (refTopo==NULL) refTopo=mggeo->get_mg_volumeid(idRefTopo);
                        MG_ELEMENT_TOPOLOGIQUE* mcTopo=NULL;
                        if (mcTopo==NULL) mcTopo=mggeo->get_mg_sommetid(idMCTopo);
                        if (mcTopo==NULL) mcTopo=mggeo->get_mg_areteid(idMCTopo);
                        if (mcTopo==NULL) mcTopo=mggeo->get_mg_faceid(idMCTopo);
                        if (mcTopo==NULL) mcTopo=mggeo->get_mg_volumeid(idMCTopo);
                        MCNode* noeud=new MCNode(id,mcTopo,refTopo,x,y,z);
                        mgmai->ajouter_mg_noeud(noeud);
                        }
                    if (entite=="CAD4FE_MCTRIANGLE")
                    {
                        parse.decode(data.c_str(),"@,@,@,@",param+2);
                        long idtopo=cid(param[2].argument[0]);
                        long idn1=cid(param[3].argument[0]);
                        long idn2=cid(param[4].argument[0]);
                        long idn3=cid(param[5].argument[0]);
                        MCNode* noeud1=(MCNode*)mgmai->get_mg_noeudid(idn1);
                        MCNode* noeud2=(MCNode*)mgmai->get_mg_noeudid(idn2);
                        MCNode* noeud3=(MCNode*)mgmai->get_mg_noeudid(idn3);
                        MCSegment* mgsegment[3];
                        MCNode *mgnoeud[3]={noeud1,noeud2,noeud3};
                        // la topo. d'un MCTriangle est obligatoirement une MCFace
                        MCFace* mcFace=(MCFace*)mggeo->get_mg_faceid(idtopo);
                        for (int i=0;i<3;i++)
                        {
                            mgsegment[i]=(MCSegment*)mgmai->get_mg_segment(mgnoeud[i]->get_id(),mgnoeud[(i+1)%3]->get_id());
                            if (mgsegment[i]==NULL)
                            {
                                mgsegment[i]=new MCSegment(mcFace,mgnoeud[i],mgnoeud[(i+1)%3]);
                                mgmai->ajouter_mg_segment(mgsegment[i]);
                            }
                        }
                        M3D_MCTriangle *triangle = new M3D_MCTriangle(id,mcFace,noeud1,noeud2,noeud3,mgsegment[0],mgsegment[1],mgsegment[2]);
                        mgmai->ajouter_mg_triangle(triangle);
                    }
                                                            
                }
        }
while (param[0].argument[0]!="FIN");

fclose(in);

}

}

MCFile::~MCFile()
{
}


Revision:	1158
Committed:	Thu Jun 13 22:18:49 2024 UTC (11 months ago) by francois
File size:	11067 byte(s)
Log Message:	compatibilité Ubuntu 22.04 Suppression des refeences à Windows Ajout d'une banière
#	User	Rev	Content
1	francois	1158	//####//------------------------------------------------------------
2			//####//------------------------------------------------------------
3			//####// MAGiC
4			//####// Jean Christophe Cuilliere et Vincent FRANCOIS
5			//####// Departement de Genie Mecanique - UQTR
6			//####//------------------------------------------------------------
7			//####// MAGIC est un projet de recherche de l equipe ERICCA
8			//####// du departement de genie mecanique de l Universite du Quebec a Trois Rivieres
9			//####// http://www.uqtr.ca/ericca
10			//####// http://www.uqtr.ca/
11			//####//------------------------------------------------------------
12			//####//------------------------------------------------------------
13			//####//
14			//####// CAD4FE_MCFile.cpp
15			//####//
16			//####//------------------------------------------------------------
17			//####//------------------------------------------------------------
18			//####// COPYRIGHT 2000-2024
19			//####// jeu 13 jun 2024 11:58:56 EDT
20			//####//------------------------------------------------------------
21			//####//------------------------------------------------------------
22	foucault	27
23
24	foucault	176	#pragma hdrstop
25	foucault	27
26			#include "gestionversion.h"
27
28			#include "CAD4FE_MCFile.h"
29
30
31			#pragma package(smart_init)
32
33			#include "mg_file.h"
34			#include "CAD4FE_MCFile.h"
35			#include "pars_argument.h"
36			#include "parse.h"
37
38			#include "CAD4FE_MCEdge.h"
39			#include "CAD4FE_MCFace.h"
40			#include "CAD4FE_MCVertex.h"
41			#include "CAD4FE_PolySurface.h"
42			#include "CAD4FE_PolyCurve.h"
43			#include "CAD4FE_MCNode.h"
44			#include "CAD4FE_MCSegment.h"
45			#include "CAD4FE_MCTriangle.h"
46	foucault	569	#include "CAD4FE_m3d_MCTriangle.h"
47	foucault	27
48
49
50			#pragma package(smart_init)
51
52			using namespace CAD4FE;
53
54			MCFile::MCFile():MG_FILE()
55			{
56			}
57
58			MCFile::MCFile(char* chemin):MG_FILE(chemin)
59			{
60			FILE *in;
61			char ligne[3000];
62			in=fopen(chemin,"rt");
63			if (in!=NULL)
64			{
65			PARS_ARGUMENT param[100];
66			PARSE parse;
67
68			MG_GEOMETRIE *mggeo;
69			MG_MAILLAGE *mgmai;
70
71
72			do
73			{
74			int ierr;
75			std::string ligne=parse.lire(in,';',&ierr);
76			parse.decode(ligne.c_str(),"@;",param);
77			const char *chaine=param[0].argument[0].c_str();
78			if ((chaine[0]=='/') && (chaine[1]=='/'))
79			{
80			if ((chaine[2]=='*')&&(chaine[3]=='i'))
81			{
82			unsigned long id;
83			sscanf(chaine,"//*i:%lu;",&id);
84			entiteidmax=id;
85			}
86			}
87			else if (param[0].argument[0]!="FIN")
88			{
89			parse.decode(ligne.c_str(),"%@=@(@);",param);
90			std::string entite=param[1].argument[0];
91			long id=atol(param[0].argument[0].c_str());
92			std::string data=param[2].argument[0] ;
93			MG_MAILLAGE *mgmai;
94			FEM_MAILLAGE *femmai;
95
96			if (entite=="MAILLAGE")
97			{
98			parse.decode(data.c_str(),"@",param+2);
99			long idgeo=cid(param[2].argument[0]);
100			mggeo=get_mg_geometrieid(idgeo);
101			mgmai=get_mg_maillageid(id);
102			}
103			if (entite=="TETRAEDRE")
104			{
105			if (mgmai->get_mg_tetraid(id) != NULL) continue;
106			parse.decode(data.c_str(),"@,@,@,@,@",param+2);
107			long idtopo=cid(param[2].argument[0]);
108			long idn1=cid(param[3].argument[0]);
109			long idn2=cid(param[4].argument[0]);
110			long idn3=cid(param[5].argument[0]);
111			long idn4=cid(param[6].argument[0]);
112			MG_NOEUD* noeud1=mgmai->get_mg_noeudid(idn1);
113			MG_NOEUD* noeud2=mgmai->get_mg_noeudid(idn2);
114			MG_NOEUD* noeud3=mgmai->get_mg_noeudid(idn3);
115			MG_NOEUD* noeud4=mgmai->get_mg_noeudid(idn4);
116			MG_ELEMENT_TOPOLOGIQUE* topo;
117			if (idtopo>-1) topo=mggeo->get_mg_volumeid(idtopo);
118			else topo=NULL;
119			if (noeud1&&noeud2&&noeud3&&noeud4) mgmai->ajouter_mg_tetra(topo,noeud1,noeud2,noeud3,noeud4,id);
120			}
121			/*if (entite=="TRIANGLE")
122			{
123			if (mgmai->get_mg_triangleid(id) != NULL) continue;
124			parse.decode(data.c_str(),"@,@,@,@",param+2);
125			long idtopo=cid(param[2].argument[0]);
126			long idn1=cid(param[3].argument[0]);
127			long idn2=cid(param[4].argument[0]);
128			long idn3=cid(param[5].argument[0]);
129			MG_NOEUD* noeud1=mgmai->get_mg_noeudid(idn1);
130			MG_NOEUD* noeud2=mgmai->get_mg_noeudid(idn2);
131			MG_NOEUD* noeud3=mgmai->get_mg_noeudid(idn3);
132			MG_ELEMENT_TOPOLOGIQUE* topo;
133			if (idtopo>-1)
134			{
135			topo=mggeo->get_mg_faceid(idtopo);
136			if (topo==NULL) topo=mggeo->get_mg_volumeid(idtopo);
137			}
138			else topo=NULL;
139			if (noeud1&&noeud2&&noeud3)mgmai->ajouter_mg_triangle(topo,noeud1,noeud2,noeud3,id);
140			}
141			if (entite=="SEGMENT")
142			{
143			if (mgmai->get_mg_segmentid(id) != NULL) continue;
144			parse.decode(data.c_str(),"@,@,@",param+2);
145			long idtopo=cid(param[2].argument[0]);
146			long idn1=cid(param[3].argument[0]);
147			long idn2=cid(param[4].argument[0]);
148			MG_NOEUD* noeud1=mgmai->get_mg_noeudid(idn1);
149			MG_NOEUD* noeud2=mgmai->get_mg_noeudid(idn2);
150			MG_ELEMENT_TOPOLOGIQUE* topo;
151			if (idtopo>-1)
152			{
153			topo=mggeo->get_mg_areteid(idtopo);
154			if (topo==NULL) topo=mggeo->get_mg_faceid(idtopo);
155			if (topo==NULL) topo=mggeo->get_mg_volumeid(idtopo);
156			}
157			else topo=NULL;
158			MG_SEGMENT* seg=new MG_SEGMENT(id,topo,noeud1,noeud2);
159			if (noeud1&&noeud2) mgmai->ajouter_mg_segment(seg);
160			} */
161			if (entite=="CAD4FE_MCSEGMENT")
162			{
163			parse.decode(data.c_str(),"@,@,@",param+2);
164			long idtopo=cid(param[2].argument[0]);
165			long idn1=cid(param[3].argument[0]);
166			long idn2=cid(param[4].argument[0]);
167			MCNode* noeud1=(MCNode*)mgmai->get_mg_noeudid(idn1);
168			MCNode* noeud2=(MCNode*)mgmai->get_mg_noeudid(idn2);
169			MG_ELEMENT_TOPOLOGIQUE* topo;
170			if (idtopo>-1)
171			{
172			topo=mggeo->get_mg_areteid(idtopo);
173			if (topo==NULL) topo=mggeo->get_mg_faceid(idtopo);
174			if (topo==NULL) topo=mggeo->get_mg_volumeid(idtopo);
175			}
176			else topo=NULL;
177			MCSegment * seg = new MCSegment(id,topo,noeud1,noeud2);
178			mgmai->ajouter_mg_segment(seg);
179			}
180			if (entite=="CAD4FE_MCNODE")
181			{
182			parse.decode(data.c_str(),"@,@,@,@,@",param+2);
183			long idRefTopo=cid(param[2].argument[0]);
184			long idMCTopo=cid(param[3].argument[0]);
185			double x=atof(param[4].argument[0].c_str());
186			double y=atof(param[5].argument[0].c_str());
187			double z=atof(param[6].argument[0].c_str());
188			MG_ELEMENT_TOPOLOGIQUE* refTopo=NULL;
189			if (refTopo==NULL) refTopo=mggeo->get_mg_sommetid(idRefTopo);
190			if (refTopo==NULL) refTopo=mggeo->get_mg_areteid(idRefTopo);
191			if (refTopo==NULL) refTopo=mggeo->get_mg_faceid(idRefTopo);
192			if (refTopo==NULL) refTopo=mggeo->get_mg_volumeid(idRefTopo);
193			MG_ELEMENT_TOPOLOGIQUE* mcTopo=NULL;
194			if (mcTopo==NULL) mcTopo=mggeo->get_mg_sommetid(idMCTopo);
195			if (mcTopo==NULL) mcTopo=mggeo->get_mg_areteid(idMCTopo);
196			if (mcTopo==NULL) mcTopo=mggeo->get_mg_faceid(idMCTopo);
197			if (mcTopo==NULL) mcTopo=mggeo->get_mg_volumeid(idMCTopo);
198			MCNode* noeud=new MCNode(id,mcTopo,refTopo,x,y,z);
199			mgmai->ajouter_mg_noeud(noeud);
200			}
201			if (entite=="CAD4FE_MCTRIANGLE")
202			{
203			parse.decode(data.c_str(),"@,@,@,@",param+2);
204			long idtopo=cid(param[2].argument[0]);
205			long idn1=cid(param[3].argument[0]);
206			long idn2=cid(param[4].argument[0]);
207			long idn3=cid(param[5].argument[0]);
208			MCNode* noeud1=(MCNode*)mgmai->get_mg_noeudid(idn1);
209			MCNode* noeud2=(MCNode*)mgmai->get_mg_noeudid(idn2);
210			MCNode* noeud3=(MCNode*)mgmai->get_mg_noeudid(idn3);
211			MCSegment* mgsegment[3];
212			MCNode *mgnoeud[3]={noeud1,noeud2,noeud3};
213			// la topo. d'un MCTriangle est obligatoirement une MCFace
214			MCFace* mcFace=(MCFace*)mggeo->get_mg_faceid(idtopo);
215			for (int i=0;i<3;i++)
216			{
217			mgsegment[i]=(MCSegment*)mgmai->get_mg_segment(mgnoeud[i]->get_id(),mgnoeud[(i+1)%3]->get_id());
218			if (mgsegment[i]==NULL)
219			{
220			mgsegment[i]=new MCSegment(mcFace,mgnoeud[i],mgnoeud[(i+1)%3]);
221			mgmai->ajouter_mg_segment(mgsegment[i]);
222			}
223			}
224			M3D_MCTriangle *triangle = new M3D_MCTriangle(id,mcFace,noeud1,noeud2,noeud3,mgsegment[0],mgsegment[1],mgsegment[2]);
225			mgmai->ajouter_mg_triangle(triangle);
226			}
227
228			}
229			}
230			while (param[0].argument[0]!="FIN");
231
232			fclose(in);
233
234			}
235
236			}
237
238			MCFile::~MCFile()
239			{
240			}
241
242
243
244