int main (int argc, char **argv) { med_err ret,lret; med_idt fid; char * fichier = NULL; char maa[MED_TAILLE_NOM+1]=""; char desc[MED_TAILLE_DESC+1]=""; char pflname[MED_TAILLE_NOM+1]="",nomlien[MED_TAILLE_NOM+1]=""; char locname[MED_TAILLE_NOM+1]=""; char * lien = NULL; char *comp= NULL, *unit= NULL; char nomcha [MED_TAILLE_NOM+1]=""; med_int mdim,ncomp,ncha,npro,nln,pflsize,*pflval,nval; med_int ngauss,nloc; int t1,t2,t3; med_type_champ typcha; med_maillage type; med_geometrie_element type_geo; med_float *refcoo, *gscoo, *wg; int i,j; if (argc != 2) { MESSAGE("Aucun nom de fichier precise, fichier test10.med utilise "); fichier = "test10.med"; } else { fichier = argv[1]; }; /* Ouverture du fichier med */ if ((fid = MEDouvrir(fichier,MED_LECTURE)) < 0){ MESSAGE("Erreur a l'ouverture du fichier : ");SSCRUTE(fichier); return -1; } ret = 0; /* infos sur le premier maillage */ if ( MEDmaaInfo(fid,1,maa,&mdim,&type,desc) < 0 ) { MESSAGE("Erreur a la lecture des informations sur le maillage : "); SSCRUTE(maa);ISCRUTE(mdim);ISCRUTE(type);SSCRUTE(desc); return -1; } printf("Maillage de nom |%s| et de dimension %d \n",maa,mdim); /* combien de champs dans le fichier */ if ((ncha = MEDnChamp(fid,0)) < 0) { MESSAGE("Impossible de lire le nombre de champs : ");ISCRUTE(ncha); return ncha; } printf("Nombre de champs : %d \n",ncha); /* lecture de tous les champs */ for (i =0;i<ncha;i++) { lret = 0; printf("\nChamp numero : %d \n",i+1); /* Lecture du nombre de composantes */ if ((ncomp = MEDnChamp(fid,i+1)) < 0) { MESSAGE("Erreur a la lecture du nombre de composantes : "); ISCRUTE(ncomp); ret = -1; continue; } /* Lecture du type du champ, des noms des composantes et du nom de l'unite*/ comp = (char*) malloc(ncomp*MED_TAILLE_PNOM+1); EXIT_IF(comp == NULL,NULL,NULL); unit = (char*) malloc(ncomp*MED_TAILLE_PNOM+1); EXIT_IF(unit == NULL,NULL,NULL); if ( MEDchampInfo(fid,i+1,nomcha,&typcha,comp,unit,ncomp) < 0 ) { MESSAGE("Erreur a la demande d'information sur les champs : "); ISCRUTE(i+1);SSCRUTE(nomcha);ISCRUTE(typcha);SSCRUTE(comp);SSCRUTE(unit); ISCRUTE(ncomp); ret = -1; continue; } printf("Nom du champ : |%s| de type %d\n",nomcha,typcha); printf("Nom des composantes : |%s|\n",comp); printf("Unites des composantes : |%s| \n",unit); free(comp); free(unit); lret = getFieldsOn(fid, nomcha, typcha, ncomp, MED_NOEUD, USER_INTERLACE ); if (lret == 0) lret = getFieldsOn(fid, nomcha, typcha, ncomp, MED_MAILLE, USER_INTERLACE ); else { MESSAGE("Erreur a la lecture des champs aux noeuds "); ret = -1; continue;} if (lret == 0) lret = getFieldsOn(fid, nomcha, typcha, ncomp, MED_FACE,USER_INTERLACE); else { MESSAGE("Erreur a la lecture des champs aux mailles "); ret = -1; continue;} if (lret == 0) lret = getFieldsOn(fid, nomcha, typcha, ncomp, MED_ARETE,USER_INTERLACE); else {MESSAGE("Erreur a la lecture des champs aux faces "); ret = -1; continue;} if (lret != 0) {MESSAGE("Erreur a la lecture des champs aux aretes "); ret = -1;}; } /* Interrogation des profils */ npro = MEDnProfil(fid); printf("\nNombre de profils stockes : %i\n\n",npro); for (i=1 ; i <= npro ; i++ ) { if ( MEDprofilInfo(fid, i, pflname, &nval) < 0) { MESSAGE("Erreur a la demande d'information sur le profil n° : "); ISCRUTE(i); ret = -1;continue; } printf("\t- Profil n°%i de nom |%s| et de taille %i\n",i,pflname,nval); pflval = (med_int*) malloc(sizeof(med_int)*nval); if ( MEDprofilLire(fid, pflval, pflname) < 0) { MESSAGE("Erreur a la lecture des valeurs du profil : "); SSCRUTE(pflname); ret = -1; } else { printf("\t"); for (j=0;j<nval;j++) printf(" %i ",*(pflval+j)); printf("\n\n"); } free(pflval); } /* Interrogation des liens */ nln = MEDnLien(fid); printf("\nNombre de liens stockes : %i\n\n",nln); for (i=1 ; i <= nln ; i++ ) { if ( MEDlienInfo(fid, i, nomlien, &nval) < 0) { MESSAGE("Erreur a la demande d'information sur le lien n° : "); ISCRUTE(i); ret = -1;continue; } printf("\t- Lien n°%i de nom |%s| et de taille %i\n",i,nomlien,nval); lien = malloc((nval+1)*sizeof(char)); EXIT_IF(lien == NULL,NULL,NULL); if ( MEDlienLire(fid, lien, nomlien) < 0 ) { MESSAGE("Erreur a la lecture du lien : "); SSCRUTE(nomlien);SSCRUTE(lien); ret = -1; } else printf("\t\t|%s|\n\n",lien); free(lien); } /* Interrogation des localisations des points de GAUSS */ nloc = MEDnGauss(fid); printf("\nNombre de localisations stockees : %i\n\n",nloc); for (i=1 ; i <= nloc ; i++ ) { if ( MEDgaussInfo(fid, i, locname, &type_geo, &ngauss) < 0) { MESSAGE("Erreur a la demande d'information sur la localisation n° : "); ISCRUTE(i); ret = -1;continue; } printf("\t- Loc. n°%i de nom |%s| et nbr. de pts de GAUSS %i\n",i,locname,ngauss); t1 = (type_geo%100)*(type_geo/100); t2 = ngauss*(type_geo/100); t3 = ngauss; refcoo = (med_float *) malloc(sizeof(med_float)*t1 ); gscoo = (med_float *) malloc(sizeof(med_float)*t2 ); wg = (med_float *) malloc(sizeof(med_float)*t3 ); if ( MEDgaussLire(fid, refcoo, gscoo, wg, USER_INTERLACE, locname ) < 0) { MESSAGE("Erreur a la lecture des valeurs de la localisation : "); SSCRUTE(locname); ret = -1; } else { printf("\t Coordonnees de l'element de reference de type %i :\n\t\t",type_geo); for (j=0;j<t1;j++) printf(" %f ",*(refcoo+j)); printf("\n"); printf("\t Localisation des points de GAUSS : \n\t\t"); for (j=0;j<t2;j++) printf(" %f ",*(gscoo+j)); printf("\n"); printf("\t Poids associes aux points de GAUSS :\n\t\t"); for (j=0;j<t3;j++) printf(" %f ",*(wg+j)); printf("\n\n"); } free(refcoo); free(gscoo); free(wg); } /* fermeture du fichier */ if ( MEDfermer(fid) < 0) return -1; return ret; }
void _MEDlocalizationRd236(int dummy, ...) { med_err _ret = -1; med_idt _lzid=0, _root=0; med_int _nentity=0,_nipoint=0,_spacedimension=0; med_int _intgeotype = -1; char _path[MED_LOCALIZATION_GRP_SIZE+MED_NAME_SIZE+1]=MED_LOCALIZATION_GRP; med_filter _filter = MED_FILTER_INIT; MED_VARGS_DECL(const, med_idt , , fid ); MED_VARGS_DECL(const, char* , const , localizationname ); MED_VARGS_DECL(const, med_switch_mode , , switchmode ); MED_VARGS_DECL(, med_float* , const , elementcoordinate ); MED_VARGS_DECL(, med_float* , const , ipointcoordinate ); MED_VARGS_DECL(, med_float* , const , weight ); MED_VARGS_DECL(, med_err * ,, fret ); va_list params; va_start(params,dummy); MED_VARGS_DEF(const, med_idt , , fid ); MED_VARGS_DEF(const, char* , const , localizationname ); MED_VARGS_DEF(const, med_switch_mode , , switchmode ); MED_VARGS_DEF(, med_float* , const , elementcoordinate ); MED_VARGS_DEF(, med_float* , const , ipointcoordinate ); MED_VARGS_DEF(, med_float* , const , weight ); MED_VARGS_DEF(, med_err * ,, fret ); /* * On inhibe le gestionnaire d'erreur HDF 5 */ _MEDmodeErreurVerrouiller(); /* * MED_GAUSS_ELNO est un mot cle reserve */ if (! strcmp(localizationname,MED_GAUSS_ELNO)) { MED_ERR_(_ret,MED_ERR_INVALID,MED_ERR_LOCALIZATION,localizationname); goto ERROR; } if ( MEDgaussLire(fid, elementcoordinate, ipointcoordinate,weight, switchmode, (char *) localizationname) < 0 ) { MED_ERR_(_ret,MED_ERR_CALL,MED_ERR_API,"MEDgaussLire"); SSCRUTE(localizationname); goto ERROR; } _ret = 0; ERROR: va_end(params); *fret = _ret; return; }