      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 : test17.c
     21  *
     22  * - Description : lecture d'elements de maillages MED ecrits par test16
     23  *                 via les routines de niveau 2
     24  *                 - equivalent a test7.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   med_int nse2;
     47   med_int *se2;
     48   char *nomse2;
     49   med_int *numse2;
     50   med_int *nufase2;
     51   med_int ntr3;
     52   med_int *tr3;
     53   char *nomtr3;
     54   med_int *numtr3;
     55   med_int *nufatr3;
     56   char maa[MED_TAILLE_NOM+1] ="maa1";
     57   med_int mdim;
     58   med_booleen inoele1,inoele2,inuele1,inuele2;
     59   med_int tse2,ttr3;
     60   med_int i;
     61   char str[MED_TAILLE_PNOM+1];
     62   char desc[MED_TAILLE_DESC+1];
     63   med_maillage type;
     64
     65   /* Ouverture du fichier test16.med en lecture seule */
     66   if ((fid = MEDouvrir("test16.med",MED_LECTURE)) < 0) {
     67     MESSAGE("Erreur a l'ouverture du fichier test16.med");
     68     return -1;
     69   }
     70
     71   /* Lecture des infos sur le 1er maillage */
     72   if (MEDmaaInfo(fid,1,maa,&mdim,&type,desc) < 0) {
     73     MESSAGE("Erreur a la lecture des infos sur le 1er maillage");
     74     return -1;
     75   }
     76   printf("Maillage de nom : %s et de dimension %d \n",maa,mdim);
     77
     78   /* Lecture du nombre de triangle et de segment */
     79   if ((nse2 = MEDnEntMaa(fid,maa,MED_CONN,MED_ARETE,MED_SEG2,MED_DESC)) < 0) {
     80     MESSAGE("Erreur a la lecture du nombre de segments ");
     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 triangles ");
     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  = (med_int*) malloc(sizeof(med_int)*tse2*nse2);
     92   nomse2 = (char*) malloc(MED_TAILLE_PNOM*nse2+1);
     93   numse2 = (med_int*) malloc(sizeof(med_int)*nse2);
     94   nufase2 = (med_int*) malloc(sizeof(med_int)*nse2);
     95
     96   ttr3 = 3;
     97   tr3 = (med_int*) malloc(sizeof(med_int)*ntr3*ttr3);
     98   nomtr3 = (char*) malloc(MED_TAILLE_PNOM*ntr3+1);
     99   numtr3 = (med_int*) malloc(sizeof(med_int)*ntr3);
    100   nufatr3 = (med_int*) malloc(sizeof(med_int)*ntr3);
    101
    102   /* Lecture des aretes segments MED_SEG2 :
    103      - Connectivite,
    104      - Noms (optionnel)
    105      - Numeros (optionnel)
    106      - Numeros de familles */
    107   if (MEDelementsLire(fid,maa,mdim,se2,MED_NO_INTERLACE,nomse2,&inoele1,numse2,&inuele1,
    108                       nufase2,nse2,MED_ARETE,MED_SEG2,MED_DESC) < 0) {
    109     MESSAGE("Erreur a la lecture des segments");
    110     ret = -1;
    111   }
    112
    113   /* Lecture des mailles triangles MED_TRIA3 :
    114      - Connectivite,
    115      - Noms (optionnel)
    116      - Numeros (optionnel)
    117      - Numeros de familles */
    118   if (MEDelementsLire(fid,maa,mdim,tr3,MED_NO_INTERLACE,nomtr3,&inoele2,numtr3,&inuele2,
    119                       nufatr3,ntr3,MED_MAILLE,MED_TRIA3,MED_DESC) < 0) {
    120     MESSAGE("Erreur a la lecture des triangles");
    121     ret = -1;
    122   }
    123
    124   /* Fermeture du fichier */
    125   if (MEDfermer(fid) < 0) {
    126     MESSAGE("Erreur a la fermeture du fichier");
    127     ret = -1;
    128   }
    129
    130   /* Affichage */
    131   if (ret == 0) {
    132     if (nse2 > 0) {
    133       printf("Connectivite des segments : \n");
    134       for (i=0;i<nse2*tse2;i++)
    135         printf("%d ",*(se2+i));
    136       if (inoele1) {
    137         printf("\nNoms des segments :\n");
    138         for (i=0;i<nse2;i++) {
    139           strncpy(str,nomse2+i*MED_TAILLE_PNOM,MED_TAILLE_PNOM);
    140           str[MED_TAILLE_PNOM] = '\0';
    141           printf("%s ",str);
    142         }
    143       }
    144       if (inuele1) {
    145         printf("\nNumeros des segments :\n");
    146         for (i=0;i<nse2;i++)
    147           printf("%d ",*(numse2+i));
    148       }
    149       printf("\nNumeros des familles des segments :\n");
    150       for (i=0;i<nse2;i++)
    151         printf("%d ",*(nufase2+i));
    152     }
    153
    154     if (ntr3 > 0) {
    155       printf("\nConnectivite des triangles : \n");
    156       for (i=0;i<ntr3*ttr3;i++)
    157         printf("%d ",*(tr3+i));
    158       if (inoele2) {
    159         printf("\nNoms des triangles :\n");
    160         for (i=0;i<ntr3;i++) {
    161           strncpy(str,nomtr3+i*MED_TAILLE_PNOM,MED_TAILLE_PNOM);
    162           str[MED_TAILLE_PNOM] = '\0';
    163           printf("%s ",str);
    164         }
    165       }
    166       if (inuele2) {
    167         printf("\nNumeros des triangles :\n");
    168         for (i=0;i<ntr3;i++)
    169           printf("%d ",*(numtr3+i));
    170       }
    171       printf("\nNumeros des familles des triangles :\n");
    172       for (i=0;i<ntr3;i++)
    173         printf("%d ",*(nufatr3+i));
    174
    175       printf("\n");
    176     }
    177   }
    178
    179   /* Nettoyage memoire */
    180   free(se2);
    181   free(nomse2);
    182   free(numse2);
    183   free(nufase2);
    184
    185   free(tr3);
    186   free(nomtr3);
    187   free(numtr3);
    188   free(nufatr3);
    189
    190   return ret;
    191 }
    192
    193
    194
    195