# test5.c

Aller à la documentation de ce fichier.

    00001 /*  This file is part of MED.
    00002  *
    00003  *  COPYRIGHT (C) 1999 - 2012  EDF R&D, CEA/DEN
    00004  *  MED is free software: you can redistribute it and/or modify
    00005  *  it under the terms of the GNU Lesser General Public License as published by
    00006  *  the Free Software Foundation, either version 3 of the License, or
    00007  *  (at your option) any later version.
    00008  *
    00009  *  MED is distributed in the hope that it will be useful,
    00010  *  but WITHOUT ANY WARRANTY; without even the implied warranty of
    00011  *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    00012  *  GNU Lesser General Public License for more details.
    00013  *
    00014  *  You should have received a copy of the GNU Lesser General Public License
    00015  *  along with MED.  If not, see <http://www.gnu.org/licenses/>.
    00016  */
    00017
    00018
    00019 /******************************************************************************
    00020  * - Nom du fichier : test5.c
    00021  *
    00022  * - Description : lecture des noeuds d'un maillage MED.
    00023  *
    00024  *****************************************************************************/
    00025
    00026 #include <med.h>
    00027 #define MESGERR 1
    00028 #include "med_utils.h"
    00029 #include <string.h>
    00030
    00031 #ifdef DEF_LECT_ECR
    00032 #define MODE_ACCES MED_ACC_RDWR
    00033 #elif DEF_LECT_AJOUT
    00034 #define MODE_ACCES MED_ACC_RDEXT
    00035 #else
    00036 #define MODE_ACCES MED_ACC_CREAT
    00037 #endif
    00038
    00039 int main (int argc, char **argv)
    00040
    00041
    00042 {
    00043   med_err ret = 0;
    00044   med_idt fid=0;
    00045   /* la dimension du maillage */
    00046   med_int mdim=0,sdim=0;
    00047   /* nom du maillage de longueur maxi MED_NAME_SIZE */
    00048   char maa[MED_NAME_SIZE+1]="";
    00049   /* le nombre de noeuds */
    00050   med_int nnoe = 0;
    00051   /* table des coordonnees */
    00052   med_float *coo1=NULL,*coo2=NULL;
    00053   /* tables des noms et des unites des coordonnees
    00054      flt : (dimension*MED_SNAME_SIZE+1) */
    00055   char nomcoo[2*MED_SNAME_SIZE+1];
    00056   char unicoo[2*MED_SNAME_SIZE+1];
    00057   /* tables des noms, numeros, numeros de familles des noeuds
    00058      autant d'elements que de noeuds - les noms ont pout longueur
    00059      MED_SNAME_SIZE */
    00060   char *nomnoe=NULL;
    00061   med_int *numnoe=NULL, *nufano=NULL;
    00062   med_axis_type rep;
    00063   med_bool inonoe=MED_FALSE,inunoe=MED_FALSE,chgt=MED_FALSE,trsf=MED_FALSE;
    00064   char str[MED_SNAME_SIZE+1];
    00065   med_int flt[2] = { 2, 3 };
    00066   char desc[MED_COMMENT_SIZE+1]="";
    00067   char dtunit[MED_SNAME_SIZE+1]="";
    00068   med_mesh_type type;
    00069   med_sorting_type sort;
    00070   med_int nstep=0,i=0;
    00071   med_filter filter=MED_FILTER_INIT;
    00072   med_int isolatednodes=0;
    00073   med_int verticesnodes=0;
    00074   med_int cellmaxnodes=0;
    00075
    00076   /* Ouverture du fichier "test4.med" en lecture seule */
    00077   fid = MEDfileOpen("test4.med",MED_ACC_RDONLY);
    00078   if (fid < 0) {
    00079     MESSAGE("Erreur a l'ouverture du fichier test4.med");
    00080     return -1;
    00081   }
    00082
    00083   if ((sdim=MEDmeshnAxis(fid, 1)) <0) {
    00084     MESSAGE("Erreur a la lecture de la dimension de l'espace du maillage :");
    00085     SSCRUTE(maa);
    00086     ret = -1;
    00087   }
    00088
    00089   /* Lecture des infos concernant le premier maillage */
    00090   if ( MEDmeshInfo( fid, 1,  maa,  &sdim, &mdim, &type, desc, dtunit, &sort,
    00091                     &nstep,  &rep, nomcoo,unicoo) < 0 ) {
    00092     MESSAGE("Erreur a la lecture des informations sur le maillage : ");SSCRUTE(maa);
    00093     return -1;
    00094   } else {
    00095     printf("Maillage de nom : |%s| , de dimension : "IFORMAT" , et de type %d\n",maa,mdim,type);
    00096     printf("\t -Dimension de l'espace : "IFORMAT"\n",sdim);
    00097     printf("\t -Description du maillage : |%s|\n",desc);
    00098     printf("\t -Noms des axes : |%s|\n",nomcoo);
    00099     printf("\t -Unités des axes : |%s|\n",unicoo);
    00100     printf("\t -Type de repère : %d\n",rep);
    00101     printf("\t -Nombre d'étape de calcul : "IFORMAT"\n",nstep);
    00102     printf("\t -Unité des dates : |%s|\n",dtunit);
    00103   }
    00104   /* Ecriture des attributs des noeuds du maillage  */
    00105   if (MEDmeshAttributeRd( fid, maa, &isolatednodes, &verticesnodes, &cellmaxnodes) < 0 ) {
    00106     MESSAGE("Aucune définition des attributs des noeuds du maillage");
    00107   } else {
    00108     printf("\t -Nombre de noeuds isolés             : "IFORMAT"\n",isolatednodes);
    00109     printf("\t -Nombre de noeuds sommets            : "IFORMAT"\n",verticesnodes);
    00110     printf("\t -Nombre maximun de noeuds par maille : "IFORMAT"\n",cellmaxnodes);
    00111   }
    00112
    00113
    00114   /* Combien de noeuds a lire ? */
    00115   nnoe = MEDmeshnEntity(fid,maa,MED_NO_DT,MED_NO_IT,
    00116                         MED_NODE,MED_NONE,MED_COORDINATE,MED_NO_CMODE,
    00117                         &chgt,&trsf);
    00118   if (nnoe < 0) {
    00119     MESSAGE("Erreur a la lecture du nombre de noeuds dans : ");
    00120     ret = -1;
    00121   } else
    00122     printf("Nombre de noeuds : "IFORMAT" \n",nnoe);
    00123
    00124   /* Allocations memoires */
    00125   if (nnoe > 0) {
    00126     /* table des coordonnees
    00127        flt : (dimension * nombre de noeuds ) */
    00128     coo1 = (med_float*) calloc(nnoe*sdim,sizeof(med_float));
    00129     coo2 = (med_float*) calloc(nnoe*sdim,sizeof(med_float));
    00130     /* table des des numeros, des numeros de familles des noeuds
    00131        flt : (nombre de noeuds) */
    00132     numnoe = (med_int*) malloc(sizeof(med_int)*nnoe);
    00133     nufano = (med_int*) malloc(sizeof(med_int)*nnoe);
    00134     /* table des noms des noeuds
    00135        flt : (nnoe*MED_SNAME_SIZE+1) */
    00136     nomnoe = (char*) malloc(MED_SNAME_SIZE*nnoe+1);
    00137   }
    00138
    00139   if ( MEDfilterEntityCr( fid, nnoe, 1, sdim, 2,
    00140                           MED_FULL_INTERLACE, MED_GLOBAL_PFLMODE,
    00141                           MED_NO_PROFILE, MED_UNDEF_SIZE,
    00142                           NULL, &filter) < 0 ) {
    00143     MESSAGE("Erreur à la création du filtre 1.");
    00144   }
    00145
    00146   /* Lecture des composantes n°2 des coordonnees des noeuds */
    00147   if (nnoe > 0) {
    00148     if ( MEDmeshNodeCoordinateAdvancedRd(fid, maa, MED_NO_DT, MED_NO_IT,
    00149                                          &filter, coo1) < 0 ) {
    00150       MESSAGE("Erreur a la lecture des coordonnees des noeuds");
    00151       ret = -1;
    00152     } else {
    00153       printf("Valeur de coo1 : ");
    00154       for (i=0;i<nnoe*sdim;i++)
    00155         printf("%4.2f ",coo1[i]);
    00156       printf("\n");
    00157     }
    00158   }
    00159
    00160   MEDfilterClose(&filter);
    00161   if ( MEDfilterEntityCr( fid, nnoe, 1, sdim, 1,
    00162                           MED_FULL_INTERLACE, MED_GLOBAL_PFLMODE,
    00163                           MED_NO_PROFILE, MED_UNDEF_SIZE,
    00164                           NULL, &filter) < 0 ) {
    00165     MESSAGE("Erreur à la création du filtre 2.");
    00166   }
    00167
    00168
    00169   /* Lecture des composantes n°1 des coordonnees des noeuds */
    00170   if (nnoe > 0) {
    00171     if ( MEDmeshNodeCoordinateAdvancedRd(fid, maa, MED_NO_DT, MED_NO_IT,
    00172                                          &filter, coo1) < 0 ) {
    00173       MESSAGE("Erreur a la lecture des coordonnees des noeuds");
    00174       ret = -1;
    00175     } else {
    00176       printf("Valeur de coo1 : ");
    00177       for (i=0;i<nnoe*sdim;i++)
    00178         printf("%4.2f ",coo1[i]);
    00179       printf("\n");
    00180     }
    00181   }
    00182
    00183   MEDfilterClose(&filter);
    00184   if ( MEDfilterEntityCr( fid, nnoe, 1, sdim, 1,
    00185                           MED_FULL_INTERLACE, MED_GLOBAL_PFLMODE,
    00186                           MED_NO_PROFILE, 2,
    00187                           flt, &filter) < 0 ) {
    00188     MESSAGE("Erreur à la création du filtre 3.");
    00189   }
    00190
    00191   /* Lecture des composantes n°1 des coordonnees des noeuds du filtre */
    00192   if (nnoe > 0) {
    00193     if ( MEDmeshNodeCoordinateAdvancedRd(fid, maa, MED_NO_DT, MED_NO_IT,
    00194                                          &filter, coo2) < 0 ) {
    00195       MESSAGE("Erreur a la lecture des coordonnees des noeuds");
    00196       ret = -1;
    00197     } else {
    00198       printf("Valeur de coo2 : ");
    00199       for (i=0;i<nnoe*sdim;i++)
    00200         printf("%4.2f ",coo2[i]);
    00201       printf("\n");
    00202     }
    00203   }
    00204
    00205   MEDfilterClose(&filter);
    00206   if ( MEDfilterEntityCr( fid, nnoe, 1, sdim, 2,
    00207                           MED_FULL_INTERLACE, MED_GLOBAL_PFLMODE,
    00208                           MED_NO_PROFILE, 2,
    00209                           flt, &filter) < 0 ) {
    00210     MESSAGE("Erreur à la création du filtre 4.");
    00211   }
    00212
    00213   /* Lecture des composantes n°2 des coordonnees des noeuds du filtre */
    00214   if (nnoe > 0) {
    00215     if ( MEDmeshNodeCoordinateAdvancedRd(fid, maa, MED_NO_DT, MED_NO_IT,
    00216                                          &filter, coo2) < 0 ) {
    00217       MESSAGE("Erreur a la lecture des coordonnees des noeuds");
    00218       ret = -1;
    00219     } else {
    00220       printf("Valeur de coo2 : ");
    00221       for (i=0;i<nnoe*sdim;i++) {
    00222         printf("%4.2f ",coo2[i]);
    00223         coo2[i] = 0.0;
    00224       }
    00225       printf("\n");
    00226     }
    00227   }
    00228
    00229   MEDfilterClose(&filter);
    00230   if ( MEDfilterEntityCr( fid, nnoe, 1, sdim, MED_ALL_CONSTITUENT,
    00231                           MED_FULL_INTERLACE, MED_GLOBAL_PFLMODE,
    00232                           MED_NO_PROFILE, 2,
    00233                           flt, &filter) < 0 ) {
    00234     MESSAGE("Erreur à la création du filtre 5.");
    00235   }
    00236
    00237   /* Lecture de toutes les composantes des coordonnees des noeuds du filtre */
    00238   if (nnoe > 0) {
    00239     if ( MEDmeshNodeCoordinateAdvancedRd(fid, maa, MED_NO_DT, MED_NO_IT,
    00240                                          &filter, coo2) < 0 ) {
    00241       MESSAGE("Erreur a la lecture des coordonnees des noeuds");
    00242       ret = -1;
    00243     } else {
    00244       printf("Valeur de coo2 : ");
    00245       for (i=0;i<nnoe*sdim;i++) {
    00246         printf("%4.2f ",coo2[i]);
    00247         coo2[i] = 0.0;
    00248       }
    00249       printf("\n");
    00250     }
    00251   }
    00252   MEDfilterClose(&filter);
    00253
    00254   /* Lecture des composantes des coordonnees des noeuds */
    00255   if (nnoe > 0) {
    00256     if ( MEDmeshNodeCoordinateRd(fid, maa, MED_NO_DT, MED_NO_IT,
    00257                                  MED_FULL_INTERLACE, coo2) < 0 ) {
    00258       MESSAGE("Erreur a la lecture des coordonnees des noeuds");
    00259       ret = -1;
    00260     } else {
    00261       printf("Valeur de coo2 : ");
    00262       for (i=0;i<nnoe*sdim;i++)
    00263         printf("%4.2f ",coo2[i]);
    00264       printf("\n");
    00265     }
    00266   }
    00267
    00268   /* Lecture des noms des noeuds (optionnel dans un maillage MED) */
    00269   if ((nnoe > 0)) {
    00270     if (MEDmeshEntityNameRd(fid,maa, MED_NO_DT, MED_NO_IT,
    00271                             MED_NODE,MED_NONE,nomnoe) < 0)
    00272       inonoe = MED_FALSE;
    00273     else
    00274       inonoe = MED_TRUE;
    00275   }
    00276
    00277   /* Lecture des numeros des noeuds (optionnel dans un maillage MED) */
    00278   if ((nnoe > 0)) {
    00279     if (MEDmeshEntityNumberRd(fid,maa, MED_NO_DT, MED_NO_IT,
    00280                               MED_NODE,MED_NONE,numnoe) < 0)
    00281       inunoe = MED_FALSE;
    00282     else
    00283       inunoe = MED_TRUE;
    00284   }
    00285
    00286   /* Lecture des numeros de familles des noeuds */
    00287   if ((nnoe > 0))
    00288     if (MEDmeshEntityFamilyNumberRd(fid,maa, MED_NO_DT, MED_NO_IT,
    00289                                     MED_NODE,MED_NONE,nufano ) < 0) {
    00290       MESSAGE("Erreur a la lecture des numeros de famille des noeuds");
    00291       ret = -1;
    00292     }
    00293
    00294   /* Fermeture du fichier */
    00295   if (MEDfileClose(fid) < 0){
    00296     MESSAGE("Erreur a la fermeture du fichier");
    00297     ret = -1;
    00298   }
    00299
    00300   /* Affichage des resulats */
    00301   if (ret == 0 && nnoe > 0)
    00302     {
    00303       printf("Type de repere : %d \n",rep);
    00304       printf("Nom des coordonnees : \n");
    00305       for (i=0;i<sdim;i++)
    00306         {
    00307           strncpy(str,nomcoo+i*MED_SNAME_SIZE,MED_SNAME_SIZE);
    00308           str[MED_SNAME_SIZE] = '\0';
    00309           printf("|%s| ",str);
    00310         }
    00311       printf("\nUnites des coordonnees : \n");
    00312       for (i=0;i<sdim;i++)
    00313         {
    00314           strncpy(str,unicoo+i*MED_SNAME_SIZE,MED_SNAME_SIZE);
    00315           str[MED_SNAME_SIZE] = '\0';
    00316           printf("|%s| ",str);
    00317         }
    00318       printf("\nCoordonnees des noeuds : \n");
    00319       for (i=0;i<nnoe*sdim;i++)
    00320         printf("%f ",*(coo2+i));
    00321       if (inonoe)
    00322         {
    00323           printf("\nNoms des noeuds : \n");
    00324           for (i=0;i<nnoe;i++)
    00325             {
    00326               strncpy(str,nomnoe+i*MED_SNAME_SIZE,MED_SNAME_SIZE);
    00327               str[MED_SNAME_SIZE] = '\0';
    00328               printf("|%s|",str);
    00329             }
    00330         }
    00331       if (inunoe)
    00332         {
    00333           printf("\nNumeros des noeuds : \n");
    00334           for (i=0;i<nnoe;i++)
    00335               printf(IFORMAT" ",*(numnoe+i));
    00336         }
    00337       printf("\nNumeros des familles des noeuds : \n");
    00338       for (i=0;i<nnoe;i++)
    00339         printf(IFORMAT" ",*(nufano+i));
    00340       printf("\n");
    00341     }
    00342
    00343   /* liberation memoire */
    00344   if (nnoe > 0) {
    00345     free(coo1);
    00346     free(coo2);
    00347     free(nomnoe);
    00348     free(numnoe);
    00349     free(nufano);
    00350   }
    00351
    00352   return ret;
    00353 }
    00354
    00355
    00356
    00357