src/MEDWrapper/V2_1/MEDequivEcr.cxx

   1 // Copyright (C) 2005  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
   2 // CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
   3 //
   4 // This library is free software; you can redistribute it and/or
   5 // modify it under the terms of the GNU Lesser General Public
   6 // License as published by the Free Software Foundation; either
   7 // version 2.1 of the License.
   8 //
   9 // This library is distributed in the hope that it will be useful
  10 // but WITHOUT ANY WARRANTY; without even the implied warranty of
  11 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  12 // Lesser General Public License for more details.
  13 //
  14 // You should have received a copy of the GNU Lesser General Public
  15 // License along with this library; if not, write to the Free Software
  16 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  17 //
  18 // See http://www.salome-platform.org/
  19 //
  20 /*************************************************************************
  21 * COPYRIGHT (C) 1999 - 2002  EDF R&D
  22 * THIS LIBRARY IS FREE SOFTWARE; YOU CAN REDISTRIBUTE IT AND/OR MODIFY
  23 * IT UNDER THE TERMS OF THE GNU LESSER GENERAL PUBLIC LICENSE
  24 * AS PUBLISHED BY THE FREE SOFTWARE FOUNDATION;
  25 * EITHER VERSION 2.1 OF THE LICENSE, OR (AT YOUR OPTION) ANY LATER VERSION.
  26 *
  27 * THIS LIBRARY IS DISTRIBUTED IN THE HOPE THAT IT WILL BE USEFUL, BUT
  28 * WITHOUT ANY WARRANTY; WITHOUT EVEN THE IMPLIED WARRANTY OF
  29 * MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. SEE THE GNU
  30 * LESSER GENERAL PUBLIC LICENSE FOR MORE DETAILS.
  31 *
  32 * YOU SHOULD HAVE RECEIVED A COPY OF THE GNU LESSER GENERAL PUBLIC LICENSE
  33 * ALONG WITH THIS LIBRARY; IF NOT, WRITE TO THE FREE SOFTWARE FOUNDATION,
  34 * INC., 59 TEMPLE PLACE, SUITE 330, BOSTON, MA 02111-1307 USA
  35 *
  36 *************************************************************************/
  37
  38 #include "med.hxx"
  39 #include "med_outils.hxx"
  40
  41 #include <string.h>
  42 #include <stdlib.h>
  43
  44 extern int mode_interlace;
  45
  46 namespace med_2_1{
  47
  48 med_err
  49 MEDequivEcr(med_idt fid, char *maa, char *eq, med_int *corr, med_int n,
  50             med_mode_acces mode, med_entite_maillage typ_ent, med_geometrie_element typ_geo)
  51 {
  52   med_idt eqid, datagroup;
  53   med_err ret;
  54   char chemin[MED_TAILLE_MAA+MED_TAILLE_EQS+2*MED_TAILLE_NOM+1];
  55   char nomdatagroup[MED_TAILLE_NOM+1];
  56   char tmp[MED_TAILLE_NOM_ENTITE+1];
  57   med_size dimd[1];
  58
  59   if (typ_geo == MED_TETRA4 || typ_geo == MED_TETRA10 ||
  60       typ_geo == MED_HEXA8  || typ_geo == MED_HEXA20  ||
  61       typ_geo == MED_PENTA6 || typ_geo == MED_PENTA15 ||
  62       typ_geo == MED_PYRA5  || typ_geo == MED_PYRA13)
  63     return -1;
  64
  65   /*
  66    * On inhibe le gestionnaire d'erreur HDF 5
  67    */
  68   _MEDmodeErreurVerrouiller();
  69
  70   /*
  71    * Si le Data Group de "eq" n'existe pas => erreur
  72    */
  73   strcpy(chemin,MED_MAA);
  74   strcat(chemin,maa);
  75   strcat(chemin,MED_EQS);
  76   strcat(chemin,eq);
  77   if ((eqid = _MEDdatagroupOuvrir(fid,chemin)) < 0)
  78     return -1;
  79
  80   /*
  81    * Ecriture de l'equivalence
  82    */
  83   if ((ret = _MEDnomEntite(nomdatagroup,typ_ent)) < 0)
  84     return -1;
  85   if ((typ_ent != MED_NOEUD))
  86     {
  87       if ((ret = _MEDnomGeometrie(tmp,typ_geo)) < 0)
  88         return -1;
  89       strcat(nomdatagroup,".");
  90       strcat(nomdatagroup,tmp);
  91     }
  92   datagroup = 0;
  93   if (((datagroup = _MEDdatagroupOuvrir(eqid,nomdatagroup)) > 0) &&
  94       (mode != MED_REMP))
  95     return -1;
  96   else
  97     if (datagroup > 0)
  98       _MEDdatagroupFermer(datagroup);
  99
 100   /* EF : verifier que çà marche si le data groupe existe déjà */
 101   if ((datagroup = _MEDdatagroupCreer(eqid,nomdatagroup)) < 0)
 102     return -1;
 103
 104   if ((ret = _MEDattrEntierEcrire(datagroup,MED_NOM_NBR,&n,mode)) < 0)
 105     return -1;
 106
 107   dimd[0] = 2*n;
 108
 109 #if defined(IRIX64)||defined(OSF1)
 110   if ((ret =  _MEDdatasetNumEcrire(datagroup,MED_NOM_COR,MED_INT64,MED_NO_INTERLACE,MED_DIM1,MED_ALL,MED_NOPF,0,MED_NOPG,dimd,
 111                                 (unsigned char*) corr,mode)) < 0)
 112     return -1;
 113 #else
 114   if ((ret =  _MEDdatasetNumEcrire(datagroup,MED_NOM_COR,MED_INT32,MED_NO_INTERLACE,MED_DIM1,MED_ALL,MED_NOPF,0,MED_NOPG,dimd,
 115                                 (unsigned char*) corr,mode)) < 0)
 116     return -1;
 117 #endif
 118
 119   /*
 120    * On ferme tout
 121    */
 122   if ((ret = _MEDdatagroupFermer(datagroup)) < 0)
 123     return -1;
 124   if ((ret = _MEDdatagroupFermer(eqid)) < 0)
 125     return -1;
 126
 127   return 0;
 128 }
 129
 130 }