src/MEDWrapper/V2_1/Core/MEDequivLire.cxx

   1 /*************************************************************************
   2 * COPYRIGHT (C) 1999 - 2002  EDF R&D
   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 #include "med.hxx"
  20 #include "med_outils.hxx"
  21
  22 #include <string.h>
  23 #include <stdlib.h>
  24
  25 extern int mode_interlace;
  26
  27 namespace med_2_1{
  28
  29 med_err
  30 MEDequivLire(med_idt fid, char *maa, char *eq, med_int *corr, med_int n,
  31              med_entite_maillage typ_ent,med_geometrie_element typ_geo)
  32 {
  33   med_idt eqid, datagroup;
  34   med_err ret;
  35   char chemin[MED_TAILLE_MAA+MED_TAILLE_EQS+2*MED_TAILLE_NOM+1];
  36   char nomdatagroup[MED_TAILLE_NOM+1];
  37   char tmp[MED_TAILLE_NOM_ENTITE+1];
  38
  39   if (typ_geo == MED_TETRA4 || typ_geo == MED_TETRA10 ||
  40       typ_geo == MED_HEXA8  || typ_geo == MED_HEXA20  ||
  41       typ_geo == MED_PENTA6 || typ_geo == MED_PENTA15 ||
  42       typ_geo == MED_PYRA5  || typ_geo == MED_PYRA13)
  43     return -1;
  44
  45   /*
  46    * On inhibe le gestionnaire d'erreur HDF 5
  47    */
  48   _MEDmodeErreurVerrouiller();
  49
  50   /*
  51    * Si le Data Group de "eq" n'existe pas => erreur
  52    */
  53   strcpy(chemin,MED_MAA);
  54   strcat(chemin,maa);
  55   strcat(chemin,MED_EQS);
  56   strcat(chemin,eq);
  57   if ((eqid = _MEDdatagroupOuvrir(fid,chemin)) < 0)
  58     return -1;
  59
  60   /*
  61    * Lecture de l'equivalence
  62    */
  63   if ((ret = _MEDnomEntite(nomdatagroup,typ_ent)) < 0)
  64     return -1;
  65   if ((typ_ent != MED_NOEUD))
  66     {
  67       if ((ret = _MEDnomGeometrie(tmp,typ_geo)) < 0)
  68         return -1;
  69       strcat(nomdatagroup,".");
  70       strcat(nomdatagroup,tmp);
  71     }
  72   if ((datagroup = _MEDdatagroupOuvrir(eqid,nomdatagroup)) < 0)
  73     return -1;
  74 #if defined(HAVE_F77INT64)
  75   if ((ret =  _MEDdatasetNumLire(datagroup,MED_NOM_COR,MED_INT64,
  76                                  MED_NO_INTERLACE,1,MED_ALL,
  77                                  MED_NOPF,0,MED_NOPG,
  78                                  (unsigned char *) corr)) < 0)
  79     return -1;
  80 #else
  81   if ((ret =  _MEDdatasetNumLire(datagroup,MED_NOM_COR,MED_INT32,
  82                                  MED_NO_INTERLACE,1,MED_ALL,
  83                                  MED_NOPF,0,MED_NOPG,
  84                                  (unsigned char *) corr)) < 0)
  85     return -1;
  86 #endif
  87
  88   /*
  89    * On ferme tout
  90    */
  91   if ((ret = _MEDdatagroupFermer(datagroup)) < 0)
  92     return -1;
  93   if ((ret = _MEDdatagroupFermer(eqid)) < 0)
  94     return -1;
  95
  96   return 0;
  97 }
  98
  99 }