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 : test15.c 21 * 22 * - Description : lecture des noeuds d'un maillage MED 23 * a l'aide des routines de niveau 2 24 * - equivalent a test5.c 25 * 26 *****************************************************************************/ 27 28 #include <med.h> 29 #define MESGERR 30 #include <med_utils.h> 31 32 #ifdef DEF_LECT_ECR 33 #define MODE_ACCES MED_LECTURE_ECRITURE 34 #elif DEF_LECT_AJOUT 35 #define MODE_ACCES MED_LECTURE_AJOUT 36 #else 37 #define MODE_ACCES MED_CREATION 38 #endif 39 40 int main (int argc, char **argv) 41 42 43 { 44 med_err ret = 0; 45 med_idt fid; 46 /* la dimension du maillage */ 47 med_int mdim; 48 /* nom du maillage de longueur maxi MED_TAILLE_NOM */ 49 char maa[MED_TAILLE_NOM+1]; 50 /* le nombre de noeuds */ 51 med_int nnoe = 0; 52 /* table des coordonnees */ 53 med_float *coo; 54 /* tables des noms et des unites des coordonnees 55 profil : (dimension*MED_TAILLE_PNOM+1) */ 56 char nomcoo[3*MED_TAILLE_PNOM+1]; 57 char unicoo[3*MED_TAILLE_PNOM+1]; 58 /* tables des noms, numeros, numeros de familles des noeuds 59 autant d'elements que de noeuds - les noms ont pout longueur 60 MED_TAILLE_PNOM */ 61 char *nomnoe; 62 med_int *numnoe; 63 med_int *nufano; 64 med_repere rep; 65 med_booleen inonoe,inunoe; 66 char str[MED_TAILLE_PNOM+1]; 67 med_int i; 68 char desc[MED_TAILLE_DESC+1]; 69 med_maillage type; 70 71 /* Ouverture du fichier passe en argument en lecture seule */ 72 if ((fid= MEDouvrir(argv[1],MED_LECTURE)) < 0) { 73 MESSAGE("Erreur a l'ouverture du fichier :"); 74 SSCRUTE(argv[1]); 75 return -1; 76 } 77 78 /* Lecture des infos concernant le premier maillage */ 79 if (MEDmaaInfo(fid,1,maa,&mdim,&type,desc) < 0) { 80 MESSAGE("Erreur a la lecture des informations du 1er maillage"); 81 return -1; 82 } 83 printf("Maillage de nom : %s et de dimension : %d \n",maa,mdim); 84 85 /* Lecture du nombre de noeuds */ 86 if ((nnoe = MEDnEntMaa(fid,maa,MED_COOR,MED_NOEUD,0,0)) < 0) { 87 MESSAGE("Erreur a la lecture du nombre de noeuds "); 88 return -1; 89 } 90 printf("Nombre de noeuds : %d \n",nnoe); 91 92 /* Allocations memoires */ 93 /* table des coordonnees 94 profil : (dimension * nombre de noeuds ) */ 95 if (nnoe > 0) { 96 coo = (med_float*) malloc(sizeof(med_float)*nnoe*mdim); 97 /* table des des numeros, des numeros de familles des noeuds 98 profil : (nombre de noeuds) */ 99 numnoe = (med_int*) malloc(sizeof(med_int)*nnoe); 100 nufano = (med_int*) malloc(sizeof(med_int)*nnoe); 101 /* table des noms des noeuds 102 profil : (nnoe*MED_TAILLE_PNOM+1) */ 103 nomnoe = (char*) malloc(MED_TAILLE_PNOM*nnoe+1); 104 105 /* Lecture des noeuds : 106 - Coordonnees 107 - Noms (optionnel dans un fichier MED) 108 - Numeros (optionnel dans un fichier MED) 109 - Numeros de familles */ 110 if (MEDnoeudsLire(fid,maa,mdim,coo,MED_FULL_INTERLACE,&rep,nomcoo,unicoo, 111 nomnoe,&inonoe,numnoe,&inunoe,nufano,nnoe) < 0) { 112 MESSAGE("Erreur a la lecture des noeuds du maillage"); 113 ret = -1; 114 } 115 116 /* Affichage */ 117 if (ret == 0) { 118 printf("Type de repere : %d \n",rep); 119 printf("Nom des coordonnees : \n"); 120 for (i=0;i<mdim;i++) { 121 strncpy(str,nomcoo+i*MED_TAILLE_PNOM,MED_TAILLE_PNOM); 122 str[MED_TAILLE_PNOM] = '\0'; 123 printf("%s ",str); 124 } 125 printf("\nUnites des coordonnees : \n"); 126 for (i=0;i<mdim;i++) { 127 strncpy(str,unicoo+i*MED_TAILLE_PNOM,MED_TAILLE_PNOM); 128 str[MED_TAILLE_PNOM] = '\0'; 129 printf("%s ",str); 130 } 131 printf("\nCoordonnees des noeuds : \n"); 132 for (i=0;i<nnoe*mdim;i++) 133 printf("%f ",*(coo+i)); 134 if (inonoe) { 135 printf("\nNoms des noeuds : \n"); 136 for (i=0;i<nnoe;i++) { 137 strncpy(str,nomnoe+i*MED_TAILLE_PNOM,MED_TAILLE_PNOM); 138 str[MED_TAILLE_PNOM] = '\0'; 139 printf(" %s ",str); 140 } 141 } 142 if (inunoe) { 143 printf("\nNumeros des noeuds : \n"); 144 for (i=0;i<nnoe;i++) 145 printf("%d ",*(numnoe+i)); 146 } 147 printf("\nNumeros des familles des noeuds : \n"); 148 for (i=0;i<nnoe;i++) 149 printf("%d ",*(nufano+i)); 150 printf("\n"); 151 } 152 153 /* Liberation memoire */ 154 free(coo); 155 free(nomnoe); 156 free(numnoe); 157 free(nufano); 158 } 159 160 /* Fermeture du fichier */ 161 if (MEDfermer(fid) < 0) { 162 MESSAGE("Erreur a la fermeture du fichier"); 163 return -1; 164 } 165 166 return ret; 167 } 168 169 170 171