src/HDFPersist/HDFdatasetWrite.c

   1 // Copyright (C) 2005  OPEN CASCADE, CEA, EDF R&D, LEG
   2 //           PRINCIPIA R&D, EADS CCR, Lip6, BV, CEDRAT
   3 // This library is free software; you can redistribute it and/or
   4 // modify it under the terms of the GNU Lesser General Public
   5 // License as published by the Free Software Foundation; either
   6 // version 2.1 of the License.
   7 //
   8 // This library is distributed in the hope that it will be useful
   9 // but 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
  14 // License along with this library; if not, write to the Free Software
  15 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  16 //
  17 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  18 //
  19 /*----------------------------------------------------------------------------
  20 SALOME HDFPersist : implementation of HDF persitent ( save/ restore )
  21   File   : HDFdatasetWrite.c
  22   Module : SALOME
  23 ----------------------------------------------------------------------------*/
  24
  25 #include "hdfi.h"
  26 #include <stdlib.h>
  27
  28 /*
  29  * - Name : _MEDdatasetWrite
  30  * - Description : writes a HDF dataset
  31  * - Parameters :
  32  *     - id  (IN)     : dataset ID
  33  *     - val  (IN)    : datset values
  34  * - Result :
  35  *     - if success : 0
  36  *     - if failure : -1
  37  */
  38 hdf_err HDFdatasetWrite(hdf_idt id, void *val)
  39 {
  40   hdf_idt datatype;
  41   hdf_err ret;
  42 #ifdef PCLINUX
  43   int isI32BE = 0;
  44   int size = 0;
  45 #endif
  46
  47   if ((datatype = H5Dget_type(id)) < 0)
  48     return -1;
  49
  50 #ifdef PCLINUX
  51   if((H5Tget_class(datatype) == H5T_INTEGER) && (H5Tget_size(datatype) == 4)) {
  52     isI32BE = 1; /* See HDFdatasetCreate */
  53
  54     /*SRN : bug IPAL9619:  replaced the method of getting the size of INT32 dataset */
  55     int i, ndim = HDFdatasetGetnDim(id);
  56     if(ndim < 0) return -1;
  57
  58     hdf_size *dim = (hdf_size *) malloc(sizeof(hdf_size)*ndim);
  59     if ((ret = HDFdatasetGetDim(id, dim)) < 0)  return -1;
  60
  61     for(i=0; i<ndim; i++) size+=dim[i];
  62     free(dim);
  63     /*SRN : end of the fix */
  64
  65     if(size == 0)
  66       return -1;
  67     if(H5Tconvert(H5T_NATIVE_INT, H5T_STD_I32BE, size, (void *)val, NULL, (hid_t)0) < 0)
  68       return -1;
  69   }
  70 #endif
  71
  72   if ((ret = H5Dwrite(id, datatype, H5S_ALL, H5S_ALL, H5P_DEFAULT, val)) < 0)
  73     return -1;
  74
  75 #ifdef PCLINUX
  76   if (isI32BE && (H5Tconvert(H5T_STD_I32BE, H5T_NATIVE_INT, size, (void *)val, NULL, (hid_t)0) < 0))
  77     return -1;
  78 #endif
  79
  80   return 0;
  81 }