src/HDFPersist/HDFdatasetWrite.c

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