mgoperation/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


#include <string.h>
#include "mg_file.h"
#include "mg_export.h"
#include "fct_export.h"
#include "mg_import.h"
#include "step_import.h"
#include "acis_import.h"
#include "step_import.h"
#include "stl_import.h"
#include "occ_import.h"
#include "mc_gestionnaire.h"
#include "vct_multi_modele.h"
#include "mgaster.h"
#include <ctype.h>
#include <math.h>
#include "mailleur_analyse.h"
#include "magic_application.h"
#include "fct_taille_fem_solution.h"
#include "fem_maillage_quadratique_outils.h"
#include "fem_maillage_outils.h"
#include "mg_table.h"
//---------------------------------------------------------------------------


void minuscule(char * mess)
{
    int nb=strlen(mess);
    for (int i=0;i<nb;i++)
        mess[i]=tolower(mess[i]);
}


int main(int argc,char **argv)
{
std::vector<MAGIC_PARAMETRE_APPLICATION> lst;
MAGIC_PARAMETRE_APPLICATION p1(2,(char*)"-gmsh",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Conversion du fichier MAGiC en gmsh");  
MAGIC_PARAMETRE_APPLICATION p2(3,(char*)"-nodeforme",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Ne convertit pas les deformes des maillages FEM");  
MAGIC_PARAMETRE_APPLICATION p3(4,(char*)"-copie",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Copie un fichier MAGiC en un autre");  
MAGIC_PARAMETRE_APPLICATION p4(5,(char*)"-out",MAGIC_PARAMETRE_APPLICATION::STRING,(char*)"Fichier de sortie");  
MAGIC_PARAMETRE_APPLICATION p5(6,(char*)"-aster",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Conversion MAGiC en format ASTER (uniquement calcul mecanique)");  
MAGIC_PARAMETRE_APPLICATION p6(7,(char*)"-coderesu",MAGIC_PARAMETRE_APPLICATION::STRING,(char*)"Choix des valeurs à calculer par EF - Defaut 11111111 - |Surface moyenne|Surface inferieure|Surface superieur|Deplacement|Contraintes|Deformations|Contraintes equivalentes|Energie|");  
MAGIC_PARAMETRE_APPLICATION p7(8,(char*)"-carte",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Conversion carte de taille ancien format en format gmsh");  
MAGIC_PARAMETRE_APPLICATION p8(9,(char*)"-step",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Conversion de fichier step en format MAGiC");  
MAGIC_PARAMETRE_APPLICATION p9(10,(char*)"-unite",MAGIC_PARAMETRE_APPLICATION::DOUBLE,(char*)"Facteur d'unité pour obtenir des métres - Defaut 1.");  
MAGIC_PARAMETRE_APPLICATION p10(11,(char*)"-stepocc",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Conversion de fichier step en format MAGiC via opencascade");  
MAGIC_PARAMETRE_APPLICATION p11(12,(char*)"-sat",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Conversion de fichier sat en format MAGiC");  
MAGIC_PARAMETRE_APPLICATION p12(13,(char*)"-eps",MAGIC_PARAMETRE_APPLICATION::DOUBLE,(char*)"Tolérance pour un double - Defaut 0.000001");  
MAGIC_PARAMETRE_APPLICATION p13(14,(char*)"-epsgeo",MAGIC_PARAMETRE_APPLICATION::DOUBLE,(char*)"Distance max entre la triangulation et la géométrie - Defaut 1.");  
MAGIC_PARAMETRE_APPLICATION p14(15,(char*)"-importstl",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Importation de la triangulation STL");  
MAGIC_PARAMETRE_APPLICATION p15(16,(char*)"-importtriangulation",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Importation de la triangulation STL avec noeuds fusionnés");  
MAGIC_PARAMETRE_APPLICATION p16(17,(char*)"-occ",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Conversion de fichier opencascade en format MAGiC");  
MAGIC_PARAMETRE_APPLICATION p17(18,(char*)"-fem",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Création d'un maillage fem");  
MAGIC_PARAMETRE_APPLICATION p18(19,(char*)"-nummai",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Numéro de maillage - Defaut le premier");  
MAGIC_PARAMETRE_APPLICATION p19(20,(char*)"-degre",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Degre du maillage à créer");  
MAGIC_PARAMETRE_APPLICATION p20(21,(char*)"-astermagic",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Importation des resultats ASTER dans un maillage");  
MAGIC_PARAMETRE_APPLICATION p21(22,(char*)"-fichieraster",MAGIC_PARAMETRE_APPLICATION::STRING,(char*)"Fichier aster");  
MAGIC_PARAMETRE_APPLICATION p22(23,(char*)"-optis",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Conversion de fichier magic en format OPTIS - Fonction esactivée");  
MAGIC_PARAMETRE_APPLICATION p23(24,(char*)"-optismagic",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Importation des résultats OPTIS");  
MAGIC_PARAMETRE_APPLICATION p24(25,(char*)"-fichieroptis",MAGIC_PARAMETRE_APPLICATION::STRING,(char*)"Fichier optis");  
MAGIC_PARAMETRE_APPLICATION p25(26,(char*)"-fichierflottant",MAGIC_PARAMETRE_APPLICATION::STRING,(char*)"Fichier flottant - Optionnel");  
MAGIC_PARAMETRE_APPLICATION p26(27,(char*)"-seuil",MAGIC_PARAMETRE_APPLICATION::DOUBLE,(char*)"Seuil de densité pour conserver les éléments");  
MAGIC_PARAMETRE_APPLICATION p27(28,(char*)"-niveau",MAGIC_PARAMETRE_APPLICATION::DOUBLE,(char*)"Facteur pour le seuil (parametre incertain)");  
MAGIC_PARAMETRE_APPLICATION p28(29,(char*)"-selectgr",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Selection de groupe topologique");  
MAGIC_PARAMETRE_APPLICATION p29(30,(char*)"-deselectgr",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Deselection de groupe topologique");  
MAGIC_PARAMETRE_APPLICATION p30(31,(char*)"-numgeo",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Numéro de geometrie - Defaut la première");  
MAGIC_PARAMETRE_APPLICATION p31(32,(char*)"-numgt",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Numéro de groupe topologique - Defaut le premier");  
MAGIC_PARAMETRE_APPLICATION p32(33,(char*)"-numentite",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Numéro d'entité à selectionner");  
MAGIC_PARAMETRE_APPLICATION p33(34,(char*)"-expand",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Prend en compte la topologie sous-jacente dans la selection");  
MAGIC_PARAMETRE_APPLICATION p34(35,(char*)"-examine",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Contenu d'un fichier MAGiC. Si un numero de maillage est présent l'action analyse ce maillage");  
MAGIC_PARAMETRE_APPLICATION p35(36,(char*)"-angle",MAGIC_PARAMETRE_APPLICATION::DOUBLE,(char*)"Angle limite pour l'analyse de retournement de triangle - Defaut 0.03");  
MAGIC_PARAMETRE_APPLICATION p36(37,(char*)"-fichierlog",MAGIC_PARAMETRE_APPLICATION::STRING,(char*)"Nom du fichier log");  
MAGIC_PARAMETRE_APPLICATION p37(38,(char*)"-deforme",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Calcul un maillage déforme");  
MAGIC_PARAMETRE_APPLICATION p38(39,(char*)"-fusiongeo",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Fusionne deux géométries");  
MAGIC_PARAMETRE_APPLICATION p39(40,(char*)"-numsol1",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Numéro de solution 1");  
MAGIC_PARAMETRE_APPLICATION p40(41,(char*)"-numsol2",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Numéro de solution 2");  
MAGIC_PARAMETRE_APPLICATION p41(42,(char*)"-numsol3",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Numéro de solution 3");  
MAGIC_PARAMETRE_APPLICATION p42(43,(char*)"-numchamp1",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Numéro de champs de la solution 1");  
MAGIC_PARAMETRE_APPLICATION p43(44,(char*)"-numchamp2",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Numéro de champs de la solution 2");  
MAGIC_PARAMETRE_APPLICATION p44(45,(char*)"-numchamp3",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Numéro de champs de la solution 3");  
MAGIC_PARAMETRE_APPLICATION p45(46,(char*)"-ccf",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Application de conditions aux limites sur une topologie");  
MAGIC_PARAMETRE_APPLICATION p46(47,(char*)"-ccft",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Application de conditions aux limites sur une topologie et sur les topologie sous jacente");  
MAGIC_PARAMETRE_APPLICATION p47(48,(char*)"-ccfhelp",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Liste des conditions aux limites disponibles");  
MAGIC_PARAMETRE_APPLICATION p48(49,(char*)"-valeur",MAGIC_PARAMETRE_APPLICATION::STRING,(char*)"Condition aux limites à appliquer sous la forme id:ccf:valeur(:NS|SU) ou id:ccf:formule:var1,var2(:NS|SU)");  
MAGIC_PARAMETRE_APPLICATION p49(50,(char*)"-oricoque",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Orientation d'une coque");  
MAGIC_PARAMETRE_APPLICATION p50(51,(char*)"-numcoque",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Numéro de la coque - Defaut la premiere");  
MAGIC_PARAMETRE_APPLICATION p51(52,(char*)"-stlmagic",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Importation de fichier STL");  
MAGIC_PARAMETRE_APPLICATION p52(53,(char*)"-epsf",MAGIC_PARAMETRE_APPLICATION::DOUBLE,(char*)"Tolérance pour un float - Defaut 0.0001");  
MAGIC_PARAMETRE_APPLICATION p53(54,(char*)"-calculaster",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Calcul ASTER"); 
MAGIC_PARAMETRE_APPLICATION p54(55,(char*)"-correspondance",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Ajout d'un fichier de correspondance entre GMSH et MAGIC");  
MAGIC_PARAMETRE_APPLICATION p55(56,(char*)"-carte2",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Conversion carte de taille nouvelle version en format gmsh");  
MAGIC_PARAMETRE_APPLICATION p56(57,(char*)"-stlout",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Conversion du fichier MAGiC en stl");  
MAGIC_PARAMETRE_APPLICATION p57(58,(char*)"-comparetaille",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Compare la taille d'un maillage avec la taille d'une carte de taille isotrope ancienne version");  
MAGIC_PARAMETRE_APPLICATION p58(59,(char*)"-cartetaille",MAGIC_PARAMETRE_APPLICATION::STRING,(char*)"Nom de la carte de taille a fournir");  
MAGIC_PARAMETRE_APPLICATION p59(60,(char*)"-comparetaille2",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Compare la taille d'un maillage avec la taille d'une carte de taille isotrope nouvelle version");  
MAGIC_PARAMETRE_APPLICATION p60(61,(char*)"-coefdef",MAGIC_PARAMETRE_APPLICATION::DOUBLE,(char*)"Coefficient de déformé - par defaut 0. = pas de creation de maillage");  
MAGIC_PARAMETRE_APPLICATION p61(62,(char*)"-detachemaillage",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Détache un maillage des ses liens avec la géométrie");  
MAGIC_PARAMETRE_APPLICATION p62(63,(char*)"-elastique",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Calcul aster en elasticité");
MAGIC_PARAMETRE_APPLICATION p63(64,(char*)"-thermique",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Calcul aster en thermique");
MAGIC_PARAMETRE_APPLICATION p64(65,(char*)"-contrainteplane",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Calcul aster en contrainte plane");
MAGIC_PARAMETRE_APPLICATION p65(66,(char*)"-deformationplane",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Calcul aster en deformation plane");
MAGIC_PARAMETRE_APPLICATION p66(67,(char*)"-sonde",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Obtenir un resultat sur une arete dans un fichier csv");
MAGIC_PARAMETRE_APPLICATION p67(68,(char*)"-numsol",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Numero de la solution (la premiere par defaut)");
MAGIC_PARAMETRE_APPLICATION p68(69,(char*)"-numchamp",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Numero du champs de la solution (le premier par defaut)");
MAGIC_PARAMETRE_APPLICATION p69(70,(char*)"-csvout",MAGIC_PARAMETRE_APPLICATION::STRING,(char*)"Nom du fichier csv");
MAGIC_PARAMETRE_APPLICATION p70(71,(char*)"-numare",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Numero de l'arete à sonder (la premiere par default)");
MAGIC_PARAMETRE_APPLICATION p71(72,(char*)"-projsolution",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Projeter une solution sur un autre maillage");
MAGIC_PARAMETRE_APPLICATION p72(73,(char*)"-fichiersolution",MAGIC_PARAMETRE_APPLICATION::STRING,(char*)"Nom du fichier contenant la solution");
MAGIC_PARAMETRE_APPLICATION p73(74,(char*)"-comparesolution",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Comparer deux champs de solution. La reference est la deuxième solution.");
MAGIC_PARAMETRE_APPLICATION p74(75,(char*)"-analyse",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Analyse le maillage quadratique. Defaut non.");
MAGIC_PARAMETRE_APPLICATION p75(76,(char*)"-table",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Extraire des tables de valeurs d'un maillage fem");
MAGIC_PARAMETRE_APPLICATION p76(77,(char*)"-noeud",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Action sur les noeuds");
MAGIC_PARAMETRE_APPLICATION p77(78,(char*)"-element1",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Action sur les element1D");
MAGIC_PARAMETRE_APPLICATION p78(79,(char*)"-element2",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Action sur les element2D");
MAGIC_PARAMETRE_APPLICATION p79(80,(char*)"-element3",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Action sur les element3D");
MAGIC_PARAMETRE_APPLICATION p80(81,(char*)"-numtopo",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Action uniquement sur le numero de topologie mentionnée");
MAGIC_PARAMETRE_APPLICATION p81(82,(char*)"-maigmsh",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Importation d'un maillage de GMSH (-in peut être omis -out peut etre omis si -in n'est pas omis -numgeo le defaut correspond à un non attechement du maillage)");
MAGIC_PARAMETRE_APPLICATION p82(83,(char*)"-fichiergmsh",MAGIC_PARAMETRE_APPLICATION::STRING,(char*)"Nom du fichier gmsh");
MAGIC_PARAMETRE_APPLICATION p83(84,(char*)"-param",MAGIC_PARAMETRE_APPLICATION::STRING,(char*)"Nom du fichier parametre");
MAGIC_PARAMETRE_APPLICATION p84(85,(char*)"-creeparam",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Création d'un fichier parametre");
MAGIC_PARAMETRE_APPLICATION p85(86,(char*)"-importefonction",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Importation d'une fonction dans une géométrie");
MAGIC_PARAMETRE_APPLICATION p86(87,(char*)"-fichier",MAGIC_PARAMETRE_APPLICATION::STRING,(char*)"Nom du fichier à importer");
MAGIC_PARAMETRE_APPLICATION p87(88,(char*)"-elastique_gd",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Calcul aster en elasticite grandes deformations");
MAGIC_PARAMETRE_APPLICATION p88(89,(char*)"-plastique",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Calcul aster en plasticite");
MAGIC_PARAMETRE_APPLICATION p89(90,(char*)"-base",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Base de numerotation des entités code aster - Defaut 10");
MAGIC_PARAMETRE_APPLICATION p90(91,(char*)"-element_gauss1",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Action aux points de gauss sur les element1D");
MAGIC_PARAMETRE_APPLICATION p91(92,(char*)"-element_gauss2",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Action aux points de gauss sur les element2D");
MAGIC_PARAMETRE_APPLICATION p92(93,(char*)"-element_gauss3",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Action aux points de gauss sur les element3D");
MAGIC_PARAMETRE_APPLICATION p93(94,(char*)"-ccfi",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Application de conditions aux limites sur une topologie par identificateur d'une autre entité MAGiC");  
MAGIC_PARAMETRE_APPLICATION p94(95,(char*)"-ccfit",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Application de conditions aux limites sur une topologie et sur les topologie sous jacente par identificateur d'une autre entité MAGiC");  
MAGIC_PARAMETRE_APPLICATION p95(96,(char*)"-ccff",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Application de conditions aux limites sur une topologie sous forme de formule");  
MAGIC_PARAMETRE_APPLICATION p96(97,(char*)"-ccfft",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Application de conditions aux limites sur une topologie et sur les topologie sous forme de formule");  
MAGIC_PARAMETRE_APPLICATION p97(98,(char*)"-opt",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Optimisation de la numerotation des noeuds");  
MAGIC_PARAMETRE_APPLICATION p98(99,(char*)"-rigidite",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Matrice de rigidite elastique par code aster");
MAGIC_PARAMETRE_APPLICATION p99(100,(char*)"-rigiditecontrainteplane",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Matrice de rigidite en contrainte plane par code aster");
MAGIC_PARAMETRE_APPLICATION p100(101,(char*)"-rigiditedeformationplane",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Matrice de rigidite en deformation plane par code aster");
MAGIC_PARAMETRE_APPLICATION p101(102,(char*)"-elementnd",MAGIC_PARAMETRE_APPLICATION::BOOL,(char*)"Action sur les tous les elements (solution de calcul mixte)");
MAGIC_PARAMETRE_APPLICATION p102(103,(char*)"-changefonction",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Modifie une fonction dans une géométrie");
MAGIC_PARAMETRE_APPLICATION p103(104,(char*)"-numfonction",MAGIC_PARAMETRE_APPLICATION::INTEGER,(char*)"Numero d'une fonction dans la géométrie (defaut la premiere)");
MAGIC_PARAMETRE_APPLICATION p104(105,(char*)"-offout",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Conversion du fichier MAGiC en off");  
MAGIC_PARAMETRE_APPLICATION p105(106,(char*)"-mince",MAGIC_PARAMETRE_APPLICATION::ACTION,(char*)"Conversion d'une coquille mince");  
MAGIC_PARAMETRE_APPLICATION p106(107,(char*)"-volcor",MAGIC_PARAMETRE_APPLICATION::STRING,(char*)"Ajout d'une face de base et d'une face en vis a vis pour un volume mince sous la forme idgeo:idvol:idfacebase:idfacevisavis");  


p1.ajouter_dependance(1);p1.ajouter_dependance(3);p1.ajouter_dependance(55);
p3.ajouter_dependance(1);p3.ajouter_dependance(5);
p5.ajouter_dependance(1);p5.ajouter_dependance(7);p5.ajouter_dependance(85);
p7.ajouter_dependance(1);
p8.ajouter_dependance(1);p8.ajouter_dependance(10);
p10.ajouter_dependance(1);p10.ajouter_dependance(13);p10.ajouter_dependance(14);p10.ajouter_dependance(15);;p10.ajouter_dependance(16);
p11.ajouter_dependance(1);p11.ajouter_dependance(10);
p16.ajouter_dependance(1);p16.ajouter_dependance(13);p16.ajouter_dependance(14);p16.ajouter_dependance(15);;p16.ajouter_dependance(16);
p17.ajouter_dependance(1);p17.ajouter_dependance(5);p17.ajouter_dependance(19);p17.ajouter_dependance(20);p17.ajouter_dependance(75);p17.ajouter_dependance(98);
p20.ajouter_dependance(1);p20.ajouter_dependance(19);p20.ajouter_dependance(22);p20.ajouter_dependance(90);
//p22.ajouter_dependance(1);p22.ajouter_dependance(7);
p23.ajouter_dependance(1);p23.ajouter_dependance(19);p23.ajouter_dependance(25);p23.ajouter_dependance(26);p23.ajouter_dependance(27);p23.ajouter_dependance(28);
p28.ajouter_dependance(1);p28.ajouter_dependance(31);p28.ajouter_dependance(32);p28.ajouter_dependance(33);p28.ajouter_dependance(34);
p29.ajouter_dependance(1);p29.ajouter_dependance(31);p29.ajouter_dependance(32);p29.ajouter_dependance(33);p29.ajouter_dependance(34);
p34.ajouter_dependance(1);p34.ajouter_dependance(19);p34.ajouter_dependance(36);p34.ajouter_dependance(37);
p37.ajouter_dependance(1);p37.ajouter_dependance(19);p37.ajouter_dependance(40);p37.ajouter_dependance(41);p37.ajouter_dependance(42);p37.ajouter_dependance(43);p37.ajouter_dependance(44);p37.ajouter_dependance(45);p37.ajouter_dependance(61);
p38.ajouter_dependance(1);p38.ajouter_dependance(31);
p45.ajouter_dependance(1);p45.ajouter_dependance(31);p45.ajouter_dependance(49);
p46.ajouter_dependance(1);p46.ajouter_dependance(31);p46.ajouter_dependance(49);
p49.ajouter_dependance(1);p49.ajouter_dependance(19);p49.ajouter_dependance(51);p49.ajouter_dependance(5);
p51.ajouter_dependance(1);p51.ajouter_dependance(10);p51.ajouter_dependance(53);
p53.ajouter_dependance(1);p53.ajouter_dependance(5);p53.ajouter_dependance(7);p53.ajouter_dependance(19);p53.ajouter_dependance(63);p53.ajouter_dependance(64);p53.ajouter_dependance(65);p53.ajouter_dependance(66);p53.ajouter_dependance(84);p53.ajouter_dependance(85);p53.ajouter_dependance(88);p53.ajouter_dependance(89);p53.ajouter_dependance(99);p53.ajouter_dependance(100);p53.ajouter_dependance(101);
p55.ajouter_dependance(1);
p56.ajouter_dependance(1);p56.ajouter_dependance(19);
p57.ajouter_dependance(1);p57.ajouter_dependance(19);p57.ajouter_dependance(59);
p59.ajouter_dependance(1);p59.ajouter_dependance(19);p59.ajouter_dependance(59);
p61.ajouter_dependance(1);p61.ajouter_dependance(5);p61.ajouter_dependance(19);
p66.ajouter_dependance(1);p66.ajouter_dependance(31);p66.ajouter_dependance(68);p66.ajouter_dependance(69);p66.ajouter_dependance(70);p66.ajouter_dependance(71);
p71.ajouter_dependance(1);p71.ajouter_dependance(73);p71.ajouter_dependance(68);p71.ajouter_dependance(69);p71.ajouter_dependance(19);p71.ajouter_dependance(5);
p73.ajouter_dependance(1);p73.ajouter_dependance(5);p73.ajouter_dependance(40);p73.ajouter_dependance(41);p73.ajouter_dependance(43);p73.ajouter_dependance(44);
p75.ajouter_dependance(1);p75.ajouter_dependance(5);p75.ajouter_dependance(19);p75.ajouter_dependance(77);p75.ajouter_dependance(78);p75.ajouter_dependance(79);p75.ajouter_dependance(80);p75.ajouter_dependance(81);p75.ajouter_dependance(91);p75.ajouter_dependance(92);p75.ajouter_dependance(93);p75.ajouter_dependance(102);
p81.ajouter_dependance(1);p81.ajouter_dependance(5);p81.ajouter_dependance(31);p81.ajouter_dependance(83);
p85.ajouter_dependance(1);p85.ajouter_dependance(31);p85.ajouter_dependance(87);
p93.ajouter_dependance(1);p93.ajouter_dependance(31);p93.ajouter_dependance(49);
p94.ajouter_dependance(1);p94.ajouter_dependance(31);p94.ajouter_dependance(49);
p95.ajouter_dependance(1);p95.ajouter_dependance(31);p95.ajouter_dependance(49);
p96.ajouter_dependance(1);p96.ajouter_dependance(31);p96.ajouter_dependance(49);
p102.ajouter_dependance(1);p102.ajouter_dependance(31);p102.ajouter_dependance(87);p102.ajouter_dependance(104);
p104.ajouter_dependance(1);p104.ajouter_dependance(19);
p105.ajouter_dependance(1);p105.ajouter_dependance(107);
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(p17);   
lst.push_back(p18);   
lst.push_back(p19);   
lst.push_back(p20);   
lst.push_back(p21);   
lst.push_back(p22);   
lst.push_back(p23);   
lst.push_back(p24);   
lst.push_back(p25);   
lst.push_back(p26);   
lst.push_back(p27);   
lst.push_back(p28);   
lst.push_back(p29);   
lst.push_back(p30);   
lst.push_back(p31);   
lst.push_back(p32);   
lst.push_back(p33);   
lst.push_back(p34);   
lst.push_back(p35);   
lst.push_back(p36);   
lst.push_back(p37);   
lst.push_back(p38);   
lst.push_back(p39);   
lst.push_back(p40);   
lst.push_back(p41);   
lst.push_back(p42);   
lst.push_back(p43);   
lst.push_back(p44);   
lst.push_back(p45);   
lst.push_back(p46);   
lst.push_back(p47);   
lst.push_back(p48);   
lst.push_back(p49);   
lst.push_back(p50);   
lst.push_back(p51);   
lst.push_back(p52);   
lst.push_back(p53);   
lst.push_back(p54);   
lst.push_back(p55);   
lst.push_back(p56);   
lst.push_back(p57);   
lst.push_back(p58);   
lst.push_back(p59);   
lst.push_back(p60);   
lst.push_back(p61);   
lst.push_back(p62);
lst.push_back(p63);
lst.push_back(p64);
lst.push_back(p65);
lst.push_back(p66);
lst.push_back(p67);
lst.push_back(p68);
lst.push_back(p69);
lst.push_back(p70);
lst.push_back(p71);
lst.push_back(p72);
lst.push_back(p73);
lst.push_back(p74);
lst.push_back(p75);
lst.push_back(p76);
lst.push_back(p77);
lst.push_back(p78);
lst.push_back(p79);
lst.push_back(p80);
lst.push_back(p81);
lst.push_back(p82);
lst.push_back(p83);
lst.push_back(p84);
lst.push_back(p85);
lst.push_back(p86);
lst.push_back(p87);
lst.push_back(p88);
lst.push_back(p89);
lst.push_back(p90);
lst.push_back(p91);
lst.push_back(p92);
lst.push_back(p93);
lst.push_back(p94);
lst.push_back(p95);
lst.push_back(p96);
lst.push_back(p97);
lst.push_back(p98);
lst.push_back(p99);
lst.push_back(p100);
lst.push_back(p101);
lst.push_back(p102);
lst.push_back(p103);
lst.push_back(p104);
lst.push_back(p105);
lst.push_back(p106);
MAGIC_APPLICATION app((char*)"Operation dans l'environnement MAGiC",argc,argv,lst,true,false);  
char nomfichier[3000];
char* p;
if (app.get_action()!=48)
  {
  if (app.get_erreur()==true) return 0;     
  if ((app.get_action()!=82) && (app.get_action()!=54)&& (app.get_action()!=6))
    {
    app.recupere_parametre_string(1,nomfichier,(char*)"Fichier entrant manquant");    
    p=strrchr(nomfichier,'.');
    }
  }  
char chaine[1000];for (int i=0;i<1000;i++) chaine[i]=0;
if (app.get_action()==2)
        {
        bool desactivedeforme;
        app.recupere_parametre_bool_avec_defaut(4,desactivedeforme,false);
        bool corres;
        app.recupere_parametre_bool_avec_defaut(55,corres,false);
        if (app.get_erreur()==true) return 0;     
        char mess[2000];
        sprintf(mess,"Conversion MAGiC vers GMSH");
        app.affiche(mess);
        MG_FILE gest(nomfichier);
        int nb=gest.get_nb_mg_maillage();
        for (int i=0;i<nb;i++)
        {
            MG_MAILLAGE* mai=gest.get_mg_maillage(i);
            unsigned long id=mai->get_id();
            strncpy(chaine,nomfichier,p-nomfichier);
            char chaine2[3000];
            sprintf(chaine2,"%s%lu",chaine,id);
            std::string namefic=chaine2;
            MG_EXPORT exp;
            sprintf(mess,"      Maillage geometrique : creation de %s.msh",namefic.c_str());
            app.affiche(mess);
            if (corres==false) exp.gmsh(mai,namefic); else exp.gmsh_avec_correspondance(mai,namefic);
            }
        nb=gest.get_nb_fem_maillage();
        for (int i=0;i<nb;i++)
        {
            FEM_MAILLAGE* mai=gest.get_fem_maillage(i);
            if (desactivedeforme) mai->desactive_deforme();
            unsigned long id=mai->get_id();
            strncpy(chaine,nomfichier,p-nomfichier);
            char chaine2[3000];
            sprintf(chaine2,"%s%lu",chaine,id);
            std::string namefic=chaine2;
            MG_EXPORT exp;
            sprintf(mess,"      Maillage EF : Creation de %s.msh",namefic.c_str());
            app.affiche(mess);
            if (corres==false) exp.gmsh(mai,namefic); else exp.gmsh_avec_correspondance(mai,namefic);
        }
    }
if (app.get_action()==4)
      {
      char fichierout[2000];
      app.recupere_parametre_string(5,fichierout,(char*)"Fichier sortant manquant");    
      if (app.get_erreur()==true) return 0;
      char mess[2000];
      sprintf(mess,"Copie de fichier MAGiC");
      app.affiche(mess);
      MG_FILE gest(nomfichier);
      gest.enregistrer(fichierout);
      }        
if (app.get_action()==6)
    {
        char mess[500];
        char coderesu[10];
        char nomparam[500];
        bool creeparam;
        app.recupere_parametre_string_avec_defaut(7,coderesu,(char*)"11111111");     
        app.recupere_parametre_bool_avec_defaut(85,creeparam,false);
        if (!creeparam) app.recupere_parametre_string(84,nomparam,(char*)"Fichier parametre manquant");
        if (!creeparam) app.recupere_parametre_string(1,nomfichier,(char*)"Fichier entrant manquant");
        app.recupere_parametre_string(84,nomparam,(char*)"Fichier parametre manquant");
        if (app.get_erreur()==true) return 0;
        if (creeparam)
        {
        char mess[500];
        sprintf(mess,"   Calcul ASTER\n\nCréation d'un fichier parametre\n");
        app.affiche(mess);
        MG_EXPORT exp;
        exp.ecrire_params_aster(nomparam);
        }
      else 
        {
        sprintf(mess,"Conversion MAGiC vers code ASTER");
        app.affiche(mess);
        MG_FILE gest(nomfichier);
        int nb=gest.get_nb_fem_maillage();
        for (int i=0;i<nb;i++)
        {
            FEM_MAILLAGE* mai=gest.get_fem_maillage(i);
            MG_GEOMETRIE* geo=mai->get_mg_geometrie();
            unsigned long id=mai->get_id();
            p=strrchr(nomfichier,'.');
            strncpy(chaine,nomfichier,p-nomfichier);
            char chaine2[3000];
            sprintf(chaine2,"%s%lu",chaine,id);
            std::string namefic=chaine2;
            MG_EXPORT exp;
            exp.lire_params_aster(nomparam);
            exp.aster(geo,mai,namefic,MAGIC::CALCUL_ASTER::ELASTIQUE,coderesu);
        }
          
        }
    }
 if (app.get_action()==8)
    {
        if (app.get_erreur()==true) return 0;
        char mess[500];
        sprintf(mess,"   Conversion carte de taille ancienne version vers GMSH");
        app.affiche(mess);
        FCT_GENERATEUR_3D<4> carte;
        carte.lire(nomfichier);
        strncpy(chaine,nomfichier,p-nomfichier);
        std::string namefic=chaine;
        FCT_EXPORT exp;
        exp.gmsh(&carte,namefic);
    }     
#ifdef BREP_STEP
if (app.get_action()==9)
    {
        double unite;
        app.recupere_parametre_double_avec_defaut(10,unite,1.);
        if (app.get_erreur()==true) return 0;
        char mess[500];
        sprintf(mess,"   Conversion STEP vers MAGiC");
        app.affiche(mess);
        STEP_IMPORT stepimport;
        MG_GESTIONNAIRE gest;
        stepimport.importer(gest,nomfichier) ;
        strncpy(chaine,nomfichier,p-nomfichier);
        std::string namefic=chaine;
        namefic=namefic  + ".magic";
        int nb=gest.get_nb_mg_geometrie();
        for (int i=0;i<nb;i++)
            gest.get_mg_geometrie(i)->change_valeur_unite(unite);
        gest.enregistrer(namefic.c_str());
    }
#endif     
#ifdef BREP_OCC
if (app.get_action()==11)
    {
        double eps,eps2;
        bool importstl,importtriangulation;
        app.recupere_parametre_double_avec_defaut(13,eps,0.000001);
        app.recupere_parametre_double_avec_defaut(14,eps2,1.);
        app.recupere_parametre_bool_avec_defaut(15,importstl,false);
        app.recupere_parametre_bool_avec_defaut(16,importtriangulation,false);
        if (app.get_erreur()==true) return 0;
        int importtri=0;
        if (importstl) importtri=1;
        if (importtriangulation) importtri=2;
        char mess[500];
        sprintf(mess,"   Conversion STEP via OpenCascade vers MAGiC avec epsilon=%lf",eps);
        app.affiche(mess);
        MG_GESTIONNAIRE gest;
        OCC_IMPORT occimport;
        MG_GEOMETRIE* mggeo=occimport.importer(gest,nomfichier,FICHIERSTEP,eps) ;
        if (importtri>0) occimport.importer(gest,mggeo,eps2,importtri);
        strncpy(chaine,nomfichier,p-nomfichier);
        std::string namefic=chaine;
        namefic=namefic  + ".magic";
        gest.enregistrer(namefic.c_str());
    }
#endif    
#ifdef BREP_SAT
if (app.get_action()==12)
    {
        double unite;
        app.recupere_parametre_double_avec_defaut(10,unite,1.);
        if (app.get_erreur()==true) return 0;
        char mess[500];
        sprintf(mess,"   Conversion SAT vers MAGiC");
        app.affiche(mess);
        MG_GESTIONNAIRE gest;
        ACIS_IMPORT acisimport;
        acisimport.importer(gest,nomfichier) ;
        strncpy(chaine,nomfichier,p-nomfichier);
        std::string namefic=chaine;
        namefic=namefic  + ".magic";
        int nb=gest.get_nb_mg_geometrie();
        for (int i=0;i<nb;i++)
            gest.get_mg_geometrie(i)->change_valeur_unite(unite);
        gest.enregistrer(namefic.c_str());
    }
#endif     
#ifdef BREP_OCC
if (app.get_action()==17)
    {
        double eps,eps2;
        bool importstl,importtriangulation;
        app.recupere_parametre_double_avec_defaut(13,eps,0.000001);
        app.recupere_parametre_double_avec_defaut(14,eps2,1.);
        app.recupere_parametre_bool_avec_defaut(15,importstl,false);
        app.recupere_parametre_bool_avec_defaut(16,importtriangulation,false);
        if (app.get_erreur()==true) return 0;
        int importtri=0;
        if (importstl) importtri=1;
        if (importtriangulation) importtri=2;
        char mess[500];
        sprintf(mess,"   Conversion OpenCascade vers MAGiC");
        app.affiche(mess);
        MG_GESTIONNAIRE gest;
        OCC_IMPORT occimport;
        MG_GEOMETRIE* mggeo=occimport.importer(gest,nomfichier,FICHIEROCC,eps) ;
        if (importtri>0) occimport.importer(gest,mggeo,eps2,importtri);
        strncpy(chaine,nomfichier,p-nomfichier);
        std::string namefic=chaine;
        namefic=namefic  + ".magic";
        gest.enregistrer(namefic.c_str());
    }
#endif
if (app.get_action()==18)
    {
        int degre,nummai;
        bool analyse;
        bool optno;
        char fichierout[2000];
        app.recupere_parametre_string_avec_defaut(5,fichierout,nomfichier);    
        app.recupere_parametre_int_avec_defaut(19,nummai,0);
        app.recupere_parametre_int(20,degre,(char*)"Degré du maillage manquant");
        app.recupere_parametre_bool_avec_defaut(75,analyse,false);
        app.recupere_parametre_bool_avec_defaut(98,optno,false);
        if (app.get_erreur()==true) return 0;
        char mess[500];
        sprintf(mess,"   Creation d'un maillage FEM");
        app.affiche(mess);
        MG_FILE gest(nomfichier);
        MG_MAILLAGE* mai;
        if (nummai==0)    mai=gest.get_mg_maillage(nummai); else mai=gest.get_mg_maillageid(nummai);
        FEM_MAILLAGE* fem=new FEM_MAILLAGE(mai->get_mg_geometrie(),mai,degre);
        gest.ajouter_fem_maillage(fem);
        if (!optno) fem->construire(0); else fem->construire(1);
        gest.enregistrer(fichierout);
        if (analyse==true) 
            {
            MAILLEUR_ANALYSE m3d(fem);
            m3d.active_affichage(app.affiche);
            m3d.analyse();
            }
            /*if (degre==2)
           { 
           //char message [5000];
           FEM_MAILLAGE_QUADRATIQUE_OUTILS ou ;
           double moy=ou.get_distortion_moy(fem);
           double min=ou.get_distortion_min (fem);
           int nb=ou.get_nombre_elements_invalides(fem);
           sprintf(mess," maillage quadratique :");
           app.affiche(mess);
            sprintf(mess," ");
           app.affiche(mess);
           sprintf(mess," distorsion moy               %f  ",moy);
           app.affiche(mess);
           sprintf(mess," distorsion min               %f  ",min);
           app.affiche(mess);
           sprintf(mess," nombre elements invalides    %d  ",nb);
           app.affiche(mess);
           }*/ 
          
          
    }      
if (app.get_action()==21)
    {
        int nummai;
        char fichieraster[500];
        int base;
        app.recupere_parametre_int_avec_defaut(19,nummai,0);
        app.recupere_parametre_string(22,fichieraster,(char*)"Fichier aster manquant");
        app.recupere_parametre_int_avec_defaut(90,base,10);
        if (app.get_erreur()==true) return 0;
        char mess[500];
        sprintf(mess,"   Importation des résultats ASTER vers MAGIC");
        app.affiche(mess);
        MG_FILE gest(nomfichier);
        FEM_MAILLAGE* mai;
        if (nummai==0)    mai=gest.get_fem_maillage(nummai); else mai=gest.get_fem_maillageid(nummai);
        MG_IMPORT imp;
        std::string res=imp.aster(base,mai,fichieraster,nomfichier);
        app.affiche((char*)("\033[1;31m"+res+"\033[1;32m").c_str());
        gest.enregistrer(nomfichier);
    }     
 
if (app.get_action()==23)
    {

    }
if (app.get_action()==24)
    {
       int nummai;
       char fichieroptis[500];
       char fichierflottant[500];
       double niveau,seuil;
       app.recupere_parametre_int_avec_defaut(19,nummai,0);
       app.recupere_parametre_string(25,fichieroptis,(char*)"Fichier OPTIS manquant");
       app.recupere_parametre_string_avec_defaut(26,fichieroptis,(char*)"");
       app.recupere_parametre_double(27,seuil,(char*)"Seuil manquant");
       app.recupere_parametre_double(28,niveau,(char*)"Niveau manquant");
       if (app.get_erreur()==true) return 0;
        char mess[500];
        sprintf(mess,"   Importation des résultats d'optimisation vers MAGIC");
        app.affiche(mess);
        MG_FILE gest(nomfichier);
        FEM_MAILLAGE* mai;
        if (nummai==0)    mai=gest.get_fem_maillage(nummai); else mai=gest.get_fem_maillageid(nummai);
        MG_IMPORT imp;
        int flottant=0;
        if (strlen(fichierflottant)>0) flottant=1;
        imp.optis(mai,fichieroptis,fichierflottant,seuil,niveau,flottant);
        gest.enregistrer(nomfichier);
    } 
 
if (app.get_action()==29)
    {
        int numgeo,numgt,numentite;
        bool expand;
        app.recupere_parametre_int_avec_defaut(31,numgeo,0);
        app.recupere_parametre_int_avec_defaut(32,numgt,0);
        app.recupere_parametre_int(33,numentite,(char*)"Numéro d'entité manquant");
        app.recupere_parametre_bool_avec_defaut(34,expand,false);
        if (app.get_erreur()==true) return 0;
        char mess[500];
        sprintf(mess,"   Selection de groupe topologique");
        app.affiche(mess);
        MG_FILE gest(nomfichier);
        MG_GEOMETRIE* mggeo;
        if (numgeo==0) mggeo=gest.get_mg_geometrie(numgeo); else mggeo=gest.get_mg_geometrieid(numgeo);
        MG_GROUPE_TOPOLOGIQUE* mggt;
        if (numgt==0)
        {
            mggt=new MG_GROUPE_TOPOLOGIQUE;
            mggeo->ajouter_mg_groupe_topologique(mggt);
        }
        else mggt=mggeo->get_mg_groupe_topologiqueid(numgt);
        MG_ELEMENT_TOPOLOGIQUE* ele;
        ele=mggeo->get_mg_sommetid(numentite);
        if (ele==NULL) ele=mggeo->get_mg_areteid(numentite);
        if (ele==NULL) ele=mggeo->get_mg_faceid(numentite);
        if (ele==NULL) ele=mggeo->get_mg_volumeid(numentite);
        if (ele!=NULL)
        {
            mggt->ajouter(ele);
            if (expand)
            {
                TPL_MAP_ENTITE<MG_ELEMENT_TOPOLOGIQUE*> lst;
                ele->get_topologie_sousjacente(&lst);
                TPL_MAP_ENTITE<MG_ELEMENT_TOPOLOGIQUE*>::ITERATEUR it;
                for (MG_ELEMENT_TOPOLOGIQUE *ele2=lst.get_premier(it);ele2!=NULL;ele2=lst.get_suivant(it))
                    mggt->ajouter(ele2);
            }
        }
        gest.enregistrer(nomfichier);
    }
if (app.get_action()==30)
    {
        int numgeo,numgt,numentite;
        bool expand;
        app.recupere_parametre_int_avec_defaut(31,numgeo,0);
        app.recupere_parametre_int_avec_defaut(32,numgt,0);
        app.recupere_parametre_int(33,numentite,(char*)"Numéro d'entité manquant");
        app.recupere_parametre_bool_avec_defaut(34,expand,false);
        if (app.get_erreur()==true) return 0;
        char mess[500];
        sprintf(mess,"   Deselection de groupe topologique");
        app.affiche(mess);
        MG_FILE gest(nomfichier);
        MG_GEOMETRIE* mggeo=gest.get_mg_geometrieid(numgeo);
        MG_GROUPE_TOPOLOGIQUE* mggt=mggeo->get_mg_groupe_topologiqueid(numgt);
        MG_ELEMENT_TOPOLOGIQUE* ele;
        ele=mggeo->get_mg_sommetid(numentite);
        if (ele==NULL) ele=mggeo->get_mg_areteid(numentite);
        if (ele==NULL) ele=mggeo->get_mg_faceid(numentite);
        if (ele==NULL) ele=mggeo->get_mg_volumeid(numentite);
        if (ele!=NULL)
        {
            mggt->supprimer(ele);
            if (expand)
            {
                TPL_MAP_ENTITE<MG_ELEMENT_TOPOLOGIQUE*> lst;
                ele->get_topologie_sousjacente(&lst);
                TPL_MAP_ENTITE<MG_ELEMENT_TOPOLOGIQUE*>::ITERATEUR it;
                for (MG_ELEMENT_TOPOLOGIQUE *ele2=lst.get_premier(it);ele2!=NULL;ele2=lst.get_suivant(it))
                    mggt->supprimer(ele2);
            }
        }
        gest.enregistrer(nomfichier);
    } 
 
 if (app.get_action()==35)
    {
        int nummai;
        double angle;
        char fichierlog[500];
        strcpy(fichierlog,"");
        app.recupere_parametre_int_avec_defaut(19,nummai,-1);
        app.recupere_parametre_double_avec_defaut(36,angle,0.03);
        app.recupere_parametre_string_avec_defaut(37,fichierlog,(char*)fichierlog);
        if (app.get_erreur()==true) return 0;
        char mess[500];
        sprintf(mess,"   Contenu d'un fichier MAGIC");
        app.affiche(mess);
        MG_FILE gest(nomfichier);
        if (nummai!=-1)
                {
                MG_MAILLAGE* mai=gest.get_mg_maillageid(nummai);
                if (mai!=NULL)
                  {
                  MAILLEUR_ANALYSE m3d(mai);
                  m3d.active_affichage(app.affiche);
                  m3d.change_eps_angle_retourne(angle);
                  if (strlen(fichierlog)==0)
                    m3d.analyse();
                  else 
                    m3d.analyse(fichierlog);
                  }
                else
                {
                  FEM_MAILLAGE* mai=gest.get_fem_maillageid(nummai);
                if (mai!=NULL)
                  {
                  MAILLEUR_ANALYSE m3d(mai);
                  m3d.active_affichage(app.affiche);
                  m3d.change_eps_angle_retourne(angle);
                  if (strlen(fichierlog)==0)
                    m3d.analyse();
                  else 
                    m3d.analyse(fichierlog);
                  }
                
                  }
                  
                }
                
        else 
        {
        char message [5000];
        int nbarbre=gest.get_nb_mg_arbre();
        sprintf(message,"           \033[1;32m%d\033[1;33m arbres :",nbarbre);
        for (int i=0;i<nbarbre;i++)
        {
            MG_ARBRE* arb=gest.get_mg_arbre(i);
            sprintf(message,"%s \033[1;31m%lu\033[1;33m ",message,arb->get_id());
        }
        app.affiche(message);
        for (int i=0;i<nbarbre;i++)
        {
            MG_ARBRE* arb=gest.get_mg_arbre(i);
            sprintf(message,"              arbre \033[1;31m%lu\033[1;33m",arb->get_id());
            app.affiche(message);
            int nbass=arb->get_nb_mg_assemblage();
            sprintf(message,"                 \033[1;32m%d\033[1;33m assemblages",nbass);
            app.affiche(message);
            int nbprim=arb->get_nb_mg_primitive();
            sprintf(message,"                 \033[1;32m%d\033[1;33m primitives",nbprim);
            //for (int i=0;i<nbprim;i++) sprintf(message,"%s \033[1;31m%lu\033[1;33m",message,arb->get_mg_primitive(i)->get_id());
            app.affiche(message);
            int nbope=arb->get_nb_mg_operateur_boolean();
            sprintf(message,"                 \033[1;32m%d\033[1;33m operateurs booleans",nbope);
            app.affiche(message);
        }
        int nbgeo=gest.get_nb_mg_geometrie();
        sprintf(message,"           \033[1;32m%d\033[1;33m geometries :",nbgeo);
        for (int i=0;i<nbgeo;i++)
        {
            MG_GEOMETRIE* geo=gest.get_mg_geometrie(i);
            sprintf(message,"%s \033[1;31m%lu\033[1;33m ",message,geo->get_id());
        }
        app.affiche(message);
        for (int i=0;i<nbgeo;i++)
        {
            MG_GEOMETRIE* geo=gest.get_mg_geometrie(i);
            sprintf(message,"              geometrie \033[1;31m%lu\033[1;33m",geo->get_id());
            app.affiche(message);
            int nbmg=geo->get_nb_mg_groupe_topologique();
            sprintf(message,"                 \033[1;32m%d\033[1;33m groupe topologique",nbmg);
            for (int j=0;j<nbmg;j++)
            {
                MG_GROUPE_TOPOLOGIQUE* mggp=geo->get_mg_groupe_topologique(j);
                sprintf(message,"%s \033[1;31m%lu\033[1;33m ",message,mggp->get_id());
            }
            app.affiche(message);
            int nbvolume=geo->get_nb_mg_volume();
            if (nbvolume==0) sprintf(message,"                 \033[1;32m%d\033[1;33m volumes",geo->get_nb_mg_volume());
            else sprintf(message,"                 \033[1;32m%d\033[1;33m volumes : ",geo->get_nb_mg_volume());
            for (int i=0;i<nbvolume;i++) sprintf(message,"%s \033[1;31m%lu\033[1;33m",message,geo->get_mg_volume(i)->get_id()); 
            app.affiche(message);
            int nbcoque=geo->get_nb_mg_coque();
            if (nbcoque==0) sprintf(message,"                 \033[1;32m%d\033[1;33m coques",geo->get_nb_mg_coque());
            else sprintf(message,"                 \033[1;32m%d\033[1;33m coques : ",geo->get_nb_mg_coque());
            for (int i=0;i<nbcoque;i++) sprintf(message,"%s \033[1;31m%lu\033[1;33m",message,geo->get_mg_coque(i)->get_id()); 
            app.affiche(message);
            int nbpoutre=geo->get_nb_mg_poutre();
            if (nbpoutre==0) sprintf(message,"                 \033[1;32m%d\033[1;33m poutres",geo->get_nb_mg_poutre());
            else sprintf(message,"                 \033[1;32m%d\033[1;33m poutres : ",geo->get_nb_mg_poutre());
            for (int i=0;i<nbpoutre;i++) sprintf(message,"%s \033[1;31m%lu\033[1;33m",message,geo->get_mg_poutre(i)->get_id()); 
            app.affiche(message);
            sprintf(message,"                 \033[1;32m%d\033[1;33m coquilles",geo->get_nb_mg_coquille());
            app.affiche(message);
            sprintf(message,"                 \033[1;32m%d\033[1;33m faces",geo->get_nb_mg_face());
            app.affiche(message);
            sprintf(message,"                 \033[1;32m%d\033[1;33m aretes",geo->get_nb_mg_arete());
            app.affiche(message);
            sprintf(message,"                 \033[1;32m%d\033[1;33m sommets",geo->get_nb_mg_sommet());
            app.affiche(message);
            int nbgef=geo->get_nb_mg_geom_fonction();
            if (nbgef==0) sprintf(message,"                 \033[1;32m%d\033[1;33m fonctions geometriques",nbgef);
            else sprintf(message,"                 \033[1;32m%d\033[1;33m fonctions geometriques : ",nbgef);
            for (int i=0;i<nbgef;i++) sprintf(message,"%s \033[1;31m%lu\033[1;33m",message,geo->get_mg_geom_fonction(i)->get_id()); 
            app.affiche(message);

        }
        int nbmai=gest.get_nb_mg_maillage();
        sprintf(message,"          \033[1;32m%d\033[1;33m maillage geometrique :",nbmai);
        app.affiche(message);
        for (int i=0;i<nbmai;i++)
        {
            MG_MAILLAGE* mai=gest.get_mg_maillage(i);
            sprintf(message,"             maillage geometrique \033[1;31m%lu\033[1;33m ",mai->get_id());
            app.affiche(message);
            sprintf(message,"                 \033[1;32m%d\033[1;33m noeuds",mai->get_nb_mg_noeud());
            app.affiche(message);
            sprintf(message,"                 \033[1;32m%d\033[1;33m segments",mai->get_nb_mg_segment());
            app.affiche(message);
            sprintf(message,"                 \033[1;32m%d\033[1;33m triangles",mai->get_nb_mg_triangle());
            app.affiche(message);
            sprintf(message,"                 \033[1;32m%d\033[1;33m quadrangles",mai->get_nb_mg_quadrangle());
            app.affiche(message);
            sprintf(message,"                 \033[1;32m%d\033[1;33m tetras",mai->get_nb_mg_tetra());
            app.affiche(message);
            sprintf(message,"                 \033[1;32m%d\033[1;33m hexas",mai->get_nb_mg_hexa());
            app.affiche(message);
            sprintf(message,"                 \033[1;32m%d\033[1;33m pentas",mai->get_nb_mg_penta());
            app.affiche(message);
        }
        int nbmaif=gest.get_nb_fem_maillage();
        sprintf(message,"          \033[1;32m%d\033[1;33m maillage FEM :",nbmaif);
        app.affiche(message);
        for (int i=0;i<nbmaif;i++)
        {
            FEM_MAILLAGE* maif=gest.get_fem_maillage(i);
            sprintf(message,"             maillage FEM \033[1;31m%lu\033[1;33m ",maif->get_id());
            app.affiche(message);
            sprintf(message,"                 \033[1;32m%d\033[1;33m noeuds",maif->get_nb_fem_noeud());
            app.affiche(message);
            sprintf(message,"                 \033[1;32m%d\033[1;33m elements 1D",maif->get_nb_fem_element1());
            app.affiche(message);
            sprintf(message,"                 \033[1;32m%d\033[1;33m elements 2D",maif->get_nb_fem_element2());
            app.affiche(message);
            sprintf(message,"                 \033[1;32m%d\033[1;33m elements 3D",maif->get_nb_fem_element3());
            app.affiche(message);
            double rx,ry,rz;
            maif->calcul_somme_reaction(rx,ry,rz);
            sprintf(message,"                 Reactions : \n                  Rx=\033[1;31m%le\033[1;33m\n                  Ry=\033[1;31m%le\033[1;33m\n                  Rz=\033[1;31m%le\033[1;33m",rx,ry,rz);
            app.affiche(message);
            
        }
        int nbsol=gest.get_nb_mg_solution();
        sprintf(message,"          \033[1;32m%d\033[1;33m solution maillage geometrique :",nbsol);
        app.affiche(message);
        message[0]=0;
        for (int i=0;i<nbsol;i++)
        {
            MG_SOLUTION* sol=gest.get_mg_solution(i);
            char legende[500];
            legende[0]=0;
            message[0]=0;
            for (int j=0;j<sol->get_nb_champ();j++)
              sprintf(legende,"%s%s_%s ",legende,(char*)sol->get_nom().c_str(),(char*)sol->get_legende(j).c_str());
            sprintf(message,"                \033[1;31m%lu(%s)\033[1;33m ",sol->get_id(),legende);
            app.affiche(message);
        }
        int nbsolf=gest.get_nb_fem_solution();
        sprintf(message,"          \033[1;32m%d\033[1;33m solution maillage FEM :",nbsolf);
        app.affiche(message);
        message[0]=0;
        for (int i=0;i<nbsolf;i++)
        {
            FEM_SOLUTION* sol=gest.get_fem_solution(i);
            char legende[500];
            legende[0]=0;
            message[0]=0;
            for (int j=0;j<sol->get_nb_champ();j++)
              sprintf(legende,"%s%s_%s ",legende,(char*)sol->get_nom().c_str(),(char*)sol->get_legende(j).c_str());
            sprintf(message,"                 \033[1;31m%lu(%s)\033[1;33m ",sol->get_id(),legende);
            app.affiche(message);
        }
        }
    }
if (app.get_action()==38)
    {
        int nummai;
        int numsol1,numsol2,numsol3;
        int numchamp1,numchamp2,numchamp3;
        double coef;
        app.recupere_parametre_int_avec_defaut(19,nummai,0);
        app.recupere_parametre_int(40,numsol1,(char*)"Numéro solution 1 manquant");
        app.recupere_parametre_int(41,numsol2,(char*)"Numéro solution 2 manquant");
        app.recupere_parametre_int(42,numsol3,(char*)"Numéro solution 3 manquant");
        app.recupere_parametre_int(43,numchamp1,(char*)"Numéro champs 1 manquant");
        app.recupere_parametre_int(44,numchamp2,(char*)"Numéro champs 2 manquant");
        app.recupere_parametre_int(45,numchamp3,(char*)"Numéro champs 3 manquant");
        app.recupere_parametre_double_avec_defaut(61,coef,0.);
        if (app.get_erreur()==true) return 0;
        char mess[500];
        sprintf(mess,"   Calcul de deforme");
        app.affiche(mess);
        MG_FILE gest(nomfichier);
        FEM_MAILLAGE* mai;
        if (nummai==0)    mai=gest.get_fem_maillage(nummai); else mai=gest.get_fem_maillageid(nummai);
        FEM_SOLUTION* sol1=gest.get_fem_solutionid(numsol1);
        FEM_SOLUTION* sol2=gest.get_fem_solutionid(numsol2);
        FEM_SOLUTION* sol3=gest.get_fem_solutionid(numsol3);
        mai->calcul_deforme(sol1,numchamp1,sol2,numchamp2,sol3,numchamp3);
        if (coef>1e-16)
          {
          sprintf(mess,"   Création du maillage deforme");
          app.affiche(mess);
          MG_MAILLAGE* mgmai=new MG_MAILLAGE(mai,coef);
          }
        sprintf(mess,"   Enregistrement");
        app.affiche(mess);
        gest.enregistrer(nomfichier);
    }
if (app.get_action()==39)
    {
        int numgeo;
        app.recupere_parametre_int_avec_defaut (31,numgeo,0);
        if (app.get_erreur()==true) return 0;
        char mess[500];
         sprintf(mess,"   Fusion de géométrie");
        app.affiche(mess);
        MG_FILE gest(nomfichier);
        MG_GEOMETRIE* mggeo;
        if (numgeo==0) mggeo=gest.get_mg_geometrie(numgeo); else mggeo=gest.get_mg_geometrieid(numgeo);
        VCT_MULTI_MODELE fusion(mggeo);
        fusion.recherche_identite();
        gest.enregistrer(nomfichier);
    }  
 if ((app.get_action()==46) || (app.get_action()==47)) 
    {
        int numgeo;
        char ccf[500];
        app.recupere_parametre_int_avec_defaut (31,numgeo,0);
        app.recupere_parametre_string(49,ccf,(char*)"Valeur de la condition aux limites manquant");
        if (app.get_erreur()==true) return 0;
        char mess[500];
        sprintf(mess,"   Selection de condition aux limites");
        unsigned long identity;
        char typeccf[3];
        char effortsuiv[3]={'N','S',0};
        typeccf[2]=0;
        double valccf=0.5;      // Valeur par défaut
        sscanf(ccf,"%lu:%c%c:%le:%c%c\n",&identity,&typeccf[0],&typeccf[1],&valccf,&effortsuiv[0],&effortsuiv[1]);
        MC_GESTIONNAIRE mcgest;
        int res=mcgest.existe(typeccf);
        if (res==0)
          {
          sprintf(mess,"   ERREUR!!!! Code conditions aux limites inconnu");
          app.affiche(mess);
          return 0;
          }
        std::string effortsuivstr=effortsuiv;
        if (!((effortsuivstr=="NS")||(effortsuivstr=="SU")))
          {
          sprintf(mess,"   ERREUR!!!! type de suivi des efforts inconnus");
          app.affiche(mess);
          return 0;
          }
        sprintf(mess,"   Ajout d'une condition aux limites : %s",mcgest.get_description(typeccf).c_str());
        app.affiche(mess);
        MG_FILE gest(nomfichier);
        MG_GEOMETRIE* geo;
        if (numgeo==0) geo=gest.get_mg_geometrie(numgeo); else geo=gest.get_mg_geometrieid(numgeo);
        TPL_MAP_ENTITE<MG_ELEMENT_TOPOLOGIQUE*> liste;
        MG_VOLUME* ele;
        ele=geo->get_mg_volumeid(identity);
        if (ele!=NULL) liste.ajouter(ele);
        else
          {
          MG_COQUE* ele;
          ele=geo->get_mg_coqueid(identity);
          if (ele!=NULL) liste.ajouter(ele);
          else 
            {
            MG_POUTRE* ele;
            ele=geo->get_mg_poutreid(identity);
            if (ele!=NULL) liste.ajouter(ele);
            else
              {
              MG_FACE* ele;
              ele=geo->get_mg_faceid(identity);
              if (ele!=NULL) liste.ajouter(ele);
              else 
                  {
                  MG_ARETE* ele;
                  ele=geo->get_mg_areteid(identity);
                  if (ele!=NULL) liste.ajouter(ele);
                  else
                    {
                    MG_SOMMET* ele;
                    ele=geo->get_mg_sommetid(identity);
                    if (ele!=NULL) liste.ajouter(ele);
                    }
                  }
              }
            }
          }
        if (app.get_action()==47)
          {
          int nb=liste.get_nb();
          for (int i=0;i<nb;i++) liste.get(i)->get_topologie_sousjacente(&liste);
          }
        int nb=liste.get_nb();
        for (int i=0;i<nb;i++) liste.get(i)->ajouter_ccf(typeccf,valccf,effortsuiv);
        gest.enregistrer(nomfichier);
    }    
 if ((app.get_action()==94) || (app.get_action()==95)) 
    {
        int numgeo;
        char ccf[500];
        app.recupere_parametre_int_avec_defaut (31,numgeo,0);
        app.recupere_parametre_string(49,ccf,(char*)"Valeur de la condition aux limites manquant");
        if (app.get_erreur()==true) return 0;
        char mess[500];
        sprintf(mess,"   Selection de condition aux limites");
        unsigned long identity;
        char typeccf[3];
        char effortsuiv[3]={'N','S',0};
        typeccf[2]=0;
        unsigned long valccf=0; // Valeur par défaut
        sscanf(ccf,"%lu:%c%c:%lu:%c%c\n",&identity,&typeccf[0],&typeccf[1],&valccf,&effortsuiv[0],&effortsuiv[1]);
        MC_GESTIONNAIRE mcgest;
        int res=mcgest.existe(typeccf);
        if (res==0)
          {
          sprintf(mess,"   ERREUR!!!! Code conditions aux limites inconnu");
          app.affiche(mess);
          return 0;
          }
        std::string effortsuivstr=effortsuiv;
        if (!((effortsuivstr=="NS")||(effortsuivstr=="SU")))
          {
          sprintf(mess,"   ERREUR!!!! type de suivi des efforts inconnus");
          app.affiche(mess);
          return 0;
          }
        sprintf(mess,"   Ajout d'une condition aux limites : %s",mcgest.get_description(typeccf).c_str());
        app.affiche(mess);
        MG_FILE gest(nomfichier);
        MG_GEOMETRIE* geo;
        if (numgeo==0) geo=gest.get_mg_geometrie(numgeo); else geo=gest.get_mg_geometrieid(numgeo);
        TPL_MAP_ENTITE<MG_ELEMENT_TOPOLOGIQUE*> liste;
        MG_VOLUME* ele;
        ele=geo->get_mg_volumeid(identity);
        if (ele!=NULL) liste.ajouter(ele);
        else
          {
          MG_COQUE* ele;
          ele=geo->get_mg_coqueid(identity);
          if (ele!=NULL) liste.ajouter(ele);
          else 
            {
            MG_POUTRE* ele;
            ele=geo->get_mg_poutreid(identity);
            if (ele!=NULL) liste.ajouter(ele);
            else
              {
              MG_FACE* ele;
              ele=geo->get_mg_faceid(identity);
              if (ele!=NULL) liste.ajouter(ele);
              else 
                  {
                  MG_ARETE* ele;
                  ele=geo->get_mg_areteid(identity);
                  if (ele!=NULL) liste.ajouter(ele);
                  else
                    {
                    MG_SOMMET* ele;
                    ele=geo->get_mg_sommetid(identity);
                    if (ele!=NULL) liste.ajouter(ele);
                    }
                  }
              }
            }
          }
        if (app.get_action()==95)
          {
          int nb=liste.get_nb();
          for (int i=0;i<nb;i++) liste.get(i)->get_topologie_sousjacente(&liste);
          }
        int nb=liste.get_nb();
        for (int i=0;i<nb;i++) liste.get(i)->ajouter_ccf(typeccf,valccf,effortsuiv);
        gest.enregistrer(nomfichier);
    }  
if ((app.get_action()==96) || (app.get_action()==97)) 
    {
        int numgeo;
        char ccf[500];
        app.recupere_parametre_int_avec_defaut (31,numgeo,0);
        app.recupere_parametre_string(49,ccf,(char*)"Valeur de la condition aux limites manquant");
        if (app.get_erreur()==true) return 0;
        char mess[500];
        sprintf(mess,"   Selection de condition aux limites");
        unsigned long identity;
        char typeccf[3];
        char effortsuiv[3]={'N','S',0};
        typeccf[2]=0;
        char  valccf[500];
        char  formule[500];
        sscanf(ccf,"%lu:%c%c:%s:%c%c\n",&identity,&typeccf[0],&typeccf[1],formule,&effortsuiv[0],&effortsuiv[1]);
        OT_CHAINE ot;
        std::vector<std::string> variable=ot.split(formule,':');
        strcpy(valccf,variable[0].c_str());
        std::vector<std::string> listvariable=ot.split(variable[1].c_str(),',');
        MC_GESTIONNAIRE mcgest;
        int res=mcgest.existe(typeccf);
        if (res==0)
          {
          sprintf(mess,"   ERREUR!!!! Code conditions aux limites inconnu");
          app.affiche(mess);
          return 0;
          }
        std::string effortsuivstr=effortsuiv;
        if (!((effortsuivstr=="NS")||(effortsuivstr=="SU")))
          {
          sprintf(mess,"   ERREUR!!!! type de suivi des efforts inconnus");
          app.affiche(mess);
          return 0;
          }
        sprintf(mess,"   Ajout d'une condition aux limites : %s",mcgest.get_description(typeccf).c_str());
        app.affiche(mess);
        MG_FILE gest(nomfichier);
        MG_GEOMETRIE* geo;
        if (numgeo==0) geo=gest.get_mg_geometrie(numgeo); else geo=gest.get_mg_geometrieid(numgeo);
        TPL_MAP_ENTITE<MG_ELEMENT_TOPOLOGIQUE*> liste;
        MG_VOLUME* ele;
        ele=geo->get_mg_volumeid(identity);
        if (ele!=NULL) liste.ajouter(ele);
        else
          {
          MG_COQUE* ele;
          ele=geo->get_mg_coqueid(identity);
          if (ele!=NULL) liste.ajouter(ele);
          else 
            {
            MG_POUTRE* ele;
            ele=geo->get_mg_poutreid(identity);
            if (ele!=NULL) liste.ajouter(ele);
            else
              {
              MG_FACE* ele;
              ele=geo->get_mg_faceid(identity);
              if (ele!=NULL) liste.ajouter(ele);
              else 
                  {
                  MG_ARETE* ele;
                  ele=geo->get_mg_areteid(identity);
                  if (ele!=NULL) liste.ajouter(ele);
                  else
                    {
                    MG_SOMMET* ele;
                    ele=geo->get_mg_sommetid(identity);
                    if (ele!=NULL) liste.ajouter(ele);
                    }
                  }
              }
            }
          }
        if (app.get_action()==97)
          {
          int nb=liste.get_nb();
          for (int i=0;i<nb;i++) liste.get(i)->get_topologie_sousjacente(&liste);
          }
        int nb=liste.get_nb();
        for (int i=0;i<nb;i++) liste.get(i)->ajouter_ccf(typeccf,valccf,listvariable,effortsuiv);
        gest.enregistrer(nomfichier);
    }    
 if (app.get_action()==48)
      {
      char mess[500];
      sprintf(mess,"   Liste des conditions aux limites disponibles");
      app.affiche(mess);
      MC_GESTIONNAIRE mcgest;
      int ok=0;
      do
        {
        std::string code;
        std::string description;
        int res=mcgest.get_description(code,description);
        if (res)
          {
          sprintf(mess,"   %s   : %s",code.c_str(),description.c_str());
          app.affiche(mess);
          }
         else ok=1;
        }
      while (ok==0);
      }    
if (app.get_action()==50)
    {
        int numgeo,nummai,numcoque;
        char fichierout[500];
        app.recupere_parametre_int_avec_defaut(31,numgeo,0);
        app.recupere_parametre_int_avec_defaut(19,nummai,0);
        app.recupere_parametre_int_avec_defaut(51,numcoque,0);
        app.recupere_parametre_string_avec_defaut(5,fichierout,nomfichier);
        if (app.get_erreur()==true) return 0;
        char mess[500];
        sprintf(mess,"   Orientation d'une coque");
        app.affiche(mess);
        MG_FILE gest(nomfichier);
        MG_GEOMETRIE* geo;
        if (numgeo==0) geo=gest.get_mg_geometrie(numgeo); else geo=gest.get_mg_geometrieid(numgeo);
        MG_MAILLAGE* mai;
        if (nummai==0) mai=gest.get_mg_maillage(nummai); else mai=gest.get_mg_maillageid(nummai);
        MG_COQUE* coque;
        if (numcoque==0) coque=mai->get_mg_geometrie()->get_mg_coque(numcoque); else coque=mai->get_mg_geometrie()->get_mg_coqueid(numcoque);
        coque->orienter(mai);
        gest.enregistrer(fichierout);     
    }   
if (app.get_action()==52)
    {
        double epsf,unite;
        app.recupere_parametre_double_avec_defaut(10,unite,1.);
        app.recupere_parametre_double_avec_defaut(53,epsf,0.0001);
        if (app.get_erreur()==true) return 0;
        char mess[500];
        sprintf(mess,"   Conversion STL vers MAGiC");
        app.affiche(mess);
        MG_GESTIONNAIRE gest;
        STL_IMPORT stlimport;
        stlimport.change_eps(epsf);
        stlimport.importer(gest,nomfichier) ;
        strncpy(chaine,nomfichier,p-nomfichier);
        std::string namefic=chaine;
        namefic=namefic  + ".magic";
        int nb=gest.get_nb_mg_geometrie();
        for (int i=0;i<nb;i++)
            gest.get_mg_geometrie(i)->change_valeur_unite(unite);
        gest.enregistrer(namefic.c_str());
    }
if (app.get_action()==54)
    {  
      int nummai;
      char coderesu[10];
      bool elastique;
      bool elastique_gd;
      bool plastique;
      bool thermique;
      bool conplane;
      bool defplane;
      bool rigiditeconplane;
      bool rigiditedefplane;
      bool creeparam;
      bool rigidite;
      char nomparam[500];
      char nomfichierout[500];
      app.recupere_parametre_int_avec_defaut(19,nummai,0);
      app.recupere_parametre_string_avec_defaut(7,coderesu,(char*)"11111111");
      app.recupere_parametre_bool_avec_defaut(63,elastique,false);
      app.recupere_parametre_bool_avec_defaut(64,thermique,false);
      app.recupere_parametre_bool_avec_defaut(65,conplane,false);
      app.recupere_parametre_bool_avec_defaut(66,defplane,false);
      app.recupere_parametre_bool_avec_defaut(85,creeparam,false);
      app.recupere_parametre_bool_avec_defaut(88,elastique_gd,false);
      app.recupere_parametre_bool_avec_defaut(89,plastique,false);
      app.recupere_parametre_bool_avec_defaut(99,rigidite,false);
      app.recupere_parametre_bool_avec_defaut(100,rigiditeconplane,false);
      app.recupere_parametre_bool_avec_defaut(101,rigiditedefplane,false);
      if (!creeparam) app.recupere_parametre_string(1,nomfichier,(char*)"Fichier entrant manquant");
      app.recupere_parametre_string(84,nomparam,(char*)"Fichier parametre manquant");
      app.recupere_parametre_string_avec_defaut(5,nomfichierout,nomfichier);
      if (app.get_erreur()==true) return 0;
      if (creeparam)
        {
        char mess[500];
        sprintf(mess,"   Calcul ASTER\n\nCréation d'un fichier parametre\n");
        app.affiche(mess);
        MG_EXPORT exp;
        exp.ecrire_params_aster(nomparam);
        }
      else 
      {
      char mess[500];
      sprintf(mess,"   Calcul ASTER");
      app.affiche(mess);
      char nometude[1000];
      char *p=strchr(nomfichierout,'.');
      strncpy(nometude,nomfichierout,p-nomfichierout);
      nometude[p-nomfichierout]=0;
      MG_FILE gest(nomfichier);
      FEM_MAILLAGE* mai;
      if (nummai==0) mai=gest.get_fem_maillage(nummai); else mai=gest.get_fem_maillageid(nummai);
      MGASTER mgaster;
      mgaster.active_affichage(app.affiche);
      if (elastique) mgaster.calcule(nomparam,mai,nometude,MAGIC::CALCUL_ASTER::ELASTIQUE,coderesu);
      if (rigidite) mgaster.calcule(nomparam,mai,nometude,MAGIC::CALCUL_ASTER::RIGIDITE,coderesu);
      if (conplane) mgaster.calcule(nomparam,mai,nometude,MAGIC::CALCUL_ASTER::CONTRAINTE_PLANE,coderesu);
      if (defplane) mgaster.calcule(nomparam,mai,nometude,MAGIC::CALCUL_ASTER::DEFORMATION_PLANE,coderesu);
      if (rigiditeconplane) mgaster.calcule(nomparam,mai,nometude,MAGIC::CALCUL_ASTER::RIGIDITE_CONTRAINTE_PLANE,coderesu);
      if (rigiditedefplane) mgaster.calcule(nomparam,mai,nometude,MAGIC::CALCUL_ASTER::RIGIDITE_DEFORMATION_PLANE,coderesu);
      if (thermique) mgaster.calcule(nomparam,mai,nometude,MAGIC::CALCUL_ASTER::THERMIQUE,(char*)"");
      if (elastique_gd) mgaster.calcule(nomparam,mai,nometude,MAGIC::CALCUL_ASTER::ELASTIQUE_GROT_GDEF,coderesu);
      if (plastique) mgaster.calcule(nomparam,mai,nometude,MAGIC::CALCUL_ASTER::PLASTIQUE,coderesu);
      
      }
    }
 if (app.get_action()==56)
    {
        if (app.get_erreur()==true) return 0;
        char mess[500];
        sprintf(mess,"   Conversion carte de taille nouvelle version vers GMSH");
        app.affiche(mess);
        FCT_TAILLE_FEM_SOLUTION carte(nomfichier);
        strncpy(chaine,nomfichier,p-nomfichier);
        std::string namefic=chaine;
        MG_EXPORT exp;
        exp.gmsh(&carte,namefic);
        sprintf(mess,"      Carte de taille FEM : Creation de %s.msh",namefic.c_str());
        app.affiche(mess);
    }
if ((app.get_action()==58) || (app.get_action()==60))
    {
    int nummai;
    app.recupere_parametre_int_avec_defaut(19,nummai,0);
    char nomcarte[500];
    if (app.get_action()==58)
    app.recupere_parametre_string(59,nomcarte,(char*)"Fichier carte manquant");
    if (app.get_action()==60)
    app.recupere_parametre_string(59,nomcarte,(char*)"Fichier carte manquant");
    if (app.get_erreur()==true) return 0;
    char mess[500];
    sprintf(mess,"   Comparaison de la taille reelle d'un maillage avec la carte de taille");
    app.affiche(mess);
    sprintf(mess,"   Lecture fichier MAGiC");
    app.affiche(mess);
    MG_FILE gest(nomfichier);
    MG_MAILLAGE* mai;
    if (nummai==0) mai=gest.get_mg_maillage(nummai); else mai=gest.get_mg_maillageid(nummai);
    sprintf(mess,"   Lecture fichier carte de taille");
    app.affiche(mess);
    FCT_TAILLE *carte;
    if (app.get_action()==58)
          carte=new FCT_GENERATEUR_3D<4>(nomcarte);
    if (app.get_action()==60)
          carte=new FCT_TAILLE_FEM_SOLUTION(nomcarte);
    char nometude[1000];
    char *p=strchr(nomfichier,'.');
    strncpy(nometude,nomfichier,p-nomfichier);
    nometude[p-nomfichier]=0;
    char nomsolution[500];
    sprintf(nomsolution,"%s_comp%lu.sol",nometude,mai->get_id());
    sprintf(mess,"   Comparaison");
    app.affiche(mess);
    MAILLEUR_ANALYSE m3d(mai);  
    sprintf(mess,"   Respect global de la taille demandee");
    app.affiche(mess);
    double tabtaille[3];
    m3d.compare_maillage_carte_isotrope(carte,nomsolution,tabtaille);
    sprintf(mess,"        Nombre de maille théorique : %lf - Nombre de maille reelle : %lf",tabtaille[0],tabtaille[1]);
    app.affiche(mess);
    sprintf(mess,"              Erreur : %lf%%",tabtaille[2]);
    app.affiche(mess);
    MG_SOLUTION *sol=gest.get_mg_solution(gest.get_nb_mg_solution()-1);
    sol->active_solution(2);
    LISTE_MG_NOEUD::iterator it;
    int tab[201];
    for (int i=0;i<201;i++) tab[i]=0;
    for (MG_NOEUD* no=mai->get_premier_noeud(it);no!=NULL;no=mai->get_suivant_noeud(it))
        {
        double val=no->get_solution();
        int valint=(int)val+100;
        if (valint<0) valint=0;
        if (valint>200) valint=200;
        tab[valint]++;
        }
    sprintf(mess,"   Repartition du respect de la taille demandee");
    app.affiche(mess);
    char nomanalyse[500];
    sprintf(nomanalyse,"%s_comp%lu.csv",nometude,mai->get_id());
    FILE *in=fopen(nomanalyse,"wt");
    for (int i=0;i<201;i++) 
      if (i<101) fprintf(in,"%d;%d;\n",i-100,tab[i]); 
      else  fprintf(in,"%d;%d;%d;\n",i-100,tab[i],tab[i]+tab[200-i]);
    fclose(in);
    sprintf(mess,"   Enregistrement");
    app.affiche(mess);
    sprintf(nomsolution,"%s_comp%lu.magic",nometude,mai->get_id());
    gest.enregistrer(nomsolution);
    delete carte;
    }
if (app.get_action()==57)
    {
        int nummai;
        app.recupere_parametre_int_avec_defaut(19,nummai,0);
        if (app.get_erreur()==true) return 0;
        char mess[500];
        sprintf(mess,"   Conversion mg maillage en stl");
        app.affiche(mess);
        MG_FILE gest(nomfichier);
        strncpy(chaine,nomfichier,p-nomfichier);
        MG_MAILLAGE* mai;
        if (nummai==0) mai=gest.get_mg_maillage(0); else mai=gest.get_mg_maillageid(nummai);
        sprintf(chaine,"%s%lu",chaine,mai->get_id());
        std::string namefic=chaine;
        MG_EXPORT exp;
        exp.stl(mai,namefic);
        sprintf(mess,"      creation de %s.stl",namefic.c_str());
        app.affiche(mess);
            
    }     
if (app.get_action()==105)
    {
        int nummai;
        app.recupere_parametre_int_avec_defaut(19,nummai,0);
        if (app.get_erreur()==true) return 0;
        char mess[500];
        sprintf(mess,"   Conversion mg maillage en off");
        app.affiche(mess);
        MG_FILE gest(nomfichier);
        strncpy(chaine,nomfichier,p-nomfichier);
        MG_MAILLAGE* mai;
        if (nummai==0) mai=gest.get_mg_maillage(0); else mai=gest.get_mg_maillageid(nummai);
        sprintf(chaine,"%s%lu",chaine,mai->get_id());
        std::string namefic=chaine;
        MG_EXPORT exp;
        exp.off(mai,namefic);
        sprintf(mess,"      creation de %s.off",namefic.c_str());
        app.affiche(mess);
            
    }     
if (app.get_action()==62)
    {
        int nummai;
        app.recupere_parametre_int_avec_defaut(19,nummai,0);
        char outfile[500];
        app.recupere_parametre_string(5,outfile,(char*)"Fichier de sortie manquant");
        if (app.get_erreur()==true) return 0;
        char mess[500];
        sprintf(mess,"   Détachement d'un mg maillage de sa geometrie");
        app.affiche(mess);
        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 gest2;
        MG_MAILLAGE *mai2=mai->detacher(&gest2);
        sprintf(mess,"   Enregistrement");
        app.affiche(mess);
        gest2.enregistrer(outfile);
    }  
if (app.get_action()==67)
    {
        int numsol;
        int numchamp;
        int numare;
        int numgeo;
        char outfile[500];
        app.recupere_parametre_int_avec_defaut(31,numgeo,0);
        app.recupere_parametre_int_avec_defaut(68,numsol,0);
        app.recupere_parametre_int_avec_defaut(69,numchamp,0);
        app.recupere_parametre_int_avec_defaut(71,numare,0);
        app.recupere_parametre_string(70,outfile,(char*)"Fichier de sortie manquant");
        if (app.get_erreur()==true) return 0;
        MG_FILE gest(nomfichier);
        FEM_SOLUTION* sol;
        if (numsol==0) sol=gest.get_fem_solution(0); else sol=gest.get_fem_solutionid(numsol);
        MG_GEOMETRIE* geo;
        if (numgeo==0) geo=gest.get_mg_geometrie(0); else geo=gest.get_mg_geometrie(numgeo);
        MG_ARETE *are;
        if (numare==0) are=geo->get_mg_arete(0); else are=geo->get_mg_areteid(numare);
        char mess[500];
        sprintf(mess,"   Sonde de %s_%s sur l'arète %lu",(char*)sol->get_nom().c_str(),(char*)sol->get_legende(numchamp).c_str(),are->get_id());
        app.affiche(mess);
        are->recupere_resultat(sol,numchamp,outfile);
        
    }    
if (app.get_action()==72)
    {
        int numsol;
        int numchamp;
        int nummai;
        char outfile[500];
        char infile[500];
        app.recupere_parametre_int_avec_defaut(68,numsol,0);
        app.recupere_parametre_int_avec_defaut(69,numchamp,0);
        app.recupere_parametre_int_avec_defaut(19,nummai,0);
        app.recupere_parametre_string(73,infile,(char*)"Fichier de solution manquant");
        app.recupere_parametre_string_avec_defaut(5,outfile,infile);
        if (app.get_erreur()==true) return 0;
        app.affiche((char*)"Lecture du fichier MAGIC");
        MG_FILE gest(nomfichier);
        app.affiche((char*)"Lecture de la solution a projeter");
        MG_FILE gestsol(infile);
        FEM_SOLUTION* sol;
        if (numsol==0) sol=gestsol.get_fem_solution(0); else sol=gestsol.get_fem_solutionid(numsol);
        FEM_MAILLAGE* fem;
        if (nummai==0) fem=gest.get_fem_maillage(0); else fem=gest.get_fem_maillageid(nummai);
        char mess[500];
        sprintf(mess,"   Projection de la solution %s_%s du maillage %lu sur le maillage %lu",(char*)sol->get_nom().c_str(),(char*)sol->get_legende(numchamp).c_str(),sol->get_maillage()->get_id(),fem->get_id());
        app.affiche(mess);
        FEM_MAILLAGE_OUTILS otfem;
        int res=otfem.projete_solution_maillage(sol,numchamp,fem);
        if (res==MAGIC::PROJECTION::SUCCES) gest.enregistrer(outfile);
        if (res==MAGIC::PROJECTION::DIMENSION_MAILLAGE_INCOMPATIBLE) app.affiche((char*)"   Erreur : Dimension des maillages incompatibles");
        if (res==MAGIC::PROJECTION::TYPE_SOLUTION_INCORRECT) app.affiche((char*)"   Erreur : Solution a projeter n'est pas exprimer aux noeuds");
        if (res==MAGIC::PROJECTION::NON_IMPLEMENTER) app.affiche((char*)"   Erreur : Cas de projection non implementer");
    }    
if (app.get_action()==74)
    {
    int numsol1;
    int numchamp1;  
    int numsol2;
    int numchamp2;
    char outfile[500];
    app.recupere_parametre_int_avec_defaut(40,numsol1,0);
    app.recupere_parametre_int_avec_defaut(41,numsol2,0);
    app.recupere_parametre_int_avec_defaut(43,numchamp1,0);
    app.recupere_parametre_int_avec_defaut(44,numchamp2,0);
    app.recupere_parametre_string_avec_defaut(5,outfile,nomfichier);
    if (app.get_erreur()==true) return 0;
    MG_FILE gest(nomfichier);
    FEM_SOLUTION* sol1;
    if (numsol1==0) sol1=gest.get_fem_solution(0); else sol1=gest.get_fem_solutionid(numsol1);
    FEM_SOLUTION* sol2;
    if (numsol2==0) sol2=gest.get_fem_solution(0); else sol2=gest.get_fem_solutionid(numsol2);
    char mess[500];
    sprintf(mess,"   Comparaison de %s_%s et de %s_%s",(char*)sol1->get_nom().c_str(),(char*)sol1->get_legende(numchamp1).c_str(),(char*)sol2->get_nom().c_str(),(char*)sol2->get_legende(numchamp2).c_str());
    app.affiche(mess);
    FEM_MAILLAGE_OUTILS otfem;
    otfem.compare_champs_solution(sol1,numchamp1,sol2,numchamp2);
    gest.enregistrer(outfile);  
    }
if (app.get_action()==76)
    {
    bool bnoeud,belement1,belement2,belement3,bgelement1,bgelement2,bgelement3,belementnd;  
    int nummai,numtopo;
    char out[500];
    app.recupere_parametre_string(5,out,(char*)"Fichier sortant manquant");
    app.recupere_parametre_int_avec_defaut(19,nummai,0);
    app.recupere_parametre_bool_avec_defaut(77,bnoeud,false);
    app.recupere_parametre_bool_avec_defaut(78,belement1,false);
    app.recupere_parametre_bool_avec_defaut(79,belement2,false);
    app.recupere_parametre_bool_avec_defaut(80,belement3,false);
    app.recupere_parametre_bool_avec_defaut(91,bgelement3,false);
    app.recupere_parametre_bool_avec_defaut(92,bgelement3,false);
    app.recupere_parametre_bool_avec_defaut(93,bgelement3,false);
    app.recupere_parametre_bool_avec_defaut(102,belementnd,false);
    app.recupere_parametre_int_avec_defaut(81,numtopo,0);
    if (app.get_erreur()==true) return 0;
    char mess[500];
    sprintf(mess,"   Extraction de tables de valeurs");
    app.affiche(mess);
    MG_FILE gest(nomfichier);
    FEM_MAILLAGE *fem;
    MG_ELEMENT_TOPOLOGIQUE *topo;
    if (nummai==0) fem=gest.get_fem_maillage(0); else fem=gest.get_fem_maillageid(nummai);
    if (numtopo==0) topo=NULL;
    else 
      {
      MG_GEOMETRIE* geo=fem->get_mg_geometrie();
      topo=geo->get_mg_faceid(numtopo);
      if (topo==NULL) topo=geo->get_mg_volumeid(numtopo);
      if (topo==NULL) topo=geo->get_mg_areteid(numtopo);
      if (topo==NULL) topo=geo->get_mg_sommetid(numtopo);
      }
    MG_TABLE tab;
    if (bnoeud)
        {
        app.affiche((char*)"     Création de la table aux noeuds");
        char nomout[500];
        sprintf(nomout,"%s_noeud.csv",out);
        tab.info_noeud(fem,nomout,topo);
        }
    if (belement1)
        {
        app.affiche((char*)"     Création de la table aux elements 1D");
        char nomout[500];
        sprintf(nomout,"%s_element1.csv",out);
        tab.info_element1(fem,nomout,topo);
        }
    if (belement2)
        {
        app.affiche((char*)"     Création de la table aux elements 2D");
        char nomout[500];
        sprintf(nomout,"%s_element2.csv",out);
        tab.info_element2(fem,nomout,topo);
        }
    if (belement3)
        {
        app.affiche((char*)"     Création de la table aux elements 3D");
        char nomout[500];
        sprintf(nomout,"%s_element3.csv",out);
        tab.info_element3(fem,nomout,topo);
        }
    if (belementnd)
        {
        app.affiche((char*)"     Création de la table aux elements de solutions mixtes");
        char nomout[500];
        sprintf(nomout,"%s_elementnd.csv",out);
        tab.info_element_nd(fem,nomout,topo);
        }
    if (bgelement1)
        {
        app.affiche((char*)"     Création de la table aux point de gauss des elements 1D");
        char nomout[500];
        sprintf(nomout,"%s_elementgauss1.csv",out);
        tab.info_element_gauss1(fem,nomout,topo);
        }
    if (bgelement2)
        {
        app.affiche((char*)"     Création de la table aux point de gauss des elements 2D");
        char nomout[500];
        sprintf(nomout,"%s_elementgauss2.csv",out);
        tab.info_element_gauss2(fem,nomout,topo);
        }
    if (bgelement3)
        {
        app.affiche((char*)"     Création de la table aux point de gauss des elements 3D");
        char nomout[500];
        sprintf(nomout,"%s_elementgauss3.csv",out);
        tab.info_element_gauss3(fem,nomout,topo);
        }
    }
if (app.get_action()==82)
    {
    int numgeo;
    char out[500];
    char in[500];
    char fgmsh[500];
    app.recupere_parametre_string_avec_defaut(1,in,(char*)"");
    if (strcmp(in,"")==0)
      app.recupere_parametre_string(5,out,(char*)"Fichier sortant manquant");
    else 
      app.recupere_parametre_string_avec_defaut(5,out,in);
    app.recupere_parametre_int_avec_defaut(31,numgeo,0);
    app.recupere_parametre_string(83,fgmsh,(char*)"Fichier GMSH manquant");
    if (app.get_erreur()==true) return 0;
    char mess[2000];
    sprintf(mess,"Importation d'un maillage GMSH");
    app.affiche(mess);
    MG_GESTIONNAIRE *gest;
    if (strcmp(in,"")==0) gest=new MG_GESTIONNAIRE; else gest=new MG_FILE(in);
    MG_GEOMETRIE *geo;
    if ((numgeo!=0) && (gest->get_nb_mg_geometrie()>0)) geo=gest->get_mg_geometrieid(numgeo);
    MG_IMPORT imp;
    imp.maillage_gmsh(gest,fgmsh,geo);
    sprintf(mess,"Enregistrement du fichier");
    app.affiche(mess);
    gest->enregistrer(out);
    delete gest;
    }
if ((app.get_action()==86)||(app.get_action()==103))
    {
    int numgeo,numfonc;
    char chemin[1000];
    app.recupere_parametre_int_avec_defaut(31,numgeo,0);
    app.recupere_parametre_string(87,chemin,(char*)"Fichier de fonction manquant");
    if (app.get_action()==103) app.recupere_parametre_int_avec_defaut(104,numfonc,0);
    if (app.get_erreur()==true) return 0;
    char mess[2000];
    if (app.get_action()==86) sprintf(mess,"Importation d'une fonction geometrique");
    if (app.get_action()==103) sprintf(mess,"Modification d'une fonction geometrique");
    app.affiche(mess);
    sprintf(mess,"\nLecture du fichier MAGiC");
    app.affiche(mess);
    MG_FILE gest(nomfichier);
    MG_GEOMETRIE *geo;
    if (numgeo==0) geo=gest.get_mg_geometrie(numgeo); else geo=gest.get_mg_geometrieid(numgeo);
    if (app.get_action()==86)
      {
      sprintf(mess,"Importation de la courbe");
      app.affiche(mess);
      MG_GEOM_FONCTION *gf=new MG_GEOM_FONCTION(chemin);
      geo->ajouter_mg_geom_fonction(gf);
      gest.enregistrer(nomfichier);
      }
    if (app.get_action()==103)
      {
      sprintf(mess,"Importation de la courbe");
      app.affiche(mess);
      MG_GEOM_FONCTION *gf;
      if (numfonc==0) gf=geo->get_mg_geom_fonction(0); else gf=geo->get_mg_geom_fonctionid(numfonc);
      gf->change(chemin);
      gest.enregistrer(nomfichier);
      }
    }
if (app.get_action()==106)
      {
      char mess[255];
      sprintf(mess,"Définition d'une coque mince");
      app.affiche(mess);
      char valeur[500];
      app.recupere_parametre_string(107,valeur,(char*)"Faces en vis a vis manquantes");
      if (app.get_erreur()==true) return 0;
      MG_FILE gest(nomfichier);
      unsigned long id,id1,id2,idgeo;
      sscanf(valeur,"%lu:%lu:%lu:%lu\n",&idgeo,&id,&id1,&id2);
      MG_GEOMETRIE* geo=gest.get_mg_geometrieid(idgeo);
      if (geo!=NULL)
        {
        MG_VOLUME* vol=geo->get_mg_volumeid(id);
        if (vol!=NULL)
            {
            MG_FACE* face1=geo->get_mg_faceid(id1);
            MG_FACE* face2=geo->get_mg_faceid(id2);
            if (face1!=NULL)
              if (face2!=NULL)
              {
              vol->ajouter_face_correspondante(face1,face2);
              sprintf(mess,"     Ajout de la face de base %lu et de la face correspondante %lu dans le volume mince %lu",id1,id2,id);
              app.affiche(mess);
              }
            }
        }
      sprintf(mess,"Enregistrement du fichier");
      gest.enregistrer(nomfichier);
      }
      
char mess[255];
sprintf(mess,"Fin");
app.affiche(mess);
} 
 
 
Revision:	881
Committed:	Wed Apr 12 23:14:07 2017 UTC (8 years, 2 months ago) by francois
File size:	78859 byte(s)
Log Message:	Changement de strategie pour le mailleur couche. Il fonctionne maintenant avec des volumes minces et non des coquilles mince est une propriète de volume et ne change pas son nom. cela necessite un changement de version de fichier -> maintenant 2.3 compatible en arrière valider sur ancien exemple plus sphere creuse