pre-pos_processing_gmsh-herezh

#!/usr/bin/perl

use strict;

##############################################################################################

# script pour passer du format de maillage genere par gmsh au format lu par herezh++ #

# version 1: version initiale #

# version 1.1: 30 avr 2008 prise en compte du fait qu'il peut ne pas y avoir de reference #

# de type chaine de caractères #

# version 1.2: 9 mai 2008 Dans le cas où on ne conserve pas les éléments noeuds, on conserve#

# quand m^eme les références associées. Au niveau des tag on différencie entre #

# 1) les physical name (c'est-à-dire les ref que l'on met #

# 2) le numéro de l'entité géométrique élémentaires #

# 3) le numéro de la partition de maillage #

# version 1.3: 13 mai 2008 #

# 1) introduction de références d'aretes et de surfaces pour les #

# éléments lignes et surfaces (uniquement) #

# 2) construction de références de noeuds correspondants aux refs d'éléments #

# que ces éléments soient concervés ou non #

# version 1.4: 23 mai 2008 #

# correction d'un bug sur l'affichage du nombre total d'element #

# version 1.5: 6 fevrier 2009 #

# correction de la numérotation pour les éléments quadratiques #

# suppression blancs et tab en debut de ligne, en lecture du .msh #

# version 1.6: 13 avril 2010 - S. Bochard #

# maillage msh au format 2.1 depuis gmsh-2.4 #

# correction lecture references pour sorties générées par gmsh-2.4 #

# version 1.7: 5 fevrier 2011 #

# re-correction de la numérotation pour les éléments quadratiques #

# hexaedre incomplet, introduction du nombre de pti #

# version 1.8: 19 octobre 2011 - S. Bochard #

# maillage msh au format 2.2 depuis gmsh-2.5 #

# correction format maillage 2.1 et supérieur pour gmsh 2.4 et 2.5 #

# version 1.9: 1 fevrier 2012 - G. Rio #

# correction bug sur la sortie des tetraedres quadratiques #

# version 2.0: 10 decembre 2012 - G. Rio #

# correction bug sur la sortie des tetraedres quadratiques #

# version 2.1: 15 decembre 2012 - G. Rio #

# correction bug sur la sortie des tetraedres quadratiques #

# version 2.2: 2 decembre 2015 - G. Rio #

# ajout du nom du type d'élément, quand question nombre de pti #

# version 2.3: 19 fevrier 2017 - G. Rio #

# sortie sytématique du nom du fichier traité #

# version 2.4: 26 novembre 2018 - G. Rio #

# modification de POUT QUADRATIQUE en pout BIE2 #

##############################################################################################

my $nom;

my $version = "2.4";

if (@ARGV[0] eq '-h')

# si on n'a pas d'argument, on donne le mode d'emploi

{

print ("\n");

print ("\n" . '#############################################################################');

print ("\n" . '# #');

print ("\n" . '#msh2her.pl #');

print ("\n" . '# transformation d\'un maillage genere par le programme gmsh #');

print ("\n" . '# en un maillage utilisable par Herezh++ #');

print ("\n" . '#############################################################################');

print ("\n" . '# licence GPL, gerard.rio@univ-ubs.fr #');

print ("\n" . '# http://www-lg2m.univ-ubs.fr/~rio #');

print ("\n" . '# version '."$version".' #');

print ("\n" . '#############################################################################');

print ("\n\n\n");

print ("\n usage: msh2her.pl ");

print ("\n ----> entree interactive des noms de fichier \n");

print ("\n usage: msh2her.pl <nom fichier.msh> <nom fichier.her> \n");

print ("\n ----> indication de la source et de la cible ");

print ("\n NB: les noms de fichier peuvent avoir ou non leur extension \n");

print ("\n usage: msh2her.pl -v ");

print ("\n ----> affiche la version du programme \n");

print ("\n usage: msh2her.pl -h ");

print ("\n ----> information sur l'utilisation \n");

print ("\n\n");

die ("arret normal");

}

elsif (@ARGV[0] eq '-v')

{ print ("\n msh2her.pl: version "."$version"."\n\n");

die ("arret normal");

};

#------------------- écriture de l'entete à l'écran -----------------

print ("\n");

print ("\n" . '#############################################################################');

print ("\n" . '# #');

print ("\n" . '#msh2her.pl #');

print ("\n" . '# transformation d\'un maillage genere par le programme gmsh #');

print ("\n" . '# en un maillage utilisable par Herezh++ #');

print ("\n" . '#############################################################################');

print ("\n" . '# licence GPL, gerard.rio@univ-ubs.fr #');

print ("\n" . '# http://www-lg2m.univ-ubs.fr/~rio #');

print ("\n" . '# version '."$version".' #');

print ("\n" . '#############################################################################');

print ("\n");

#------------------- fin écriture de l'entete à l'écran -----------------

#------------------- on s'occupe tout d'abord des fichiers --------------------

# tout d'abord on récupère le nom du fichier à traiter

my $nombase;

if (@ARGV[0])

# si on a un argument, c'est le nom du fichier

{ $nombase = @ARGV[0];}

else

# sinon on demande interactivement le nom du fichier

{ print ("\n nom du fichier ? ");

$nombase = <STDIN>; chomp($nombase); # on supprime le retour à la ligne

};

# on supprime l'extension .msh si elle existe

$nombase =~ s/\.msh//;

# on ouvre le fichier msh en lecture tout d'abord, il doit donc exister

my $fichier_msh = $nombase . '.msh';

open (fic_msh,$fichier_msh) || die (" \n erreur en ouverture du fichier " . $fichier_msh . " \n");

print (" --> traitement du fichier $fichier_msh ---\n ");

# on ouvre le fichier .her en écriture, on écrase donc le fichier existant

my $fichier_her;

if (@ARGV[1])

# si on a un argument 2, c'est le nom du fichier de sortie

{ $fichier_her = @ARGV[1];}

else

# sinon on construit le nom du fichier

{$fichier_her = $nombase . '.her';

# on regarde s'il existe déjà

foreach (<*.her>)

{ if ($_ eq $fichier_her)

{ print ("\n le fichier " . $fichier_her . " existe deja ! ");

print ("\n voulez-vous l\' ecraser ? (0 ou N ) ");

my $rep = <STDIN>; chomp($rep); $rep=lc($rep); # lc passage en minuscule

if ($rep eq 'n')

{print ("\n nom du fichier de sortie ? ");

$fichier_her = <STDIN>; chomp($fichier_her);

};

};

}; # fin du foreach

};

# on supprime l'extension .her si elle existe

$fichier_her =~ s/\.her//;

# et on la rajoute

$fichier_her .= '.her';

# maintenant on tente l'ouverture de fichier

open (fic_her,">$fichier_her") || die (" \n erreur en ouverture du fichier " . $fichier_her . " \n");

#------------------------ fin des ouvertures de fichiers ---------------------

#------------------ lecture des tag qui vont servir pour fabriquer les références -------------

print ("\n lecture des tags (preparation a la creation des references) ");

# on lit dans les lignes jusqu'au mot clé $PhysicalNames

my $ligne;my @champ; my $existe_physicalNames = 0;

my $position_sur_node = 0; # pour signaler que l'on est à la position des noeuds

#

# debut des modifications - Stéphane

#

my $mesh_format; # ajout Stéphane

while (<fic_msh>)

{ $ligne = $_; chomp($ligne); # sauvegarde de la ligne

$ligne =~ s/^[ \t]+//mg; # suppression des blancs et tabulation en début de ligne

@champ = split (/\s+/,$ligne); # le /s+/ permet de séparer sur des espaces

# on lit dans le fichier jusqu'au mot clé $PhysicalNames

if ( $champ[0] eq '$MeshFormat' ) # ajout Stéphane

{$ligne = <fic_msh>; chomp($ligne); # sauvegarde de la ligne suivante # ajout Stéphane

$ligne =~ s/^[ \t]+//mg; # suppression des blancs et tabulation en début de ligne # ajout Stéphane

@champ = split (/\s+/,$ligne); # le /s+/ permet de séparer sur des espaces # ajout Stéphane

$mesh_format = $champ[0]; # on recupere le format du maillage # ajout Stéphane

print ("\n format de sauvegarde du maillage : msh " . $mesh_format . "\n" ); # ajout Stéphane

} # ajout Stéphane

if ( $champ[0] eq '$PhysicalNames' )

{$existe_physicalNames = 1;# on a trouvé le premier enreg,

last;

}

elsif ( $champ[0] eq '$Nodes' )

{$existe_physicalNames = 0;# on a pas trouvé de $PhysicalNames

print "\n **** NB: pas de chaine de caractere pour references a lire *****\n";

$position_sur_node = 1; # drapeau qui signale que l'on est au début des noeuds

last;

};

};

# --- on lit les tag si il en existe

my %les_references = () ; # déclaration de la table des reference

if ($existe_physicalNames)

{ $ligne = <fic_msh>; chomp($ligne);

$ligne =~ s/^[ \t]+//mg; # suppression des blancs et tabulation en début de ligne

@champ = split (/\s+/,$ligne); # le /s+/ permet de séparer sur des espaces

print (" nombre total de references lus : " . $champ[0] );

my $nb_reference = $champ[0] ;

# --- on lit chaque ligne de reference que l'on recopie dans un tableau jusqu'au mot clé $EndPhysicalNames

while (<fic_msh>)

{ $ligne = $_; chomp($ligne); # sauvegarde de la ligne

$ligne =~ s/^[ \t]+//mg; # suppression des blancs et tabulation en début de ligne

@champ = split (/\s+/,$ligne); # le /s+/ permet de séparer sur des espaces

# on lit dans le fichier jusqu'au mot clé $EndPhysicalNames

if ( $champ[0] eq '$EndPhysicalNames' )

{last;}; # on a trouvé le mot cle de fin de reference,

# ce n'est pas le mot clé, donc on enregistre la reference

if ($mesh_format eq "2") # ajout Stéphane

{$champ[1] =~ s/\"+//g; # on retire les guillemets éventuelles # ajout Stéphane

$les_references{"$champ[0]"}= "$champ[1]"; # ajout Stéphane

} # ajout Stéphane

elsif ($mesh_format ge "2.1") # ajout Stéphane # modif version 1.8

{$champ[2] =~ s/\"+//g; # on retire les guillemets éventuelles # ajout Stéphane

$les_references{"$champ[1]"}= "$champ[2]"; # ajout Stéphane

}; # ajout Stéphane

#

# fin des modifications - Stéphane

#

};

print "\n les references ";

while ( (my $cle, my $valeur) = each %les_references)

{ print "\n ".$cle . ' => ' .$valeur;};

};

#------------------ fin de la lecture des tag qui vont servir pour fabriquer les références -------------

# ++++++++++++++++++++++ lecture des coordonnées des noeuds ++++++++++++++++++++++++++

print ("\n lecture des coordonnees des noeuds ");

# on lit dans les lignes jusqu'au mot clé $Nodes

my $ligne;my @champ;

if ($position_sur_node != 1)

{ while (<fic_msh>)

{ $ligne = $_; chomp($ligne); # sauvegarde de la ligne

$ligne =~ s/^[ \t]+//mg; # suppression des blancs et tabulation en début de ligne

@champ = split (/\s+/,$ligne); # le /s+/ permet de séparer sur des espaces

# on lit dans le fichier jusqu'au mot clé $Nodes

if ( $champ[0] eq '$Nodes' )

{$position_sur_node = 1;

last;}; # on a trouvé le premier enreg,

};

if ($position_sur_node != 1)

{ die ( "\n on n'a pas trouve le mot cle \$Nodes ?? , on ne peut pas continuer \n");

};

};

# --- on lit le nombre de noeud

$ligne = <fic_msh>; chomp($ligne);

$ligne =~ s/^[ \t]+//mg; # suppression des blancs et tabulation en début de ligne

@champ = split (/\s+/,$ligne); # le /s+/ permet de séparer sur des espaces

print (" nombre total de noeuds lus : " . $champ[0] );

my $nb_noeud = $champ[0] ;

# --- on lit chaque ligne de coordonnée que l'on recopie dans un tableau jusqu'au mot clé $EndNodes

my @les_noeuds = () ; # déclaration du tableau des coordonnées des noeuds

while (<fic_msh>)

{ $ligne = $_; chomp($ligne); # sauvegarde de la ligne

$ligne =~ s/^[ \t]+//mg; # suppression des blancs et tabulation en début de ligne

@champ = split (/\s+/,$ligne); # le /s+/ permet de séparer sur des espaces

# on lit dans le fichier jusqu'au mot clé $EndNodes

if ( $champ[0] eq '$EndNodes' )

{last;}; # on a trouvé le premier enreg,

# ce n'est pas le mot clé, donc on enregistre la ligne

push @les_noeuds, $ligne; # sauvegarde des coordonnées et du numéro initial du noeud

};

# ++++++++++++++++++++++ fin de lecture des coordonnees noeuds ++++++++++++++++++++++++++++++++++++

#-------- def de tableaux de conversion entre les numérotations gmsh et celles d'herezh++ ----------

# définition des types d'éléments

my @type_element = ( 'rien pour le numero 0', # '0'

'segment lineaire a 2 noeuds' , # '1'

'triangle lineaire a 3 noeuds' , # '2'

'quadrangle bilineaire a 4 noeuds' , # '3'

'tetrahedre lineaire a 4 noeuds' , # '4'

'hexahedre trilineaire a 8 noeuds' , # '5'

'pentaedre bilineaire a 6 noeuds' , # '6'

'pyramide a 5 noeuds' , # '7'

'segment quadratique a 3 noeuds' , # '8'

'triangle quadratique a 6 noeuds' , # '9'

'quadrangle quadratique complet a 9 noeuds' , #'10'

'tetraedre quadratique a 10 noeuds' , #'11'

'hexaedre quadratique complet a 27 noeuds' , #'12'

'pentaedre quadratique complet a 18 noeuds' , #'13'

'pyramide quadratique a 14 noeuds' , #'14'

' point a 1 noeud ' , #'15'

'quadrangle quadratique incomplet a 8 noeuds' , #'16'

'hexaedre quadratique incomplet a 20 noeuds' , #'17'

'pentaedre quadratique incomplet a 15 noeuds' , #'18'

'pyramide quadratique incomplete a 13 noeuds' #'19'

);

my @type_herezh_element = ( 'rien pour le numero 0', # '0'

' POUT BIE1 ' , # '1'

' TRIANGLE LINEAIRE ' , # '2'

' QUADRANGLE LINEAIRE ' , # '3'

' TETRAEDRE LINEAIRE ' , # '4'

' HEXAEDRE LINEAIRE ' , # '5'

' PENTAEDRE LINEAIRE ' , # '6'

' PYRAMIDE LINEAIRE ' , # '7'

' POUT BIE2 ' , # '8'

' TRIANGLE QUADRACOMPL ' , # '9'

' QUADRANGLE QUADRACOMPL ' , #'10'

' TETRAEDRE QUADRACOMPL ' , #'11'

' HEXAEDRE QUADRACOMPL ' , #'12'

' PENTAEDRE QUADRACOMPL ' , #'13'

' PYRAMIDE QUADRACOMPL ' , #'14'

' POINT LINEAIRE ' , #'15'

' QUADRANGLE QUADRATIQUE ' , #'16'

' HEXAEDRE QUADRATIQUE ' , #'17'

' PENTAEDRE QUADRATIQUE ' , #'18'

' PYRAMIDE QUADRATIQUE ' #'19'

);

# les types et nombre de points d'intégration pour chaque éléments:

# le premier chiffre indique le nombre de pti par défaut, si 0 : il n'y a pas plusieurs choix: un seul cas par défaut

# puis les différents cas

my @type_herezh_integration = ( [0], #'rien pour le numero 0', # '0'

[0], # ' POUT BIE1 ' , # '1'

[0], # ' TRIANGLE LINEAIRE ' , # '2'

[0], # ' QUADRANGLE LINEAIRE ' , # '3'

[0], # ' TETRAEDRE LINEAIRE ' , # '4'

[8,1,27,64], # ' HEXAEDRE LINEAIRE ' , # '5'

[2,1,6], # ' PENTAEDRE LINEAIRE ' , # '6'

[0], # ' PYRAMIDE LINEAIRE ' , # '7'

[0], # ' POUT BIE2 -> QUADRATIQUE ' , # '8'

[0], # ' TRIANGLE QUADRACOMPL ' , # '9'

[0], # ' QUADRANGLE QUADRACOMPL ' , #'10'

[0], # ' TETRAEDRE QUADRACOMPL ' , #'11'

[8,1,27,64], # ' HEXAEDRE QUADRACOMPL ' , #'12'

[6,9,12,18], #' PENTAEDRE QUADRACOMPL ' , #'13'

[0], # ' PYRAMIDE QUADRACOMPL ' , #'14'

[0], # ' POINT LINEAIRE ' , #'15'

[0], # ' QUADRANGLE QUADRATIQUE ' , #'16'

[8,1,27,64], # ' HEXAEDRE QUADRATIQUE ' , #'17'

[6,9,12,18], # ' PENTAEDRE QUADRATIQUE ' , #'18'

[0], # ' PYRAMIDE QUADRATIQUE ' #'19'

);

my @nbn_par_elem;

$nbn_par_elem[0]=0; $nbn_par_elem[1]=2; $nbn_par_elem[2]=3; $nbn_par_elem[3]=4;

$nbn_par_elem[4]=4; $nbn_par_elem[5]=8; $nbn_par_elem[6]=6; $nbn_par_elem[7]=5;

$nbn_par_elem[8]=3; $nbn_par_elem[9]=6; $nbn_par_elem[10]=9; $nbn_par_elem[11]=10;

$nbn_par_elem[12]=27; $nbn_par_elem[13]=18; $nbn_par_elem[14]=14; $nbn_par_elem[15]=1;

$nbn_par_elem[16]=8; $nbn_par_elem[17]=20; $nbn_par_elem[18]=15; $nbn_par_elem[19]=13;

# $connex[i][j] donne pour l'element i, le num de gmsh pour le numéro j d'herezh

my @connex = ( [], # le premier avec rien, pour que la numérotation commence à 1

[1,2], # 1: ' POUT BIE1 ' ,

[1,2,3], # 2: ' TRIANGLE LINEAIRE ' ,

[1,2,3,4], # 3: ' QUADRANGLE LINEAIRE ' ,

[1,2,3,4], # 4: ' TETRAEDRE LINEAIRE ' ,

[1,2,3,4,5,6,7,8], # 5: ' HEXAEDRE LINEAIRE ' ,

[1,2,3,4,5,6], # 6: ' PENTAEDRE LINEAIRE ' ,

[1,2,3,4,5], # 7: ' PYRAMIDE LINEAIRE ' ,

[1,3,2], # 8: ' POUT BIE2 ' # pour QUADRATIQUE ' ,

[1,2,3,4,5,6], # 9: ' TRIANGLE QUADRACOMPL ' ,

[1,2,3,4,5,6,7,8,9], # 10: ' QUADRANGLE QUADRACOMPL ' ,

[1,2,3,4,5,6,7,8,10,9], # 11: ' TETRAEDRE QUADRACOMPL ' ,

[1,2,3,4,5,6,7,8,9,12,14,10,11,13,15,16,17,19,20,18,21,22,24,25,23,26,27], # 12: ' HEXAEDRE QUADRACOMPL ' ,

# [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27] uniquement pour voir les numéros

[1,2,3,4,5,6,7,10,8,9, 11,12,13,15,14,16,18,17], # 13:' PENTAEDRE QUADRACOMPL ' ,

[1,2,3,4,5,6,7,8,9,10,11,12,13,14], # 14:' PYRAMIDE QUADRACOMPL ' ,

[1], # 15: ' POINT LINEAIRE ' ,

[1,2,3,4,5,6,7,8], # 16: ' QUADRANGLE QUADRATIQUE ' ,

[1,2,3,4,5,6,7,8,9,12,14,10,11,13,15,16,17,19,20,18], # 17: ' HEXAEDRE QUADRATIQUE ' ,

# [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20], # 17: ' HEXAEDRE QUADRATIQUE ' ,

[1,2,3,4,5,6,7,10,8,9, 11,12,13,15,14], # 18:' PENTAEDRE QUADRATIQUE ' ,

[1,2,3,4,5,6,7,8,9,10,11,12,13], # 19: ' PYRAMIDE QUADRATIQUE '

);

my @type_geom =

( 'RIEN_TYPE_GEOM', # '0'

'LIGNE' , # '1'

'SURFACE' , # '2'

'SURFACE' , # '3'

'VOLUME' , # '4'

'VOLUME' , # '5'

'VOLUME' , # '6'

'VOLUME' , # '7'

'LIGNE' , # '8'

'SURFACE' , # '9'

'SURFACE' , #'10'

'VOLUME' , #'11'

'VOLUME' , #'12'

'VOLUME' , #'13'

'VOLUME' , #'14'

'POINT_G' , #'15'

'SURFACE' , #'16'

'VOLUME' , #'17'

'VOLUME' , #'18'

'VOLUME' #'19'

);

#-------- fin def de tableaux de conversion entre les numérotations gmsh et celles d'herezh++ ----------

# ++++++++++++++++++++++ lecture des éléments +++++++++++++++++++++++++++++++++++++++++++++++++++++++++

print ("\n lecture des elements ");

while (<fic_msh>)

{ $ligne = $_; chomp($ligne); # sauvegarde de la ligne

$ligne =~ s/^[ \t]+//mg; # suppression des blancs et tabulation en début de ligne

@champ = split (/\s+/,$ligne); # le /s+/ permet de séparer sur des espaces

# on lit dans le fichier jusqu'au mot clé $Elements

if ( $champ[0] eq '$Elements' )

{last;}; # on a trouvé le premier enreg,

};

# ----- on lit le nombre d'éléments

$ligne = <fic_msh>; chomp($ligne);

$ligne =~ s/^[ \t]+//mg; # suppression des blancs et tabulation en début de ligne

@champ = split (/\s+/,$ligne); # le /s+/ permet de séparer sur des espaces

print (" nombre total d'elements lus : " . $champ[0] );

my $nb_element = $champ[0] ;

# --- on lit chaque ligne que l'on recopie dans des tableaux jusqu'au mot clé $EndElements

my @liste_element = (); # ce sera un tableau multidimentionel

# chaque élément de @liste_element est une liste d'éléments d'un m^eme type

# $liste_element[i] = la liste des éléments de type i

# $liste_element[i][j] un élément = une liste d'entiers donc un tableau

# $liste_element [i][j][k] = un scalaire (entier) ouf !!

my @number_type = (); # les nombres d'élément de chaque type

while (<fic_msh>)

{ $ligne = $_; chomp($ligne); # sauvegarde de la ligne

$ligne =~ s/^[ \t]+//mg; # suppression des blancs et tabulation en début de ligne

@champ = split (/\s+/,$ligne); # le /s+/ permet de séparer sur des espaces

# on lit dans le fichier jusqu'au mot clé $EndElements

if ( $champ[0] eq '$EndElements' )

{last;}; # on a trouvé le premier enreg,

# traitement et enregistrement

push @{$liste_element[$champ[1]]}, [@champ ]; # tableau multidimentionnel, syntaxe particulière

#if ( $champ[1] == 15)

# {print "\n trouve " . @champ; print " ligne = " . $ligne;

# print ",,, " .$champ[0]." ".$champ[1]." ".$champ[2]." ".$champ[3]." ".$champ[4]." ".$champ[5]." ";};

$number_type[$champ[1]]++;

};

# test débug

#print "\n essai= " . $liste_element[15][1][3] ;

#print "\n une tranche

# on recherche le nombre d'élément différent

my $nb_type_elem = 0;

my $taille_list_elem = $#liste_element;

#print "\n number_type = " . @number_type;

for (my $j=1;$j<=$taille_list_elem;$j++)

{ if ($number_type[$j] != 0)

{ print "\n nombre d'elements du type ".$type_element[$j] ." = " . $number_type[$j] ;

$nb_type_elem++;

};

};

print "\n nombre initial d'elements differents " . $nb_type_elem ;

# maintenant on demande si l'on veut tous les éléments ou certain uniquement

my @number_type_initial = @number_type; # sauvegarde pour le traitement des références

print "\n\n === choix d'elements: voulez-vous ? ";

for (my $j=1;$j<=$taille_list_elem;$j++)

{ if ($number_type[$j] != 0)

{ print "\n conserver les elements du type ".$type_element[$j] ." (rep o ou n) " ;

my $rep = <STDIN>; chomp($rep); $rep=lc($rep); # lc passage en minuscule

if ($rep eq 'o')

{print (' ===> ok, on conserve les elements '.$type_element[$j] . "\n");}

else

{print (' *** on supprime les elements '.$type_element[$j] .'*** '. "\n");

$number_type[$j] = 0;

}

};

};

# on demande éventuellement le nombre de point d'intégration

my @type_out_integration; # le tableau contenant les types d'integration

for (my $j=1;$j<=$taille_list_elem;$j++)

{ if ($number_type[$j] != 0) # on ne traite que les éléments que l'on a gardé

{if ($type_herezh_integration[$j] == 0)

{$type_out_integration[$j] = ' '; }#push @type_out_integration, ' ';} # rien pour la sortie

else

{ # cas où plusieurs type de point d'intégration sont possible

print ($type_element[$j].": nombre de points d'integration par defaut ".$type_herezh_integration[$j][0]." ? (rep o/n) ");

my $rep = <STDIN>; chomp($rep); $rep=lc($rep); # lc passage en minuscule

if ($rep eq 'o')

{$type_out_integration[$j] = ' '; }

else

# cas où plusieurs choix existent

{ print (" valeur preconisee parmi: ");

foreach my $nbpti ( @{$type_herezh_integration[$j]})

{ print ( $nbpti . " ");};

print (" : ? ");

my $rep = <STDIN>; chomp($rep);

if ($rep != $type_herezh_integration[$j][0])

{$type_out_integration[$j] = ' _cm'.$rep.'pti ';}

else {$type_out_integration[$j] = ' '; };

};

print ("\n");

}

};

};

# ++++++++++++++++++++++ fin lecture des éléments ++++++++++++++++++++++++++

# ++++++++++++++++++++++ création des références ++++++++++++++++++++++++++++++

# les références de noeuds se trouvent par l'intermédiaire des éléments "points"

# on passe donc en revue ces éléments

my %ref_noeud; # def de conteneur des références de noeuds

my %ref_element; # def du conteneur des références d'éléments

my %ref_arete; # def du conteneur des références d'arretes

my %ref_surface; # def du conteneur des références de surfaces

#print "\n taille de la liste d'elem = " . $taille_list_elem;

my $numero = 1;

for (my $j=1;$j<=$taille_list_elem;$j++)

{ #print "\n taille des elem = " . $number_type[$j] ;

# !!! la sélection est maintenant faite à l'intérieur de la boucle

# c'est moins performant mais cela permet d'avoir des ref de points de ligne, surface, volume

# m^eme si ces derniers ne sont pas sélectionné comme éléments à retenir !!

# if (($number_type[$j] != 0) || ($liste_element[$j][0][1] == 15))

{ my $nb_el = $number_type_initial[$j];

#print "\n nb_el= ".$nb_el;

for (my $el = 0;$el < $nb_el;$el++)

{

# ------- cas des références de noeuds -----

if ($liste_element[$j][$el][1] == 15) # 15 c'est le type point

{ # on balaie les tag pour ajouter

my $nbtag = $liste_element[$j][$el][2]; # récup du nombre de tag

# on différencie les deux derniers tags qui sont :

# 1) le numéro de l'entité géométrique élémentaires

# 2) le numéro de la partition de maillage

# en fait pour les points il semble qu'il n'y ait que un tag sup, donc on va jusqu'à $nbtag-1

# --- tout d'abord les tag des noms physiques (ceux qui ne sont pas généré automatiquement)

for (my $ic = 0; $ic < $nbtag-1; $ic++)

{ my $indic_ref = $liste_element[$j][$el][3+$ic]; # l'indicateur de la référence

#print "\n indicateur de la reference : ".$indic_ref;

my $nbnoe = $liste_element[$j][$el][3+ $nbtag]; # le numéro du noeud

# $les_references{"$indic_ref"} = le nom de la référence, si ce n'est pas directement le nombre

my $nom = $les_references{"$indic_ref"};

#print "\n l equivalent d'indicateur de ref ".$nom;

if ($nom ne undef)

{ # cas ou le nom de la référence existe

if ( $ref_noeud{'N_' . $nom} == undef)

# si la référence n'est pas définie on la définie

{ my @tab = ($nbnoe);

#print "\n la cle = " . 'N_' . $nom;

$ref_noeud{'N_' . $nom} = [@tab]; # on sauvegarde le tableau

}

else

{ #print "\n la cle = " . 'N_' . $nom;

push @{$ref_noeud{'N_' . $nom}}, $nbnoe;

};

}

else

{ # cas ou le nom de la référence n'existe pas, on utilise directement le nombre

if ( $ref_noeud{'N_' ."$indic_ref"} == undef)

# si la référence n'est pas définie on la définie

{ my @tab = ($nbnoe);

$ref_noeud{'N_' ."$indic_ref"} = [@tab]; # on sauvegarde le tableau

}

else

{ push @{$ref_noeud{'N_' ."$indic_ref"}}, $nbnoe;

};

};

};

# --- puis les 2 tags automatique:

# 1) le numéro de l'entité géométrique élémentaires

# 2) le numéro de la partition de maillage

# en fait pour les points il semble qu'il n'y ait que un tag sup, donc on va jusqu'à $nbtag-1

my $geoParti = 1; # pour choisir entre les deux

for (my $ic = $nbtag-1 ; $ic < $nbtag; $ic++,$geoParti++)

{ my $indic_ref = $liste_element[$j][$el][3+$ic]; # l'indicateur de la référence

#print "\n indicateur de la reference : ".$indic_ref;

my $nbnoe = $liste_element[$j][$el][3+ $nbtag]; # le numéro du noeud

# ici il s'agit exclusiment de nombre, on l'utilise directement

# on signale la différence avec $geoParti

if ( $ref_noeud{'N_'."$geoParti".'_'."$indic_ref"} == undef)

# si la référence n'est pas définie on la définie

{ my @tab = ($nbnoe);

$ref_noeud{'N_'."$geoParti".'_'."$indic_ref"} = [@tab]; # on sauvegarde le tableau

}

else

{ push @{$ref_noeud{'N_'."$geoParti".'_'."$indic_ref"}}, $nbnoe;

};

};

}

else # diff de 15 donc éléments normal

# ------- cas des références d'éléments -----

# on va également créer des références de noeuds correspondants à chaques références particulières

# d'éléments ce qui permettra d'avoir des ref des noeuds de lignes, de surfaces ou de volumes particulières

{ # on balaie les tag pour ajouter

my $tyyp = $liste_element[$j][$el][1]; # le type de l'élément

my $nbtag = $liste_element[$j][$el][2]; # récup du nombre de tag

# on différencie les deux derniers tags qui sont :

# 1) le numéro de l'entité géométrique élémentaires

# 2) le numéro de la partition de maillage

# --- tout d'abord les tag des noms physiques (ceux qui ne sont pas généré automatiquement)

my %tag_deja_traite;

my $fin_tag = $nbtag-2;# ajout Stéphane version 1.8

if ($mesh_format ge "2.2") {$fin_tag++;};# ajout Stéphane version 1.8

for (my $ic = 0; $ic < $fin_tag; $ic++) # ajout Stéphane version 1.8

# for (my $ic = 0; $ic < $nbtag-2; $ic++) # modif Stéphane version 1.8

{ my $indic_ref = $liste_element[$j][$el][3+$ic]; # l'indicateur de la référence

# dans le cas où il y a un nom de référence, l'indicateur de référence est doublé

# ce qui n'est pas normal, on supprime ce doublement

if ($tag_deja_traite{"$indic_ref"} ne undef)

{next;}

$tag_deja_traite{"$indic_ref"} = '1'; # de maniere à pas avoir de undef

#print "\n indicateur de la reference : ".$indic_ref;

# $les_references{"$indic_ref"} = le nom de la référence, si ce n'est pas directement le nombre

my $nom = $les_references{"$indic_ref"};

#print "\n l equivalent d'indicateur de ref ".$nom;

if ($nom ne undef)

{ # cas ou le nom de la référence d'élément existe (on applique les m^eme règles pour les N associés)

if ( $ref_element{'E_' . $nom} == undef)

# si la référence n'est pas définie on la définie

{

#print "\n la cle = " . 'E_' . $nom;

# on n'intervient pour les ref d'éléments que si ces éléments sont à retenir

if ($number_type[$j] != 0)

{ my @tab = ($numero);

$ref_element{'E_' . $nom} = [@tab]; # on sauvegarde le tableau

# idem pour la ref de surfaces associés, mais ici dans le cas ou c'est un élément de surface

#print "\n le type de geom = " . 'E_' . $nom;

if ($type_geom[$tyyp] eq "SURFACE")

{ my @tab_S = ($numero);

$ref_surface{'F_' . $nom} = [@tab_S]; # on sauvegarde le tableau

};

# idem pour la ref d'arete associés, mais ici dans le cas ou c'est un élément linéique

if ($type_geom[$tyyp] eq "LIGNE")

{ my @tab_L = ($numero);

$ref_arete{'A_' . $nom} = [@tab_L]; # on sauvegarde le tableau

};

};

}

else

{ #print "\n la cle = " . 'E_' . $nom;

# on n'intervient pour les ref d'éléments que si ces éléments sont à retenir

if ($number_type[$j] != 0)

{ push @{$ref_element{'E_' . $nom}}, $numero;

# idem pour la ref de surfaces associés, mais ici dans le cas ou c'est un élément de surface

if ($type_geom[$tyyp] eq "SURFACE")

{ push @{$ref_surface{'F_' . $nom}}, $numero;};

# idem pour la ref d'arete associés, mais ici dans le cas ou c'est un élément linéique

if ($type_geom[$tyyp] eq "LIGNE")

{ push @{$ref_arete{'A_' . $nom}}, $numero;};

};

};

# idem pour la ref de noeuds associés, mais ici dans tous les cas

if ( $ref_noeud{'N_'."$type_geom[$tyyp]".'_'.$nom} == undef)

{ my @t_noeud_ref_elem; # on sauvegarde le tableau (mais ici pour l'instant sans rien dedans)

$ref_noeud{'N_'."$type_geom[$tyyp]".'_'.$nom} = [@t_noeud_ref_elem];

};

# on rempli le tableau de num de noeud pour la ref de noeuds associé

for (my $ia =0;$ia < $nbn_par_elem[$tyyp]; $ia++)

{ my $bb = $liste_element[$j][$el][3+ $nbtag + $ia ]; # le numero de noeud

push @{$ref_noeud{'N_'."$type_geom[$tyyp]".'_'.$nom}}, $bb;

};

}

else

{ # cas ou le nom de la référence n'existe pas, on utilise directement le nombre

if ( $ref_element{'E_' ."$indic_ref"} == undef)

# si la référence n'est pas définie on la définie

{ # on n'intervient pour les ref d'éléments que si ces éléments sont à retenir

if ($number_type[$j] != 0)

{ my @tab = ($numero);

$ref_element{'E_' ."$indic_ref"} = [@tab]; # on sauvegarde le tableau

# idem pour la ref de surfaces associés, mais ici dans le cas ou c'est un élément de surface

if ($type_geom[$tyyp] eq "SURFACE")

{ my @tab_S = ($numero);

$ref_surface{'F_' ."$indic_ref"} = [@tab_S]; # on sauvegarde le tableau

};

# idem pour la ref d'arete associés, mais ici dans le cas ou c'est un élément linéique

if ($type_geom[$tyyp] eq "LIGNE")

{ my @tab_L = ($numero);

$ref_arete{'A_' . "$indic_ref"} = [@tab_L]; # on sauvegarde le tableau

};

};

}

else

{ # on n'intervient pour les ref d'éléments que si ces éléments sont à retenir

if ($number_type[$j] != 0)

{ push @{$ref_element{'E_' ."$indic_ref"}}, $numero;

# idem pour la ref de surfaces associés, mais ici dans le cas ou c'est un élément de surface

if ($type_geom[$tyyp] eq "SURFACE")

{ push @{$ref_surface{'F_' . "$indic_ref"}}, $numero;};

# idem pour la ref d'arete associés, mais ici dans le cas ou c'est un élément linéique

if ($type_geom[$tyyp] eq "LIGNE")

{ push @{$ref_arete{'A_' . "$indic_ref"}}, $numero;};

};

};

# idem pour la ref de noeuds associés, mais ici dans tous les cas

if ( $ref_noeud{'N_'."$type_geom[$tyyp]".'_'."$indic_ref"} == undef)

{ my @t_noeud_ref_elem; # on sauvegarde le tableau (mais ici pour l'instant sans rien dedans)

$ref_noeud{'N_'."$type_geom[$tyyp]".'_'."$indic_ref"} = [@t_noeud_ref_elem];

};

# on rempli le tableau de num de noeud pour la ref de noeuds associé

for (my $ia =0;$ia < $nbn_par_elem[$tyyp]; $ia++)

{ my $bb = $liste_element[$j][$el][3+ $nbtag + $ia ]; # le numero de noeud

push @{$ref_noeud{'N_'."$type_geom[$tyyp]".'_'."$indic_ref"}}, $bb;

};

};

};

# --- puis les 2 tags automatique:

# 1) le numéro de l'entité géométrique élémentaires

# 2) le numéro de la partition de maillage

# on n'intervient pour les ref d'éléments que si ces éléments sont à retenir

if ($number_type[$j] != 0)

{ my $geoParti = 1; # pour choisir entre les deux

for (my $ic = $nbtag-2 ; $ic < $nbtag; $ic++,$geoParti++)

{ my $indic_ref = $liste_element[$j][$el][3+$ic]; # l'indicateur de la référence

# dans le cas où il y a un nom de référence, l'indicateur de référence est doublé

# ce qui n'est pas normal, on supprime ce doublement

if ($tag_deja_traite{"$indic_ref"} ne undef)

{next;}

$tag_deja_traite{"$indic_ref"} = '1'; # de maniere à pas avoir de undef

{ # cas ou le nom de la référence n'existe pas, on utilise directement le nombre

if ( $ref_element{'E_'."$geoParti".'_'."$indic_ref"} == undef)

# si la référence n'est pas définie on la définie

{ my @tab = ($numero);

$ref_element{'E_'."$geoParti".'_'."$indic_ref"} = [@tab]; # on sauvegarde le tableau

}

else

{ push @{$ref_element{'E_'."$geoParti".'_'."$indic_ref"}}, $numero;

};

};

};

};

};

# on suit le cas de l'écriture finale des éléments donc on n'incrémente le numéro

# d'élément que s'il sera écrit

if ($number_type[$j] != 0)

{ $numero++;};

};

};

};

# --- maintenant, pour les références de noeuds on va supprimer les doublons, car le mode

# de construction des références d'éléments entra^ine systématiquement des doublons

foreach my $refno(keys %ref_noeud)

{ # on crée une table de hachage de num de noeuds

my %les_noeuds;

# on parcours les numéros que l'on stocke dans la table (mais ici le mode de stockage

# conduit à pas de doublon

foreach my $nbnoeu ( @{$ref_noeud{$refno}})

{ $les_noeuds{$nbnoeu} = $nbnoeu;};

# maintenant on crée un tableau intermédiaire qui va contenir les noeuds définitif

my @noeuds_def;

foreach my $cle (sort keys %les_noeuds)

{ push @noeuds_def,$les_noeuds{$cle};

};

# enfin on remplace le tableau associé à la référence de noeud par celui que l'on vient de créer

$ref_noeud{$refno} = [@noeuds_def];

};