| 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 |
#include "magic_application.h"
|
| 25 |
#include <string.h>
|
| 26 |
#include "mg_file.h"
|
| 27 |
#include "fct_generateur_3d.h"
|
| 28 |
#include "fct_taille_fem_solution.h"
|
| 29 |
#include "mailleur3d_peau.h"
|
| 30 |
|
| 31 |
int main(int argc,char **argv)
|
| 32 |
{
|
| 33 |
std::vector<MAGIC_PARAMETRE_APPLICATION> lst;
|
| 34 |
MAGIC_PARAMETRE_APPLICATION p1(3,(char*)"-carte",MAGIC_PARAMETRE_APPLICATION::STRING,(char*)"Fichier de carte de taille");lst.push_back(p1);
|
| 35 |
MAGIC_PARAMETRE_APPLICATION p2(4,(char*)"-niveauopt",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Limite supérieure des mailles à optimiser - Défaut 2");lst.push_back(p2);
|
| 36 |
MAGIC_PARAMETRE_APPLICATION p3(5,(char*)"-nummai",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Numero de maillage à mailler - Défaut le premier");lst.push_back(p3);
|
| 37 |
MAGIC_PARAMETRE_APPLICATION p4(6,(char*)"-analyse",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Ajout d'une analyse de la qualité du maillage en fin de programme");lst.push_back(p4);
|
| 38 |
MAGIC_PARAMETRE_APPLICATION p5(7,(char*)"-enregistrement",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Enregistrement de la dernière étape (pour debbugage) - Defaut false");lst.push_back(p5);
|
| 39 |
MAGIC_PARAMETRE_APPLICATION p6(8,(char*)"-numvol",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Numero de volume à mailler - Défaut le premier ou NULL s il n y a en pas");lst.push_back(p6);
|
| 40 |
MAGIC_PARAMETRE_APPLICATION p7(9,(char*)"-carte2",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Utilisation de la nouvelle version de carte de taille");lst.push_back(p7);
|
| 41 |
MAGIC_PARAMETRE_APPLICATION p8(10,(char*)"-pas",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Valeur du pas d'integration - Defaut 32");lst.push_back(p8);
|
| 42 |
|
| 43 |
MAGIC_APPLICATION app((char*)"Mailleur 3D d'une enveloppe triangulee",argc,argv,lst,true,true);
|
| 44 |
if (app.get_erreur()==true) return 0;
|
| 45 |
|
| 46 |
char fichier[500];
|
| 47 |
char fichiersortie[500];
|
| 48 |
char fichiersortie2[500];
|
| 49 |
char fichiercarte[500];
|
| 50 |
int niveau;
|
| 51 |
int nummai;
|
| 52 |
int numvol;
|
| 53 |
bool analyse;
|
| 54 |
bool save;
|
| 55 |
bool version;
|
| 56 |
int pas;
|
| 57 |
|
| 58 |
app.recupere_parametre_string(1,fichier,(char*)"Fichier entrant manquant");
|
| 59 |
app.recupere_parametre_string_avec_defaut(2,fichiersortie,fichier);
|
| 60 |
app.recupere_parametre_string(3,fichiercarte,(char*)"Fichier carte manquant");
|
| 61 |
app.recupere_parametre_int_avec_defaut(4,niveau,2);
|
| 62 |
app.recupere_parametre_int_avec_defaut(5,nummai,0);
|
| 63 |
app.recupere_parametre_bool_avec_defaut(6,analyse,false);
|
| 64 |
app.recupere_parametre_bool_avec_defaut(7,save,false);
|
| 65 |
app.recupere_parametre_int_avec_defaut(8,numvol,0);
|
| 66 |
app.recupere_parametre_bool_avec_defaut(9,version,false);
|
| 67 |
app.recupere_parametre_int_avec_defaut(10,pas,32);
|
| 68 |
|
| 69 |
if (app.get_erreur()==true) return 0;
|
| 70 |
app.affiche((char*)"Lecture fichier magic");
|
| 71 |
MG_FILE gest(fichier);
|
| 72 |
app.affiche((char*)"Lecture carte de taille");
|
| 73 |
FCT_TAILLE *metrique;
|
| 74 |
if (!version)
|
| 75 |
{
|
| 76 |
FCT_GENERATEUR_3D<4> *carte=new FCT_GENERATEUR_3D<4>;
|
| 77 |
carte->lire(fichiercarte);
|
| 78 |
metrique=carte;
|
| 79 |
}
|
| 80 |
else
|
| 81 |
{
|
| 82 |
FCT_TAILLE_FEM_SOLUTION *carte=new FCT_TAILLE_FEM_SOLUTION(fichiercarte);
|
| 83 |
metrique=carte;
|
| 84 |
}
|
| 85 |
MG_MAILLAGE* mai;
|
| 86 |
if (nummai==0) mai=gest.get_mg_maillage(0);
|
| 87 |
else mai=gest.get_mg_maillageid(nummai);
|
| 88 |
MG_GEOMETRIE* mggeo=mai->get_mg_geometrie();
|
| 89 |
MG_VOLUME *mgvol=NULL;
|
| 90 |
if (mggeo!=NULL)
|
| 91 |
{
|
| 92 |
if (numvol==0) mgvol=mggeo->get_mg_volume(0);
|
| 93 |
else mgvol=mggeo->get_mg_volumeid(numvol);
|
| 94 |
}
|
| 95 |
MAILLEUR3D_PEAU m3d(mai,NULL,NULL,metrique,save);
|
| 96 |
if (save)
|
| 97 |
{
|
| 98 |
char fichiersortie2[500];
|
| 99 |
strcpy(fichiersortie2,fichiersortie);
|
| 100 |
strcat(fichiersortie2,"dbg");
|
| 101 |
m3d.change_nom_fichierdbg(fichiersortie2);
|
| 102 |
}
|
| 103 |
m3d.active_affichage(app.affiche);
|
| 104 |
m3d.change_niveau_optimisation(niveau);
|
| 105 |
m3d.change_pas_integration(pas);
|
| 106 |
app.affiche((char*)"Maillage 3D");
|
| 107 |
m3d.maille(mgvol,NULL,NULL);
|
| 108 |
if (analyse==true)
|
| 109 |
{
|
| 110 |
app.affiche((char*)"Analyse maillage");
|
| 111 |
double vol;
|
| 112 |
int tab[4];
|
| 113 |
m3d.analyse_maillage_obtenu(vol,tab);
|
| 114 |
char mess[100];
|
| 115 |
sprintf(mess," 1) Respect de la carte de taille");
|
| 116 |
app.affiche(mess);
|
| 117 |
sprintf(mess," nombre de tetra theorique %f, nombre de tetra pratique %d ",vol,mai->get_nb_mg_tetra());
|
| 118 |
app.affiche(mess);
|
| 119 |
sprintf(mess," erreur : %.2f%%",(mai->get_nb_mg_tetra()-vol)*100./vol);
|
| 120 |
app.affiche(mess);
|
| 121 |
sprintf(mess," 2) Qualite du maillage");
|
| 122 |
app.affiche(mess);
|
| 123 |
sprintf(mess," nombre de tetra bon (qualite>0.5) : %d (%.2f%%)",tab[3],tab[3]*100./mai->get_nb_mg_tetra());
|
| 124 |
app.affiche(mess);
|
| 125 |
sprintf(mess," nombre de tetra acceptable (qualite>0.2) : %d (%.2f%%) ",tab[2],tab[2]*100./mai->get_nb_mg_tetra());
|
| 126 |
app.affiche(mess);
|
| 127 |
sprintf(mess," nombre de tetra moyen (qualite>0.1) : %d (%.2f%%) ",tab[1],tab[1]*100./mai->get_nb_mg_tetra());
|
| 128 |
app.affiche(mess);
|
| 129 |
sprintf(mess," nombre de tetra mauvais (qualite>0.) : %d (%.2f%%) ",tab[0],tab[0]*100./mai->get_nb_mg_tetra());
|
| 130 |
app.affiche(mess);
|
| 131 |
}
|
| 132 |
app.affiche((char*)"Enregistrement");
|
| 133 |
MG_MAILLAGE::DIMENSIONMAILLAGESANSTOPO=3;
|
| 134 |
gest.enregistrer(fichiersortie);
|
| 135 |
delete metrique;
|
| 136 |
app.affiche((char*)"Fin");
|
| 137 |
return 0;
|
| 138 |
}
|
| 139 |
|
| 140 |
|
| 141 |
|
| 142 |
|
