src/HDFPersist/HDFdatasetCreate.c

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