src/HDFPersist/HDFdatasetCreate.c

   1 //  Copyright (C) 2007-2008  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.
  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 SALOME HDFPersist : implementation of HDF persitent ( save/ restore )
  24   File   : HDFdatasetCreate.c
  25   Module : SALOME
  26 ----------------------------------------------------------------------------*/
  27
  28 #include "hdfi.h"
  29
  30 /*
  31  * - Name : HDFdatasetCreate
  32  * - Description : creates a HDF dataset
  33  * - Parameters :
  34  *     - pid  (IN)     : father ID
  35  *     - name (IN)     : dataset name
  36  *     - type (IN)     : dataset type (HDF_STRING,HDF_INT32,HDF_INT64,HDF_FLOAT64)
  37  *     - dimd (IN)     : dataset size
  38  *     - order(IN)     : byte order (H5T_ORDER_NONE, H5T_ORDER_LE, H5T_ORDER_BE)
  39  * - Result :
  40  *     - if success : returns dataset ID
  41  *     - if failure : -1
  42  */
  43
  44 hdf_idt HDFdatasetCreate(hdf_idt pid,char *name,hdf_type type,
  45                          hdf_size *dimd, int ndim, hdf_byte_order order)
  46 {
  47   hdf_idt dataset, dataspace = 0;
  48   hdf_err ret;
  49   hdf_idt type_hdf, new_type_hdf = -1;
  50
  51   switch(type)
  52     {
  53     case HDF_FLOAT64 :
  54 #if defined (PCLINUX) || defined (PCLINUX64)
  55       type_hdf = H5T_IEEE_F64BE;
  56 #else
  57       type_hdf = H5T_IEEE_F64LE;
  58 #endif
  59       break;
  60
  61     case HDF_INT32 :
  62 #if defined (PCLINUX) || defined (PCLINUX64)
  63       type_hdf = H5T_STD_I32BE;
  64 #else
  65       type_hdf = H5T_NATIVE_INT;
  66 #endif
  67       break;
  68
  69     case HDF_INT64 :
  70       type_hdf = H5T_NATIVE_LONG;
  71       break;
  72
  73     case HDF_STRING :
  74       if((new_type_hdf = H5Tcopy(H5T_C_S1)) < 0)
  75         return -1;
  76       if((ret = H5Tset_size(new_type_hdf,1)) < 0)
  77         return -1;
  78       break;
  79
  80     default :
  81       return -1;
  82     }
  83
  84   /* set order */
  85   if ( order != H5T_ORDER_ERROR &&
  86        order != H5T_ORDER_NONE &&
  87        type  != HDF_STRING )
  88     {
  89       if (( new_type_hdf = H5Tcopy( type_hdf )) < 0 )
  90         return -1;
  91       if (( ret = H5Tset_order (new_type_hdf, order )) < 0 )
  92         return -1;
  93     }
  94
  95   if ((dataset = H5Dopen(pid,name)) < 0)
  96     {
  97       if ((dataspace = H5Screate_simple(ndim, dimd, NULL)) < 0)
  98         return -1;
  99       if ((dataset = H5Dcreate(pid,name,
 100                                new_type_hdf < 0 ? type_hdf : new_type_hdf,
 101                                dataspace, H5P_DEFAULT)) < 0)
 102         return -1;
 103     }
 104   else
 105     return -1;
 106
 107   if ( ! (new_type_hdf < 0) && (ret = H5Tclose(new_type_hdf)) < 0)
 108     return -1;
 109   if ((ret = H5Sclose(dataspace)) < 0)
 110     return -1;
 111
 112   return dataset;
 113 }