| 1 |
//------------------------------------------------------------
|
| 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 "tpl_octree.h"
|
| 33 |
#include "tpl_grille.h"
|
| 34 |
#include "ot_decalage_parametre.h"
|
| 35 |
#include "ot_cpu.h"
|
| 36 |
#include "fct_generateur_3d.h"
|
| 37 |
#include "fct_taille_fem_solution.h"
|
| 38 |
#include "mailleur_delaunay.h"
|
| 39 |
#include <fstream>
|
| 40 |
#include <string.h>
|
| 41 |
#include "mg_maillage_algos.h"
|
| 42 |
#include "magic_application.h"
|
| 43 |
#include "mailleur_fem.h"
|
| 44 |
|
| 45 |
|
| 46 |
//---------------------------------------------------------------------------
|
| 47 |
int main(int argc,char **argv)
|
| 48 |
|
| 49 |
{
|
| 50 |
std::vector<MAGIC_PARAMETRE_APPLICATION> lst;
|
| 51 |
MAGIC_PARAMETRE_APPLICATION p1(2,(char*)"-carte",MAGIC_PARAMETRE_APPLICATION::STRING,(char*)"Fichier de carte de taille");
|
| 52 |
MAGIC_PARAMETRE_APPLICATION p2(3,(char*)"-niveauopt",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Limite supérieure des mailles à optimiser - Défaut 2");
|
| 53 |
MAGIC_PARAMETRE_APPLICATION p3(4,(char*)"-criterearret",MAGIC_PARAMETRE_APPLICATION::DOUBLE,(char*)"Valeur maximale de la longueur des arêtes par rapport à la longueur initiale - Défaut 1.5");
|
| 54 |
MAGIC_PARAMETRE_APPLICATION p4(5,(char*)"-analyse",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Ajout d'une analyse de la qualité du maillage en fin de programme");
|
| 55 |
MAGIC_PARAMETRE_APPLICATION p5(6,(char*)"-delaunay",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Maillage d'une boite par la méthode de Delaunay avec ancienne carte de taille");
|
| 56 |
MAGIC_PARAMETRE_APPLICATION p6(7,(char*)"-octal",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Maillage d'une boite par la méthode octal");
|
| 57 |
MAGIC_PARAMETRE_APPLICATION p7(8,(char*)"-tailleoctree",MAGIC_PARAMETRE_APPLICATION::DOUBLE,(char*)"Taille des octants par rapport à la taille souhaitée - Défaut 2.");
|
| 58 |
MAGIC_PARAMETRE_APPLICATION p8(9,(char*)"-fem",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Création du maillage fem");
|
| 59 |
MAGIC_PARAMETRE_APPLICATION p9(10,(char*)"-degre",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Degre du maillage fem - Défaut 1");
|
| 60 |
MAGIC_PARAMETRE_APPLICATION p10(11,(char*)"-pas",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Valeur du pas d'integration - Defaut 5");
|
| 61 |
MAGIC_PARAMETRE_APPLICATION p11(12,(char*)"-delaunay2",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Maillage d'une boite par la méthode de Delaunay avec nouvelle carte de taille");
|
| 62 |
p5.ajouter_dependance(1);p5.ajouter_dependance(2);p5.ajouter_dependance(3);p5.ajouter_dependance(4),p5.ajouter_dependance(5);p5.ajouter_dependance(11);
|
| 63 |
p6.ajouter_dependance(1);p6.ajouter_dependance(2);p6.ajouter_dependance(8);p6.ajouter_dependance(9);p6.ajouter_dependance(10);
|
| 64 |
p11.ajouter_dependance(1);p11.ajouter_dependance(2);p11.ajouter_dependance(3);p11.ajouter_dependance(4),p11.ajouter_dependance(5);p11.ajouter_dependance(11);
|
| 65 |
lst.push_back(p1);
|
| 66 |
lst.push_back(p2);
|
| 67 |
lst.push_back(p3);
|
| 68 |
lst.push_back(p4);
|
| 69 |
lst.push_back(p5);
|
| 70 |
lst.push_back(p6);
|
| 71 |
lst.push_back(p7);
|
| 72 |
lst.push_back(p8);
|
| 73 |
lst.push_back(p9);
|
| 74 |
lst.push_back(p10);
|
| 75 |
lst.push_back(p11);
|
| 76 |
|
| 77 |
MAGIC_APPLICATION app((char*)"Mailleur d'un espace défini par une carte de taille",argc,argv,lst,false,true);
|
| 78 |
if (app.get_erreur()==true) return 0;
|
| 79 |
|
| 80 |
|
| 81 |
|
| 82 |
if ((app.get_action()==6) || (app.get_action()==12))
|
| 83 |
{
|
| 84 |
char fichier[500];
|
| 85 |
char fichiercarte[500];
|
| 86 |
double critere;
|
| 87 |
bool analyse;
|
| 88 |
char message[500];
|
| 89 |
int niveau;
|
| 90 |
int pas;
|
| 91 |
app.recupere_parametre_string(1,fichier,(char*)"Fichier de sortie manquant");
|
| 92 |
app.recupere_parametre_string(2,fichiercarte,(char*)"Fichier de carte de taille manquant",MAGIC_APPLICATION::VERIFFICHIEREXISTE);
|
| 93 |
app.recupere_parametre_int_avec_defaut(3,niveau,2);
|
| 94 |
app.recupere_parametre_double_avec_defaut(4,critere,1.5);
|
| 95 |
app.recupere_parametre_bool_avec_defaut(5,analyse,false);
|
| 96 |
app.recupere_parametre_int_avec_defaut(11,pas,5);
|
| 97 |
if (app.get_erreur()==true) return 0;
|
| 98 |
app.affiche((char*)"Méthode de Delaunay");
|
| 99 |
sprintf(message,"Fichier carte : %s",fichiercarte);
|
| 100 |
app.affiche(message);
|
| 101 |
app.affiche((char*)"");
|
| 102 |
app.affiche((char*)"Lecture carte de taille");
|
| 103 |
OT_CPU compteur;
|
| 104 |
FCT_TAILLE* metrique;
|
| 105 |
if (app.get_action()==6)
|
| 106 |
{
|
| 107 |
FCT_GENERATEUR_3D<4> *carte=new FCT_GENERATEUR_3D<4>;
|
| 108 |
carte->lire(fichiercarte);
|
| 109 |
metrique=carte;
|
| 110 |
}
|
| 111 |
//MG_FILE *gestcarte;
|
| 112 |
if (app.get_action()==12)
|
| 113 |
{
|
| 114 |
FCT_TAILLE_FEM_SOLUTION *carte=new FCT_TAILLE_FEM_SOLUTION(fichiercarte);
|
| 115 |
metrique=carte;
|
| 116 |
}
|
| 117 |
MG_GESTIONNAIRE* gest;
|
| 118 |
MG_FILE* file=new MG_FILE(fichier);
|
| 119 |
if (file->get_code_de_lecture()==0)
|
| 120 |
{
|
| 121 |
delete file;
|
| 122 |
MG_GESTIONNAIRE* gest2=new MG_GESTIONNAIRE;
|
| 123 |
gest=gest2;
|
| 124 |
}
|
| 125 |
else gest=file;
|
| 126 |
MG_MAILLAGE* mai=new MG_MAILLAGE(NULL);
|
| 127 |
gest->ajouter_mg_maillage(mai);
|
| 128 |
MAILLEUR_DELAUNAY m(mai,metrique,critere,&compteur);
|
| 129 |
m.active_affichage(app.affiche);
|
| 130 |
m.change_niveau_optimisation(niveau);
|
| 131 |
m.change_pas_integration(pas);
|
| 132 |
m.maille();
|
| 133 |
char fich[500];
|
| 134 |
strcpy(fich,fichier);
|
| 135 |
gest->enregistrer(fich);
|
| 136 |
if (analyse)
|
| 137 |
{
|
| 138 |
MG_MAILLAGE_ALGOS::calcule_caracteristiques(mai,fich,metrique);
|
| 139 |
strcat(fich,".magic");
|
| 140 |
gest->enregistrer(fich);
|
| 141 |
}
|
| 142 |
delete metrique;
|
| 143 |
compteur.ajouter_etape("Enregistrement");
|
| 144 |
std::vector<std::string> tab;
|
| 145 |
compteur.get_tabfinal(tab);
|
| 146 |
for (int i=0;i<tab.size();i++)
|
| 147 |
app.affiche((char*)tab[i].c_str());
|
| 148 |
}
|
| 149 |
if (app.get_action()==7)
|
| 150 |
{
|
| 151 |
char fichier[500];
|
| 152 |
char fichiercarte[500];
|
| 153 |
double coef;
|
| 154 |
bool fem;
|
| 155 |
int degre;
|
| 156 |
app.recupere_parametre_string(1,fichier,(char*)"Fichier de sortie manquant");
|
| 157 |
app.recupere_parametre_string(2,fichiercarte,(char*)"Fichier de carte de taille manquant",MAGIC_APPLICATION::VERIFFICHIEREXISTE);
|
| 158 |
app.recupere_parametre_double_avec_defaut(8,coef,2.);
|
| 159 |
app.recupere_parametre_bool_avec_defaut(9,fem,false);
|
| 160 |
app.recupere_parametre_int_avec_defaut(10,degre,1);
|
| 161 |
if (app.get_erreur()==true) return 0;
|
| 162 |
char message[5000];
|
| 163 |
app.affiche((char*)"Méthode octale");
|
| 164 |
sprintf(message," Fichier carte : %s",fichiercarte);
|
| 165 |
app.affiche(message);
|
| 166 |
app.affiche((char*)"Lecture carte de taille");
|
| 167 |
FCT_GENERATEUR_3D<4> carte;
|
| 168 |
carte.lire(fichiercarte);
|
| 169 |
app.affiche((char*)"Maillage");
|
| 170 |
MG_GESTIONNAIRE* gest;
|
| 171 |
MG_FILE* file=new MG_FILE(fichier);
|
| 172 |
if (file->get_code_de_lecture()==0)
|
| 173 |
{
|
| 174 |
delete file;
|
| 175 |
MG_GESTIONNAIRE* gest2=new MG_GESTIONNAIRE;
|
| 176 |
gest=gest2;
|
| 177 |
}
|
| 178 |
else gest=file;
|
| 179 |
MG_MAILLAGE* mgmai=new MG_MAILLAGE(NULL);
|
| 180 |
gest->ajouter_mg_maillage(mgmai);
|
| 181 |
MAILLEUR3D mailleur(&carte,mgmai);
|
| 182 |
mailleur.change_coef(coef);
|
| 183 |
mailleur.maille();
|
| 184 |
if (fem)
|
| 185 |
{
|
| 186 |
app.affiche((char*)"Maillage FEM");
|
| 187 |
FEM_MAILLAGE* femmai=new FEM_MAILLAGE(mgmai->get_mg_geometrie(),mgmai,degre);
|
| 188 |
gest->ajouter_fem_maillage(femmai);
|
| 189 |
MAILLEUR_FEM mf;
|
| 190 |
mf.maille(femmai,0);
|
| 191 |
}
|
| 192 |
app.affiche((char*)"Enregistrement");
|
| 193 |
gest->enregistrer(fichier);
|
| 194 |
delete gest;
|
| 195 |
app.affiche((char*)"Fin");
|
| 196 |
|
| 197 |
|
| 198 |
|
| 199 |
}
|
| 200 |
|
| 201 |
|
| 202 |
|
| 203 |
return 0;
|
| 204 |
}
|
| 205 |
|
| 206 |
|
| 207 |
|
| 208 |
#pragma package(smart_init)
|
