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 : test7.c 21 * 22 * - Description : lecture des elements du maillage MED crees par test6 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 int main (int argc, char **argv) 39 40 41 { 42 med_err ret = 0; 43 med_idt fid; 44 med_int nse2; 45 med_int *se2_1; 46 med_int *se2_2; 47 char *nomse2; 48 med_int *numse2; 49 med_int *nufase2; 50 med_int ntr3; 51 med_int *tr3; 52 char *nomtr3; 53 med_int *numtr3; 54 med_int *nufatr3; 55 char maa[MED_TAILLE_NOM+1] ="maa1"; 56 med_int mdim = 2; 57 med_booleen inoele,inuele; 58 med_int tse2,ttr3; 59 med_int i; 60 char str[MED_TAILLE_PNOM+1]; 61 med_int profil[2] = { 2, 3 }; 62 char desc[MED_TAILLE_DESC+1]; 63 med_maillage type; 64 65 /* Ouverture du fichier en mode lecture seule */ 66 if ((fid = MEDouvrir("test6.med",MED_LECTURE)) < 0) { 67 MESSAGE("Erreur a l'ouverture du fichier test6.med"); 68 return -1; 69 } 70 71 /* Lecture des informations sur le premier maillage */ 72 if (MEDmaaInfo(fid,1,maa,&mdim,&type,desc) < 0) { 73 MESSAGE("Erreur a la lecture des information sur le 1er maillage"); 74 return -1; 75 } else 76 printf("Maillage de nom : %s et de dimension %d \n",maa,mdim); 77 78 /* Combien de triangles et de segments */ 79 if ((nse2 = MEDnEntMaa(fid,maa,MED_CONN,MED_ARETE,MED_SEG2,MED_DESC)) < 0) { 80 MESSAGE("Erreur a la lecture du nombre de faces MED_SEG2"); 81 return -1; 82 } 83 if ((ntr3 = MEDnEntMaa(fid,maa,MED_CONN,MED_MAILLE,MED_TRIA3,MED_DESC))<0) { 84 MESSAGE("Erreur a la lecture du nombre de mailles MED_TRIA3"); 85 return -1; 86 } 87 printf("Nombre de MED_SEG2 : %d - nombre de MED_TRIA3 : %d\n",nse2,ntr3); 88 89 /* Allocations memoire */ 90 tse2 = 2; 91 se2_1 = (med_int*) calloc(tse2*nse2,sizeof(med_int)); 92 se2_2 = (med_int*) malloc(sizeof(med_int)*tse2*nse2); 93 nomse2 = (char*) malloc(MED_TAILLE_PNOM*nse2+1); 94 numse2 = (med_int*) malloc(sizeof(med_int)*nse2); 95 nufase2 = (med_int*) malloc(sizeof(med_int)*nse2); 96 97 ttr3 = 3; 98 tr3 = (med_int*) malloc(sizeof(med_int)*ntr3*ttr3); 99 nomtr3 = (char*) malloc(MED_TAILLE_PNOM*ntr3+1); 100 numtr3 = (med_int*) malloc(sizeof(med_int)*ntr3); 101 nufatr3 = (med_int*) malloc(sizeof(med_int)*ntr3); 102 103 /* Lecture des connectivites des segments avec profil */ 104 if (MEDconnLire(fid,maa,mdim,se2_1,MED_FULL_INTERLACE,profil,2, 105 MED_ARETE,MED_SEG2,MED_DESC) < 0) { 106 MESSAGE("Erreur a la lecture de la connectivite des segments"); 107 return -1; 108 } 109 110 /* Lecture de la connectivite des segments */ 111 if (MEDconnLire(fid,maa,mdim,se2_2,MED_FULL_INTERLACE,NULL,0, 112 MED_ARETE ,MED_SEG2,MED_DESC) < 0) { 113 MESSAGE("Erreur a la lecture de la connectivite des segments"); 114 return -1; 115 } 116 117 /* Lecture (optionnelle) des noms des segments */ 118 if (MEDnomLire(fid,maa,nomse2,nse2,MED_ARETE,MED_SEG2) < 0) 119 inoele = MED_FAUX; 120 else 121 inoele = MED_VRAI; 122 123 /* Lecture (optionnelle) des numeros des segments */ 124 if (MEDnumLire(fid,maa,numse2,nse2,MED_ARETE,MED_SEG2) < 0) 125 inuele = MED_FAUX; 126 else 127 inuele = MED_VRAI; 128 129 /* Lecture des numeros des familles des segments */ 130 if (MEDfamLire(fid,maa,nufase2,nse2,MED_ARETE,MED_SEG2) < 0) { 131 MESSAGE("Erreur a la lecture des numéros de famille des segments"); 132 return -1; 133 } 134 135 /* Lecture de la connectivite des triangles */ 136 if (MEDconnLire(fid,maa,mdim,tr3,MED_NO_INTERLACE,NULL,0,MED_MAILLE,MED_TRIA3, 137 MED_DESC) < 0) { 138 MESSAGE("Erreur a la lecture de la connectivite des triangles"); 139 return -1; 140 } 141 142 /* Lecture (optionnelle) des noms des triangles */ 143 if (MEDnomLire(fid,maa,nomtr3,ntr3,MED_MAILLE,MED_TRIA3) < 0) 144 inoele = MED_FAUX; 145 else 146 inoele = MED_VRAI; 147 148 /* Lecture (optionnelle) des numeros des triangles */ 149 if (MEDnumLire(fid,maa,numtr3,ntr3,MED_MAILLE,MED_TRIA3) < 0) 150 inuele = MED_FAUX; 151 else 152 inuele = MED_VRAI; 153 154 /* Lecture des numeros des familles des triangles */ 155 if (ret = MEDfamLire(fid,maa,nufatr3,ntr3,MED_MAILLE,MED_TRIA3) < 0) { 156 MESSAGE("Erreur a la lecture des numeros de famille des segments"); 157 return -1; 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 /* Affichage */ 167 if (ret == 0) { 168 printf("Connectivite des segments (1): \n"); 169 for (i=0;i<nse2*tse2;i++) 170 printf("%d ",*(se2_1+i)); 171 printf("\n"); 172 printf("Connectivite des segments (2): \n"); 173 for (i=0;i<nse2*tse2;i++) 174 printf("%d ",*(se2_2+i)); 175 if (inoele) { 176 printf("\nNoms des segments :\n"); 177 for (i=0;i<nse2;i++) { 178 strncpy(str,nomse2+i*MED_TAILLE_PNOM,MED_TAILLE_PNOM); 179 str[MED_TAILLE_PNOM] = '\0'; 180 printf("%s ",str); 181 } 182 } 183 if (inuele) { 184 printf("\nNumeros des segments :\n"); 185 for (i=0;i<nse2;i++) 186 printf("%d ",*(numse2+i)); 187 } 188 printf("\nNumeros des familles des segments :\n"); 189 for (i=0;i<nse2;i++) 190 printf("%d ",*(nufase2+i)); 191 192 printf("\nConnectivite des triangles : \n"); 193 for (i=0;i<ntr3*ttr3;i++) 194 printf("%d ",*(tr3+i)); 195 if (inoele) { 196 printf("\nNoms des triangles :\n"); 197 for (i=0;i<ntr3;i++) { 198 strncpy(str,nomtr3+i*MED_TAILLE_PNOM,MED_TAILLE_PNOM); 199 str[MED_TAILLE_PNOM] = '\0'; 200 printf("%s ",str); 201 } 202 } 203 if (inuele) { 204 printf("\nNumeros des triangles :\n"); 205 for (i=0;i<ntr3;i++) 206 printf("%d ",*(numtr3+i)); 207 } 208 printf("\nNumeros des familles des triangles :\n"); 209 for (i=0;i<ntr3;i++) 210 printf("%d ",*(nufatr3+i)); 211 212 printf("\n"); 213 } 214 215 /* Nettoyage memoire */ 216 free(se2_1); 217 free(se2_2); 218 free(nomse2); 219 free(numse2); 220 free(nufase2); 221 222 free(tr3); 223 free(nomtr3); 224 free(numtr3); 225 free(nufatr3); 226 227 return ret; 228 } 229