      1 /*************************************************************************
      2  * COPYRIGHT (C) 1999 - 2007  EDF R&D, CEA/DEN
      3  * THIS LIBRARY IS FREE SOFTWARE; YOU CAN REDISTRIBUTE IT AND/OR MODIFY
      4  * IT UNDER THE TERMS OF THE GNU LESSER GENERAL PUBLIC LICENSE
      5  * AS PUBLISHED BY THE FREE SOFTWARE FOUNDATION;
      6  * EITHER VERSION 2.1 OF THE LICENSE, OR (AT YOUR OPTION) ANY LATER VERSION.
      7  *
      8  * THIS LIBRARY IS DISTRIBUTED IN THE HOPE THAT IT WILL BE USEFUL, BUT
      9  * WITHOUT ANY WARRANTY; WITHOUT EVEN THE IMPLIED WARRANTY OF
     10  * MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. SEE THE GNU
     11  * LESSER GENERAL PUBLIC LICENSE FOR MORE DETAILS.
     12  *
     13  * YOU SHOULD HAVE RECEIVED A COPY OF THE GNU LESSER GENERAL PUBLIC LICENSE
     14  * ALONG WITH THIS LIBRARY; IF NOT, WRITE TO THE FREE SOFTWARE FOUNDATION,
     15  * INC., 59 TEMPLE PLACE, SUITE 330, BOSTON, MA 02111-1307 USA
     16  *
     17  *************************************************************************/
     18
     19 /******************************************************************************
     20  * - Nom du fichier : test11.c
     21  *
     22  * - Description : lecture de champs de resultats MED
     23  *
     24  *****************************************************************************/
     25
     26 #include <med.h>
     27 #define MESGERR
     28 #include <med_utils.h>
     29
     30 #ifdef DEF_LECT_ECR
     31 #define MODE_ACCES MED_LECTURE_ECRITURE
     32 #elif DEF_LECT_AJOUT
     33 #define MODE_ACCES MED_LECTURE_AJOUT
     34 #else
     35 #define MODE_ACCES MED_CREATION
     36 #endif
     37
     38 #ifndef USER_INTERLACE
     39 #define USER_INTERLACE MED_FULL_INTERLACE
     40 #endif
     41
     42 #define USER_MODE MED_COMPACT
     43
     44 med_err getFieldsOn(med_idt fid, char * nomcha, med_type_champ typcha, med_int ncomp,
     45                     med_entite_maillage entite, med_mode_switch stockage);
     46
     47 int main (int argc, char **argv)
     48
     49
     50 {
     51   med_err ret,lret;
     52   med_idt fid;
     53   char * fichier = NULL;
     54   char maa[MED_TAILLE_NOM+1]="";
     55   char desc[MED_TAILLE_DESC+1]="";
     56   char pflname[MED_TAILLE_NOM+1]="",nomlien[MED_TAILLE_NOM+1]="";
     57   char locname[MED_TAILLE_NOM+1]="";
     58   char * lien = NULL;
     59   char *comp= NULL, *unit= NULL;
     60   char nomcha  [MED_TAILLE_NOM+1]="";
     61   med_int mdim,ncomp,ncha,npro,nln,pflsize,*pflval,nval;
     62   med_int ngauss,nloc;
     63   int t1,t2,t3;
     64   med_type_champ typcha;
     65   med_maillage type;
     66   med_geometrie_element type_geo;
     67   med_float *refcoo, *gscoo, *wg;
     68   int i,j;
     69
     70   if (argc != 2) {
     71     MESSAGE("Aucun nom de fichier precise, fichier test10.med utilise ");
     72     fichier = "test10.med";
     73   } else {
     74     fichier = argv[1];
     75   };
     76
     77
     78   /* Ouverture du fichier med */
     79   if ((fid = MEDouvrir(fichier,MED_LECTURE)) < 0){
     80     MESSAGE("Erreur a l'ouverture du fichier : ");SSCRUTE(fichier);
     81     return -1;
     82   }
     83
     84   ret = 0;
     85
     86   /* infos sur le premier maillage */
     87   if (  MEDmaaInfo(fid,1,maa,&mdim,&type,desc) < 0 ) {
     88     MESSAGE("Erreur a la lecture des informations sur le maillage : ");
     89     SSCRUTE(maa);ISCRUTE(mdim);ISCRUTE(type);SSCRUTE(desc);
     90     return -1;
     91   }
     92
     93   printf("Maillage de nom |%s| et de dimension %d \n",maa,mdim);
     94
     95   /* combien de champs dans le fichier */
     96   if ((ncha = MEDnChamp(fid,0)) < 0) {
     97     MESSAGE("Impossible de lire le nombre de champs : ");ISCRUTE(ncha);
     98     return ncha;
     99   }
    100
    101   printf("Nombre de champs : %d \n",ncha);
    102
    103   /* lecture de tous les champs  */
    104   for (i =0;i<ncha;i++) {
    105     lret = 0;
    106     printf("\nChamp numero : %d \n",i+1);
    107
    108     /* Lecture du nombre de composantes */
    109     if ((ncomp = MEDnChamp(fid,i+1)) < 0) {
    110       MESSAGE("Erreur a la lecture du nombre de composantes : "); ISCRUTE(ncomp);
    111       ret = -1; continue;
    112     }
    113
    114     /* Lecture du type du champ, des noms des composantes et du nom de l'unite*/
    115     comp = (char*) malloc(ncomp*MED_TAILLE_PNOM+1);
    116     EXIT_IF(comp == NULL,NULL,NULL);
    117     unit = (char*) malloc(ncomp*MED_TAILLE_PNOM+1);
    118     EXIT_IF(unit == NULL,NULL,NULL);
    119
    120     if ( MEDchampInfo(fid,i+1,nomcha,&typcha,comp,unit,ncomp) < 0 ) {
    121       MESSAGE("Erreur a la demande d'information sur les champs : ");
    122       ISCRUTE(i+1);SSCRUTE(nomcha);ISCRUTE(typcha);SSCRUTE(comp);SSCRUTE(unit);
    123       ISCRUTE(ncomp);
    124       ret = -1; continue;
    125     }
    126
    127     printf("Nom du champ : |%s| de type %d\n",nomcha,typcha);
    128     printf("Nom des composantes : |%s|\n",comp);
    129     printf("Unites des composantes : |%s| \n",unit);
    130
    131     free(comp);
    132     free(unit);
    133
    134
    135     lret = getFieldsOn(fid, nomcha, typcha, ncomp, MED_NOEUD, USER_INTERLACE );
    136
    137     if (lret == 0) lret = getFieldsOn(fid, nomcha, typcha, ncomp, MED_MAILLE, USER_INTERLACE );
    138     else { MESSAGE("Erreur a la lecture des champs aux noeuds "); ret = -1; continue;}
    139
    140     if (lret == 0) lret = getFieldsOn(fid, nomcha, typcha, ncomp, MED_FACE,USER_INTERLACE);
    141     else { MESSAGE("Erreur a la lecture des champs aux mailles "); ret = -1; continue;}
    142
    143     if (lret == 0) lret = getFieldsOn(fid, nomcha, typcha, ncomp, MED_ARETE,USER_INTERLACE);
    144     else {MESSAGE("Erreur a la lecture des champs aux faces "); ret = -1; continue;}
    145
    146     if (lret == 0) lret = getFieldsOn(fid, nomcha, typcha, ncomp, MED_NOEUD_MAILLE,USER_INTERLACE);
    147     else {MESSAGE("Erreur a la lecture des champs aux aretes"); ret = -1; continue;}
    148
    149     if  (lret != 0) {MESSAGE("Erreur a la lecture des champs aux noeuds des mailles "); ret = -1;};
    150   }
    151
    152
    153   /* Interrogation des profils */
    154   npro = MEDnProfil(fid);
    155
    156   printf("\nNombre de profils stockes : %i\n\n",npro);
    157   for (i=1 ; i <= npro ; i++ ) {
    158     if ( MEDprofilInfo(fid, i, pflname, &nval) < 0)  {
    159       MESSAGE("Erreur a la demande d'information sur le profil n° : "); ISCRUTE(i);
    160       ret = -1;continue;
    161     }
    162     printf("\t- Profil n°%i de nom |%s| et de taille %i\n",i,pflname,nval);
    163     pflval = (med_int*) malloc(sizeof(med_int)*nval);
    164     if ( MEDprofilLire(fid, pflval, pflname) < 0) {
    165       MESSAGE("Erreur a la lecture des valeurs du profil : ");
    166       SSCRUTE(pflname);
    167       ret = -1;
    168     } else {
    169       printf("\t");
    170       for (j=0;j<nval;j++) printf(" %i ",*(pflval+j));
    171       printf("\n\n");
    172     }
    173     free(pflval);
    174   }
    175
    176   /* Interrogation des liens */
    177   nln = MEDnLien(fid);
    178
    179   printf("\nNombre de liens stockes : %i\n\n",nln);
    180   for (i=1 ; i <= nln ; i++ ) {
    181     if ( MEDlienInfo(fid, i, nomlien, &nval) < 0)  {
    182       MESSAGE("Erreur a la demande d'information sur le lien n° : "); ISCRUTE(i);
    183       ret = -1;continue;
    184     }
    185     printf("\t- Lien n°%i de nom |%s| et de taille %i\n",i,nomlien,nval);
    186
    187     lien = malloc((nval+1)*sizeof(char));
    188     EXIT_IF(lien == NULL,NULL,NULL);
    189
    190     if ( MEDlienLire(fid, lien, nomlien) < 0 )  {
    191       MESSAGE("Erreur a la lecture du lien : ");
    192       SSCRUTE(nomlien);SSCRUTE(lien);
    193       ret = -1;
    194     } else {
    195       lien[nval] = '\0';
    196       printf("\t\t|%s|\n\n",lien);
    197     }
    198     free(lien);
    199   }
    200
    201   /* Interrogation des localisations des points de GAUSS */
    202   nloc = MEDnGauss(fid);
    203
    204   printf("\nNombre de localisations stockees : %i\n\n",nloc);
    205   for (i=1 ; i <= nloc ; i++ ) {
    206     if ( MEDgaussInfo(fid, i, locname, &type_geo, &ngauss) < 0)  {
    207       MESSAGE("Erreur a la demande d'information sur la localisation n° : "); ISCRUTE(i);
    208       ret = -1;continue;
    209     }
    210     printf("\t- Loc. n°%i de nom |%s| et nbr. de pts de GAUSS %i\n",i,locname,ngauss);
    211     t1 = (type_geo%100)*(type_geo/100);
    212     t2 = ngauss*(type_geo/100);
    213     t3 = ngauss;
    214     refcoo = (med_float *) malloc(sizeof(med_float)*t1 );
    215     gscoo  = (med_float *) malloc(sizeof(med_float)*t2 );
    216     wg     = (med_float *) malloc(sizeof(med_float)*t3 );
    217
    218     if ( MEDgaussLire(fid, refcoo, gscoo, wg, USER_INTERLACE, locname ) < 0) {
    219       MESSAGE("Erreur a la lecture des valeurs de la localisation : ");
    220       SSCRUTE(locname);
    221       ret = -1;
    222     } else {
    223       printf("\t  Coordonnees de l'element de reference de type %i :\n\t\t",type_geo);
    224       for (j=0;j<t1;j++) printf(" %f ",*(refcoo+j));
    225       printf("\n");
    226       printf("\t  Localisation des points de GAUSS : \n\t\t");
    227       for (j=0;j<t2;j++) printf(" %f ",*(gscoo+j));
    228       printf("\n");
    229       printf("\t  Poids associes aux points de GAUSS :\n\t\t");
    230       for (j=0;j<t3;j++) printf(" %f ",*(wg+j));
    231       printf("\n\n");
    232     }
    233     free(refcoo);
    234     free(gscoo);
    235     free(wg);
    236   }
    237
    238
    239   /* fermeture du fichier */
    240   if ( MEDfermer(fid) < 0) return -1;
    241
    242   return ret;
    243 }
    244
    245 med_err getFieldsOn(med_idt fid, char * nomcha, med_type_champ typcha, med_int ncomp,
    246                     med_entite_maillage entite, med_mode_switch stockage) {
    247
    248   int j,k,l,m,n,nb_geo;
    249   med_int nbpdtnor=0,pflsize,*pflval,ngauss=0,ngroup,*vale=NULL,nval;
    250   med_int numdt=0,numo=0,lnsize,nbrefmaa;
    251   med_float *valr=NULL,dt=0.0;
    252   med_err ret=0;
    253   med_booleen local;
    254   char pflname [MED_TAILLE_NOM+1]="";
    255   char locname [MED_TAILLE_NOM+1]="";
    256   char * lien = NULL;
    257   char maa_ass [MED_TAILLE_NOM+1]="";
    258   char dt_unit [MED_TAILLE_PNOM+1]="";
    259
    260
    261   med_geometrie_element * type_geo;
    262   med_geometrie_element typ_noeud[1] = { MED_NONE };
    263   med_geometrie_element typmai[MED_NBR_GEOMETRIE_MAILLE+2] = {MED_POINT1, MED_SEG2, MED_SEG3, MED_TRIA3,
    264                                                               MED_QUAD4, MED_TRIA6,MED_QUAD8, MED_TETRA4,
    265                                                               MED_PYRA5, MED_PENTA6, MED_HEXA8, MED_TETRA10,
    266                                                               MED_PYRA13, MED_PENTA15, MED_HEXA20,
    267                                                               MED_POLYGONE, MED_POLYEDRE};
    268   med_geometrie_element typfac[MED_NBR_GEOMETRIE_FACE+1] = {MED_TRIA3,MED_TRIA6,
    269                                                             MED_QUAD4,MED_QUAD8,
    270                                                             MED_POLYGONE};
    271   med_geometrie_element typare[MED_NBR_GEOMETRIE_ARETE] = {MED_SEG2,MED_SEG3};
    272
    273   char ** AFF;
    274
    275   switch (entite) {
    276   case MED_NOEUD :
    277     type_geo = typ_noeud;
    278     nb_geo   = 1;
    279     AFF      = MED_GEOMETRIE_NOEUD_AFF;
    280     break;
    281   case  MED_MAILLE :
    282   case  MED_NOEUD_MAILLE :
    283     type_geo = typmai;
    284     nb_geo   = MED_NBR_GEOMETRIE_MAILLE+2;
    285     AFF      = MED_GEOMETRIE_MAILLE_AFF;
    286     break;
    287   case  MED_FACE :
    288     type_geo = typfac;
    289     nb_geo   = MED_NBR_GEOMETRIE_FACE+1;
    290     AFF      = MED_GEOMETRIE_FACE_AFF;
    291     break;
    292   case  MED_ARETE :
    293     type_geo = typare;
    294     nb_geo   = MED_NBR_GEOMETRIE_ARETE;
    295     AFF      = MED_GEOMETRIE_ARETE_AFF;
    296     break;
    297   }
    298
    299
    300   for (k=0;k<nb_geo;k++) {
    301
    302     /* Combien de (PDT,NOR) a lire */
    303     nbpdtnor = MEDnPasdetemps(fid,nomcha,entite,type_geo[k]);
    304     if (nbpdtnor < 1 ) continue;
    305
    306     for (j=0;j<nbpdtnor;j++) {
    307
    308       if ( MEDpasdetempsInfo(fid,nomcha,entite,type_geo[k],
    309                              j+1, &ngauss, &numdt, &numo, dt_unit,
    310                              &dt, maa_ass, &local, &nbrefmaa) <0) {
    311         MESSAGE("Erreur a la demande d'information sur (pdt,nor) : ");
    312         ISCRUTE(numdt); ISCRUTE(numo);
    313         ret = -1; continue;
    314       };
    315
    316       printf("\n  +Pas de Temps n.%i (%f) [%s], n. d'ordre %i, avec %i pts de gauss sur le maillage par defaut.\n",numdt,dt,dt_unit,numo,ngauss);
    317
    318       printf("\tLe maillage par defaut est : |%s|, sur un total de : %i maillages associes\n",
    319              maa_ass, nbrefmaa);
    320
    321       /* Le maillage reference est-il porte par un autre fichier */
    322       if ( !local ) {
    323
    324         if ( (lnsize=MEDnValLien(fid,maa_ass) ) < 0 )  {
    325           MESSAGE("Erreur a la lecture de la taille du lien : ");
    326           SSCRUTE(maa_ass);
    327           ret = -1;
    328         } else {
    329
    330           lien = malloc(lnsize*sizeof(char));
    331           EXIT_IF(lien == NULL,NULL,NULL);
    332
    333           if ( MEDlienLire(fid, lien, maa_ass) < 0 )  {
    334             MESSAGE("Erreur a la lecture du lien : ");
    335             SSCRUTE(maa_ass);SSCRUTE(lien);
    336             ret = -1;
    337           } else {
    338             printf("\tLe maillage |%s| est porte par un fichier distant |%s|\n",maa_ass,lien);
    339           }
    340           free(lien);
    341         }
    342       }
    343
    344       /* Combien de maillages lies aux (nomcha,ent,geo,numdt,numo)  */
    345       /* Notons que cette information est egalement disponible a partir de MEDpasdetempsInfo */
    346       if ( (nbrefmaa = MEDnChampRef(fid,nomcha,entite,type_geo[k],numdt,numo) ) < 0 ) {
    347         MESSAGE("Erreur a la demande du nombre de maillages references par le champ : ");
    348         SSCRUTE(nomcha);
    349         ISCRUTE(entite);ISCRUTE(type_geo[k]);
    350         SSCRUTE(MED_ENTITE_MAILLAGE_AFF[(int)entite]);SSCRUTE(AFF[k]);
    351         ISCRUTE(numdt); ISCRUTE(numo);
    352         ret = -1; continue;
    353       };
    354
    355       for (l=0;l<nbrefmaa;l++) {
    356
    357         if ( MEDchampRefInfo(fid,nomcha,entite,type_geo[k],
    358                              l+1,numdt, numo, maa_ass, &local, &ngauss) <0 ) {
    359           MESSAGE("Erreur a la demande d'information sur le maillage utilise par le champ n° : ");
    360           SSCRUTE(nomcha);ISCRUTE(entite);ISCRUTE(type_geo[k]);
    361           ISCRUTE(l+1);ISCRUTE(numdt); ISCRUTE(numo);SSCRUTE(maa_ass);
    362           ret = -1; continue;
    363         };
    364
    365         /* Prend en compte le nbre de pt de gauss automatiquement */
    366         if ((nval = MEDnVal(fid,nomcha,entite,type_geo[k],numdt,numo,maa_ass,USER_MODE)) <= 0)   {
    367           MESSAGE("Erreur a la lecture du nombre de valeurs du champ : ");
    368           SSCRUTE(nomcha);ISCRUTE(entite);ISCRUTE(type_geo[k]);
    369           ISCRUTE(numdt);ISCRUTE(numo);SSCRUTE(maa_ass);ISCRUTE(USER_MODE);
    370           ret = -1; continue;
    371         };
    372
    373
    374         printf("\t- Il y a %d valeurs en mode %i. Chaque entite %s\
    375  de type geometrique %s associes au maillage |%s| a %i pts de gauss \n",
    376                nval,USER_MODE,MED_ENTITE_MAILLAGE_AFF[(int)entite],AFF[k],maa_ass,ngauss);
    377
    378         /* Le maillage reference est-il porte par un autre fichier */
    379         if ( !local ) {
    380
    381           if ( (lnsize=MEDnValLien(fid,maa_ass) ) < 0 )  {
    382             MESSAGE("Erreur a la lecture de la taille du lien : ");
    383             SSCRUTE(maa_ass);
    384             ret = -1;
    385           } else {
    386
    387             lien = malloc(lnsize*sizeof(char) + 1);
    388             EXIT_IF(lien == NULL,NULL,NULL);
    389
    390             if ( MEDlienLire(fid, lien, maa_ass) < 0 )  {
    391               MESSAGE("Erreur a la lecture du lien : ");
    392               SSCRUTE(maa_ass);SSCRUTE(lien);
    393               ret = -1;
    394             } else {
    395               lien[lnsize] = '\0';
    396               printf("\tLe maillage |%s| est porte par un fichier distant |%s|\n",maa_ass,lien);
    397             }
    398             free(lien);
    399           }
    400         }
    401
    402         /*Lecture des valeurs du champ */
    403         if (typcha == MED_FLOAT64) {
    404
    405           valr = (med_float*) calloc(ncomp*nval,sizeof(med_float));
    406           EXIT_IF(valr == NULL,NULL,NULL);
    407           if ( MEDchampLire(fid,maa_ass,nomcha,(unsigned char*)valr,stockage,MED_ALL,locname,
    408                             pflname,USER_MODE,entite,type_geo[k],numdt,numo) < 0 ) {
    409             MESSAGE("Erreur a la lecture du nombre de valeurs du champ : ");
    410             SSCRUTE(nomcha);ISCRUTE(entite);ISCRUTE(type_geo[k]);
    411             ISCRUTE(numdt);ISCRUTE(numo);SSCRUTE(maa_ass);
    412             ret = -1;
    413           };
    414
    415         } else {
    416
    417           vale = (med_int*) calloc(ncomp*nval,sizeof(med_int));
    418           EXIT_IF(vale == NULL,NULL,NULL);
    419           if ( MEDchampLire(fid,maa_ass,nomcha,(unsigned char*)vale,stockage,MED_ALL,locname,
    420                             pflname,USER_MODE,entite,type_geo[k],numdt,numo) < 0 ) {
    421             MESSAGE("Erreur a la lecture des valeurs du champ : ");
    422             SSCRUTE(nomcha);ISCRUTE(entite);ISCRUTE(type_geo[k]);
    423             ISCRUTE(numdt);ISCRUTE(numo);SSCRUTE(maa_ass);
    424             ret = -1;
    425           };
    426
    427         }
    428
    429         if ( ngauss > 1 )
    430           printf("\t- Modèle de localisation des points de Gauss de nom |%s|\n",locname);
    431
    432         if (entite == MED_NOEUD_MAILLE)
    433           ngroup = type_geo[k] % 100;
    434         else
    435           ngroup = ngauss;
    436
    437         switch (stockage) {
    438
    439         case MED_FULL_INTERLACE :
    440           printf("\t- Valeurs :\n\t");
    441           for (m=0;m<nval/ngroup;m++) {
    442             printf("|");
    443             for (n=0;n<ngroup*ncomp;n++)
    444               if (typcha == MED_FLOAT64)
    445                 printf(" %f ",*(valr+(m*ngroup*ncomp)+n));
    446               else
    447                 printf(" %d ",*(vale+(m*ngroup*ncomp)+n));
    448
    449           }
    450           break;
    451
    452           /*Affichage en fonction du profil à traiter*/
    453         case MED_NO_INTERLACE :
    454           printf("\t- Valeurs :\n\t");
    455           for (m=0;m<ncomp;m++) {
    456             printf("|");
    457             for (n=0;n<nval;n++)
    458               if (typcha == MED_FLOAT64)
    459                 printf(" %f ",*(valr+(m*nval)+n));
    460               else
    461                 printf(" %d ",*(vale+(m*nval)+n));
    462           }
    463           break;
    464         }
    465
    466         printf("|\n");
    467         if (typcha == MED_FLOAT64) {
    468           if ( valr ) {free(valr);valr = NULL;}}
    469         else
    470           if (vale) { free(vale);vale = NULL; }
    471
    472         /*Lecture du profil associe */
    473         if (strcmp(pflname,MED_NOPFL) == 0 )
    474           printf("\t- Profil : MED_NOPFL\n");
    475         else {
    476
    477           if ( (pflsize = MEDnValProfil(fid,pflname)) <0 )  {
    478             MESSAGE("Erreur a la lecture du nombre de valeurs du profil : ");
    479             SSCRUTE(pflname);
    480             ret = -1; continue;
    481           }
    482
    483           printf("\t- Profil : |%s| de taille %i\n",pflname,pflsize);
    484
    485           pflval = (med_int*) malloc(sizeof(med_int)*pflsize);
    486           EXIT_IF(pflval == NULL,NULL,NULL);
    487           if ( MEDprofilLire(fid,pflval,pflname) <0) {
    488             MESSAGE("Erreur a la lecture des valeurs du profil : ");
    489             SSCRUTE(pflname);
    490             ret = -1;
    491           }
    492           printf("\t");
    493           for (m=0;m<pflsize;m++) printf(" %i ",*(pflval+m));
    494           printf("\n");
    495           free(pflval);
    496
    497         }
    498
    499       }
    500     }
    501   } /* fin for sur les mailles*/
    502
    503   return ret;
    504 }
    505