carte/src/main.cpp

//------------------------------------------------------------
//------------------------------------------------------------
//  MAGiC
//  Jean Christophe Cuilli�re et Vincent FRANCOIS
//  D�partement de G�nie M�canique - UQTR
//------------------------------------------------------------
//  Le projet MAGIC est un projet de recherche du d�partement
//  de g�nie m�canique de l'Universit� du Qu�bec �
//  Trois Rivi�res
//  Les librairies ne peuvent �tre utilis�es 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 "mg_file.h"
#include "mg_export.h"
#include "fct_generateur_frontiere.h"
#include "fct_generateur_constante.h"
#include "fct_generateur_fichier.h"
#include "fct_taille_fem_solution_generateur_constant.h"
#include "fct_taille_fem_solution_generateur_echantillon_frontiere.h"
#include "fct_taille_fem_solution_generateur_fichier.h"
#include "magic_application.h"
#include <string.h>
//---------------------------------------------------------------------------


int main(int argc,char **argv)
{
std::vector<MAGIC_PARAMETRE_APPLICATION> lst;
MAGIC_PARAMETRE_APPLICATION p1(3,(char*)"-fichierpoint",MAGIC_PARAMETRE_APPLICATION::STRING,(char*)"Fichier d'échantillonage");
MAGIC_PARAMETRE_APPLICATION p2(4,(char*)"-dg",MAGIC_PARAMETRE_APPLICATION::DOUBLE,(char*)"Densité globale");
MAGIC_PARAMETRE_APPLICATION p3(5,(char*)"-eps",MAGIC_PARAMETRE_APPLICATION::DOUBLE,(char*)"Écart entre le maillage et la géométrie");
MAGIC_PARAMETRE_APPLICATION p4(6,(char*)"-taillezone",MAGIC_PARAMETRE_APPLICATION::DOUBLE,(char*)"Facteur pour la distance de variation de la carte autour de l'échantillon -  Defaut 4");
MAGIC_PARAMETRE_APPLICATION p5(7,(char*)"-tauxtransfert",MAGIC_PARAMETRE_APPLICATION::DOUBLE,(char*)"Poid de l'échantillon virtuel place en dehors de la zone de raffinement - Defaut 0.25");
MAGIC_PARAMETRE_APPLICATION p6(8,(char*)"-fechantillonnage",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Nombre d'échantillon par entité - Defaut 20");
MAGIC_PARAMETRE_APPLICATION p7(9,(char*)"-festimation",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Nombre de cellule dans une direction pour la grille virtuelle - Defaut 10");
MAGIC_PARAMETRE_APPLICATION p8(10,(char*)"-constante",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Carte de taille constante");
MAGIC_PARAMETRE_APPLICATION p9(11,(char*)"-fichierechan",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Carte de taille par un echantillonnage dans un fichier exterieur");
MAGIC_PARAMETRE_APPLICATION p10(12,(char*)"-courbure",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Carte de taille par approximation des courbures");
MAGIC_PARAMETRE_APPLICATION p11(13,(char*)"-groupe",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Groupe topologique utilise pour faire le calcul - Defaut 0");
MAGIC_PARAMETRE_APPLICATION p12(14,(char*)"-numgeo",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Géométrie du groupe topologique - Defaut 0");
MAGIC_PARAMETRE_APPLICATION p13(15,(char*)"-constante2",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Carte de taille constante définie par une FEM_SOLUTION");
MAGIC_PARAMETRE_APPLICATION p14(16,(char*)"-nx",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Nombre d'intervals sur x - Defaut 10");
MAGIC_PARAMETRE_APPLICATION p15(17,(char*)"-ny",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Nombre d'intervals sur y - Defaut 10");
MAGIC_PARAMETRE_APPLICATION p16(18,(char*)"-nz",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Nombre d'intervals sur z - Defaut 10");
MAGIC_PARAMETRE_APPLICATION p19(19,(char*)"-courbure2",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Carte de taille par approximation des courbures définie par une FEM_SOLUTION");
MAGIC_PARAMETRE_APPLICATION p18(20,(char*)"-fichierechan2",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Carte de taille par un echantillonnage dans un fichier exterieur définie par une FEM_SOLUTION");
MAGIC_PARAMETRE_APPLICATION p20(21,(char*)"-constante2mai",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Carte de taille constante définie par une FEM_SOLUTION basé sur maillage");
MAGIC_PARAMETRE_APPLICATION p21(22,(char*)"-nummai",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Numéro de maillage - Defaut le premier");
p8.ajouter_dependance(1);p8.ajouter_dependance(2);p8.ajouter_dependance(4);p8.ajouter_dependance(8);p8.ajouter_dependance(9);p8.ajouter_dependance(13);p8.ajouter_dependance(14);
p9.ajouter_dependance(1);p9.ajouter_dependance(2);p9.ajouter_dependance(3);p9.ajouter_dependance(4);p9.ajouter_dependance(6);p9.ajouter_dependance(7);p9.ajouter_dependance(8);p9.ajouter_dependance(9);
p10.ajouter_dependance(1);p10.ajouter_dependance(2);p10.ajouter_dependance(3);p10.ajouter_dependance(4);p10.ajouter_dependance(5);p10.ajouter_dependance(6);p10.ajouter_dependance(7);p10.ajouter_dependance(8);p10.ajouter_dependance(9);p10.ajouter_dependance(13);p10.ajouter_dependance(14);
p13.ajouter_dependance(1);p13.ajouter_dependance(2);p13.ajouter_dependance(4);p13.ajouter_dependance(8);p13.ajouter_dependance(14);p13.ajouter_dependance(16);p13.ajouter_dependance(17);p13.ajouter_dependance(18);
p18.ajouter_dependance(1);p18.ajouter_dependance(2);p18.ajouter_dependance(3);p18.ajouter_dependance(4);p18.ajouter_dependance(6);p18.ajouter_dependance(7);p18.ajouter_dependance(8);p18.ajouter_dependance(14);p18.ajouter_dependance(16);p18.ajouter_dependance(17);p18.ajouter_dependance(18);
p19.ajouter_dependance(1);p19.ajouter_dependance(2);p19.ajouter_dependance(4);p19.ajouter_dependance(5);p19.ajouter_dependance(6);p19.ajouter_dependance(7);p19.ajouter_dependance(8);p19.ajouter_dependance(14);p19.ajouter_dependance(16);p19.ajouter_dependance(17);p19.ajouter_dependance(18);
p20.ajouter_dependance(1);p20.ajouter_dependance(2);p20.ajouter_dependance(4);p20.ajouter_dependance(8);p20.ajouter_dependance(22);p20.ajouter_dependance(16);p20.ajouter_dependance(17);p20.ajouter_dependance(18);

lst.push_back(p1);
lst.push_back(p2);
lst.push_back(p3);
lst.push_back(p4);
lst.push_back(p5);
lst.push_back(p6);
lst.push_back(p7);
lst.push_back(p8);
lst.push_back(p9);
lst.push_back(p10);
lst.push_back(p11);
lst.push_back(p12);
lst.push_back(p13);
lst.push_back(p14);
lst.push_back(p15);
lst.push_back(p16);
lst.push_back(p19);
lst.push_back(p18);
lst.push_back(p20);
lst.push_back(p21);
MAGIC_APPLICATION app((char*)"Generation de carte de taille",argc,argv,lst,true,true);
if (app.get_erreur()==true) return 0;  

    
if (app.get_action()==12)
    {
        char nomfichier[500];
        char nomfichiersortie[500];
        char nomfichierpoint[500];
        double eng;
        double eps;
        double ind1,ind2;
        int fechan,festim;
        int numgroupe,numgeo;
        app.recupere_parametre_string(1,nomfichier,(char*)"Fichier entrant manquant");
        app.recupere_parametre_string(2,nomfichiersortie,(char*)"Fichier de sortie manquant");
        app.recupere_parametre_string_avec_defaut(3,nomfichierpoint,(char*)"");
        app.recupere_parametre_double(4,eng,(char*)"Valeur de l'écart nodal global");
        app.recupere_parametre_double(5,eps,(char*)"Valeur de l'écart entre la géométrie et le maillage");
        app.recupere_parametre_double_avec_defaut(6,ind1,4.);
        app.recupere_parametre_double_avec_defaut(7,ind2,0.25);
        app.recupere_parametre_int_avec_defaut(8,fechan,20);
        app.recupere_parametre_int_avec_defaut(9,festim,10);
        app.recupere_parametre_int_avec_defaut(13,numgroupe,0);
        app.recupere_parametre_int_avec_defaut(14,numgeo,0);
        if (app.get_erreur()==true) return 0;  
        app.affiche((char*)"Generation de carte de taille basee sur la courbure des frontieres V1");  
        app.affiche((char*)"Lecture du fichier");
        MG_FILE gest(nomfichier);
        int nbptctl=0;
        double *ptctl=NULL;
        char mess[255];
        if (strlen(nomfichierpoint)>0)
          {
            FILE* in=fopen(nomfichierpoint,"rt");
            if (in!=NULL)
              {
              app.affiche((char*)"Lecture du fichier de points de controles");
              char mess[255];
              char* res=fgets(mess,255,in);
              sscanf(mess,"%d",&nbptctl);
              ptctl=new double[4*nbptctl];
              for (int i=0;i<nbptctl;i++)
              {
                char* res=fgets(mess,255,in);
                sscanf(mess,"%lf %lf %lf %lf",ptctl+4*i,ptctl+4*i+1,ptctl+4*i+2,ptctl+4*i+3);
              }
              fclose(in);
              }
           }
        FCT_GENERATEUR_FRONTIERE carte(gest,eng,eps,ind1,ind2,ptctl,nbptctl);
        app.affiche((char*)"Construction de la carte");
        MG_GROUPE_TOPOLOGIQUE* mggt=NULL;
        if (numgroupe!=0)
        {
            MG_GEOMETRIE* geo=gest.get_mg_geometrieid(numgeo);
            mggt=geo->get_mg_groupe_topologiqueid(numgroupe);
        }
        carte.construit(fechan,festim,mggt);
        app.affiche((char*)"Enregistrement de la carte");
        carte.enregistrer(nomfichiersortie);
        if (nbptctl!=0) delete [] ptctl;
        app.affiche((char*)"Fin");
    }
if (app.get_action()==19)
    {
        char nomfichier[500];
        char nomfichiersortie[500];
        char nomfichierpoint[500];
        double eng;
        double eps;
        double ind1,ind2;
        int fechan;
        int numgeo;
        int nx,ny,nz;
        app.recupere_parametre_string(1,nomfichier,(char*)"Fichier entrant manquant");
        app.recupere_parametre_string(2,nomfichiersortie,(char*)"Fichier de sortie manquant");
        app.recupere_parametre_double(4,eng,(char*)"Valeur de l'écart nodal global");
        app.recupere_parametre_double(5,eps,(char*)"Valeur de l'écart entre la géométrie et le maillage");
        app.recupere_parametre_double_avec_defaut(6,ind1,4.);
        app.recupere_parametre_double_avec_defaut(7,ind2,0.25);
        app.recupere_parametre_int_avec_defaut(8,fechan,20);
        app.recupere_parametre_int_avec_defaut(14,numgeo,0);
        app.recupere_parametre_int_avec_defaut(16,nx,10);
        app.recupere_parametre_int_avec_defaut(17,ny,10);
        app.recupere_parametre_int_avec_defaut(18,nz,10);
        if (app.get_erreur()==true) return 0;  
        app.affiche((char*)"Generation de carte de taille basee sur la courbure des frontieres V2");  
        app.affiche((char*)"Lecture du fichier");
        MG_FILE gest(nomfichier);
        char mess[255];
        MG_GESTIONNAIRE gest1;
        char nomsolution[500];
        sprintf(nomsolution,"%sctt.sol",nomfichiersortie);
        app.affiche((char*)"Construction de la carte");
        MG_GEOMETRIE *geo;
        if (numgeo==0) geo=gest.get_mg_geometrie(0); else geo=gest.get_mg_geometrieid(numgeo);
        FCT_TAILLE_FEM_SOLUTION_GENERATEUR_ECHANTILLON_FRONTIERE carte(&gest1,geo,eng,fechan,nx,ny,nz,nomsolution,eps,ind1,ind2);
        carte.construit();
        app.affiche((char*)"Enregistrement de la carte");
        carte.enregistrer(nomfichiersortie);
        app.affiche((char*)"Fin");
    }
if (app.get_action()==10)
    {
        char nomfichier[500];
        char nomfichiersortie[500];
        double eng;
        int fechan,festim;
        int numgroupe,numgeo;
        app.recupere_parametre_string(1,nomfichier,(char*)"Fichier entrant manquant");
        app.recupere_parametre_string(2,nomfichiersortie,(char*)"Fichier de sortie manquant");
        app.recupere_parametre_double(4,eng,(char*)"Valeur de l'écart nodal global");
        app.recupere_parametre_int_avec_defaut(8,fechan,20);
        app.recupere_parametre_int_avec_defaut(9,festim,10);
        app.recupere_parametre_int_avec_defaut(13,numgroupe,0);
        app.recupere_parametre_int_avec_defaut(14,numgeo,0);
        if (app.get_erreur()==true) return 0;  
        app.affiche((char*)"Generation de carte de taille constante V1");
        app.affiche((char*)"Lecture du fichier");
        MG_FILE gest(nomfichier);
        FCT_GENERATEUR_CONSTANTE carte(gest,eng);
        app.affiche((char*)"Construction de la carte");
        MG_GROUPE_TOPOLOGIQUE* mggt=NULL;
        if (numgroupe!=0)
        {
            MG_GEOMETRIE* geo=gest.get_mg_geometrieid(numgeo);
            mggt=geo->get_mg_groupe_topologiqueid(numgroupe);
        }
        carte.construit(fechan,festim,mggt);
        app.affiche((char*)"Enregistrement de la carte");
        carte.enregistrer(nomfichiersortie);
        app.affiche((char*)"Fin");
    }
  if (app.get_action()==11)
    {
        char nomfichier[500];
        char nomfichiersortie[500];
        char nomfichierpoint[500];
        double eng;
        int fechan,festim;
        double ind1,ind2;
        app.recupere_parametre_string(1,nomfichier,(char*)"Fichier entrant manquant");
        app.recupere_parametre_string(2,nomfichiersortie,(char*)"Fichier de sortie manquant");
        app.recupere_parametre_string(3,nomfichierpoint,(char*)"Fichier de points manquant");
        app.recupere_parametre_double(4,eng,(char*)"Valeur de l'écart nodal global");
        app.recupere_parametre_double_avec_defaut(6,ind1,4.);
        app.recupere_parametre_double_avec_defaut(7,ind2,0.25);
        app.recupere_parametre_int_avec_defaut(8,fechan,20);
        app.recupere_parametre_int_avec_defaut(9,festim,10);
        if (app.get_erreur()==true) return 0;  
        app.affiche((char*)"Generation de carte de taille");
        app.affiche((char*)"   a partir du fichier d'estimation V1");
        FCT_GENERATEUR_FICHIER carte(nomfichierpoint,eng,ind1,ind2);
        app.affiche((char*)"Construction de la carte");
        carte.construit(fechan,festim);
        app.affiche((char*)"Enregistrement de la carte");
        carte.enregistrer(nomfichiersortie);
        app.affiche((char*)"Fin");
    }
  if (app.get_action()==20)
    {
        char nomfichier[500];
        char nomfichiersortie[500];
        char nomfichierpoint[500];
        double eng;
        int fechan,nx,ny,nz,numgeo;
        double ind1,ind2;
        app.recupere_parametre_string(1,nomfichier,(char*)"Fichier entrant manquant");
        app.recupere_parametre_string(2,nomfichiersortie,(char*)"Fichier de sortie manquant");
        app.recupere_parametre_string(3,nomfichierpoint,(char*)"Fichier de points manquant");
        app.recupere_parametre_double(4,eng,(char*)"Valeur de l'écart nodal global");
        app.recupere_parametre_double_avec_defaut(6,ind1,4.);
        app.recupere_parametre_double_avec_defaut(7,ind2,0.25);
        app.recupere_parametre_int_avec_defaut(8,fechan,20);
        app.recupere_parametre_int_avec_defaut(14,numgeo,0);
        app.recupere_parametre_int_avec_defaut(16,nx,10);
        app.recupere_parametre_int_avec_defaut(17,ny,10);
        app.recupere_parametre_int_avec_defaut(18,nz,10);
        if (app.get_erreur()==true) return 0;  
        app.affiche((char*)"Generation de carte de taille");
        app.affiche((char*)"   a partir du fichier d'estimation V2");
        app.affiche((char*)"Lecture du fichier");
        MG_FILE gest(nomfichier);
        MG_GESTIONNAIRE gest1;
        char nomsolution[500];
        sprintf(nomsolution,"%sctt.sol",nomfichiersortie);
        MG_GEOMETRIE *geo;
        if (numgeo==0) geo=gest.get_mg_geometrie(0); else geo=gest.get_mg_geometrieid(numgeo);
        FCT_TAILLE_FEM_SOLUTION_GENERATEUR_FICHIER carte(&gest1,geo,eng,fechan,nx,ny,nz,nomsolution,nomfichierpoint,ind1,ind2);
        app.affiche((char*)"Construction de la carte");
        carte.construit();
        app.affiche((char*)"Enregistrement de la carte");
        carte.enregistrer(nomfichiersortie);
        app.affiche((char*)"Fin");
    }
if (app.get_action()==15)
    {
        char nomfichier[500];
        char nomfichiersortie[500];
        double eng;
        int fechan,nx,ny,nz;
        int numgeo;
        app.recupere_parametre_string(1,nomfichier,(char*)"Fichier entrant manquant");
        app.recupere_parametre_string(2,nomfichiersortie,(char*)"Fichier de sortie manquant");
        app.recupere_parametre_double(4,eng,(char*)"Valeur de l'écart nodal global");
        app.recupere_parametre_int_avec_defaut(8,fechan,20);
        app.recupere_parametre_int_avec_defaut(14,numgeo,0);
        app.recupere_parametre_int_avec_defaut(16,nx,10);
        app.recupere_parametre_int_avec_defaut(17,ny,10);
        app.recupere_parametre_int_avec_defaut(18,nz,10);
        if (app.get_erreur()==true) return 0;  
        app.affiche((char*)"Generation de carte de taille constante V2");
        app.affiche((char*)"Lecture du fichier");
        MG_FILE gest(nomfichier);
        MG_GEOMETRIE* geo;
        if (numgeo==0) geo=gest.get_mg_geometrie(0); else geo=gest.get_mg_geometrieid(numgeo);
        MG_GESTIONNAIRE gest1;
        char nomsolution[500];
        sprintf(nomsolution,"%sctt.sol",nomfichiersortie);
        app.affiche((char*)"Construction de la carte");
        FCT_TAILLE_FEM_SOLUTION_GENERATEUR_CONSTANT carte(&gest1,geo,eng,fechan,nx,ny,nz,nomsolution);
        carte.construit();
        app.affiche((char*)"Enregistrement de la carte");
        carte.enregistrer(nomfichiersortie);
        app.affiche((char*)"Fin");
    }
if (app.get_action()==21)
    {
        char nomfichier[500];
        char nomfichiersortie[500];
        double eng;
        int fechan,nx,ny,nz;
        int nummai;
        app.recupere_parametre_string(1,nomfichier,(char*)"Fichier entrant manquant");
        app.recupere_parametre_string(2,nomfichiersortie,(char*)"Fichier de sortie manquant");
        app.recupere_parametre_double(4,eng,(char*)"Valeur de l'écart nodal global");
        app.recupere_parametre_int_avec_defaut(8,fechan,20);
        app.recupere_parametre_int_avec_defaut(22,nummai,0);
        app.recupere_parametre_int_avec_defaut(16,nx,10);
        app.recupere_parametre_int_avec_defaut(17,ny,10);
        app.recupere_parametre_int_avec_defaut(18,nz,10);
        if (app.get_erreur()==true) return 0;  
        app.affiche((char*)"Generation de carte de taille constante V2");
        app.affiche((char*)"Lecture du fichier");
        MG_FILE gest(nomfichier);
        MG_MAILLAGE* mai;
        if (nummai==0) mai=gest.get_mg_maillage(0); else mai=gest.get_mg_maillageid(nummai);
        MG_GESTIONNAIRE gest1;
        char nomsolution[500];
        sprintf(nomsolution,"%sctt.sol",nomfichiersortie);
        app.affiche((char*)"Construction de la carte");
        FCT_TAILLE_FEM_SOLUTION_GENERATEUR_CONSTANT carte(&gest1,mai,eng,fechan,nx,ny,nz,nomsolution);
        carte.construit();
        app.affiche((char*)"Enregistrement de la carte");
        carte.enregistrer(nomfichiersortie);
        app.affiche((char*)"Fin");
    }

    return 0;
}


#pragma package(smart_init)
Revision:	588
Committed:	Mon Nov 3 22:20:27 2014 UTC (10 years, 9 months ago) by sattarpa
File size:	18102 byte(s)
Log Message:	resolve a bug in reading a mesh without geometry, and also resolve bugs in point insertion and stl refine.
#	User	Rev	Content
1	francois	283	//------------------------------------------------------------
2			//------------------------------------------------------------
3			// MAGiC
4			// Jean Christophe Cuilli�re et Vincent FRANCOIS
5			// D�partement de G�nie M�canique - UQTR
6			//------------------------------------------------------------
7			// Le projet MAGIC est un projet de recherche du d�partement
8			// de g�nie m�canique de l'Universit� du Qu�bec �
9			// Trois Rivi�res
10			// Les librairies ne peuvent �tre utilis�es 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
24
25
26
27			#include "gestionversion.h"
28			#ifdef WINDOWS_VERSION
29			#include "fenetre.h"
30			#endif
31
32			#pragma hdrstop
33			#include "mg_file.h"
34			#include "mg_export.h"
35			#include "fct_generateur_frontiere.h"
36			#include "fct_generateur_constante.h"
37			#include "fct_generateur_fichier.h"
38	francois	469	#include "fct_taille_fem_solution_generateur_constant.h"
39	francois	523	#include "fct_taille_fem_solution_generateur_echantillon_frontiere.h"
40			#include "fct_taille_fem_solution_generateur_fichier.h"
41	francois	432	#include "magic_application.h"
42	francois	283	#include <string.h>
43			//---------------------------------------------------------------------------
44
45
46
47
48			int main(int argc,char **argv)
49			{
50	francois	432	std::vector<MAGIC_PARAMETRE_APPLICATION> lst;
51	francois	523	MAGIC_PARAMETRE_APPLICATION p1(3,(char)"-fichierpoint",MAGIC_PARAMETRE_APPLICATION::STRING,(char)"Fichier d'échantillonage");
52	francois	432	MAGIC_PARAMETRE_APPLICATION p2(4,(char)"-dg",MAGIC_PARAMETRE_APPLICATION::DOUBLE,(char)"Densité globale");
53			MAGIC_PARAMETRE_APPLICATION p3(5,(char)"-eps",MAGIC_PARAMETRE_APPLICATION::DOUBLE,(char)"Écart entre le maillage et la géométrie");
54			MAGIC_PARAMETRE_APPLICATION p4(6,(char)"-taillezone",MAGIC_PARAMETRE_APPLICATION::DOUBLE,(char)"Facteur pour la distance de variation de la carte autour de l'échantillon - Defaut 4");
55			MAGIC_PARAMETRE_APPLICATION p5(7,(char)"-tauxtransfert",MAGIC_PARAMETRE_APPLICATION::DOUBLE,(char)"Poid de l'échantillon virtuel place en dehors de la zone de raffinement - Defaut 0.25");
56			MAGIC_PARAMETRE_APPLICATION p6(8,(char)"-fechantillonnage",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char)"Nombre d'échantillon par entité - Defaut 20");
57			MAGIC_PARAMETRE_APPLICATION p7(9,(char)"-festimation",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char)"Nombre de cellule dans une direction pour la grille virtuelle - Defaut 10");
58			MAGIC_PARAMETRE_APPLICATION p8(10,(char)"-constante",MAGIC_PARAMETRE_APPLICATION::ACTION,(char)"Carte de taille constante");
59			MAGIC_PARAMETRE_APPLICATION p9(11,(char)"-fichierechan",MAGIC_PARAMETRE_APPLICATION::ACTION,(char)"Carte de taille par un echantillonnage dans un fichier exterieur");
60			MAGIC_PARAMETRE_APPLICATION p10(12,(char)"-courbure",MAGIC_PARAMETRE_APPLICATION::ACTION,(char)"Carte de taille par approximation des courbures");
61			MAGIC_PARAMETRE_APPLICATION p11(13,(char)"-groupe",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char)"Groupe topologique utilise pour faire le calcul - Defaut 0");
62			MAGIC_PARAMETRE_APPLICATION p12(14,(char)"-numgeo",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char)"Géométrie du groupe topologique - Defaut 0");
63	francois	469	MAGIC_PARAMETRE_APPLICATION p13(15,(char)"-constante2",MAGIC_PARAMETRE_APPLICATION::ACTION,(char)"Carte de taille constante définie par une FEM_SOLUTION");
64			MAGIC_PARAMETRE_APPLICATION p14(16,(char)"-nx",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char)"Nombre d'intervals sur x - Defaut 10");
65			MAGIC_PARAMETRE_APPLICATION p15(17,(char)"-ny",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char)"Nombre d'intervals sur y - Defaut 10");
66			MAGIC_PARAMETRE_APPLICATION p16(18,(char)"-nz",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char)"Nombre d'intervals sur z - Defaut 10");
67	francois	523	MAGIC_PARAMETRE_APPLICATION p19(19,(char)"-courbure2",MAGIC_PARAMETRE_APPLICATION::ACTION,(char)"Carte de taille par approximation des courbures définie par une FEM_SOLUTION");
68			MAGIC_PARAMETRE_APPLICATION p18(20,(char)"-fichierechan2",MAGIC_PARAMETRE_APPLICATION::ACTION,(char)"Carte de taille par un echantillonnage dans un fichier exterieur définie par une FEM_SOLUTION");
69	sattarpa	545	MAGIC_PARAMETRE_APPLICATION p20(21,(char)"-constante2mai",MAGIC_PARAMETRE_APPLICATION::ACTION,(char)"Carte de taille constante définie par une FEM_SOLUTION basé sur maillage");
70			MAGIC_PARAMETRE_APPLICATION p21(22,(char)"-nummai",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char)"Numéro de maillage - Defaut le premier");
71	francois	432	p8.ajouter_dependance(1);p8.ajouter_dependance(2);p8.ajouter_dependance(4);p8.ajouter_dependance(8);p8.ajouter_dependance(9);p8.ajouter_dependance(13);p8.ajouter_dependance(14);
72			p9.ajouter_dependance(1);p9.ajouter_dependance(2);p9.ajouter_dependance(3);p9.ajouter_dependance(4);p9.ajouter_dependance(6);p9.ajouter_dependance(7);p9.ajouter_dependance(8);p9.ajouter_dependance(9);
73	francois	523	p10.ajouter_dependance(1);p10.ajouter_dependance(2);p10.ajouter_dependance(3);p10.ajouter_dependance(4);p10.ajouter_dependance(5);p10.ajouter_dependance(6);p10.ajouter_dependance(7);p10.ajouter_dependance(8);p10.ajouter_dependance(9);p10.ajouter_dependance(13);p10.ajouter_dependance(14);
74			p13.ajouter_dependance(1);p13.ajouter_dependance(2);p13.ajouter_dependance(4);p13.ajouter_dependance(8);p13.ajouter_dependance(14);p13.ajouter_dependance(16);p13.ajouter_dependance(17);p13.ajouter_dependance(18);
75			p18.ajouter_dependance(1);p18.ajouter_dependance(2);p18.ajouter_dependance(3);p18.ajouter_dependance(4);p18.ajouter_dependance(6);p18.ajouter_dependance(7);p18.ajouter_dependance(8);p18.ajouter_dependance(14);p18.ajouter_dependance(16);p18.ajouter_dependance(17);p18.ajouter_dependance(18);
76			p19.ajouter_dependance(1);p19.ajouter_dependance(2);p19.ajouter_dependance(4);p19.ajouter_dependance(5);p19.ajouter_dependance(6);p19.ajouter_dependance(7);p19.ajouter_dependance(8);p19.ajouter_dependance(14);p19.ajouter_dependance(16);p19.ajouter_dependance(17);p19.ajouter_dependance(18);
77	sattarpa	545	p20.ajouter_dependance(1);p20.ajouter_dependance(2);p20.ajouter_dependance(4);p20.ajouter_dependance(8);p20.ajouter_dependance(22);p20.ajouter_dependance(16);p20.ajouter_dependance(17);p20.ajouter_dependance(18);
78	francois	523
79	francois	432	lst.push_back(p1);
80			lst.push_back(p2);
81			lst.push_back(p3);
82			lst.push_back(p4);
83			lst.push_back(p5);
84			lst.push_back(p6);
85			lst.push_back(p7);
86			lst.push_back(p8);
87			lst.push_back(p9);
88			lst.push_back(p10);
89			lst.push_back(p11);
90			lst.push_back(p12);
91	francois	469	lst.push_back(p13);
92			lst.push_back(p14);
93			lst.push_back(p15);
94			lst.push_back(p16);
95	francois	523	lst.push_back(p19);
96			lst.push_back(p18);
97	sattarpa	545	lst.push_back(p20);
98			lst.push_back(p21);
99	francois	432	MAGIC_APPLICATION app((char*)"Generation de carte de taille",argc,argv,lst,true,true);
100			if (app.get_erreur()==true) return 0;
101
102
103			if (app.get_action()==12)
104	francois	283	{
105	francois	432	char nomfichier[500];
106			char nomfichiersortie[500];
107			char nomfichierpoint[500];
108			double eng;
109			double eps;
110			double ind1,ind2;
111			int fechan,festim;
112			int numgroupe,numgeo;
113			app.recupere_parametre_string(1,nomfichier,(char*)"Fichier entrant manquant");
114			app.recupere_parametre_string(2,nomfichiersortie,(char*)"Fichier de sortie manquant");
115			app.recupere_parametre_string_avec_defaut(3,nomfichierpoint,(char*)"");
116			app.recupere_parametre_double(4,eng,(char*)"Valeur de l'écart nodal global");
117			app.recupere_parametre_double(5,eps,(char*)"Valeur de l'écart entre la géométrie et le maillage");
118			app.recupere_parametre_double_avec_defaut(6,ind1,4.);
119			app.recupere_parametre_double_avec_defaut(7,ind2,0.25);
120			app.recupere_parametre_int_avec_defaut(8,fechan,20);
121			app.recupere_parametre_int_avec_defaut(9,festim,10);
122			app.recupere_parametre_int_avec_defaut(13,numgroupe,0);
123			app.recupere_parametre_int_avec_defaut(14,numgeo,0);
124			if (app.get_erreur()==true) return 0;
125	francois	529	app.affiche((char*)"Generation de carte de taille basee sur la courbure des frontieres V1");
126	francois	432	app.affiche((char*)"Lecture du fichier");
127	francois	283	MG_FILE gest(nomfichier);
128			int nbptctl=0;
129			double *ptctl=NULL;
130			char mess[255];
131	francois	432	if (strlen(nomfichierpoint)>0)
132			{
133			FILE* in=fopen(nomfichierpoint,"rt");
134			if (in!=NULL)
135			{
136			app.affiche((char*)"Lecture du fichier de points de controles");
137			char mess[255];
138			char* res=fgets(mess,255,in);
139			sscanf(mess,"%d",&nbptctl);
140			ptctl=new double[4*nbptctl];
141			for (int i=0;i<nbptctl;i++)
142			{
143	francois	283	char* res=fgets(mess,255,in);
144			sscanf(mess,"%lf %lf %lf %lf",ptctl+4i,ptctl+4i+1,ptctl+4i+2,ptctl+4i+3);
145	francois	432	}
146			fclose(in);
147			}
148			}
149	francois	283	FCT_GENERATEUR_FRONTIERE carte(gest,eng,eps,ind1,ind2,ptctl,nbptctl);
150	francois	432	app.affiche((char*)"Construction de la carte");
151	francois	283	MG_GROUPE_TOPOLOGIQUE* mggt=NULL;
152			if (numgroupe!=0)
153			{
154			MG_GEOMETRIE* geo=gest.get_mg_geometrieid(numgeo);
155			mggt=geo->get_mg_groupe_topologiqueid(numgroupe);
156			}
157			carte.construit(fechan,festim,mggt);
158	francois	432	app.affiche((char*)"Enregistrement de la carte");
159	francois	283	carte.enregistrer(nomfichiersortie);
160			if (nbptctl!=0) delete [] ptctl;
161	francois	432	app.affiche((char*)"Fin");
162	francois	283	}
163	francois	523	if (app.get_action()==19)
164			{
165			char nomfichier[500];
166			char nomfichiersortie[500];
167			char nomfichierpoint[500];
168			double eng;
169			double eps;
170			double ind1,ind2;
171			int fechan;
172			int numgeo;
173			int nx,ny,nz;
174			app.recupere_parametre_string(1,nomfichier,(char*)"Fichier entrant manquant");
175			app.recupere_parametre_string(2,nomfichiersortie,(char*)"Fichier de sortie manquant");
176			app.recupere_parametre_double(4,eng,(char*)"Valeur de l'écart nodal global");
177			app.recupere_parametre_double(5,eps,(char*)"Valeur de l'écart entre la géométrie et le maillage");
178			app.recupere_parametre_double_avec_defaut(6,ind1,4.);
179			app.recupere_parametre_double_avec_defaut(7,ind2,0.25);
180			app.recupere_parametre_int_avec_defaut(8,fechan,20);
181			app.recupere_parametre_int_avec_defaut(14,numgeo,0);
182			app.recupere_parametre_int_avec_defaut(16,nx,10);
183			app.recupere_parametre_int_avec_defaut(17,ny,10);
184			app.recupere_parametre_int_avec_defaut(18,nz,10);
185			if (app.get_erreur()==true) return 0;
186	francois	529	app.affiche((char*)"Generation de carte de taille basee sur la courbure des frontieres V2");
187	francois	523	app.affiche((char*)"Lecture du fichier");
188			MG_FILE gest(nomfichier);
189			char mess[255];
190			MG_GESTIONNAIRE gest1;
191			char nomsolution[500];
192			sprintf(nomsolution,"%sctt.sol",nomfichiersortie);
193			app.affiche((char*)"Construction de la carte");
194			MG_GEOMETRIE *geo;
195			if (numgeo==0) geo=gest.get_mg_geometrie(0); else geo=gest.get_mg_geometrieid(numgeo);
196			FCT_TAILLE_FEM_SOLUTION_GENERATEUR_ECHANTILLON_FRONTIERE carte(&gest1,geo,eng,fechan,nx,ny,nz,nomsolution,eps,ind1,ind2);
197			carte.construit();
198			app.affiche((char*)"Enregistrement de la carte");
199			carte.enregistrer(nomfichiersortie);
200			app.affiche((char*)"Fin");
201			}
202	francois	469	if (app.get_action()==10)
203	francois	283	{
204	francois	432	char nomfichier[500];
205			char nomfichiersortie[500];
206			double eng;
207			int fechan,festim;
208			int numgroupe,numgeo;
209			app.recupere_parametre_string(1,nomfichier,(char*)"Fichier entrant manquant");
210			app.recupere_parametre_string(2,nomfichiersortie,(char*)"Fichier de sortie manquant");
211			app.recupere_parametre_double(4,eng,(char*)"Valeur de l'écart nodal global");
212			app.recupere_parametre_int_avec_defaut(8,fechan,20);
213			app.recupere_parametre_int_avec_defaut(9,festim,10);
214			app.recupere_parametre_int_avec_defaut(13,numgroupe,0);
215			app.recupere_parametre_int_avec_defaut(14,numgeo,0);
216			if (app.get_erreur()==true) return 0;
217	francois	529	app.affiche((char*)"Generation de carte de taille constante V1");
218	francois	432	app.affiche((char*)"Lecture du fichier");
219	francois	283	MG_FILE gest(nomfichier);
220			FCT_GENERATEUR_CONSTANTE carte(gest,eng);
221	francois	432	app.affiche((char*)"Construction de la carte");
222	francois	283	MG_GROUPE_TOPOLOGIQUE* mggt=NULL;
223			if (numgroupe!=0)
224			{
225			MG_GEOMETRIE* geo=gest.get_mg_geometrieid(numgeo);
226			mggt=geo->get_mg_groupe_topologiqueid(numgroupe);
227			}
228			carte.construit(fechan,festim,mggt);
229	francois	432	app.affiche((char*)"Enregistrement de la carte");
230	francois	283	carte.enregistrer(nomfichiersortie);
231	francois	432	app.affiche((char*)"Fin");
232	francois	283	}
233	francois	432	if (app.get_action()==11)
234	francois	283	{
235	francois	432	char nomfichier[500];
236			char nomfichiersortie[500];
237			char nomfichierpoint[500];
238			double eng;
239			int fechan,festim;
240			double ind1,ind2;
241			app.recupere_parametre_string(1,nomfichier,(char*)"Fichier entrant manquant");
242			app.recupere_parametre_string(2,nomfichiersortie,(char*)"Fichier de sortie manquant");
243			app.recupere_parametre_string(3,nomfichierpoint,(char*)"Fichier de points manquant");
244			app.recupere_parametre_double(4,eng,(char*)"Valeur de l'écart nodal global");
245			app.recupere_parametre_double_avec_defaut(6,ind1,4.);
246			app.recupere_parametre_double_avec_defaut(7,ind2,0.25);
247			app.recupere_parametre_int_avec_defaut(8,fechan,20);
248			app.recupere_parametre_int_avec_defaut(9,festim,10);
249			if (app.get_erreur()==true) return 0;
250			app.affiche((char*)"Generation de carte de taille");
251	francois	529	app.affiche((char*)" a partir du fichier d'estimation V1");
252	francois	394	FCT_GENERATEUR_FICHIER carte(nomfichierpoint,eng,ind1,ind2);
253	francois	432	app.affiche((char*)"Construction de la carte");
254	francois	283	carte.construit(fechan,festim);
255	francois	432	app.affiche((char*)"Enregistrement de la carte");
256	francois	283	carte.enregistrer(nomfichiersortie);
257	francois	432	app.affiche((char*)"Fin");
258	francois	283	}
259	francois	523	if (app.get_action()==20)
260			{
261			char nomfichier[500];
262			char nomfichiersortie[500];
263			char nomfichierpoint[500];
264			double eng;
265			int fechan,nx,ny,nz,numgeo;
266			double ind1,ind2;
267			app.recupere_parametre_string(1,nomfichier,(char*)"Fichier entrant manquant");
268			app.recupere_parametre_string(2,nomfichiersortie,(char*)"Fichier de sortie manquant");
269			app.recupere_parametre_string(3,nomfichierpoint,(char*)"Fichier de points manquant");
270			app.recupere_parametre_double(4,eng,(char*)"Valeur de l'écart nodal global");
271			app.recupere_parametre_double_avec_defaut(6,ind1,4.);
272			app.recupere_parametre_double_avec_defaut(7,ind2,0.25);
273			app.recupere_parametre_int_avec_defaut(8,fechan,20);
274			app.recupere_parametre_int_avec_defaut(14,numgeo,0);
275			app.recupere_parametre_int_avec_defaut(16,nx,10);
276			app.recupere_parametre_int_avec_defaut(17,ny,10);
277			app.recupere_parametre_int_avec_defaut(18,nz,10);
278			if (app.get_erreur()==true) return 0;
279			app.affiche((char*)"Generation de carte de taille");
280	francois	529	app.affiche((char*)" a partir du fichier d'estimation V2");
281	francois	523	app.affiche((char*)"Lecture du fichier");
282			MG_FILE gest(nomfichier);
283			MG_GESTIONNAIRE gest1;
284			char nomsolution[500];
285			sprintf(nomsolution,"%sctt.sol",nomfichiersortie);
286			MG_GEOMETRIE *geo;
287			if (numgeo==0) geo=gest.get_mg_geometrie(0); else geo=gest.get_mg_geometrieid(numgeo);
288			FCT_TAILLE_FEM_SOLUTION_GENERATEUR_FICHIER carte(&gest1,geo,eng,fechan,nx,ny,nz,nomsolution,nomfichierpoint,ind1,ind2);
289			app.affiche((char*)"Construction de la carte");
290			carte.construit();
291			app.affiche((char*)"Enregistrement de la carte");
292			carte.enregistrer(nomfichiersortie);
293			app.affiche((char*)"Fin");
294			}
295	francois	469	if (app.get_action()==15)
296			{
297			char nomfichier[500];
298			char nomfichiersortie[500];
299			double eng;
300			int fechan,nx,ny,nz;
301			int numgeo;
302			app.recupere_parametre_string(1,nomfichier,(char*)"Fichier entrant manquant");
303			app.recupere_parametre_string(2,nomfichiersortie,(char*)"Fichier de sortie manquant");
304			app.recupere_parametre_double(4,eng,(char*)"Valeur de l'écart nodal global");
305			app.recupere_parametre_int_avec_defaut(8,fechan,20);
306			app.recupere_parametre_int_avec_defaut(14,numgeo,0);
307			app.recupere_parametre_int_avec_defaut(16,nx,10);
308			app.recupere_parametre_int_avec_defaut(17,ny,10);
309			app.recupere_parametre_int_avec_defaut(18,nz,10);
310			if (app.get_erreur()==true) return 0;
311	francois	529	app.affiche((char*)"Generation de carte de taille constante V2");
312	francois	469	app.affiche((char*)"Lecture du fichier");
313			MG_FILE gest(nomfichier);
314			MG_GEOMETRIE* geo;
315			if (numgeo==0) geo=gest.get_mg_geometrie(0); else geo=gest.get_mg_geometrieid(numgeo);
316			MG_GESTIONNAIRE gest1;
317	francois	523	char nomsolution[500];
318			sprintf(nomsolution,"%sctt.sol",nomfichiersortie);
319	francois	469	app.affiche((char*)"Construction de la carte");
320	francois	523	FCT_TAILLE_FEM_SOLUTION_GENERATEUR_CONSTANT carte(&gest1,geo,eng,fechan,nx,ny,nz,nomsolution);
321			carte.construit();
322	francois	469	app.affiche((char*)"Enregistrement de la carte");
323			carte.enregistrer(nomfichiersortie);
324			app.affiche((char*)"Fin");
325			}
326	sattarpa	545	if (app.get_action()==21)
327			{
328			char nomfichier[500];
329			char nomfichiersortie[500];
330			double eng;
331			int fechan,nx,ny,nz;
332			int nummai;
333			app.recupere_parametre_string(1,nomfichier,(char*)"Fichier entrant manquant");
334			app.recupere_parametre_string(2,nomfichiersortie,(char*)"Fichier de sortie manquant");
335			app.recupere_parametre_double(4,eng,(char*)"Valeur de l'écart nodal global");
336			app.recupere_parametre_int_avec_defaut(8,fechan,20);
337			app.recupere_parametre_int_avec_defaut(22,nummai,0);
338			app.recupere_parametre_int_avec_defaut(16,nx,10);
339			app.recupere_parametre_int_avec_defaut(17,ny,10);
340			app.recupere_parametre_int_avec_defaut(18,nz,10);
341			if (app.get_erreur()==true) return 0;
342			app.affiche((char*)"Generation de carte de taille constante V2");
343			app.affiche((char*)"Lecture du fichier");
344			MG_FILE gest(nomfichier);
345			MG_MAILLAGE* mai;
346	sattarpa	588	if (nummai==0) mai=gest.get_mg_maillage(0); else mai=gest.get_mg_maillageid(nummai);
347	sattarpa	545	MG_GESTIONNAIRE gest1;
348			char nomsolution[500];
349			sprintf(nomsolution,"%sctt.sol",nomfichiersortie);
350			app.affiche((char*)"Construction de la carte");
351			FCT_TAILLE_FEM_SOLUTION_GENERATEUR_CONSTANT carte(&gest1,mai,eng,fechan,nx,ny,nz,nomsolution);
352			carte.construit();
353			app.affiche((char*)"Enregistrement de la carte");
354			carte.enregistrer(nomfichiersortie);
355			app.affiche((char*)"Fin");
356			}
357	francois	283
358			return 0;
359			}
360
361
362
363			#pragma package(smart_init)