1 // Copyright (C) 2007-2008 CEA/DEN, EDF R&D, OPEN CASCADE
3 // Copyright (C) 2003-2007 OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
4 // CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
6 // This library is free software; you can redistribute it and/or
7 // modify it under the terms of the GNU Lesser General Public
8 // License as published by the Free Software Foundation; either
9 // version 2.1 of the License.
11 // This library is distributed in the hope that it will be useful,
12 // but WITHOUT ANY WARRANTY; without even the implied warranty of
13 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
14 // Lesser General Public License for more details.
16 // You should have received a copy of the GNU Lesser General Public
17 // License along with this library; if not, write to the Free Software
18 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
20 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
22 #include "MEDMEM_MEDMEMchampLire.hxx"
24 * En attendant une correction de la gestion du mode d'accès au fichier dans MEDfichier
25 * on intègre la correction ici.
33 fprintf(stderr, "%s [%d] : " , __FILE__ , __LINE__ ) ; \
37 # define ISCRUTE_MED(entier) { \
39 fprintf(stderr,"%s = %d\n",#entier,entier) ; \
43 # define SSCRUTE_MED(chaine) { \
45 fprintf(stderr,"%s = \"%s\"\n",#chaine,chaine) ; \
48 # define MESSAGE_MED(chaine) { \
50 fprintf(stderr,"%s\n",chaine) ; \
54 extern void _MEDmodeErreurVerrouiller(void);
57 #include <med_outils.h>
63 * - Nom de la fonction : MEDchampLire
64 * - Description : Lecture d'un Champ Résultat
66 * - fid (IN) : ID du fichier HDF courant
67 * - maa (IN) : le nom du maillage sur lequel s'applique le champ
68 * - cha (IN) : le nom du champ
69 * - val (OUT) : valeurs du champ à lire
70 * - interlace(IN) : entrelacement voulu en mémoire {MED_FULL_INTERLACE,MED_NO_INTERLACE}
71 * - numco (IN) : n° de la composante à lire (MED_ALL si toutes)
72 * - profil (OUT) : nom du profil utilisé (MED_NOPFL si inutilisé)
73 * - pflmod (IN) : Indique comment lire les informations en mémoire { MED_COMPACT, MED_GLOBALE }.
74 * - type_ent (IN) : entité concerné par le champ {MED_NOEUD,MED_ARETE,MED_FACE,MED_MAILLE}
75 * - type_geo (IN) : type géométrique de l'entité concerné {MED_POINT,MED_SEG2 ......}
76 * - numdt (IN) : n° du pas de temps (MED_NOPDT si aucun)
77 * - numo (IN) : n° d'ordre utilisé MED_NONOR si inutile
78 * - Resultat : 0 en cas de succes, -1 sinon
81 /* REM : La taille de val allouée par l'utilisateur doit prendre en compte le nbre de points de gauss et le nbre de composantes*/
84 MEDMEMchampLire(med_idt fid,char *maa, char *cha, unsigned char *val,med_mode_switch interlace,med_int numco,
85 char * locname, char *profil, med_mode_profil pflmod,
86 med_entite_maillage type_ent, med_geometrie_element type_geo,
87 med_int numdt, med_int numo)
91 med_idt gid=0, datagroup1=0, datagroup2=0,datagroup3=0;
92 med_int ncomp=0, chtype=0, ngauss=0, i=0, pfluse=0;
93 char nomdatagroup1[2*MED_TAILLE_NOM_ENTITE+2]="",nomdatagroup2[2*MED_MAX_PARA+1]="";
94 char tmp1[MED_TAILLE_NOM_ENTITE+1]="", pfltmp[MED_TAILLE_NOM+1]="";
95 char chemin[MED_TAILLE_CHA+MED_TAILLE_NOM+1]="";
101 * On inhibe le gestionnaire d'erreur HDF 5
103 _MEDmodeErreurVerrouiller();
106 * Si le Data Group cha n'existe pas => erreur
108 strcpy(chemin,MED_CHA);
110 if ((gid = _MEDdatagroupOuvrir(fid,chemin)) < 0)
113 /* Lecture du nbre de composantes du champ */
115 if (_MEDattrEntierLire(gid,MED_NOM_NCO,&ncomp) < 0)
119 * Si le Data Group de niveau 1 <type_ent>[.<type_geo>] n'existe pas => erreur
121 if (_MEDnomEntite(nomdatagroup1,type_ent) < 0)
123 if ((type_ent != MED_NOEUD))
125 if (_MEDnomGeometrie(tmp1,type_geo) < 0)
127 strcat(nomdatagroup1,".");
128 strcat(nomdatagroup1,tmp1);
131 if ( (datagroup1 = _MEDdatagroupOuvrir(gid,nomdatagroup1)) < 0 )
135 * Si le Data Group de niveau 2 <numdt>.<numoo> n'existe pas => erreur
137 sprintf(nomdatagroup2,"%*li%*li",MED_MAX_PARA,(long ) numdt,MED_MAX_PARA,(long ) numo);
140 if ( (datagroup2 = _MEDdatagroupOuvrir(datagroup1,nomdatagroup2)) < 0)
144 * Ouvre le datagroup de niveau 3 <maa>
147 if ( ! strcmp(maa,MED_NOREF) )
148 if (_MEDattrStringLire(datagroup2,MED_NOM_MAI,MED_TAILLE_NOM,maa) < 0)
152 if ( (datagroup3 = _MEDdatagroupOuvrir(datagroup2,maa)) < 0 )
155 /* Gestion des profils*/
161 if (_MEDattrStringLire(datagroup3,MED_NOM_PFL,MED_TAILLE_NOM,pfltmp) < 0)
164 if ( pfluse = (strcmp(pfltmp,MED_NOPFLi) && strcmp(pfltmp,"")) ) /* le test MED_NOPFLi pour des raisons de compatibilité */
166 strcpy(profil,pfltmp);
167 if ( (i = MEDnValProfil(fid,profil)) < 0 )
172 pfltabtmp = (med_int *) malloc (sizeof(med_int)*psize);
173 pfltab = (med_ssize *) malloc (sizeof(med_ssize)*psize);
174 if (MEDprofilLire(fid,pfltabtmp,profil) < 0)
176 for (i=0;i<psize;i++)
177 pfltab[i] = (med_ssize) pfltabtmp[i];
182 strcpy(profil,MED_NOPFL);
186 /* Lire le nbre des points de GAUSS*/
187 if (_MEDattrEntierLire(datagroup3,MED_NOM_NGA,&ngauss) < 0) {
188 MESSAGE_MED("Erreur à la lecture de l'attribut MED_NOM_NGA : ");
189 ISCRUTE_MED(ngauss);goto ERROR;
192 /* Lire l'identificateur de localisation des points de GAUSS*/
193 if ( _MEDattrStringLire(datagroup3,MED_NOM_GAU,MED_TAILLE_NOM,locname) < 0) {
194 MESSAGE_MED("Erreur à la lecture de l'attribut MED_NOM_GAU : ");
195 SSCRUTE_MED(locname); goto ERROR;
202 if (_MEDattrEntierLire(gid,MED_NOM_TYP,&chtype) < 0)
208 if ( _MEDdatasetNumLire(datagroup3,MED_NOM_CO,MED_FLOAT64,
209 interlace,ncomp,numco,
210 psize,pflmod,pfltab,ngauss,val)< 0)
215 #if defined(F77INT64)
216 if ( _MEDdatasetNumLire(datagroup3,MED_NOM_CO,MED_INT64,
217 interlace,ncomp,numco,
218 psize,pflmod,pfltab,ngauss,val)< 0)
221 if ( _MEDdatasetNumLire(datagroup3,MED_NOM_CO,MED_INT32,
222 interlace,ncomp,numco,
223 psize,pflmod,pfltab,ngauss,val)< 0)
229 #if defined(F77INT64)
230 if ( _MEDdatasetNumLire(datagroup3,MED_NOM_CO,MED_INT64,
231 interlace,ncomp,numco,
232 psize,pflmod,pfltab,ngauss,val)< 0)
251 if ( pfluse ) { free(pfltab); free(pfltabtmp);}
253 if (datagroup3>0) if (_MEDdatagroupFermer(datagroup3) < 0) {
254 MESSAGE_MED("Impossible de fermer le datagroup : ");
255 ISCRUTE_MED(datagroup3); ret = -1;
258 if (datagroup2>0) if (_MEDdatagroupFermer(datagroup2) < 0) {
259 MESSAGE_MED("Impossible de fermer le datagroup : ");
260 ISCRUTE_MED(datagroup2); ret = -1;
263 if (datagroup1>0) if (_MEDdatagroupFermer(datagroup1) < 0) {
264 MESSAGE_MED("Impossible de fermer le datagroup : ");
265 ISCRUTE_MED(datagroup1); ret = -1;
268 if (gid>0) if (_MEDdatagroupFermer(gid) < 0) {
269 MESSAGE_MED("Impossible de fermer le datagroup : ");
270 ISCRUTE_MED(gid); ret = -1;