remailleur/src/main.cpp

//------------------------------------------------------------
//------------------------------------------------------------
//  MAGiC
//  Jean Christophe Cuilli�e et Vincent FRANCOIS
//  D�artement de G�ie M�anique - UQTR
//------------------------------------------------------------
//  Le projet MAGIC est un projet de recherche du d�artement
//  de g�ie m�anique de l'Universit�du Qu�ec �
//  Trois Rivi�es
//  Les librairies ne peuvent �re utilis�s sans l'accord
//  des auteurs (contact : francois@uqtr.ca)
//------------------------------------------------------------
//------------------------------------------------------------
//
//       main.cpp
//
//------------------------------------------------------------
//------------------------------------------------------------
//  COPYRIGHT 2000
//  Version du 02/03/2006 �11H25
//------------------------------------------------------------
//------------------------------------------------------------
#include "gestionversion.h"
#ifdef WINDOWS_VERSION
#include "fenetre.h"
#endif

#pragma hdrstop
#include <string.h>
#include "vct_comparaison_resultat.h"
#include "vct_comparaison.h"
#include "mg_file.h"
#include "remailleur.h"
#include <iostream>
//---------------------------------------------------------------------------

void affiche(char* message)
{
#ifdef WINDOWS_VERSION
Form1->Memo1->Lines->Add(message);
#else
std::cout << message << std::endl;
#endif
}


#ifdef WINDOWS_VERSION
int amain(int argc,char **argv)
#else
int main(int argc,char **argv)
#endif
{
if (argc<3) return 1;
char fichier[1000],fichiers[1000],fichier2[1000],fichiercarte1[1000],fichiercarte2[1000],fichierres[1000];
int nummaiorigine=0,numgeomodifie=0,rescmp=0;
int diff=0;
for (int i=0;i<argc;i++)
        {
        if (strcmp(argv[i],"-origine")==0)        strcpy(fichier,argv[i+1]);
        if (strcmp(argv[i],"-modifie")==0)        strcpy(fichier2,argv[i+1]);
        if (strcmp(argv[i],"-carteorigine")==0)        strcpy(fichiercarte1,argv[i+1]);
        if (strcmp(argv[i],"-cartemodifie")==0)       strcpy(fichiercarte2,argv[i+1]);
        if (strcmp(argv[i],"-nummaiorigine")==0)       nummaiorigine=atoi(argv[i+1]);
        if (strcmp(argv[i],"-numgeomodifie")==0)       numgeomodifie=atoi(argv[i+1]);
        if (strcmp(argv[i],"-out")==0)       {strcpy(fichiers,argv[i+1]);diff=1;}
        if (strcmp(argv[i],"-rescmp")==0)       {strcpy(fichierres,argv[i+1]);rescmp=1;}

        }
if (!diff) strcpy(fichiers,fichier2);
affiche("*******************");
affiche("     REMAILLEUR    ");
affiche("*******************");
affiche("");
affiche("");
affiche("Lecture fichier 1");
MG_FILE *gest1=new MG_FILE(fichier);
affiche("Lecture fichier 2");
MG_FILE *gest2=new MG_FILE(fichier2);
MG_MAILLAGE* mai1=gest1->get_mg_maillage(nummaiorigine);
MG_GEOMETRIE* geo2=gest2->get_mg_geometrie(numgeomodifie);
affiche("Comparaison");
VCT_COMPARAISON_RESULTAT cmp(gest1,gest2,mai1->get_mg_geometrie(),geo2);
VCT_COMPARAISON compare(gest1,mai1->get_mg_geometrie(),gest2,geo2);
clock_t cpu=clock();
compare.trouver_les_identites(1);
clock_t cpu2=clock();
cout << "    identification en " << (cpu2-cpu)*1.0/CLOCKS_PER_SEC << " secondes" << endl;
affiche("Resultat comparaison");
compare.identifier_les_modifications_appliquees(cmp);
clock_t cpu3=clock();
cout << "    Localisation en " << (cpu3-cpu2)*1.0/CLOCKS_PER_SEC << " secondes" << endl;
if (rescmp) cmp.enregistrer(fichierres,fichier,fichiers);
FCT_GENERATEUR_3D<4> carte1;
carte1.lire(fichiercarte1);
FCT_GENERATEUR_3D<4> carte2;
carte2.lire(fichiercarte2);
MG_MAILLAGE* mmodi=new MG_MAILLAGE(geo2);
gest2->ajouter_mg_maillage(mmodi);
REMAILLEUR rm(gest1,gest2,mai1->get_mg_geometrie(),geo2,&carte1,&carte1,mai1,mmodi,cmp);
rm.active_affichage(affiche);
//rm.maille();
affiche("Enregistrer");
gest2->enregistrer(fichiers);
if (rescmp) cmp.enregistrer(fichierres,fichier,fichiers);
affiche("Fin");
{
int nbpoint=0;
int nbvecteur=0;
int taillemin=1000000;
int taillemax=0;

int nbface1=mai1->get_mg_geometrie()->get_nb_mg_face();
for (int i=0;i<nbface1;i++)
        {
        MG_FACE* face=mai1->get_mg_geometrie()->get_mg_face(i);
        int nb1=face->get_vectorisation().get_nb_points();
        nbpoint=nbpoint+nb1;    
        int nb2=face->get_vectorisation().get_vecteurs().size();
        nbvecteur=nbvecteur+nb2;
        if (nb2<taillemin) taillemin=nb2;
        if (nb2>taillemax) taillemax=nb2;
        int nb3=face->get_surface()->get_vectorisation().get_nb_points();
        nbpoint=nbpoint+nb3;    
        int nb4=face->get_surface()->get_vectorisation().get_vecteurs().size();
        nbvecteur=nbvecteur+nb4;
        if (nb4<taillemin) taillemin=nb4;
        if (nb4>taillemax) taillemax=nb4;
        }
cout << "Modele initial" << endl;
cout << nbface1 << " faces" << endl;
cout << nbpoint << " points" << endl;
cout << nbvecteur << " vecteurs" << endl;
cout << taillemin << " taille tenseur min" << endl;
cout << taillemax << " taille tenseur max" << endl;
}
{
int nbpoint=0;
int nbvecteur=0;
int taillemin=1000000;
int taillemax=0;
cout << "Modele modifie" << endl;
int nbface1=geo2->get_nb_mg_face();
for (int i=0;i<nbface1;i++)
        {
        MG_FACE* face=geo2->get_mg_face(i);
        int nb1=face->get_vectorisation().get_nb_points();
        nbpoint=nbpoint+nb1;    
        int nb2=face->get_vectorisation().get_vecteurs().size();
        nbvecteur=nbvecteur+nb2;
        if (nb2<taillemin) taillemin=nb2;
        if (nb2>taillemax) taillemax=nb2;
        int nb3=face->get_surface()->get_vectorisation().get_nb_points();
        nbpoint=nbpoint+nb3;    
        int nb4=face->get_surface()->get_vectorisation().get_vecteurs().size();
        nbvecteur=nbvecteur+nb4;
        if (nb4<taillemin) taillemin=nb4;
        if (nb4>taillemax) taillemax=nb4;
        }

cout << nbface1 << " faces" << endl;
cout << nbpoint << " points4" << endl;
cout << nbvecteur << " vecteurs" << endl;
cout << taillemin << " taille tenseur min" << endl;
cout << taillemax << " taille tenseur max" << endl;
}


delete gest1;
delete gest2;


return 0;
}


#pragma package(smart_init)
Revision:	253
Committed:	Tue Jul 13 19:40:46 2010 UTC (14 years, 11 months ago) by francois
File size:	6006 byte(s)
Log Message:	changement de hiearchie et utilisation de ccmake + mise a jour
#	User	Rev	Content
1	francois	147	//------------------------------------------------------------
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	francois	169	#include <string.h>
30	francois	147	#include "vct_comparaison_resultat.h"
31			#include "vct_comparaison.h"
32			#include "mg_file.h"
33			#include "remailleur.h"
34			#include <iostream>
35			//---------------------------------------------------------------------------
36
37			void affiche(char* message)
38			{
39			#ifdef WINDOWS_VERSION
40			Form1->Memo1->Lines->Add(message);
41			#else
42			std::cout << message << std::endl;
43			#endif
44			}
45
46
47			#ifdef WINDOWS_VERSION
48			int amain(int argc,char **argv)
49			#else
50			int main(int argc,char **argv)
51			#endif
52			{
53			if (argc<3) return 1;
54	francois	164	char fichier[1000],fichiers[1000],fichier2[1000],fichiercarte1[1000],fichiercarte2[1000],fichierres[1000];
55			int nummaiorigine=0,numgeomodifie=0,rescmp=0;
56	francois	147	int diff=0;
57			for (int i=0;i<argc;i++)
58			{
59			if (strcmp(argv[i],"-origine")==0) strcpy(fichier,argv[i+1]);
60			if (strcmp(argv[i],"-modifie")==0) strcpy(fichier2,argv[i+1]);
61			if (strcmp(argv[i],"-carteorigine")==0) strcpy(fichiercarte1,argv[i+1]);
62			if (strcmp(argv[i],"-cartemodifie")==0) strcpy(fichiercarte2,argv[i+1]);
63			if (strcmp(argv[i],"-nummaiorigine")==0) nummaiorigine=atoi(argv[i+1]);
64			if (strcmp(argv[i],"-numgeomodifie")==0) numgeomodifie=atoi(argv[i+1]);
65			if (strcmp(argv[i],"-out")==0) {strcpy(fichiers,argv[i+1]);diff=1;}
66	francois	164	if (strcmp(argv[i],"-rescmp")==0) {strcpy(fichierres,argv[i+1]);rescmp=1;}
67
68	francois	147	}
69			if (!diff) strcpy(fichiers,fichier2);
70			affiche("*******************");
71			affiche(" REMAILLEUR ");
72			affiche("*******************");
73			affiche("");
74			affiche("");
75			affiche("Lecture fichier 1");
76			MG_FILE *gest1=new MG_FILE(fichier);
77			affiche("Lecture fichier 2");
78			MG_FILE *gest2=new MG_FILE(fichier2);
79	francois	155	MG_MAILLAGE* mai1=gest1->get_mg_maillage(nummaiorigine);
80			MG_GEOMETRIE* geo2=gest2->get_mg_geometrie(numgeomodifie);
81			affiche("Comparaison");
82	francois	147	VCT_COMPARAISON_RESULTAT cmp(gest1,gest2,mai1->get_mg_geometrie(),geo2);
83			VCT_COMPARAISON compare(gest1,mai1->get_mg_geometrie(),gest2,geo2);
84	francois	226	clock_t cpu=clock();
85	francois	147	compare.trouver_les_identites(1);
86	francois	253	clock_t cpu2=clock();
87			cout << " identification en " << (cpu2-cpu)*1.0/CLOCKS_PER_SEC << " secondes" << endl;
88	francois	147	affiche("Resultat comparaison");
89			compare.identifier_les_modifications_appliquees(cmp);
90	francois	253	clock_t cpu3=clock();
91			cout << " Localisation en " << (cpu3-cpu2)*1.0/CLOCKS_PER_SEC << " secondes" << endl;
92	francois	164	if (rescmp) cmp.enregistrer(fichierres,fichier,fichiers);
93	francois	147	FCT_GENERATEUR_3D<4> carte1;
94			carte1.lire(fichiercarte1);
95			FCT_GENERATEUR_3D<4> carte2;
96			carte2.lire(fichiercarte2);
97			MG_MAILLAGE* mmodi=new MG_MAILLAGE(geo2);
98			gest2->ajouter_mg_maillage(mmodi);
99			REMAILLEUR rm(gest1,gest2,mai1->get_mg_geometrie(),geo2,&carte1,&carte1,mai1,mmodi,cmp);
100			rm.active_affichage(affiche);
101	francois	253	//rm.maille();
102	francois	147	affiche("Enregistrer");
103			gest2->enregistrer(fichiers);
104	francois	164	if (rescmp) cmp.enregistrer(fichierres,fichier,fichiers);
105	francois	147	affiche("Fin");
106	francois	253	{
107			int nbpoint=0;
108			int nbvecteur=0;
109			int taillemin=1000000;
110			int taillemax=0;
111
112			int nbface1=mai1->get_mg_geometrie()->get_nb_mg_face();
113			for (int i=0;i<nbface1;i++)
114			{
115			MG_FACE* face=mai1->get_mg_geometrie()->get_mg_face(i);
116			int nb1=face->get_vectorisation().get_nb_points();
117			nbpoint=nbpoint+nb1;
118			int nb2=face->get_vectorisation().get_vecteurs().size();
119			nbvecteur=nbvecteur+nb2;
120			if (nb2<taillemin) taillemin=nb2;
121			if (nb2>taillemax) taillemax=nb2;
122			int nb3=face->get_surface()->get_vectorisation().get_nb_points();
123			nbpoint=nbpoint+nb3;
124			int nb4=face->get_surface()->get_vectorisation().get_vecteurs().size();
125			nbvecteur=nbvecteur+nb4;
126			if (nb4<taillemin) taillemin=nb4;
127			if (nb4>taillemax) taillemax=nb4;
128			}
129			cout << "Modele initial" << endl;
130			cout << nbface1 << " faces" << endl;
131			cout << nbpoint << " points" << endl;
132			cout << nbvecteur << " vecteurs" << endl;
133			cout << taillemin << " taille tenseur min" << endl;
134			cout << taillemax << " taille tenseur max" << endl;
135			}
136			{
137			int nbpoint=0;
138			int nbvecteur=0;
139			int taillemin=1000000;
140			int taillemax=0;
141			cout << "Modele modifie" << endl;
142			int nbface1=geo2->get_nb_mg_face();
143			for (int i=0;i<nbface1;i++)
144			{
145			MG_FACE* face=geo2->get_mg_face(i);
146			int nb1=face->get_vectorisation().get_nb_points();
147			nbpoint=nbpoint+nb1;
148			int nb2=face->get_vectorisation().get_vecteurs().size();
149			nbvecteur=nbvecteur+nb2;
150			if (nb2<taillemin) taillemin=nb2;
151			if (nb2>taillemax) taillemax=nb2;
152			int nb3=face->get_surface()->get_vectorisation().get_nb_points();
153			nbpoint=nbpoint+nb3;
154			int nb4=face->get_surface()->get_vectorisation().get_vecteurs().size();
155			nbvecteur=nbvecteur+nb4;
156			if (nb4<taillemin) taillemin=nb4;
157			if (nb4>taillemax) taillemax=nb4;
158			}
159
160			cout << nbface1 << " faces" << endl;
161			cout << nbpoint << " points4" << endl;
162			cout << nbvecteur << " vecteurs" << endl;
163			cout << taillemin << " taille tenseur min" << endl;
164			cout << taillemax << " taille tenseur max" << endl;
165			}
166
167
168
169
170	francois	147	delete gest1;
171	francois	226	delete gest2;
172	francois	253
173
174
175
176
177
178	francois	147	return 0;
179			}
180
181
182
183			#pragma package(smart_init)