/* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * Copyright by The HDF Group. * * All rights reserved. * * * * This file is part of HDF5. The full HDF5 copyright notice, including * * terms governing use, modification, and redistribution, is contained in * * the COPYING file, which can be found at the root of the source code * * distribution tree, or in https://www.hdfgroup.org/licenses. * * If you do not have access to either file, you may request a copy from * * help@hdfgroup.org. * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * */ /* * File/dataset generator for VDS use case 5 * * See the header file for a description. */ #include #include #include "UC_common.h" #include "UC_5.h" int main(void) { hid_t src_sid = H5I_INVALID_HID; /* source dataset's dataspace ID */ hid_t src_dcplid = H5I_INVALID_HID; /* source dataset property list ID */ hid_t vds_sid = H5I_INVALID_HID; /* VDS dataspace ID */ hid_t vds_dcplid = H5I_INVALID_HID; /* VDS dataset property list ID */ hid_t fid = H5I_INVALID_HID; /* HDF5 file ID */ hid_t did = H5I_INVALID_HID; /* dataset ID */ hid_t msid = H5I_INVALID_HID; /* memory dataspace ID */ hid_t fsid = H5I_INVALID_HID; /* file dataspace ID */ hsize_t extent[RANK]; /* source dataset extents */ hsize_t start[RANK]; /* starting point for hyperslab */ hsize_t stride[RANK]; /* hypserslab stride */ hsize_t count[RANK]; /* hypserslab count */ hsize_t map_start = 0; /* starting point in the VDS map */ int *buffer = NULL; /* data buffer */ int value = -1; /* value written to datasets */ hsize_t n = 0; /* number of elements in a plane */ int i; /* iterator */ int j; /* iterator */ hsize_t k; /* iterator */ /* Start by creating the virtual dataset (VDS) dataspace and creation * property list. The individual source datasets are then created * and the VDS map (stored in the VDS property list) is updated. */ /* Create VDS dcpl */ if ((vds_dcplid = H5Pcreate(H5P_DATASET_CREATE)) < 0) UC_ERROR if (H5Pset_fill_value(vds_dcplid, UC_5_VDS_DATATYPE, &UC_5_VDS_FILL_VALUE) < 0) UC_ERROR /* Create VDS dataspace */ if ((vds_sid = H5Screate_simple(RANK, UC_5_VDS_DIMS, UC_5_VDS_MAX_DIMS)) < 0) UC_ERROR /********************************* * Map source files and datasets * *********************************/ /* Hyperslab array setup */ start[0] = 0; start[1] = 0; start[2] = 0; map_start = 0; stride[0] = UC_5_N_SOURCES; stride[1] = 1; stride[2] = 1; count[0] = H5S_UNLIMITED; count[1] = 1; count[2] = 1; extent[0] = UC_5_SRC_PLANES; extent[1] = UC_5_HEIGHT; extent[2] = UC_5_WIDTH; for (i = 0; i < UC_5_N_SOURCES; i++) { /* source dataset dcpl */ if ((src_dcplid = H5Pcreate(H5P_DATASET_CREATE)) < 0) UC_ERROR if (H5Pset_chunk(src_dcplid, RANK, UC_5_PLANE) < 0) UC_ERROR if (H5Pset_fill_value(src_dcplid, UC_5_SOURCE_DATATYPE, &UC_5_FILL_VALUES[i]) < 0) UC_ERROR if (H5Pset_deflate(src_dcplid, COMPRESSION_LEVEL) < 0) UC_ERROR /* Create source file, dataspace, and dataset */ if ((fid = H5Fcreate(UC_5_FILE_NAMES[i], H5F_ACC_TRUNC, H5P_DEFAULT, H5P_DEFAULT)) < 0) UC_ERROR if ((src_sid = H5Screate_simple(RANK, UC_5_SOURCE_DIMS, UC_5_SOURCE_MAX_DIMS)) < 0) UC_ERROR if ((did = H5Dcreate2(fid, UC_5_SOURCE_DSET_NAME, UC_5_SOURCE_DATATYPE, src_sid, H5P_DEFAULT, src_dcplid, H5P_DEFAULT)) < 0) UC_ERROR /* Set the dataset's extent */ if (H5Dset_extent(did, extent) < 0) UC_ERROR /* Create a data buffer that represents a plane */ n = UC_5_PLANE[1] * UC_5_PLANE[2]; if (NULL == (buffer = (int *)malloc(n * sizeof(int)))) UC_ERROR /* Create the memory dataspace */ if ((msid = H5Screate_simple(RANK, UC_5_PLANE, NULL)) < 0) UC_ERROR /* Get the file dataspace */ if ((fsid = H5Dget_space(did)) < 0) UC_ERROR /* Write planes to the dataset */ for (j = 0; j < UC_5_SRC_PLANES; j++) { value = ((i + 1) * 10) + j; for (k = 0; k < n; k++) buffer[k] = value; start[0] = (hsize_t)j; start[1] = 0; start[2] = 0; if (H5Sselect_hyperslab(fsid, H5S_SELECT_SET, start, NULL, UC_5_PLANE, NULL) < 0) UC_ERROR if (H5Dwrite(did, H5T_NATIVE_INT, msid, fsid, H5P_DEFAULT, buffer) < 0) UC_ERROR } /* end for */ /* set up hyperslabs for source and destination datasets */ start[0] = 0; start[1] = 0; start[2] = 0; if (H5Sselect_hyperslab(src_sid, H5S_SELECT_SET, start, NULL, UC_5_SOURCE_MAX_DIMS, NULL) < 0) UC_ERROR start[0] = map_start; if (H5Sselect_hyperslab(vds_sid, H5S_SELECT_SET, start, stride, count, UC_5_PLANE) < 0) UC_ERROR map_start += 1; /* Add VDS mapping */ if (H5Pset_virtual(vds_dcplid, vds_sid, UC_5_FILE_NAMES[i], UC_5_SOURCE_DSET_PATH, src_sid) < 0) UC_ERROR /* close */ if (H5Sclose(msid) < 0) UC_ERROR if (H5Sclose(fsid) < 0) UC_ERROR if (H5Sclose(src_sid) < 0) UC_ERROR if (H5Pclose(src_dcplid) < 0) UC_ERROR if (H5Dclose(did) < 0) UC_ERROR if (H5Fclose(fid) < 0) UC_ERROR free(buffer); } /* end for */ /******************* * Create VDS file * *******************/ /* file */ if ((fid = H5Fcreate(UC_5_VDS_FILE_NAME, H5F_ACC_TRUNC, H5P_DEFAULT, H5P_DEFAULT)) < 0) UC_ERROR /* dataset */ if ((did = H5Dcreate2(fid, UC_5_VDS_DSET_NAME, UC_5_VDS_DATATYPE, vds_sid, H5P_DEFAULT, vds_dcplid, H5P_DEFAULT)) < 0) UC_ERROR /* close */ if (H5Pclose(vds_dcplid) < 0) UC_ERROR if (H5Sclose(vds_sid) < 0) UC_ERROR if (H5Dclose(did) < 0) UC_ERROR if (H5Fclose(fid) < 0) UC_ERROR return EXIT_SUCCESS; error: H5E_BEGIN_TRY { if (src_sid >= 0) (void)H5Sclose(src_sid); if (src_dcplid >= 0) (void)H5Pclose(src_dcplid); if (vds_sid >= 0) (void)H5Sclose(vds_sid); if (vds_dcplid >= 0) (void)H5Pclose(vds_dcplid); if (fid >= 0) (void)H5Fclose(fid); if (did >= 0) (void)H5Dclose(did); if (msid >= 0) (void)H5Sclose(msid); if (fsid >= 0) (void)H5Sclose(fsid); if (buffer != NULL) free(buffer); } H5E_END_TRY return EXIT_FAILURE; } /* end main() */