test5.c

Aller à la documentation de ce fichier.
00001 /*  This file is part of MED.
00002  *
00003  *  COPYRIGHT (C) 1999 - 2011  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 

Généré le Mon May 16 17:10:24 2011 pour MED fichier par  doxygen 1.6.1