ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/REPOS_ERICCA/magic/app/toIbrep/src/main.cpp
Revision: 375
Committed: Fri Nov 9 15:27:14 2012 UTC (12 years, 9 months ago) by francois
File size: 5235 byte(s)
Log Message:
Solution FEM sont vectorielles et tensorielles + affichage des options de compilations dans les exe

File Contents

# User Rev Content
1 francois 249 //------------------------------------------------------------
2     //------------------------------------------------------------
3     // MAGiC
4     // Jean Christophe Cuilli�e et Vincent FRANCOIS
5     // D�artement de G�ie M�anique - UQTR
6     //------------------------------------------------------------
7     // Le projet MAGIC est un projet de recherche du d�artement
8     // de g�ie m�anique de l'Universit�du Qu�ec �
9     // Trois Rivi�es
10     // Les librairies ne peuvent �re utilis�s sans l'accord
11     // des auteurs (contact : francois@uqtr.ca)
12     //------------------------------------------------------------
13     //------------------------------------------------------------
14     //
15     // main.cpp
16     //
17     //------------------------------------------------------------
18     //------------------------------------------------------------
19     // COPYRIGHT 2000
20     // Version du 02/03/2006 �11H25
21     //------------------------------------------------------------
22     //------------------------------------------------------------
23     #include "gestionversion.h"
24     #ifdef WINDOWS_VERSION
25     #include "fenetre.h"
26     #endif
27    
28     #pragma hdrstop
29    
30     #include "mg_file.h"
31     #include "vct_multi_modele.h"
32     #include "vct_comparaison.h"
33     #include "tpl_octree.h"
34     #include "tpl_grille.h"
35     #include "ot_decalage_parametre.h"
36 francois 259 #include "toibrep.h"
37 francois 276 #include "ot_cpu.h"
38     #include "fct_generateur_3D.h"
39     #include "mailleur3d.h"
40 francois 249 #include <fstream>
41     #include <string.h>
42     //---------------------------------------------------------------------------
43    
44     void affiche(char* message)
45     {
46     #ifdef WINDOWS_VERSION
47     Form1->Memo1->Lines->Add(message);
48     #else
49     std::cout << message << std::endl;
50     #endif
51     }
52    
53    
54     #ifdef WINDOWS_VERSION
55     int amain(int argc,char **argv)
56     #else
57     int main(int argc,char **argv)
58     #endif
59     {
60     char nomfichier[500];
61 francois 276 char fichiercarte[500];
62 francois 249 char nomfichiersortie[500];
63 francois 258 char nomfichieribrep[500];
64 francois 276 int nummai=0;
65     int numgeo=0;
66     int numgroupe=0;
67     int pas=50;
68     int degre=1;
69     double coef=2.;
70 francois 375 nomfichiersortie[0]='0';
71     bool versioncompilee=false;
72 francois 249 for (int i=0;i<argc;i++)
73     {
74     if (strcmp(argv[i],"-fichier")==0) strcpy(nomfichier,argv[i+1]);
75     if (strcmp(argv[i],"-out")==0) strcpy(nomfichiersortie,argv[i+1]);
76 francois 258 if (strcmp(argv[i],"-ibrep")==0) strcpy(nomfichieribrep,argv[i+1]);
77 francois 249 if (strcmp(argv[i],"-nummaillage")==0) nummai=atoi(argv[i+1]);
78     if (strcmp(argv[i],"-numgeometrie")==0) numgeo=atoi(argv[i+1]);
79     if (strcmp(argv[i],"-groupe")==0) numgroupe=atoi(argv[i+1]);
80 francois 276 if (strcmp(argv[i],"-pas")==0) pas=atoi(argv[i+1]);
81     if (strcmp(argv[i],"-carte")==0) strcpy(fichiercarte,argv[i+1]);
82     if (strcmp(argv[i],"-degre")==0) degre=atoi(argv[i+1]);
83     if (strcmp(argv[i],"-coef")==0) coef=atof(argv[i+1]);
84 francois 375 if (strcmp(argv[i],"-version")==0) versioncompilee=true;
85 francois 249 }
86 francois 276 if (nomfichiersortie[0]=='0') strcpy(nomfichiersortie,nomfichier);
87     affiche((char*)"***************************************");
88     affiche((char*)" MODELE XFEM ");
89     affiche((char*)" UQTR - ERICCA ");
90     affiche((char*)"***************************************");
91 francois 375 if (versioncompilee)
92     {
93     GESTIONVERSION v;
94     char version[500];
95     v.print(version);
96     affiche(version);
97     }
98 francois 276 affiche((char*)"");
99     affiche((char*)"");
100     OT_CPU compteur;
101 francois 272 affiche((char*)"Debut");
102     affiche((char*)"Lecture fichier magic");
103 francois 249 MG_FILE gest(nomfichier);
104     MG_GEOMETRIE* geo=gest.get_mg_geometrieid(numgeo);
105     MG_GROUPE_TOPOLOGIQUE* mggt=geo->get_mg_groupe_topologiqueid(numgroupe);
106 francois 276 compteur.ajouter_etape("Lecture fichier");
107     FEM_MAILLAGE* mai;
108     if (nummai!=0) mai=gest.get_fem_maillageid(nummai);
109     else
110     {
111     affiche((char*)"Lecture carte de taille");
112     FCT_GENERATEUR_3D<4> carte;
113     carte.lire(fichiercarte);
114     compteur.ajouter_etape("Lecture carte de taille");
115     affiche((char*)"Mailleur octal");
116     MG_MAILLAGE* mgmai=new MG_MAILLAGE(NULL);
117     gest.ajouter_mg_maillage(mgmai);
118     MAILLEUR3D mailleur(&carte,mgmai);
119     mailleur.change_coef(coef);
120     mailleur.maille();
121     compteur.ajouter_etape("Mailleur octal");
122     affiche((char*)"Maillage FEM");
123     mai=new FEM_MAILLAGE(mgmai->get_mg_geometrie(),mgmai,degre);
124     gest.ajouter_fem_maillage(mai);
125     mai->construire(0);
126     compteur.ajouter_etape("Mailleur FEM");
127     }
128     TOIBREP t(&gest,geo,mai,pas,&compteur);
129 francois 258 //t.exporter_IBrep(nomfichieribrep,mggt);
130 francois 276 t.active_affichage(affiche);
131    
132    
133     IBrep brep=t.importer(nomfichiersortie,nomfichieribrep,mggt);
134    
135    
136    
137    
138     affiche((char*)"enregistrement magic");
139 francois 249 gest.enregistrer(nomfichiersortie);
140 francois 276 affiche((char*)"enregistrement gmsh");
141     int nlong=strlen(nomfichieribrep);
142     for (int i=0;i<nlong;i++)
143     if (nomfichieribrep[i]=='.') nomfichieribrep[i]='_';
144     strcat(nomfichieribrep,".pos");
145     std::ofstream output2((char*)nomfichieribrep);
146     brep.togmshpos(output2);
147     output2.close();
148 francois 272 affiche((char*)"fin totale");
149 francois 276 compteur.ajouter_etape("Enregistrement");
150 francois 249
151 francois 276
152     std::vector<std::string> tab;
153     compteur.get_tabfinal(tab);
154     for (int i=0;i<tab.size();i++)
155     affiche((char*)tab[i].c_str());
156     affiche((char*)"");
157     compteur.get_tabfinalcompresse(tab);
158     for (int i=0;i<tab.size();i++)
159     affiche((char*)tab[i].c_str());
160    
161    
162    
163    
164    
165    
166 francois 249 return 0;
167     }
168    
169    
170    
171     #pragma package(smart_init)