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
  27 /*
  28  * - Name : _MEDdatasetWrite
  29  * - Description : writes a HDF dataset
  30  * - Parameters :
  31  *     - id  (IN)     : dataset ID
  32  *     - val  (IN)    : datset values
  33  * - Result :
  34  *     - if success : 0
  35  *     - if failure : -1
  36  */
  37 hdf_err HDFdatasetWrite(hdf_idt id, void *val)
  38 {
  39   hdf_idt datatype;
  40   hdf_err ret;
  41 #ifdef PCLINUX
  42   int isI32BE = 0;
  43   int size = 0;
  44 #endif
  45
  46   if ((datatype = H5Dget_type(id)) < 0)
  47     return -1;
  48
  49 #ifdef PCLINUX
  50   if((H5Tget_class(datatype) == H5T_INTEGER) && (H5Tget_size(datatype) == 4)) {
  51     isI32BE = 1; /* See HDFdatasetCreate */
  52
  53     /*SRN : bug IPAL9619:  replaced the method of getting the size of INT32 dataset */
  54     int i, ndim = HDFdatasetGetnDim(id);
  55     if(ndim < 0) return -1;
  56
  57     hdf_size *dim = (hdf_size *) malloc(sizeof(hdf_size)*ndim);
  58     if ((ret == HDFdatasetGetDim(id, dim)) < 0)  return -1;
  59
  60     for(i=0; i<ndim; i++) size+=dim[i];
  61     free(dim);
  62     /*SRN : end of the fix */
  63
  64     if(size == 0)
  65       return -1;
  66     if(H5Tconvert(H5T_NATIVE_INT, H5T_STD_I32BE, size, (void *)val, NULL, (hid_t)0) < 0)
  67       return -1;
  68   }
  69 #endif
  70
  71   if ((ret = H5Dwrite(id, datatype, H5S_ALL, H5S_ALL, H5P_DEFAULT, val)) < 0)
  72     return -1;
  73
  74 #ifdef PCLINUX
  75   if (isI32BE && (H5Tconvert(H5T_STD_I32BE, H5T_NATIVE_INT, size, (void *)val, NULL, (hid_t)0) < 0))
  76     return -1;
  77 #endif
  78
  79   return 0;
  80 }