ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/REPOS_ERICCA/magic/app/mailleur/src/main3.cpp
Revision: 424
Committed: Fri Sep 13 14:39:52 2013 UTC (11 years, 10 months ago) by francois
File size: 4797 byte(s)
Log Message:
analyseur de qualité de maillage plus complet apres les mailleurs 3D

File Contents

# User Rev Content
1 francois 422 //------------------------------------------------------------
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     #include "ot_cpu.h"
37     #include "fct_generateur_3D.h"
38     #include "mailleur3d_peau.h"
39     #include <fstream>
40     #include <string.h>
41     #include "mg_maillage_algos.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 fichier[500];
61     char fichiersortie[500];
62 francois 423 strcpy(fichiersortie,"");
63 francois 422 char fichiercarte[500];
64     double critere=1.56;
65     int analyse=0;
66     char message[500];
67     int niveau=2;
68     bool versioncompilee=false;
69     int nummai=0;
70     if (argc<3)
71     {
72     sprintf(message,"mailleur_peau.exe -magic fich.magic -out fich2.magic -carte fich.ctt -mesh num");
73     affiche(message);
74     return 0;
75     }
76     for (int i=0;i<argc;i++)
77     {
78     if (strcmp(argv[i],"-carte")==0) strcpy(fichiercarte,argv[i+1]);
79     if (strcmp(argv[i],"-magic")==0) strcpy(fichier,argv[i+1]);
80     if (strcmp(argv[i],"-out")==0) strcpy(fichiersortie,argv[i+1]);
81     if (strcmp(argv[i],"-niveauopt")==0) niveau=atoi(argv[i+1]);
82     if (strcmp(argv[i],"-version")==0) versioncompilee=true;
83     if (strcmp(argv[i],"-mesh")==0) nummai=atoi(argv[i+1]);
84 francois 424 if (strcmp(argv[i],"-analyse")==0) analyse=1;
85 francois 422 }
86 francois 423 if (strcmp(fichiersortie,"")==0) strcpy(fichiersortie,fichier);
87 francois 422 sprintf(message," Mailleur 3D d'une enveloppe triangulee");
88     affiche(message);
89     if (versioncompilee)
90     {
91     GESTIONVERSION v;
92     char version[500];
93     v.print(version);
94     affiche(version);
95     }
96     sprintf(message," Fichier carte : %s",fichiercarte);
97     affiche(message);
98     affiche((char*)"Lecture fichier magic");
99     MG_FILE gest(fichier);
100     affiche((char*)"Lecture carte de taille");
101     FCT_GENERATEUR_3D<4> carte;
102     carte.lire(fichiercarte);
103     MG_MAILLAGE* mai;
104     if (nummai==0) mai=gest.get_mg_maillage(0);
105     else mai=gest.get_mg_maillageid(nummai);
106     MAILLEUR3D_PEAU m3d(mai,NULL,NULL,&carte);
107     m3d.active_affichage(affiche);
108     m3d.change_niveau_optimisation(niveau);
109     affiche((char*)"Maillage 3D");
110     m3d.maille(NULL,NULL,NULL);
111 francois 424 if (analyse)
112     {
113     affiche((char*)"Analyse maillage");
114     double vol;
115     int tab[4];
116     m3d.analyse_maillage_obtenu(vol,tab);
117     char mess[100];
118     sprintf(mess," 1) Respect de la carte de taille");
119     affiche(mess);
120     sprintf(mess," nombre de tetra theorique %f, nombre de tetra pratique %d ",vol,mai->get_nb_mg_tetra());
121     affiche(mess);
122     sprintf(mess," erreur : %.2f%%",(mai->get_nb_mg_tetra()-vol)*100./vol);
123     affiche(mess);
124     sprintf(mess," 2) Qualite du maillage");
125     affiche(mess);
126     sprintf(mess," nombre de tetra bon (qualite>0.5) : %d (%.2f%%)",tab[3],tab[3]*100./mai->get_nb_mg_tetra());
127     affiche(mess);
128     sprintf(mess," nombre de tetra acceptable (qualite>0.2) : %d (%.2f%%) ",tab[2],tab[2]*100./mai->get_nb_mg_tetra());
129     affiche(mess);
130     sprintf(mess," nombre de tetra moyen (qualite>0.1) : %d (%.2f%%) ",tab[1],tab[1]*100./mai->get_nb_mg_tetra());
131     affiche(mess);
132     sprintf(mess," nombre de tetra mauvais (qualite>0.) : %d (%.2f%%) ",tab[0],tab[0]*100./mai->get_nb_mg_tetra());
133     affiche(mess);
134     }
135 francois 422 affiche((char*)"Enregistrement");
136     MG_MAILLAGE::DIMENSIONMAILLAGESANSTOPO=3;
137     gest.enregistrer(fichiersortie);
138     affiche((char*)"Fin");
139     return 0;
140     }
141    
142    
143    
144     #pragma package(smart_init)