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.
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.
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
17 *************************************************************************/
20 #include "med_outils.hxx"
25 extern int mode_interlace;
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)
33 med_idt eqid, datagroup;
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];
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)
46 * On inhibe le gestionnaire d'erreur HDF 5
48 _MEDmodeErreurVerrouiller();
51 * Si le Data Group de "eq" n'existe pas => erreur
53 strcpy(chemin,MED_MAA);
55 strcat(chemin,MED_EQS);
57 if ((eqid = _MEDdatagroupOuvrir(fid,chemin)) < 0)
61 * Lecture de l'equivalence
63 if ((ret = _MEDnomEntite(nomdatagroup,typ_ent)) < 0)
65 if ((typ_ent != MED_NOEUD))
67 if ((ret = _MEDnomGeometrie(tmp,typ_geo)) < 0)
69 strcat(nomdatagroup,".");
70 strcat(nomdatagroup,tmp);
72 if ((datagroup = _MEDdatagroupOuvrir(eqid,nomdatagroup)) < 0)
74 #if defined(HAVE_F77INT64)
75 if ((ret = _MEDdatasetNumLire(datagroup,MED_NOM_COR,MED_INT64,
76 MED_NO_INTERLACE,1,MED_ALL,
78 (unsigned char *) corr)) < 0)
81 if ((ret = _MEDdatasetNumLire(datagroup,MED_NOM_COR,MED_INT32,
82 MED_NO_INTERLACE,1,MED_ALL,
84 (unsigned char *) corr)) < 0)
91 if ((ret = _MEDdatagroupFermer(datagroup)) < 0)
93 if ((ret = _MEDdatagroupFermer(eqid)) < 0)