/* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 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. * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * */ /* * Programmer: Robb Matzke * Thursday, October 1, 1998 * * Purpose: Tests dataset fill values. */ #include "h5test.h" #include "H5srcdir.h" /* * Define NO_FILLING if you want to compare how this test works when there is * no fill value (that is, when the fill value is zero). */ /* #define NO_FILLING */ static const char *FILENAME[] = {"fillval_1", "fillval_2", "fillval_3", "fillval_4", "fillval_5", "fillval_6", "fillval_7", "fillval_8", "fillval_9", NULL}; /* Common type for compound datatype operations */ typedef struct { float a; int x; double y; char z; } comp_datatype; /* Common type for compound+vl datatype operations */ typedef struct { int x; char *a; char *b; int y; } comp_vl_datatype; /* The fill_old.h5 is generated from gen_old_fill.c in HDF5 'test' directory * for version 1.4(after 1.4.3). To get this data file, simply compile * gen_old_fill.c with HDF5 library (before v1.5) and run it. */ #define FILE_COMPATIBLE "fill_old.h5" #define FILE_NAME_RAW "fillval.raw" /*------------------------------------------------------------------------- * Function: create_compound_type * * Purpose: create a compound datatype * * Return: Success: datatype ID * * Failure: -1 * * Programmer: Raymond Lu * Monday, Jan 26, 2001 * *------------------------------------------------------------------------- */ static hid_t create_compound_type(void) { hid_t ret_value = -1; if ((ret_value = H5Tcreate(H5T_COMPOUND, sizeof(comp_datatype))) < 0) goto error; if (H5Tinsert(ret_value, "a", HOFFSET(comp_datatype, a), H5T_NATIVE_FLOAT) < 0) goto error; if (H5Tinsert(ret_value, "x", HOFFSET(comp_datatype, x), H5T_NATIVE_INT) < 0) goto error; if (H5Tinsert(ret_value, "y", HOFFSET(comp_datatype, y), H5T_NATIVE_DOUBLE) < 0) goto error; if (H5Tinsert(ret_value, "z", HOFFSET(comp_datatype, z), H5T_NATIVE_CHAR) < 0) goto error; return ret_value; error: H5E_BEGIN_TRY { H5Tclose(ret_value); } H5E_END_TRY return -1; } /*------------------------------------------------------------------------- * Function: create_compound_vl_type * * Purpose: create a compound+vl datatype * * Return: Success: datatype ID * * Failure: -1 * * Programmer: Quincey Koziol * Tuesday, July 3, 2007 * *------------------------------------------------------------------------- */ static hid_t create_compound_vl_type(void) { hid_t str_id = -1; /* Datatype for VL-string fields */ hid_t ret_value = -1; /* Create a string datatype */ if ((str_id = H5Tcopy(H5T_C_S1)) < 0) TEST_ERROR; if (H5Tset_size(str_id, H5T_VARIABLE) < 0) TEST_ERROR; if ((ret_value = H5Tcreate(H5T_COMPOUND, sizeof(comp_vl_datatype))) < 0) TEST_ERROR; if (H5Tinsert(ret_value, "x", HOFFSET(comp_vl_datatype, x), H5T_NATIVE_INT) < 0) TEST_ERROR; if (H5Tinsert(ret_value, "a", HOFFSET(comp_vl_datatype, a), str_id) < 0) TEST_ERROR; if (H5Tinsert(ret_value, "b", HOFFSET(comp_vl_datatype, b), str_id) < 0) TEST_ERROR; if (H5Tinsert(ret_value, "y", HOFFSET(comp_vl_datatype, y), H5T_NATIVE_INT) < 0) TEST_ERROR; /* Close string datatype */ if (H5Tclose(str_id) < 0) TEST_ERROR; return ret_value; error: H5E_BEGIN_TRY { H5Tclose(str_id); H5Tclose(ret_value); } H5E_END_TRY return -1; } /* end create_compound_vl_type() */ /*------------------------------------------------------------------------- * Function: test_getset * * Purpose: Tests the H5Pget_fill_value() and H5Pset_fill_value() * functions. * * Return: Success: 0 * * Failure: number of errors * * Programmer: Robb Matzke * Thursday, October 1, 1998 * *------------------------------------------------------------------------- */ static int test_getset(void) { hid_t dcpl = -1; int fill_i; hid_t type_ss = -1, type_si = -1; struct fill_si { int v1, v2; } fill_si; struct fill_ss { short v1, v2; } fill_ss, fill_ss_rd; TESTING("property lists"); /* * Create the dataset creation property list and the data types that will * be used during this test. */ if ((dcpl = H5Pcreate(H5P_DATASET_CREATE)) < 0) goto error; if ((type_ss = H5Tcreate(H5T_COMPOUND, sizeof fill_ss)) < 0 || H5Tinsert(type_ss, "v1", HOFFSET(struct fill_ss, v1), H5T_NATIVE_SHORT) < 0 || H5Tinsert(type_ss, "v2", HOFFSET(struct fill_ss, v2), H5T_NATIVE_SHORT) < 0) { goto error; } if ((type_si = H5Tcreate(H5T_COMPOUND, sizeof fill_si)) < 0 || H5Tinsert(type_si, "v1", HOFFSET(struct fill_si, v1), H5T_NATIVE_INT) < 0 || H5Tinsert(type_si, "v2", HOFFSET(struct fill_si, v2), H5T_NATIVE_INT) < 0) { goto error; } /* * Reading the fill value from a dataset creation property list that has * no fill value should result in a failure. */ H5E_BEGIN_TRY { H5Pget_fill_value(dcpl, H5T_NATIVE_INT, &fill_i); } H5E_END_TRY if (fill_i != 0) { H5_FAILED(); HDputs(" H5Pget_fill_value() should return default 0"); goto error; } /* * Set the fill value using a struct as the data type. */ fill_ss.v1 = 1111; fill_ss.v2 = 2222; if (H5Pset_fill_value(dcpl, type_ss, &fill_ss) < 0) goto error; /* * Get the fill value using the same data type that was used to set it. */ if (H5Pget_fill_value(dcpl, type_ss, &fill_ss_rd) < 0) goto error; if (fill_ss.v1 != fill_ss_rd.v1 || fill_ss.v2 != fill_ss_rd.v2) { H5_FAILED(); HDputs(" Failed to get fill value using same data type that was "); HDputs(" used to set the fill value."); goto error; } /* * Get the fill value using some other data type. */ if (H5Pget_fill_value(dcpl, type_si, &fill_si) < 0) goto error; if (fill_ss.v1 != fill_si.v1 || fill_ss.v2 != fill_si.v2) { H5_FAILED(); HDputs(" Failed to get fill value using a data type other than what"); HDputs(" was used to set the fill value."); goto error; } /* * Reset the fill value */ if (H5Pset_fill_value(dcpl, type_si, &fill_si) < 0) goto error; if (H5Pget_fill_value(dcpl, type_ss, &fill_ss) < 0) goto error; if (fill_si.v1 != fill_ss.v1 || fill_si.v2 != fill_ss.v2) { H5_FAILED(); HDputs(" Resetting the fill value was unsuccessful."); goto error; } /* Success */ if (H5Pclose(dcpl) < 0) goto error; if (H5Tclose(type_si) < 0) goto error; if (H5Tclose(type_ss) < 0) goto error; PASSED(); return 0; error: H5E_BEGIN_TRY { H5Pclose(dcpl); H5Tclose(type_si); H5Tclose(type_ss); } H5E_END_TRY return 1; } /*------------------------------------------------------------------------- * Function: test_getset_vl * * Purpose: Tests the H5Pget_fill_value() and H5Pset_fill_value() * functions, using variable-length datatype. * * Return: Success: 0 * Failure: number of errors * * Programmer: Quincey Koziol * Thursday, May 31, 2007 * *------------------------------------------------------------------------- */ static int test_getset_vl(hid_t fapl) { hsize_t dims[1] = {2}; hid_t fileid = (-1), spaceid = (-1), dtypeid = (-1), datasetid = (-1), plistid = (-1); char fill_value[] = "aaaa"; char orig_fill_value[] = "aaaa"; char *f1 = fill_value; char *f2; char filename[1024]; TESTING("property lists, with variable-length datatype"); /* Create string type. */ if ((dtypeid = H5Tcopy(H5T_C_S1)) < 0) TEST_ERROR; if (H5Tset_size(dtypeid, H5T_VARIABLE) < 0) TEST_ERROR; /* Set up dataset creation property list, with fill value */ if ((plistid = H5Pcreate(H5P_DATASET_CREATE)) < 0) TEST_ERROR; if (H5Pset_fill_value(plistid, dtypeid, &f1) < 0) TEST_ERROR; /* Modify original fill value string */ fill_value[0] = 'b'; /* Retrieve fill value from property */ if (H5Pget_fill_value(plistid, dtypeid, &f2) < 0) TEST_ERROR; /* Verify that the fill value is the original value */ if (HDstrcmp(f2, orig_fill_value) != 0) TEST_ERROR; /* Release the fill value retrieved */ HDfree(f2); /* Open file. */ h5_fixname(FILENAME[0], fapl, filename, sizeof filename); if ((fileid = H5Fcreate(filename, H5F_ACC_TRUNC, H5P_DEFAULT, H5P_DEFAULT)) < 0) TEST_ERROR; /* Write an dataset of this type. */ if ((spaceid = H5Screate_simple(1, dims, NULL)) < 0) TEST_ERROR; if ((datasetid = H5Dcreate2(fileid, "Dataset", dtypeid, spaceid, H5P_DEFAULT, plistid, H5P_DEFAULT)) < 0) TEST_ERROR; /* Close IDs (except datatype) */ if (H5Dclose(datasetid) < 0) TEST_ERROR; if (H5Pclose(plistid) < 0) TEST_ERROR; if (H5Sclose(spaceid) < 0) TEST_ERROR; if (H5Fclose(fileid) < 0) TEST_ERROR; /* Re-open file, group & dataset */ if ((fileid = H5Fopen(filename, H5F_ACC_RDONLY, H5P_DEFAULT)) < 0) TEST_ERROR; if ((datasetid = H5Dopen2(fileid, "Dataset", H5P_DEFAULT)) < 0) TEST_ERROR; /* Get dataset's creation property list */ if ((plistid = H5Dget_create_plist(datasetid)) < 0) TEST_ERROR; /* Query fill value */ if (H5Pget_fill_value(plistid, dtypeid, &f2) < 0) TEST_ERROR; /* Verify that the fill value is the original value */ if (HDstrcmp(f2, orig_fill_value) != 0) TEST_ERROR; /* Release the fill value retrieved */ HDfree(f2); /* Close IDs */ if (H5Dclose(datasetid) < 0) TEST_ERROR; if (H5Fclose(fileid) < 0) TEST_ERROR; if (H5Pclose(plistid) < 0) TEST_ERROR; if (H5Tclose(dtypeid) < 0) TEST_ERROR; PASSED(); return 0; error: H5E_BEGIN_TRY { } H5E_END_TRY return 1; } /* end test_getset_vl() */ /*------------------------------------------------------------------------- * Function: test_create * * Purpose: Tests creating datasets that have fill values. * * Return: Success: 0 * * Failure: number of errors * * Programmer: Robb Matzke * Thursday, October 1, 1998 * *------------------------------------------------------------------------- */ static int test_create(hid_t fapl, const char *base_name, H5D_layout_t layout) { hid_t file = -1, space = -1, dcpl = -1, comp_type_id = -1; hid_t dset1 = -1, dset2 = -1, dset3 = -1, dset4 = -1, dset5 = -1, dset6 = -1, /* dset7=-1, */ dset8 = -1, dset9 = -1; hsize_t cur_size[5] = {2, 8, 8, 4, 2}; hsize_t ch_size[5] = {1, 1, 1, 4, 1}; short rd_s, fill_s = 0x1234; long rd_l, fill_l = 0x4321; char filename[1024]; H5D_space_status_t allocation; H5D_alloc_time_t alloc_time; H5D_fill_time_t fill_time; comp_datatype rd_c, fill_ctype; if (H5D_CHUNKED == layout) { TESTING("chunked dataset creation"); } else if (H5D_COMPACT == layout) { TESTING("compact dataset creation"); } else { TESTING("contiguous dataset creation"); } /* * Create a file. */ h5_fixname(base_name, fapl, filename, sizeof filename); if ((file = H5Fcreate(filename, H5F_ACC_TRUNC, H5P_DEFAULT, fapl)) < 0) goto error; if ((space = H5Screate_simple(5, cur_size, cur_size)) < 0) goto error; if ((dcpl = H5Pcreate(H5P_DATASET_CREATE)) < 0) goto error; if (H5D_CHUNKED == layout) { if (H5Pset_chunk(dcpl, 5, ch_size) < 0) goto error; } else if (H5D_COMPACT == layout) { if (H5Pset_layout(dcpl, H5D_COMPACT) < 0) goto error; } /* Create a compound datatype */ if ((comp_type_id = create_compound_type()) < 0) goto error; /* I. Test cases for late space allocation except compact dataset */ if (H5D_COMPACT != layout) { if (H5Pset_alloc_time(dcpl, H5D_ALLOC_TIME_LATE) < 0) goto error; if (H5Pset_fill_time(dcpl, H5D_FILL_TIME_ALLOC) < 0) goto error; /* 1. Compound datatype test */ if (H5Pget_fill_value(dcpl, comp_type_id, &fill_ctype) < 0) goto error; fill_ctype.y = 4444; if (H5Pset_fill_value(dcpl, comp_type_id, &fill_ctype) < 0) goto error; if ((dset9 = H5Dcreate2(file, "dset9", comp_type_id, space, H5P_DEFAULT, dcpl, H5P_DEFAULT)) < 0) goto error; /* The three datasets test three fill * conversion paths: small to large, large to small, and no conversion. * They depend on `short' being smaller than `long'. */ /* 2. Small to large fill conversion */ #ifndef NO_FILLING if (H5Pset_fill_value(dcpl, H5T_NATIVE_SHORT, &fill_s) < 0) goto error; #endif if ((dset1 = H5Dcreate2(file, "dset1", H5T_NATIVE_LONG, space, H5P_DEFAULT, dcpl, H5P_DEFAULT)) < 0) goto error; /* 3. Large to small fill conversion */ #ifndef NO_FILLING if (H5Pset_fill_value(dcpl, H5T_NATIVE_LONG, &fill_l) < 0) goto error; #endif if ((dset2 = H5Dcreate2(file, "dset2", H5T_NATIVE_SHORT, space, H5P_DEFAULT, dcpl, H5P_DEFAULT)) < 0) goto error; /* 4. No conversion */ #ifndef NO_FILLING if (H5Pset_fill_value(dcpl, H5T_NATIVE_LONG, &fill_l) < 0) goto error; #endif if ((dset3 = H5Dcreate2(file, "dset3", H5T_NATIVE_LONG, space, H5P_DEFAULT, dcpl, H5P_DEFAULT)) < 0) goto error; /* 5. late space allocation and never write fill value */ if (H5Pset_fill_time(dcpl, H5D_FILL_TIME_NEVER) < 0) goto error; if ((dset4 = H5Dcreate2(file, "dset4", H5T_NATIVE_LONG, space, H5P_DEFAULT, dcpl, H5P_DEFAULT)) < 0) goto error; /* 6. fill value is undefined while fill write time is H5D_FILL_TIME_ALLOC. * Supposed to fail. */ if (H5Pset_fill_value(dcpl, (hid_t)-1, NULL) < 0) goto error; if (H5Pset_fill_time(dcpl, H5D_FILL_TIME_ALLOC) < 0) goto error; H5E_BEGIN_TRY { if (H5Dcreate2(file, "dset7", H5T_NATIVE_LONG, space, H5P_DEFAULT, dcpl, H5P_DEFAULT) != FAIL) goto error; } H5E_END_TRY } /* II. Test early space allocation cases */ if (H5Pclose(dcpl) < 0) goto error; if ((dcpl = H5Pcreate(H5P_DATASET_CREATE)) < 0) goto error; if (H5D_CHUNKED == layout) { if (H5Pset_chunk(dcpl, 5, ch_size) < 0) goto error; } else if (H5D_COMPACT == layout) { if (H5Pset_layout(dcpl, H5D_COMPACT) < 0) goto error; } if (H5Pset_alloc_time(dcpl, H5D_ALLOC_TIME_EARLY) < 0) goto error; /* 1. Compound datatype test */ if (H5Pget_fill_value(dcpl, comp_type_id, &fill_ctype) < 0) goto error; fill_ctype.y = 4444; if (H5Pset_fill_value(dcpl, comp_type_id, &fill_ctype) < 0) goto error; if ((dset8 = H5Dcreate2(file, "dset8", comp_type_id, space, H5P_DEFAULT, dcpl, H5P_DEFAULT)) < 0) goto error; if (H5Pset_fill_value(dcpl, H5T_NATIVE_LONG, &fill_l) < 0) goto error; /* 2. Never write fill value */ if (H5Pset_fill_time(dcpl, H5D_FILL_TIME_NEVER) < 0) goto error; if ((dset5 = H5Dcreate2(file, "dset5", H5T_NATIVE_INT, space, H5P_DEFAULT, dcpl, H5P_DEFAULT)) < 0) goto error; /* 3. Write fill value at space allocation time */ if (H5Pset_fill_time(dcpl, H5D_FILL_TIME_ALLOC) < 0) goto error; if ((dset6 = H5Dcreate2(file, "dset6", H5T_NATIVE_LONG, space, H5P_DEFAULT, dcpl, H5P_DEFAULT)) < 0) goto error; /* 4. fill value is undefined while fill write time is H5D_FILL_TIME_ALLOC. * Supposed to fail. */ if (H5Pset_fill_value(dcpl, (hid_t)-1, NULL) < 0) goto error; H5E_BEGIN_TRY { if (H5Dcreate2(file, "dset7", H5T_NATIVE_LONG, space, H5P_DEFAULT, dcpl, H5P_DEFAULT) != FAIL) goto error; } H5E_END_TRY /* Close everything */ if (H5D_COMPACT != layout) { if (H5Dclose(dset1) < 0) goto error; if (H5Dclose(dset2) < 0) goto error; if (H5Dclose(dset3) < 0) goto error; if (H5Dclose(dset4) < 0) goto error; if (H5Dclose(dset9) < 0) goto error; } if (H5Dclose(dset5) < 0) goto error; if (H5Dclose(dset6) < 0) goto error; if (H5Dclose(dset8) < 0) goto error; if (H5Sclose(space) < 0) goto error; if (H5Pclose(dcpl) < 0) goto error; if (H5Fclose(file) < 0) goto error; /* Open the file and get the dataset fill value from each dataset */ if ((file = H5Fopen(filename, H5F_ACC_RDONLY, fapl)) < 0) goto error; /* I. Check cases for late space allocation except compact dataset */ if (H5D_COMPACT != layout) { /* 1. Large to small conversion */ if ((dset1 = H5Dopen2(file, "dset1", H5P_DEFAULT)) < 0) goto error; if ((dcpl = H5Dget_create_plist(dset1)) < 0) goto error; #ifndef NO_FILLING if (H5Pget_fill_value(dcpl, H5T_NATIVE_SHORT, &rd_s) < 0) goto error; if (rd_s != fill_s) { H5_FAILED(); HDprintf(" %d: Got a different fill value than what was set.", __LINE__); HDprintf(" Got %d, set %d\n", rd_s, fill_s); goto error; } #endif if (H5Dclose(dset1) < 0) goto error; if (H5Pclose(dcpl) < 0) goto error; /* 2. Small to large conversion */ if ((dset2 = H5Dopen2(file, "dset2", H5P_DEFAULT)) < 0) goto error; if ((dcpl = H5Dget_create_plist(dset2)) < 0) goto error; #ifndef NO_FILLING if (H5Pget_fill_value(dcpl, H5T_NATIVE_LONG, &rd_l) < 0) goto error; if (rd_l != fill_l) { H5_FAILED(); HDprintf(" %d: Got a different fill value than what was set.", __LINE__); HDprintf(" Got %ld, set %ld\n", rd_l, fill_l); goto error; } #endif if (H5Dclose(dset2) < 0) goto error; if (H5Pclose(dcpl) < 0) goto error; /* 3. No conversion */ if ((dset3 = H5Dopen2(file, "dset3", H5P_DEFAULT)) < 0) goto error; if ((dcpl = H5Dget_create_plist(dset3)) < 0) goto error; #ifndef NO_FILLING if (H5Pget_fill_value(dcpl, H5T_NATIVE_LONG, &rd_l) < 0) goto error; if (rd_l != fill_l) { H5_FAILED(); HDprintf(" %d: Got a different fill value than what was set.", __LINE__); HDprintf(" Got %ld, set %ld\n", rd_l, fill_l); goto error; } #endif if (H5Pget_alloc_time(dcpl, &alloc_time) < 0) goto error; if (H5Pget_fill_time(dcpl, &fill_time) < 0) goto error; if (alloc_time != H5D_ALLOC_TIME_LATE) { H5_FAILED(); HDputs(" Got non-H5D_ALLOC_TIME_LATE space allocation time."); HDprintf(" Got %d\n", alloc_time); } if (fill_time != H5D_FILL_TIME_ALLOC) { H5_FAILED(); HDputs(" Got non-H5D_FILL_TIME_ALLOC fill value write time."); HDprintf(" Got %d\n", fill_time); } if (H5Dclose(dset3) < 0) goto error; if (H5Pclose(dcpl) < 0) goto error; /* 4. late space allocation and never write fill value */ if ((dset4 = H5Dopen2(file, "dset4", H5P_DEFAULT)) < 0) goto error; if (H5Dget_space_status(dset4, &allocation) < 0) goto error; if (layout == H5D_CONTIGUOUS && allocation != H5D_SPACE_STATUS_NOT_ALLOCATED) { H5_FAILED(); HDputs(" Got allocated space instead of unallocated."); HDprintf(" Got %d\n", allocation); goto error; } if ((dcpl = H5Dget_create_plist(dset4)) < 0) goto error; if (H5Pget_alloc_time(dcpl, &alloc_time) < 0) goto error; if (H5Pget_fill_time(dcpl, &fill_time) < 0) goto error; if (alloc_time != H5D_ALLOC_TIME_LATE) { H5_FAILED(); HDputs(" Got non-H5D_ALLOC_TIME_LATE space allocation time."); HDprintf(" Got %d\n", alloc_time); } if (fill_time != H5D_FILL_TIME_NEVER) { H5_FAILED(); HDputs(" Got non-H5D_FILL_TIME_NEVER fill value write time."); HDprintf(" Got %d\n", fill_time); } if (H5Dclose(dset4) < 0) goto error; if (H5Pclose(dcpl) < 0) goto error; /* 5. Compound datatype test */ if ((dset9 = H5Dopen2(file, "dset9", H5P_DEFAULT)) < 0) goto error; if ((dcpl = H5Dget_create_plist(dset9)) < 0) goto error; if (H5Pget_fill_value(dcpl, comp_type_id, &rd_c) < 0) goto error; if (!H5_FLT_ABS_EQUAL(rd_c.a, 0) || !H5_DBL_ABS_EQUAL(rd_c.y, fill_ctype.y) || rd_c.x != 0 || rd_c.z != '\0') { H5_FAILED(); HDputs(" Got wrong fill value"); HDprintf(" Got rd_c.a=%f, rd_c.y=%f and rd_c.x=%d, rd_c.z=%c\n", (double)rd_c.a, rd_c.y, rd_c.x, rd_c.z); } if (H5Dclose(dset9) < 0) goto error; if (H5Pclose(dcpl) < 0) goto error; } /* II. Check early space allocation cases */ /* 1. Never write fill value */ if ((dset5 = H5Dopen2(file, "dset5", H5P_DEFAULT)) < 0) goto error; if ((dcpl = H5Dget_create_plist(dset5)) < 0) goto error; if (H5Dget_space_status(dset5, &allocation) < 0) goto error; if (layout == H5D_CONTIGUOUS && allocation != H5D_SPACE_STATUS_ALLOCATED) { H5_FAILED(); HDprintf(" %d: Got unallocated space instead of allocated.\n", __LINE__); HDprintf(" Got %d\n", allocation); goto error; } if (H5Pget_alloc_time(dcpl, &alloc_time) < 0) goto error; if (alloc_time != H5D_ALLOC_TIME_EARLY) { H5_FAILED(); HDputs(" Got non-H5D_ALLOC_TIME_EARLY space allocation time."); HDprintf(" Got %d\n", alloc_time); } if (H5Pget_fill_time(dcpl, &fill_time) < 0) goto error; if (fill_time != H5D_FILL_TIME_NEVER) { H5_FAILED(); HDputs(" Got non-H5D_FILL_TIME_NEVER fill value write time."); HDprintf(" Got %d\n", fill_time); } if (H5Dclose(dset5) < 0) goto error; if (H5Pclose(dcpl) < 0) goto error; /* 2. test writing fill value at space allocation time */ if ((dset6 = H5Dopen2(file, "dset6", H5P_DEFAULT)) < 0) goto error; if ((dcpl = H5Dget_create_plist(dset6)) < 0) goto error; if (H5Dget_space_status(dset6, &allocation) < 0) goto error; if (layout == H5D_CONTIGUOUS && allocation != H5D_SPACE_STATUS_ALLOCATED) { H5_FAILED(); HDprintf(" %d: Got unallocated space instead of allocated.\n", __LINE__); HDprintf(" Got %d\n", allocation); goto error; } if (H5Pget_fill_value(dcpl, H5T_NATIVE_LONG, &rd_l) < 0) goto error; if (rd_l != fill_l) { H5_FAILED(); HDprintf(" %d: Got a different fill value than what was set.", __LINE__); HDprintf(" Got %ld, set %ld\n", rd_l, fill_l); goto error; } if (H5Pget_alloc_time(dcpl, &alloc_time) < 0) goto error; if (alloc_time != H5D_ALLOC_TIME_EARLY) { H5_FAILED(); HDputs(" Got non-H5D_ALLOC_TIME_EARLY space allocation time."); HDprintf(" Got %d\n", alloc_time); } if (H5Pget_fill_time(dcpl, &fill_time) < 0) goto error; if (fill_time != H5D_FILL_TIME_ALLOC) { H5_FAILED(); HDputs(" Got non-H5D_FILL_TIME_ALLOC fill value write time."); HDprintf(" Got %d\n", fill_time); } if (H5Dclose(dset6) < 0) goto error; if (H5Pclose(dcpl) < 0) goto error; /* 3. Compound datatype test */ if ((dset8 = H5Dopen2(file, "dset8", H5P_DEFAULT)) < 0) goto error; if ((dcpl = H5Dget_create_plist(dset8)) < 0) goto error; if (H5Pget_fill_value(dcpl, comp_type_id, &rd_c) < 0) goto error; if (!H5_FLT_ABS_EQUAL(rd_c.a, 0) || !H5_DBL_ABS_EQUAL(rd_c.y, fill_ctype.y) || rd_c.x != 0 || rd_c.z != '\0') { H5_FAILED(); HDputs(" Got wrong fill value"); HDprintf(" Got rd_c.a=%f, rd_c.y=%f and rd_c.x=%d, rd_c.z=%c\n", (double)rd_c.a, rd_c.y, rd_c.x, rd_c.z); } if (H5Dclose(dset8) < 0) goto error; if (H5Pclose(dcpl) < 0) goto error; if (H5Fclose(file) < 0) goto error; PASSED(); return 0; error: H5E_BEGIN_TRY { H5Pclose(dcpl); H5Sclose(space); if (H5D_COMPACT != layout) { H5Dclose(dset1); H5Dclose(dset2); H5Dclose(dset3); H5Dclose(dset4); H5Dclose(dset9); } H5Dclose(dset5); H5Dclose(dset6); H5Dclose(dset8); H5Fclose(file); } H5E_END_TRY return 1; } /*------------------------------------------------------------------------- * Function: test_rdwr_cases * * Purpose: Tests fill values read and write for datasets. * * Return: Success: 0 * * Failure: 1 * * Programmer: Robb Matzke * Thursday, October 1, 1998 * *------------------------------------------------------------------------- */ static int test_rdwr_cases(hid_t file, hid_t dcpl, const char *dname, void *_fillval, H5D_fill_time_t fill_time, H5D_layout_t layout, H5T_class_t datatype, hid_t ctype_id) { hid_t fspace = -1, mspace = -1, dset1 = -1, dset2 = -1; hsize_t cur_size[5] = {2, 8, 8, 4, 2}; hsize_t one[5] = {1, 1, 1, 1, 1}; hsize_t hs_size[5], hs_stride[5]; hsize_t hs_offset[5], nelmts; int fillval = (-1), val_rd, should_be; int i, j, *buf = NULL, odd; unsigned u; comp_datatype rd_c, fill_c, should_be_c; comp_datatype *buf_c = NULL; H5D_space_status_t allocation = H5D_SPACE_STATUS_ERROR; fill_c.a = 0; fill_c.x = 0; fill_c.y = 0; fill_c.z = 0; if (datatype == H5T_INTEGER) { fillval = *(int *)_fillval; } else if (datatype == H5T_COMPOUND) { fill_c.a = ((comp_datatype *)_fillval)->a; fill_c.x = ((comp_datatype *)_fillval)->x; fill_c.y = ((comp_datatype *)_fillval)->y; fill_c.z = ((comp_datatype *)_fillval)->z; } else { HDputs("Invalid type for test"); goto error; } /* Create dataset */ if ((fspace = H5Screate_simple(5, cur_size, cur_size)) < 0) goto error; if (datatype == H5T_INTEGER && (dset1 = H5Dcreate2(file, dname, H5T_NATIVE_INT, fspace, H5P_DEFAULT, dcpl, H5P_DEFAULT)) < 0) goto error; if (datatype == H5T_COMPOUND && (dset2 = H5Dcreate2(file, dname, ctype_id, fspace, H5P_DEFAULT, dcpl, H5P_DEFAULT)) < 0) goto error; /* Read some data and make sure it's the fill value */ if ((mspace = H5Screate_simple(5, one, NULL)) < 0) goto error; for (i = 0; i < 1000; i++) { for (j = 0; j < 5; j++) hs_offset[j] = (hsize_t)HDrand() % cur_size[j]; if (H5Sselect_hyperslab(fspace, H5S_SELECT_SET, hs_offset, NULL, one, NULL) < 0) goto error; /* case for atomic datatype */ if (datatype == H5T_INTEGER) { if (H5Dread(dset1, H5T_NATIVE_INT, mspace, fspace, H5P_DEFAULT, &val_rd) < 0) goto error; if (fill_time != H5D_FILL_TIME_NEVER && val_rd != fillval) { H5_FAILED(); HDfprintf(stdout, "%u: Value read was not a fill value.\n", (unsigned)__LINE__); HDfprintf(stdout, " Elmt={%" PRIuHSIZE ",%" PRIuHSIZE ",%" PRIuHSIZE ",%" PRIuHSIZE ",%" PRIuHSIZE "}, read: %u, " "Fill value: %u\n", hs_offset[0], hs_offset[1], hs_offset[2], hs_offset[3], hs_offset[4], val_rd, fillval); goto error; } /* case for compound datatype */ } else if (datatype == H5T_COMPOUND) { if (H5Dread(dset2, ctype_id, mspace, fspace, H5P_DEFAULT, &rd_c) < 0) goto error; if (fill_time != H5D_FILL_TIME_NEVER && (!H5_FLT_ABS_EQUAL(rd_c.a, fill_c.a) || rd_c.x != fill_c.x || !H5_DBL_ABS_EQUAL(rd_c.y, fill_c.y) || rd_c.z != fill_c.z)) { H5_FAILED(); HDfprintf(stdout, "%u: Value read was not a fill value.\n", (unsigned)__LINE__); HDfprintf(stdout, " Elmt={%" PRIuHSIZE ",%" PRIuHSIZE ",%" PRIuHSIZE ",%" PRIuHSIZE ",%" PRIuHSIZE "}, read: %f, %d, %f, %c" "Fill value: %f, %d, %f, %c\n", hs_offset[0], hs_offset[1], hs_offset[2], hs_offset[3], hs_offset[4], (double)rd_c.a, rd_c.x, rd_c.y, rd_c.z, (double)fill_c.a, fill_c.x, fill_c.y, fill_c.z); goto error; } } } if (H5Sclose(mspace) < 0) goto error; /* Select all odd data locations in the file dataset */ for (i = 0, nelmts = 1; i < 5; i++) { hs_size[i] = cur_size[i] / 2; hs_offset[i] = 0; hs_stride[i] = 2; nelmts *= hs_size[i]; } if ((mspace = H5Screate_simple(5, hs_size, hs_size)) < 0) goto error; if (H5Sselect_hyperslab(fspace, H5S_SELECT_SET, hs_offset, hs_stride, hs_size, NULL) < 0) goto error; /* Read non-contiguous selection from empty dataset */ /* case for atomic datatype */ if (datatype == H5T_INTEGER) { /*check for overflow*/ assert((nelmts * sizeof(int)) == (hsize_t)((size_t)(nelmts * sizeof(int)))); buf = (int *)HDmalloc((size_t)(nelmts * sizeof(int))); if (H5Dread(dset1, H5T_NATIVE_INT, mspace, fspace, H5P_DEFAULT, buf) < 0) goto error; /* Verify values, except if no fill value written */ if (fill_time != H5D_FILL_TIME_NEVER) { for (u = 0; u < nelmts; u++) { if (buf[u] != fillval) { H5_FAILED(); HDfprintf(stdout, "%u: Value read was not a fill value.\n", (unsigned)__LINE__); HDfprintf(stdout, " Elmt={%" PRIuHSIZE ", %" PRIuHSIZE ", %" PRIuHSIZE ", %" PRIuHSIZE ", %" PRIuHSIZE "}, read: %u, " "Fill value: %u\n", hs_offset[0], hs_offset[1], hs_offset[2], hs_offset[3], hs_offset[4], buf[u], fillval); goto error; } /* end if */ } /* end for */ } /* end if */ } /* case for compound datatype */ else if (datatype == H5T_COMPOUND) { /*check for overflow*/ assert((nelmts * sizeof(comp_datatype)) == (hsize_t)((size_t)(nelmts * sizeof(comp_datatype)))); buf_c = (comp_datatype *)HDmalloc((size_t)nelmts * sizeof(comp_datatype)); if (H5Dread(dset2, ctype_id, mspace, fspace, H5P_DEFAULT, buf_c) < 0) goto error; /* Verify values, except if no fill value written */ if (fill_time != H5D_FILL_TIME_NEVER) { for (u = 0; u < nelmts; u++) { if (!H5_FLT_ABS_EQUAL(buf_c[u].a, fill_c.a) || buf_c[u].x != fill_c.x || !H5_DBL_ABS_EQUAL(buf_c[u].y, fill_c.y) || buf_c[u].z != fill_c.z) { H5_FAILED(); HDfprintf(stdout, "%u: Value read was not a fill value.\n", (unsigned)__LINE__); HDfprintf(stdout, " Elmt={%" PRIuHSIZE ", %" PRIuHSIZE ", %" PRIuHSIZE ", %" PRIuHSIZE ", %" PRIuHSIZE "}, read: %f, %d, %f, %c" "Fill value: %f, %d, %f, %c\n", hs_offset[0], hs_offset[1], hs_offset[2], hs_offset[3], hs_offset[4], (double)buf_c[u].a, buf_c[u].x, buf_c[u].y, buf_c[u].z, (double)fill_c.a, fill_c.x, fill_c.y, fill_c.z); goto error; } /* end if */ } /* end for */ } /* end if */ } /* Write to all odd data locations */ /* case for atomic datatype */ if (datatype == H5T_INTEGER) { for (u = 0; u < nelmts; u++) buf[u] = 9999; if (H5Dwrite(dset1, H5T_NATIVE_INT, mspace, fspace, H5P_DEFAULT, buf) < 0) goto error; } /* case for compound datatype */ else if (datatype == H5T_COMPOUND) { HDmemset(buf_c, 0, ((size_t)nelmts * sizeof(comp_datatype))); for (u = 0; u < nelmts; u++) { buf_c[u].a = 1111.11F; buf_c[u].x = 2222; buf_c[u].y = 3333.3333; buf_c[u].z = 'd'; } if (H5Dwrite(dset2, ctype_id, mspace, fspace, H5P_DEFAULT, buf_c) < 0) goto error; } /* Check if space is allocated */ if (datatype == H5T_INTEGER && H5Dget_space_status(dset1, &allocation) < 0) goto error; if (datatype == H5T_COMPOUND && H5Dget_space_status(dset2, &allocation) < 0) goto error; if (layout == H5D_CONTIGUOUS && allocation != H5D_SPACE_STATUS_ALLOCATED) { H5_FAILED(); HDprintf(" %d: Got unallocated space instead of allocated.\n", __LINE__); HDprintf(" Got %d\n", allocation); goto error; } HDfree(buf); buf = NULL; H5Sclose(mspace); /* Read some data and make sure it's the right value */ if ((mspace = H5Screate_simple(5, one, NULL)) < 0) goto error; for (i = 0; i < 1000; i++) { for (j = 0, odd = 0; j < 5; j++) { hs_offset[j] = (hsize_t)HDrand() % cur_size[j]; odd += (int)(hs_offset[j] % 2); } /* end for */ if (H5Sselect_hyperslab(fspace, H5S_SELECT_SET, hs_offset, NULL, one, NULL) < 0) goto error; /* case for atomic datatype */ if (datatype == H5T_INTEGER) { if (H5Dread(dset1, H5T_NATIVE_INT, mspace, fspace, H5P_DEFAULT, &val_rd) < 0) goto error; if (fill_time == H5D_FILL_TIME_ALLOC) { should_be = odd ? fillval : 9999; if (val_rd != should_be) { H5_FAILED(); HDfprintf(stdout, "%u: Value read was not correct.\n", (unsigned)__LINE__); HDprintf(" Elmt={%ld,%ld,%ld,%ld,%ld}, read: %u, " "should be: %u\n", (long)hs_offset[0], (long)hs_offset[1], (long)hs_offset[2], (long)hs_offset[3], (long)hs_offset[4], val_rd, should_be); goto error; } } else if (fill_time == H5D_FILL_TIME_NEVER && !odd) { should_be = 9999; if (val_rd != should_be) { H5_FAILED(); HDfprintf(stdout, "%u: Value read was not correct.\n", (unsigned)__LINE__); HDprintf(" Elmt={%ld,%ld,%ld,%ld,%ld}, read: %u, " "should be: %u\n", (long)hs_offset[0], (long)hs_offset[1], (long)hs_offset[2], (long)hs_offset[3], (long)hs_offset[4], val_rd, should_be); goto error; } } else if (fill_time == H5D_FILL_TIME_NEVER && odd) { /*Trash data. Don't compare*/ } } /* end for datatype==H5T_INTEGER */ /* case for compound datatype */ else if (datatype == H5T_COMPOUND) { if (H5Dread(dset2, ctype_id, mspace, fspace, H5P_DEFAULT, &rd_c) < 0) goto error; if (fill_time == H5D_FILL_TIME_ALLOC) { if (odd) { should_be_c.a = fill_c.a; should_be_c.x = fill_c.x; should_be_c.y = fill_c.y; should_be_c.z = fill_c.z; } else { should_be_c.a = buf_c[0].a; should_be_c.x = buf_c[0].x; should_be_c.y = buf_c[0].y; should_be_c.z = buf_c[0].z; } if (!H5_FLT_ABS_EQUAL(rd_c.a, should_be_c.a) || rd_c.x != should_be_c.x || !H5_DBL_ABS_EQUAL(rd_c.y, should_be_c.y) || rd_c.z != should_be_c.z) { H5_FAILED(); HDfprintf(stdout, "%u: Value read was not correct.\n", (unsigned)__LINE__); HDprintf(" Elmt={%ld,%ld,%ld,%ld,%ld}, read: %f,%d,%f,%c " "should be: %f,%d,%f,%c\n", (long)hs_offset[0], (long)hs_offset[1], (long)hs_offset[2], (long)hs_offset[3], (long)hs_offset[4], (double)rd_c.a, rd_c.x, rd_c.y, rd_c.z, (double)should_be_c.a, should_be_c.x, should_be_c.y, should_be_c.z); goto error; } } /* end for fill_time == H5D_FILL_TIME_ALLOC */ else if (fill_time == H5D_FILL_TIME_NEVER && !odd) { should_be_c.a = buf_c[0].a; should_be_c.x = buf_c[0].x; should_be_c.y = buf_c[0].y; should_be_c.z = buf_c[0].z; if (!H5_FLT_ABS_EQUAL(rd_c.a, should_be_c.a) || rd_c.x != should_be_c.x || !H5_DBL_ABS_EQUAL(rd_c.y, should_be_c.y) || rd_c.z != should_be_c.z) { H5_FAILED(); HDfprintf(stdout, "%u: Value read was not correct.\n", (unsigned)__LINE__); HDprintf(" Elmt={%ld,%ld,%ld,%ld,%ld}, read: %f,%d,%f,%c " "should be: %f,%d,%f,%c\n", (long)hs_offset[0], (long)hs_offset[1], (long)hs_offset[2], (long)hs_offset[3], (long)hs_offset[4], (double)rd_c.a, rd_c.x, rd_c.y, rd_c.z, (double)should_be_c.a, should_be_c.x, should_be_c.y, should_be_c.z); goto error; } } /* end for fill_time == H5D_FILL_TIME_NEVER */ else if (fill_time == H5D_FILL_TIME_NEVER && odd) { /*Trash data. Don't compare*/ } } /* end for datatype==H5T_COMPOUND */ } if (datatype == H5T_COMPOUND) { HDfree(buf_c); buf_c = NULL; } /* end if */ if (H5Sclose(mspace) < 0) goto error; if (datatype == H5T_INTEGER && H5Dclose(dset1) < 0) goto error; if (datatype == H5T_COMPOUND && H5Dclose(dset2) < 0) goto error; if (H5Sclose(fspace) < 0) goto error; return 0; error: H5E_BEGIN_TRY { if (datatype == H5T_INTEGER) H5Dclose(dset1); if (datatype == H5T_COMPOUND) H5Dclose(dset2); H5Sclose(fspace); H5Sclose(mspace); } H5E_END_TRY return 1; } /*------------------------------------------------------------------------- * Function: test_rdwr * * Purpose: Tests fill values for datasets. * * Return: Success: 0 * * Failure: number of errors * * Programmer: Robb Matzke * Thursday, October 1, 1998 * *------------------------------------------------------------------------- */ static int test_rdwr(hid_t fapl, const char *base_name, H5D_layout_t layout) { char filename[1024]; hid_t file = -1, dcpl = -1, ctype_id = -1; hsize_t ch_size[5] = {2, 8, 8, 4, 2}; int nerrors = 0; int fillval = 0x4c70f1cd; comp_datatype fill_ctype = {0, 0, 0, 0}; if (H5D_CHUNKED == layout) { TESTING("chunked dataset I/O"); } else if (H5D_COMPACT == layout) { TESTING("compact dataset I/O"); } else { TESTING("contiguous dataset I/O"); } h5_fixname(base_name, fapl, filename, sizeof filename); if ((file = H5Fcreate(filename, H5F_ACC_TRUNC, H5P_DEFAULT, fapl)) < 0) goto error; if ((dcpl = H5Pcreate(H5P_DATASET_CREATE)) < 0) goto error; if (H5D_CHUNKED == layout) { if (H5Pset_chunk(dcpl, 5, ch_size) < 0) goto error; } else if (H5D_COMPACT == layout) { if (H5Pset_layout(dcpl, H5D_COMPACT) < 0) goto error; } if ((ctype_id = create_compound_type()) < 0) goto error; /* I. Test H5D_ALLOC_TIME_LATE space allocation cases */ if (H5D_COMPACT != layout) { if (H5Pset_alloc_time(dcpl, H5D_ALLOC_TIME_LATE) < 0) goto error; /* case for H5D_FILL_TIME_ALLOC as fill write time and fill value to be default */ if (H5Pset_fill_time(dcpl, H5D_FILL_TIME_ALLOC) < 0) goto error; fillval = 0; nerrors += test_rdwr_cases(file, dcpl, "dset1", &fillval, H5D_FILL_TIME_ALLOC, layout, H5T_INTEGER, (hid_t)-1); /* case for H5D_FILL_TIME_NEVER as fill write time and fill value to be default */ if (H5Pset_fill_time(dcpl, H5D_FILL_TIME_NEVER) < 0) goto error; nerrors += test_rdwr_cases(file, dcpl, "dset2", &fillval, H5D_FILL_TIME_NEVER, layout, H5T_INTEGER, (hid_t)-1); /* case for H5D_FILL_TIME_ALLOC as fill write time and fill value is user-defined */ if (H5Pset_fill_time(dcpl, H5D_FILL_TIME_ALLOC) < 0) goto error; fillval = 0x4c70f1cd; if (H5Pset_fill_value(dcpl, H5T_NATIVE_INT, &fillval) < 0) goto error; nerrors += test_rdwr_cases(file, dcpl, "dset3", &fillval, H5D_FILL_TIME_ALLOC, layout, H5T_INTEGER, (hid_t)-1); /* case for H5D_FILL_TIME_NEVER as fill write time and fill value is user-defined */ if (H5Pset_fill_time(dcpl, H5D_FILL_TIME_NEVER) < 0) goto error; if (H5Pset_fill_value(dcpl, H5T_NATIVE_INT, &fillval) < 0) goto error; nerrors += test_rdwr_cases(file, dcpl, "dset4", &fillval, H5D_FILL_TIME_NEVER, layout, H5T_INTEGER, (hid_t)-1); /* case for H5D_FILL_TIME_ALLOC as fill write time and fill value is undefined */ /* This case has been tested in test_create() function */ /* case for H5D_FILL_TIME_NEVER as fill write time and fill value is undefined */ if (H5Pset_fill_time(dcpl, H5D_FILL_TIME_NEVER) < 0) goto error; if (H5Pset_fill_value(dcpl, (hid_t)-1, NULL) < 0) goto error; nerrors += test_rdwr_cases(file, dcpl, "dset5", &fillval, H5D_FILL_TIME_NEVER, layout, H5T_INTEGER, (hid_t)-1); /* case for H5D_FILL_TIME_ALLOC as fill write time and fill value is user-defined * as compound type */ if (H5Pset_fill_time(dcpl, H5D_FILL_TIME_ALLOC) < 0) goto error; HDmemset(&fill_ctype, 0, sizeof(fill_ctype)); fill_ctype.y = 4444.4444; if (H5Pset_fill_value(dcpl, ctype_id, &fill_ctype) < 0) goto error; nerrors += test_rdwr_cases(file, dcpl, "dset11", &fill_ctype, H5D_FILL_TIME_ALLOC, layout, H5T_COMPOUND, ctype_id); if (H5Pclose(dcpl) < 0) goto error; if ((dcpl = H5Pcreate(H5P_DATASET_CREATE)) < 0) goto error; if (H5D_CHUNKED == layout) { if (H5Pset_chunk(dcpl, 5, ch_size) < 0) goto error; } } /* II. Test H5D_ALLOC_TIME_EARLY space allocation cases */ if (H5Pset_alloc_time(dcpl, H5D_ALLOC_TIME_EARLY) < 0) goto error; /* case for H5D_FILL_TIME_ALLOC as fill write time and fill value to be default */ if (H5Pset_fill_time(dcpl, H5D_FILL_TIME_ALLOC) < 0) goto error; fillval = 0; nerrors += test_rdwr_cases(file, dcpl, "dset6", &fillval, H5D_FILL_TIME_ALLOC, layout, H5T_INTEGER, (hid_t)-1); /* case for H5D_FILL_TIME_NEVER as fill write time and fill value to be default */ if (H5Pset_fill_time(dcpl, H5D_FILL_TIME_NEVER) < 0) goto error; nerrors += test_rdwr_cases(file, dcpl, "dset7", &fillval, H5D_FILL_TIME_NEVER, layout, H5T_INTEGER, (hid_t)-1); /* case for H5D_FILL_TIME_ALLOC as fill write time and fill value is user-defined */ if (H5Pset_fill_time(dcpl, H5D_FILL_TIME_ALLOC) < 0) goto error; fillval = 0x4c70f1cd; if (H5Pset_fill_value(dcpl, H5T_NATIVE_INT, &fillval) < 0) goto error; nerrors += test_rdwr_cases(file, dcpl, "dset8", &fillval, H5D_FILL_TIME_ALLOC, layout, H5T_INTEGER, (hid_t)-1); /* case for H5D_FILL_TIME_NEVER as fill write time and fill value is user-defined */ if (H5Pset_fill_time(dcpl, H5D_FILL_TIME_NEVER) < 0) goto error; if (H5Pset_fill_value(dcpl, H5T_NATIVE_INT, &fillval) < 0) goto error; nerrors += test_rdwr_cases(file, dcpl, "dset9", &fillval, H5D_FILL_TIME_NEVER, layout, H5T_INTEGER, (hid_t)-1); /* case for H5D_FILL_TIME_ALLOC as fill write time and fill value is undefined */ /* This case has been tested in test_create() function */ /* case for H5D_FILL_TIME_NEVER as fill write time and fill value is undefined */ if (H5Pset_fill_time(dcpl, H5D_FILL_TIME_NEVER) < 0) goto error; if (H5Pset_fill_value(dcpl, (hid_t)-1, NULL) < 0) goto error; nerrors += test_rdwr_cases(file, dcpl, "dset10", &fillval, H5D_FILL_TIME_NEVER, layout, H5T_INTEGER, (hid_t)-1); /* case for H5D_FILL_TIME_ALLOC as fill write time and fill value is user-defined * as compound type */ if (H5Pset_fill_time(dcpl, H5D_FILL_TIME_ALLOC) < 0) goto error; HDmemset(&fill_ctype, 0, sizeof(fill_ctype)); fill_ctype.y = 4444.4444; if (H5Pset_fill_value(dcpl, ctype_id, &fill_ctype) < 0) goto error; nerrors += test_rdwr_cases(file, dcpl, "dset12", &fill_ctype, H5D_FILL_TIME_ALLOC, layout, H5T_COMPOUND, ctype_id); if (nerrors) goto error; if (H5Pclose(dcpl) < 0) goto error; if (H5Tclose(ctype_id) < 0) goto error; if (H5Fclose(file) < 0) goto error; PASSED(); return 0; error: H5E_BEGIN_TRY { H5Pclose(dcpl); H5Tclose(ctype_id); H5Fclose(file); } H5E_END_TRY return nerrors; } /*------------------------------------------------------------------------- * Function: test_extend_init_integer * * Purpose: Initializes integer values * * Return: Success: 0 * Failure: < 0 * * Programmer: Quincey Koziol * Tuesday, July 3, 2007 * *------------------------------------------------------------------------- */ static int test_extend_init_integer(void *_buf, size_t nelmts, const void *_val) { int *buf = (int *)_buf; /* Buffer to initialize */ const int *val = (const int *)_val; /* Value to use */ while (nelmts) { *buf++ = *val; nelmts--; } /* end while */ return 0; } /* end test_extend_init_integer() */ /*------------------------------------------------------------------------- * Function: test_extend_verify_integer * * Purpose: Verifies integer values * * Return: Success: 0 * Failure: < 0 * * Programmer: Quincey Koziol * Tuesday, July 3, 2007 * *------------------------------------------------------------------------- */ static int test_extend_verify_integer(unsigned lineno, const hsize_t *offset, const void *_test_val, const void *_compare_val) { const int *test_val = (const int *)_test_val; /* Value to test */ const int *compare_val = (const int *)_compare_val; /* Value to compare against */ /* Verify value */ if (*test_val != *compare_val) { HDfprintf(stdout, "%u: Value read was not expected.\n", lineno); HDfprintf(stdout, " Elmt = {%" PRIuHSIZE ", %" PRIuHSIZE ", %" PRIuHSIZE ", %" PRIuHSIZE ", %" PRIuHSIZE "}, read: %d, " "expected: %d\n", offset[0], offset[1], offset[2], offset[3], offset[4], *test_val, *compare_val); goto error; } /* end if */ return 0; error: return -1; } /* end test_extend_verify_integer() */ /*------------------------------------------------------------------------- * Function: test_extend_release_integer * * Purpose: Release element of integer value * * Return: Success: 0 * Failure: < 0 * * Programmer: Quincey Koziol * Tuesday, July 3, 2007 * *------------------------------------------------------------------------- */ static int test_extend_release_integer(void H5_ATTR_UNUSED *_elmt) { return 0; } /* end test_extend_release_integer() */ /*------------------------------------------------------------------------- * Function: test_extend_init_cmpd_vl * * Purpose: Initializes compound+vl values * * Return: Success: 0 * Failure: < 0 * * Programmer: Quincey Koziol * Tuesday, July 3, 2007 * *------------------------------------------------------------------------- */ static int test_extend_init_cmpd_vl(void *_buf, size_t nelmts, const void *_val) { comp_vl_datatype *buf = (comp_vl_datatype *)_buf; /* Buffer to initialize */ const comp_vl_datatype *val = (const comp_vl_datatype *)_val; /* Value to use */ while (nelmts) { /* Shallow copy all fields */ *buf = *val; /* Deep copy string fields */ buf->a = HDstrdup(val->a); buf->b = HDstrdup(val->b); buf++; nelmts--; } /* end while */ return 0; } /* end test_extend_init_cmpd_vl() */ /*------------------------------------------------------------------------- * Function: test_extend_verify_cmpd_vl * * Purpose: Verifies compound+vl values * * Return: Success: 0 * Failure: < 0 * * Programmer: Quincey Koziol * Tuesday, July 3, 2007 * *------------------------------------------------------------------------- */ static int test_extend_verify_cmpd_vl(unsigned lineno, const hsize_t *offset, const void *_test_val, const void *_compare_val) { const comp_vl_datatype *test_val = (const comp_vl_datatype *)_test_val; /* Value to test */ const comp_vl_datatype *compare_val = (const comp_vl_datatype *)_compare_val; /* Value to compare against */ /* Verify value */ if ((test_val->x != compare_val->x) || HDstrcmp(test_val->a, compare_val->a) != 0 || HDstrcmp(test_val->b, compare_val->b) != 0 || (test_val->y != compare_val->y)) { HDfprintf(stdout, "%u: Value read was not expected.\n", lineno); HDfprintf(stdout, " Elmt = {%" PRIuHSIZE ", %" PRIuHSIZE ", %" PRIuHSIZE ", %" PRIuHSIZE ", %" PRIuHSIZE "}, read: {%d, '%s', '%s', %d} " "expected: {%d, '%s', '%s', %d}\n", offset[0], offset[1], offset[2], offset[3], offset[4], test_val->x, test_val->a, test_val->b, test_val->y, compare_val->x, compare_val->a, compare_val->b, compare_val->y); goto error; } /* end if */ return 0; error: return -1; } /* end test_extend_verify_cmpd_vl() */ /*------------------------------------------------------------------------- * Function: test_extend_release_cmpd_vl * * Purpose: Release element of compound+vl value * * Return: Success: 0 * Failure: < 0 * * Programmer: Quincey Koziol * Tuesday, July 3, 2007 * *------------------------------------------------------------------------- */ static int test_extend_release_cmpd_vl(void *_elmt) { comp_vl_datatype *elmt = (comp_vl_datatype *)_elmt; /* Element to free */ /* Free memory for string fields */ HDfree(elmt->a); HDfree(elmt->b); return 0; } /* end test_extend_release_integer() */ /*------------------------------------------------------------------------- * Function: test_extend_cases * * Purpose: Called to test fill values with various different values * * Return: Success: 0 * Failure: number of errors * * Programmer: Quincey Koziol * Tuesday, July 3, 2007 * *------------------------------------------------------------------------- */ static int test_extend_cases(hid_t file, hid_t _dcpl, const char *dset_name, const hsize_t *ch_size, hsize_t *start_size, hsize_t *max_size, hid_t dtype, void *fillval) { hid_t fspace = -1, mspace = -1; /* File & memory dataspaces */ hid_t dset = -1; /* Dataset ID */ hid_t dcpl = -1; /* Dataset creation property list */ hsize_t extend_size[5]; /* Dimensions to extend to */ hsize_t one[5] = {1, 1, 1, 1, 1}; /* Dimensions of single element dataspace */ hsize_t hs_size[5], hs_stride[5], hs_offset[5]; size_t nelmts; H5T_class_t dtype_class; /* Class of datatype */ int (*init_rtn)(void *, size_t, const void *); int (*verify_rtn)(unsigned, const hsize_t *, const void *, const void *); int (*release_rtn)(void *); size_t val_size; /* Size of element */ void *val_rd, *odd_val; const void *init_val, *should_be, *even_val; int val_rd_i, init_val_i = 9999; comp_vl_datatype init_val_c = {87, NULL, NULL, 129}; comp_vl_datatype val_rd_c; void *buf = NULL; unsigned odd; /* Whether an odd or even coord. was read */ unsigned i, j; /* Local index variables */ /* Set vl datatype init value strings */ init_val_c.a = HDstrdup("baz"); init_val_c.b = HDstrdup("mumble"); /* Make copy of dataset creation property list */ if ((dcpl = H5Pcopy(_dcpl)) < 0) TEST_ERROR; #ifndef NO_FILLING if (H5Pset_fill_value(dcpl, dtype, fillval) < 0) TEST_ERROR; #endif /* Get class of datatype */ if ((dtype_class = H5Tget_class(dtype)) < 0) TEST_ERROR; /* Create a dataspace */ if ((fspace = H5Screate_simple(5, start_size, max_size)) < 0) TEST_ERROR; /* Create dataset */ if ((dset = H5Dcreate2(file, dset_name, dtype, fspace, H5P_DEFAULT, dcpl, H5P_DEFAULT)) < 0) TEST_ERROR; /* Set up pointers to elements */ if (dtype_class == H5T_INTEGER) { /* Initialize specific values for this datatype */ val_size = sizeof(int); init_val = &init_val_i; init_rtn = test_extend_init_integer; verify_rtn = test_extend_verify_integer; release_rtn = test_extend_release_integer; val_rd = &val_rd_i; odd_val = fillval; even_val = &init_val_i; } /* end if */ else { /* Sanity check */ assert(dtype_class == H5T_COMPOUND); /* Initialize specific values for this datatype */ val_size = sizeof(comp_vl_datatype); init_val = &init_val_c; init_rtn = test_extend_init_cmpd_vl; verify_rtn = test_extend_verify_cmpd_vl; release_rtn = test_extend_release_cmpd_vl; val_rd = &val_rd_c; odd_val = fillval; even_val = &init_val_c; } /* end else */ /* Read some data and make sure it's the fill value */ if ((mspace = H5Screate_simple(5, one, NULL)) < 0) TEST_ERROR; for (i = 0; i < 1000; i++) { /* Set offset for random element */ for (j = 0; j < 5; j++) hs_offset[j] = (hsize_t)HDrand() % start_size[j]; /* Select the random element */ if (H5Sselect_hyperslab(fspace, H5S_SELECT_SET, hs_offset, NULL, one, NULL) < 0) TEST_ERROR; /* Read the random element */ if (H5Dread(dset, dtype, mspace, fspace, H5P_DEFAULT, val_rd) < 0) TEST_ERROR; /* Verify the element read in */ if (verify_rtn((unsigned)__LINE__, hs_offset, val_rd, fillval) < 0) TEST_ERROR; /* Release any VL components */ if (H5Treclaim(dtype, mspace, H5P_DEFAULT, val_rd) < 0) TEST_ERROR; /* Clear the read buffer */ HDmemset(val_rd, 0, val_size); } /* end for */ if (H5Sclose(mspace) < 0) TEST_ERROR; /* Initialize dataspace & hyperslab info */ for (i = 0, nelmts = 1; i < 5; i++) { hs_size[i] = (start_size[i] + 1) / 2; hs_offset[i] = 0; hs_stride[i] = 2; nelmts *= hs_size[i]; } /* end for */ /* Check for overflow */ assert((nelmts * val_size) == (hsize_t)((size_t)(nelmts * val_size))); /* Allocate & initialize buffer */ buf = HDmalloc((size_t)(nelmts * val_size)); init_rtn(buf, nelmts, init_val); /* Create dataspace describing memory buffer */ if ((mspace = H5Screate_simple(5, hs_size, hs_size)) < 0) TEST_ERROR; /* Select elements within file dataspace */ if (H5Sselect_hyperslab(fspace, H5S_SELECT_SET, hs_offset, hs_stride, hs_size, NULL) < 0) TEST_ERROR; /* Write to all even data locations */ if (H5Dwrite(dset, dtype, mspace, fspace, H5P_DEFAULT, buf) < 0) TEST_ERROR; /* Close memory dataspace */ if (H5Sclose(mspace) < 0) TEST_ERROR; /* Read some data and make sure it's the right value */ if ((mspace = H5Screate_simple(5, one, NULL)) < 0) TEST_ERROR; for (i = 0; i < 1000; i++) { /* Set offset for random element */ for (j = 0, odd = 0; j < 5; j++) { hs_offset[j] = (hsize_t)HDrand() % start_size[j]; odd += (unsigned)(hs_offset[j] % 2); } /* end for */ /* Select the random element */ if (H5Sselect_hyperslab(fspace, H5S_SELECT_SET, hs_offset, NULL, one, NULL) < 0) TEST_ERROR; /* Read the random element */ if (H5Dread(dset, dtype, mspace, fspace, H5P_DEFAULT, val_rd) < 0) TEST_ERROR; /* Verify the element read in */ should_be = odd ? odd_val : even_val; if (verify_rtn((unsigned)__LINE__, hs_offset, val_rd, should_be) < 0) TEST_ERROR; /* Release any VL components */ if (H5Treclaim(dtype, mspace, H5P_DEFAULT, val_rd) < 0) TEST_ERROR; /* Clear the read buffer */ HDmemset(val_rd, 0, val_size); } /* end for */ if (H5Sclose(mspace) < 0) TEST_ERROR; /* Extend the dataset one element in each dimension */ for (i = 0; i < 5; i++) extend_size[i] = start_size[i] + 1; if (H5Dset_extent(dset, extend_size) < 0) TEST_ERROR; /* Re-open file dataspace */ if (H5Sclose(fspace) < 0) TEST_ERROR; if ((fspace = H5Dget_space(dset)) < 0) TEST_ERROR; /* Read some data and make sure it's the right value */ if ((mspace = H5Screate_simple(5, one, NULL)) < 0) TEST_ERROR; for (i = 0; i < 1000; i++) { /* Set offset for random element */ for (j = 0, odd = 0; j < 5; j++) { hs_offset[j] = (hsize_t)HDrand() % extend_size[j]; if (hs_offset[j] >= start_size[j]) odd = 1; else odd += (unsigned)(hs_offset[j] % 2); } /* end for */ /* Select the random element */ if (H5Sselect_hyperslab(fspace, H5S_SELECT_SET, hs_offset, NULL, one, NULL) < 0) TEST_ERROR; /* Read the random element */ if (H5Dread(dset, dtype, mspace, fspace, H5P_DEFAULT, val_rd) < 0) TEST_ERROR; /* Verify the element read in */ should_be = odd ? odd_val : even_val; if (verify_rtn((unsigned)__LINE__, hs_offset, val_rd, should_be) < 0) TEST_ERROR; /* Release any VL components */ if (H5Treclaim(dtype, mspace, H5P_DEFAULT, val_rd) < 0) TEST_ERROR; /* Clear the read buffer */ HDmemset(val_rd, 0, val_size); } /* end for */ if (H5Sclose(mspace) < 0) TEST_ERROR; /* Extend the dataset to the maximum dimension sizes */ if (H5Dset_extent(dset, max_size) < 0) TEST_ERROR; /* Re-open file dataspace */ if (H5Sclose(fspace) < 0) TEST_ERROR; if ((fspace = H5Dget_space(dset)) < 0) TEST_ERROR; /* Read some data and make sure it's the right value */ if ((mspace = H5Screate_simple(5, one, NULL)) < 0) TEST_ERROR; for (i = 0; i < 1000; i++) { /* Set offset for random element */ for (j = 0, odd = 0; j < 5; j++) { hs_offset[j] = (hsize_t)HDrand() % max_size[j]; if (hs_offset[j] >= start_size[j]) odd = 1; else odd += (unsigned)(hs_offset[j] % 2); } /* end for */ /* Select the random element */ if (H5Sselect_hyperslab(fspace, H5S_SELECT_SET, hs_offset, NULL, one, NULL) < 0) TEST_ERROR; /* Read the random element */ if (H5Dread(dset, dtype, mspace, fspace, H5P_DEFAULT, val_rd) < 0) TEST_ERROR; /* Verify the element read in */ should_be = odd ? odd_val : even_val; if (verify_rtn((unsigned)__LINE__, hs_offset, val_rd, should_be) < 0) TEST_ERROR; /* Release any VL components */ if (H5Treclaim(dtype, mspace, H5P_DEFAULT, val_rd) < 0) TEST_ERROR; /* Clear the read buffer */ HDmemset(val_rd, 0, val_size); } /* end for */ if (H5Sclose(mspace) < 0) TEST_ERROR; /* Shrink the dataset to half of it's maximum size, plus 1/2 of a chunk */ for (i = 0; i < 5; i++) extend_size[i] = (max_size[i] / 2) + (ch_size[i] / 2); if (H5Dset_extent(dset, extend_size) < 0) TEST_ERROR; /* Re-open file dataspace */ if (H5Sclose(fspace) < 0) TEST_ERROR; if ((fspace = H5Dget_space(dset)) < 0) TEST_ERROR; /* Read some data and make sure it's the right value */ if ((mspace = H5Screate_simple(5, one, NULL)) < 0) TEST_ERROR; for (i = 0; i < 1000; i++) { /* Set offset for random element */ for (j = 0, odd = 0; j < 5; j++) { hs_offset[j] = (hsize_t)HDrand() % extend_size[j]; if (hs_offset[j] >= start_size[j]) odd = 1; else odd += (unsigned)(hs_offset[j] % 2); } /* end for */ /* Select the random element */ if (H5Sselect_hyperslab(fspace, H5S_SELECT_SET, hs_offset, NULL, one, NULL) < 0) TEST_ERROR; /* Read the random element */ if (H5Dread(dset, dtype, mspace, fspace, H5P_DEFAULT, val_rd) < 0) TEST_ERROR; /* Verify the element read in */ should_be = odd ? odd_val : even_val; if (verify_rtn((unsigned)__LINE__, hs_offset, val_rd, should_be) < 0) TEST_ERROR; /* Release any VL components */ if (H5Treclaim(dtype, mspace, H5P_DEFAULT, val_rd) < 0) TEST_ERROR; /* Clear the read buffer */ HDmemset(val_rd, 0, val_size); } /* end for */ if (H5Sclose(mspace) < 0) TEST_ERROR; /* Extend the dataset's size by one element, in a dimension that won't * cause additional chunks to be needed */ extend_size[2] += 1; if (H5Dset_extent(dset, extend_size) < 0) TEST_ERROR; /* Re-open file dataspace */ if (H5Sclose(fspace) < 0) TEST_ERROR; if ((fspace = H5Dget_space(dset)) < 0) TEST_ERROR; /* Create dataspace for single element sized buffer */ if ((mspace = H5Screate_simple(5, one, NULL)) < 0) TEST_ERROR; /* Set location for "top-most" element in dataset to write */ for (i = 0; i < 5; i++) hs_offset[i] = extend_size[i] - 1; /* Select the element */ if (H5Sselect_hyperslab(fspace, H5S_SELECT_SET, hs_offset, NULL, one, NULL) < 0) TEST_ERROR; /* Write one element in a chunk that's 'partial' and overwrite a fill * value that was initialized in the H5Dset_extent() routine. This will * overwrite a fill-value element, which must be de-allocated properly or * next read of another fill-value initialized element in this chunk will * fail. */ if (H5Dwrite(dset, dtype, mspace, fspace, H5P_DEFAULT, buf) < 0) TEST_ERROR; /* Read value back in */ if (H5Dread(dset, dtype, mspace, fspace, H5P_DEFAULT, val_rd) < 0) TEST_ERROR; /* Verify the element read in is the value written out */ if (verify_rtn((unsigned)__LINE__, hs_offset, val_rd, buf) < 0) TEST_ERROR; /* Set the element back to fillval */ if (H5Dwrite(dset, dtype, mspace, fspace, H5P_DEFAULT, fillval) < 0) TEST_ERROR; /* Release any VL components */ if (H5Treclaim(dtype, mspace, H5P_DEFAULT, val_rd) < 0) TEST_ERROR; /* Clear the read buffer */ HDmemset(val_rd, 0, val_size); /* Set location for another element initialized by H5Dset_extent() */ hs_offset[3] -= 1; /* Select the element */ if (H5Sselect_hyperslab(fspace, H5S_SELECT_SET, hs_offset, NULL, one, NULL) < 0) TEST_ERROR; /* Read value back in */ if (H5Dread(dset, dtype, mspace, fspace, H5P_DEFAULT, val_rd) < 0) TEST_ERROR; /* Verify the element read in is the fill-value */ if (verify_rtn((unsigned)__LINE__, hs_offset, val_rd, fillval) < 0) TEST_ERROR; /* Release any VL components */ if (H5Treclaim(dtype, mspace, H5P_DEFAULT, val_rd) < 0) TEST_ERROR; /* Clear the read buffer */ HDmemset(val_rd, 0, val_size); /* Read some data and make sure it's the right value */ for (i = 0; i < 1000; i++) { /* Set offset for random element */ for (j = 0, odd = 0; j < 5; j++) { hs_offset[j] = (hsize_t)HDrand() % extend_size[j]; if (hs_offset[j] >= start_size[j]) odd = 1; else odd += (unsigned)(hs_offset[j] % 2); } /* end for */ /* Select the random element */ if (H5Sselect_hyperslab(fspace, H5S_SELECT_SET, hs_offset, NULL, one, NULL) < 0) TEST_ERROR; /* Read the random element */ if (H5Dread(dset, dtype, mspace, fspace, H5P_DEFAULT, val_rd) < 0) TEST_ERROR; /* Verify the element read in */ should_be = odd ? odd_val : even_val; if (verify_rtn((unsigned)__LINE__, hs_offset, val_rd, should_be) < 0) TEST_ERROR; /* Release any VL components */ if (H5Treclaim(dtype, mspace, H5P_DEFAULT, val_rd) < 0) TEST_ERROR; /* Clear the read buffer */ HDmemset(val_rd, 0, val_size); } /* end for */ if (H5Sclose(mspace) < 0) TEST_ERROR; /* Release elements & memory buffer */ for (i = 0; i < nelmts; i++) release_rtn((void *)((char *)buf + (val_size * i))); HDfree(init_val_c.a); HDfree(init_val_c.b); HDfree(buf); /* Cleanup IDs */ if (H5Pclose(dcpl) < 0) TEST_ERROR; if (H5Dclose(dset) < 0) TEST_ERROR; if (H5Sclose(fspace) < 0) TEST_ERROR; return 0; error: HDfree(init_val_c.a); HDfree(init_val_c.b); HDfree(buf); H5E_BEGIN_TRY { H5Pclose(dcpl); H5Dclose(dset); H5Sclose(fspace); H5Sclose(mspace); } H5E_END_TRY return -1; } /* end test_extend_cases() */ /*------------------------------------------------------------------------- * Function: test_extend * * Purpose: Test that filling works okay when a dataset is extended. * * Return: Success: 0 * * Failure: number of errors * * Programmer: Robb Matzke * Monday, October 5, 1998 * *------------------------------------------------------------------------- */ static int test_extend(hid_t fapl, const char *base_name, H5D_layout_t layout) { hid_t file = -1; /* File ID */ hid_t dcpl = -1; /* Dataset creation property list ID */ hid_t cmpd_vl_tid = -1; /* Compound+vl datatype ID */ hsize_t start_size[5] = {8, 8, 8, 4, 2}; hsize_t max_size[5] = {32, 32, 32, 16, 8}; hsize_t ch_size[5] = {1, 8, 8, 4, 2}; #ifdef NO_FILLING int fillval_i = 0; #else int fillval_i = 0x4c70f1cd; #endif comp_vl_datatype fillval_c = {32, NULL, NULL, 64}; /* Fill value for compound+vl datatype tests */ char filename[1024]; /* Print testing message */ if (H5D_CHUNKED == layout) TESTING("chunked dataset extend"); else TESTING("contiguous dataset extend"); /* Set vl datatype fill value strings */ fillval_c.a = HDstrdup("foo"); fillval_c.b = HDstrdup("bar"); /* Create dataset creation property list */ if ((dcpl = H5Pcreate(H5P_DATASET_CREATE)) < 0) TEST_ERROR; if (H5D_CHUNKED == layout) if (H5Pset_chunk(dcpl, 5, ch_size) < 0) TEST_ERROR; #if 1 /* * Remove this once contiguous datasets can support extensions in other * than the slowest varying dimension. The purpose of this block is to * make only the slowest varying dimension extendible and yet have the * same total number of elements as originally. * * If this is removed prematurely then you will get an error `only the * first dimension can be extendible' as long as the test isn't skipped * below. */ if (H5D_CONTIGUOUS == layout) { max_size[0] = (max_size[0] * max_size[1] * max_size[2] * max_size[3] * max_size[4]) / (start_size[1] * start_size[2] * start_size[3] * start_size[4]); max_size[1] = start_size[1]; max_size[2] = start_size[2]; max_size[3] = start_size[3]; max_size[4] = start_size[4]; } #endif #if 1 /* * Remove this once internal contiguous datasets can support * extending. If it's removed prematurely you will get an error * `extendible contiguous non-external dataset' as long as the test isn't * skipped below. */ if (H5D_CONTIGUOUS == layout) { int fd; hsize_t nelmts; nelmts = max_size[0] * max_size[1] * max_size[2] * max_size[3] * max_size[4]; if ((fd = HDopen(FILE_NAME_RAW, O_RDWR | O_CREAT | O_TRUNC, H5_POSIX_CREATE_MODE_RW)) < 0 || HDclose(fd) < 0) goto error; if (H5Pset_external(dcpl, FILE_NAME_RAW, (off_t)0, (hsize_t)nelmts * sizeof(int)) < 0) goto error; } #endif #if 1 /* * Remove this when contiguous datasets can be extended to some * predetermined fininte size, even if it's just in the slowest varying * dimension. If it's removed prematurely then you'll get one of the * errors described above or `unable to select fill value region'. */ if (H5D_CONTIGUOUS == layout) { SKIPPED(); HDputs(" Not implemented yet -- needs H5S_SELECT_DIFF operator"); goto skip; } #endif /* Create a file and dataset */ h5_fixname(base_name, fapl, filename, sizeof filename); if ((file = H5Fcreate(filename, H5F_ACC_TRUNC, H5P_DEFAULT, fapl)) < 0) TEST_ERROR; /* Get the compound+vl datatype */ if ((cmpd_vl_tid = create_compound_vl_type()) < 0) TEST_ERROR; /* Test integer datatype case */ if (test_extend_cases(file, dcpl, "dset1", ch_size, start_size, max_size, H5T_NATIVE_INT, &fillval_i) < 0) TEST_ERROR; /* Test compound+vl datatype datatype case */ if (test_extend_cases(file, dcpl, "dset2", ch_size, start_size, max_size, cmpd_vl_tid, &fillval_c) < 0) TEST_ERROR; /* Cleanup */ if (H5Tclose(cmpd_vl_tid) < 0) TEST_ERROR; if (H5Pclose(dcpl) < 0) TEST_ERROR; if (H5Fclose(file) < 0) TEST_ERROR; HDfree(fillval_c.a); HDfree(fillval_c.b); PASSED(); return 0; error: HDfree(fillval_c.a); HDfree(fillval_c.b); H5E_BEGIN_TRY { H5Tclose(cmpd_vl_tid); H5Pclose(dcpl); H5Fclose(file); } H5E_END_TRY return 1; skip: HDfree(fillval_c.a); HDfree(fillval_c.b); H5E_BEGIN_TRY { H5Pclose(dcpl); H5Fclose(file); } H5E_END_TRY return 0; } /* end test_extend() */ /*------------------------------------------------------------------------- * Function: test_compatible * * Purpose: Tests fill value and dataspace for datasets created by v1.4 * library. * * Return: Success: 0 * * Failure: number of errors * * Programmer: Raymond Lu * Feb 27, 2002 * *------------------------------------------------------------------------- */ static int test_compatible(void) { hid_t file = -1, dset1 = -1, dset2 = -1; hid_t dcpl1 = -1, dcpl2 = -1, fspace = -1, mspace = -1; int rd_fill = 0, fill_val = 4444, val_rd = 0; hsize_t dims[2], one[2] = {1, 1}; hsize_t hs_offset[2] = {3, 4}; H5D_fill_value_t status; const char *testfile = H5_get_srcdir_filename(FILE_COMPATIBLE); /* Corrected test file name */ TESTING("contiguous dataset compatibility with v. 1.4"); if ((file = H5Fopen(testfile, H5F_ACC_RDONLY, H5P_DEFAULT)) < 0) { HDprintf(" Could not open file %s. Try set $srcdir to point at the " "source directory of test\n", testfile); goto error; } if ((dset1 = H5Dopen2(file, "dset1", H5P_DEFAULT)) < 0) goto error; if ((dcpl1 = H5Dget_create_plist(dset1)) < 0) goto error; if (H5Pfill_value_defined(dcpl1, &status) < 0) goto error; if (status != H5D_FILL_VALUE_UNDEFINED) { H5_FAILED(); HDprintf(" %d: Got a different fill value than what was set.", __LINE__); HDprintf(" Got status=%ld, suppose to be H5D_FILL_VALUE_UNDEFINED\n", (long)status); goto error; } if ((fspace = H5Dget_space(dset1)) < 0) goto error; if (H5Sget_simple_extent_dims(fspace, dims, NULL) < 0) goto error; if (dims[0] != 8 || dims[1] != 8) { H5_FAILED(); HDputs(" Got a different dimension size than what was set."); HDprintf(" Got dims[0]=%ld, dims[1]=%ld, set 8x8\n", (long)dims[0], (long)dims[1]); goto error; } if ((mspace = H5Screate_simple(2, one, NULL)) < 0) goto error; if (H5Sselect_hyperslab(fspace, H5S_SELECT_SET, hs_offset, NULL, one, NULL) < 0) goto error; if (H5Dread(dset1, H5T_NATIVE_INT, mspace, fspace, H5P_DEFAULT, &val_rd) < 0) goto error; if (val_rd != 0) { H5_FAILED(); HDputs(" Got a different value than what was set."); HDprintf(" Got %ld, set 0\n", (long)val_rd); goto error; } if (H5Pclose(dcpl1) < 0) goto error; if (H5Sclose(fspace) < 0) goto error; if (H5Sclose(mspace) < 0) goto error; if (H5Dclose(dset1) < 0) goto error; if ((dset2 = H5Dopen2(file, "dset2", H5P_DEFAULT)) < 0) goto error; if ((dcpl2 = H5Dget_create_plist(dset2)) < 0) goto error; if (H5Pfill_value_defined(dcpl2, &status) < 0) goto error; if (status != H5D_FILL_VALUE_USER_DEFINED) { H5_FAILED(); HDprintf(" %d: Got a different fill value than what was set.", __LINE__); HDprintf(" Got status=%ld, suppose to be H5D_FILL_VALUE_USER_DEFINED\n", (long)status); goto error; } if (H5Pget_fill_value(dcpl2, H5T_NATIVE_INT, &rd_fill) < 0) goto error; if (rd_fill != fill_val) { H5_FAILED(); HDprintf(" %d: Got a different fill value than what was set.", __LINE__); HDprintf(" Got %ld, set %ld\n", (long)rd_fill, (long)fill_val); goto error; } fspace = -1; if ((fspace = H5Dget_space(dset2)) < 0) goto error; dims[0] = dims[1] = (hsize_t)-1; if (H5Sget_simple_extent_dims(fspace, dims, NULL) < 0) goto error; if (dims[0] != 8 || dims[1] != 8) { H5_FAILED(); HDputs(" Got a different dimension size than what was set."); HDprintf(" Got dims[0]=%ld, dims[1]=%ld, set 8x8\n", (long)dims[0], (long)dims[1]); goto error; } if ((mspace = H5Screate_simple(2, one, NULL)) < 0) goto error; if (H5Sselect_hyperslab(fspace, H5S_SELECT_SET, hs_offset, NULL, one, NULL) < 0) goto error; if (H5Dread(dset2, H5T_NATIVE_INT, mspace, fspace, H5P_DEFAULT, &val_rd) < 0) goto error; if (val_rd != fill_val) { H5_FAILED(); HDputs(" Got a different value than what was set."); HDprintf(" Got %ld, set %ld\n", (long)val_rd, (long)fill_val); goto error; } if (H5Pclose(dcpl2) < 0) goto error; if (H5Sclose(fspace) < 0) goto error; if (H5Sclose(mspace) < 0) goto error; if (H5Dclose(dset2) < 0) goto error; if (H5Fclose(file) < 0) goto error; PASSED(); return 0; error: H5E_BEGIN_TRY { H5Pclose(dcpl1); H5Sclose(fspace); H5Sclose(mspace); H5Dclose(dset1); H5Pclose(dcpl2); H5Sclose(fspace); H5Dclose(dset2); H5Fclose(file); } H5E_END_TRY return 1; } /*------------------------------------------------------------------------- * Function: test_partalloc_cases * * Purpose: Tests fill values read and write for datasets. * * Return: Success: 0 * * Failure: 1 * * Programmer: Joel Plutchak * April 15, 2013 * *------------------------------------------------------------------------- */ static int test_partalloc_cases(hid_t file, hid_t dcpl, const char *dname, H5D_fill_time_t fill_time) { hid_t fspace = -1, dset1 = -1, rspace = -1; herr_t ret; hsize_t ds_size[2] = {4, 4}; hsize_t max_size[2] = {H5S_UNLIMITED, 4}; hsize_t chunk_size[2] = {1, 4}; int fillval = (-1); int w_values[] = {42}; /* New value to be written */ int f_values[4] = {88, 88, 88, 88}; /* pre-seed read buffer with known values */ int r_values[4] = {88, 88, 88, 88}; /* pre-seed read buffer with known values */ hsize_t coord[1][2]; /* coordinate(s) of point to write */ hsize_t start[2], count[2]; fillval = 0; /* default fill value is zero */ /* Create dataset with 4x4 integer dataset */ if ((fspace = H5Screate_simple(2, ds_size, max_size)) < 0) goto error; if ((dset1 = H5Dcreate2(file, dname, H5T_NATIVE_INT, fspace, H5P_DEFAULT, dcpl, H5P_DEFAULT)) < 0) goto error; /* * Select a point in the file dataspace. */ coord[0][0] = 0; coord[0][1] = 0; if (H5Sselect_elements(fspace, H5S_SELECT_SET, (size_t)1, (const hsize_t *)coord)) goto error; /* * Write single data point to the dataset. */ if ((ret = H5Dwrite(dset1, H5T_NATIVE_INT, H5S_ALL, fspace, H5P_DEFAULT, w_values)) < 0) { goto error; } /* Read a line/chunk and make sure values are right */ rspace = H5Screate_simple(2, chunk_size, NULL); /* Read the first row of elements: one known plus three fill */ start[0] = 0; start[1] = 0; count[0] = 1; count[1] = 4; if ((ret = H5Sselect_hyperslab(fspace, H5S_SELECT_SET, start, NULL, count, NULL)) < 0) goto error; if ((ret = H5Sselect_all(rspace)) < 0) goto error; if (H5Dread(dset1, H5T_NATIVE_INT, rspace, fspace, H5P_DEFAULT, &r_values) < 0) goto error; /* Read the third row of elements: all fill */ start[0] = 2; start[1] = 0; count[0] = 1; count[1] = 4; if ((ret = H5Sselect_hyperslab(fspace, H5S_SELECT_SET, start, NULL, count, NULL)) < 0) goto error; if (H5Dread(dset1, H5T_NATIVE_INT, rspace, fspace, H5P_DEFAULT, &f_values) < 0) goto error; if (fill_time != H5D_FILL_TIME_NEVER) { /* check allocated chunk */ if ((r_values[0] != w_values[0]) || (r_values[1] != fillval) || (r_values[2] != fillval) || (r_values[3] != fillval)) { H5_FAILED(); HDfprintf(stdout, "%u: Allocated chunk value read was not correct.\n", (unsigned)__LINE__); HDprintf(" {%ld,%ld,%ld,%ld} should be {%ld,%ld,%ld,%ld}\n", (long)r_values[0], (long)r_values[1], (long)r_values[2], (long)r_values[3], (long)w_values[0], (long)fillval, (long)fillval, (long)fillval); goto error; } /* check unallocated chunk */ if ((f_values[0] != fillval) || (f_values[1] != fillval) || (f_values[2] != fillval) || (f_values[3] != fillval)) { H5_FAILED(); HDfprintf(stdout, "%u: Unallocated chunk value read was not correct.\n", (unsigned)__LINE__); HDprintf(" {%ld,%ld,%ld,%ld} should be {%ld,%ld,%ld,%ld}\n", (long)f_values[0], (long)f_values[1], (long)f_values[2], (long)f_values[3], (long)fillval, (long)fillval, (long)fillval, (long)fillval); goto error; } /* for the "never fill" case expect to get trash values, so skip */ } else if (fill_time == H5D_FILL_TIME_NEVER) { } if (H5Sclose(rspace) < 0) goto error; if (H5Dclose(dset1) < 0) goto error; if (H5Sclose(fspace) < 0) goto error; return 0; error: H5E_BEGIN_TRY { H5Dclose(dset1); H5Sclose(fspace); H5Sclose(rspace); } H5E_END_TRY return 1; } /*------------------------------------------------------------------------- * Function: test_partalloc * * Purpose: Tests fill values for chunked, partially-allocated datasets. * Regression test for HDFFV-8247. * * Return: Success: 0 * * Failure: number of errors * * Programmer: Joel Plutchak * April 15, 2013 * *------------------------------------------------------------------------- */ static int test_partalloc(hid_t fapl, const char *base_name) { char filename[1024]; hid_t file = -1, dcpl = -1; hsize_t ch_size[2] = {1, 4}; int nerrors = 0; TESTING("chunked dataset partially allocated I/O"); h5_fixname(base_name, fapl, filename, sizeof filename); if ((file = H5Fcreate(filename, H5F_ACC_TRUNC, H5P_DEFAULT, fapl)) < 0) goto error; if ((dcpl = H5Pcreate(H5P_DATASET_CREATE)) < 0) goto error; if (H5Pset_chunk(dcpl, 2, ch_size) < 0) goto error; /* I. Test H5D_ALLOC_TIME_LATE space allocation cases */ if (H5Pset_alloc_time(dcpl, H5D_ALLOC_TIME_LATE) < 0) goto error; #ifdef DEBUG HDfprintf(stdout, "\nALLOC_TIME_LATE\n"); #endif /* case for H5D_FILL_TIME_ALLOC as fill write time and fill value to be default */ if (H5Pset_fill_time(dcpl, H5D_FILL_TIME_ALLOC) < 0) goto error; #ifdef DEBUG HDfprintf(stdout, " FILL_TIME_ALLOC\n"); #endif nerrors += test_partalloc_cases(file, dcpl, "dset1", H5D_FILL_TIME_ALLOC); /* case for H5D_FILL_TIME_NEVER as fill write time and fill value to be default */ if (H5Pset_fill_time(dcpl, H5D_FILL_TIME_NEVER) < 0) goto error; #ifdef DEBUG HDfprintf(stdout, " FILL_TIME_NEVER\n"); #endif nerrors += test_partalloc_cases(file, dcpl, "dset2", H5D_FILL_TIME_NEVER); /* case for H5D_FILL_TIME_IFSET as fill write time and fill value to be default */ if (H5Pset_fill_time(dcpl, H5D_FILL_TIME_IFSET) < 0) goto error; #ifdef DEBUG HDfprintf(stdout, " FILL_TIME_IFSET\n"); #endif nerrors += test_partalloc_cases(file, dcpl, "dset3", H5D_FILL_TIME_IFSET); /* II. Test H5D_ALLOC_TIME_INCR space allocation cases */ if (H5Pset_alloc_time(dcpl, H5D_ALLOC_TIME_INCR) < 0) goto error; #ifdef DEBUG HDfprintf(stdout, "\nALLOC_TIME_INCR\n"); #endif /* case for H5D_FILL_TIME_ALLOC as fill write time and fill value to be default */ if (H5Pset_fill_time(dcpl, H5D_FILL_TIME_ALLOC) < 0) goto error; #ifdef DEBUG HDfprintf(stdout, " FILL_TIME_ALLOC\n"); #endif nerrors += test_partalloc_cases(file, dcpl, "dset4", H5D_FILL_TIME_ALLOC); /* case for H5D_FILL_TIME_NEVER as fill write time and fill value to be default */ if (H5Pset_fill_time(dcpl, H5D_FILL_TIME_NEVER) < 0) goto error; #ifdef DEBUG HDfprintf(stdout, " FILL_TIME_NEVER\n"); #endif nerrors += test_partalloc_cases(file, dcpl, "dset5", H5D_FILL_TIME_NEVER); /* case for H5D_FILL_TIME_IFSET as fill write time and fill value to be default */ if (H5Pset_fill_time(dcpl, H5D_FILL_TIME_IFSET) < 0) goto error; #ifdef DEBUG HDfprintf(stdout, " FILL_TIME_IFSET\n"); #endif nerrors += test_partalloc_cases(file, dcpl, "dset6", H5D_FILL_TIME_IFSET); /* III. Test H5D_ALLOC_TIME_EARLY space allocation cases */ if (H5Pset_alloc_time(dcpl, H5D_ALLOC_TIME_EARLY) < 0) goto error; #ifdef DEBUG HDfprintf(stdout, "\nALLOC_TIME_EARLY\n"); #endif /* case for H5D_FILL_TIME_ALLOC as fill write time and fill value to be default */ if (H5Pset_fill_time(dcpl, H5D_FILL_TIME_ALLOC) < 0) goto error; #ifdef DEBUG HDfprintf(stdout, " FILL_TIME_ALLOC\n"); #endif nerrors += test_partalloc_cases(file, dcpl, "dset7", H5D_FILL_TIME_ALLOC); /* case for H5D_FILL_TIME_NEVER as fill write time and fill value to be default */ if (H5Pset_fill_time(dcpl, H5D_FILL_TIME_NEVER) < 0) goto error; #ifdef DEBUG HDfprintf(stdout, " FILL_TIME_NEVER\n"); #endif nerrors += test_partalloc_cases(file, dcpl, "dset8", H5D_FILL_TIME_NEVER); /* case for H5D_FILL_TIME_IFSET as fill write time and fill value to be default */ if (H5Pset_fill_time(dcpl, H5D_FILL_TIME_IFSET) < 0) goto error; #ifdef DEBUG HDfprintf(stdout, " FILL_TIME_IFSET\n"); #endif nerrors += test_partalloc_cases(file, dcpl, "dset9", H5D_FILL_TIME_IFSET); if (nerrors) goto error; if (H5Pclose(dcpl) < 0) goto error; if (H5Fclose(file) < 0) goto error; PASSED(); return 0; error: H5E_BEGIN_TRY { H5Pclose(dcpl); H5Fclose(file); } H5E_END_TRY return nerrors; } /*------------------------------------------------------------------------- * Function: main * * Purpose: Tests fill values * * Return: Success: * * Failure: * * Programmer: Robb Matzke * Thursday, October 1, 1998 * *------------------------------------------------------------------------- */ int main(int argc, char *argv[]) { int nerrors = 0, argno, test_contig = 1, test_chunk = 1, test_compact = 1; hid_t fapl = (-1), fapl2 = (-1); /* File access property lists */ unsigned new_format; /* Whether to use the new format or not */ hbool_t driver_is_default_compatible; if (argc >= 2) { test_contig = test_chunk = test_compact = 0; for (argno = 1; argno < argc; argno++) { if (!HDstrcmp(argv[argno], "contiguous")) test_contig = 1; else if (!HDstrcmp(argv[argno], "chunked")) test_chunk = 1; else if (!HDstrcmp(argv[argno], "compact")) test_compact = 1; else { HDfprintf(stderr, "usage: %s [contiguous] [chunked] [compact]\n", argv[0]); HDexit(EXIT_FAILURE); } } /* end for */ } /* end if */ h5_reset(); fapl = h5_fileaccess(); if (h5_driver_is_default_vfd_compatible(fapl, &driver_is_default_compatible) < 0) TEST_ERROR; /* Property list tests */ nerrors += test_getset(); nerrors += test_getset_vl(fapl); /* Copy the file access property list */ if ((fapl2 = H5Pcopy(fapl)) < 0) TEST_ERROR; /* Set the "use the latest version of the format" bounds for creating objects in the file */ if (H5Pset_libver_bounds(fapl2, H5F_LIBVER_LATEST, H5F_LIBVER_LATEST) < 0) TEST_ERROR; /* Loop over using new group format */ for (new_format = FALSE; new_format <= TRUE; new_format++) { hid_t my_fapl; /* Set the FAPL for the type of format */ if (new_format) { HDputs("\nTesting with new file format:"); my_fapl = fapl2; } /* end if */ else { HDputs("Testing with old file format:"); my_fapl = fapl; } /* end else */ /* Chunked storage layout tests */ if (test_chunk) { nerrors += test_create(my_fapl, FILENAME[0], H5D_CHUNKED); nerrors += test_rdwr(my_fapl, FILENAME[2], H5D_CHUNKED); nerrors += test_extend(my_fapl, FILENAME[4], H5D_CHUNKED); nerrors += test_partalloc(my_fapl, FILENAME[8]); } /* end if */ /* Contiguous storage layout tests */ if (test_contig) { nerrors += test_create(my_fapl, FILENAME[1], H5D_CONTIGUOUS); nerrors += test_rdwr(my_fapl, FILENAME[3], H5D_CONTIGUOUS); nerrors += test_extend(my_fapl, FILENAME[5], H5D_CONTIGUOUS); if (driver_is_default_compatible) { nerrors += test_compatible(); } } /* end if */ /* Compact dataset storage tests */ if (test_compact) { nerrors += test_create(my_fapl, FILENAME[6], H5D_COMPACT); nerrors += test_rdwr(my_fapl, FILENAME[7], H5D_COMPACT); } /* end if */ } /* end for */ /* Close 2nd FAPL */ H5Pclose(fapl2); /* Verify symbol table messages are cached */ nerrors += (h5_verify_cached_stabs(FILENAME, fapl) < 0 ? 1 : 0); if (nerrors) goto error; HDputs("All fill value tests passed."); if (h5_cleanup(FILENAME, fapl)) HDremove(FILE_NAME_RAW); HDexit(EXIT_SUCCESS); error: HDputs("***** FILL VALUE TESTS FAILED *****"); HDexit(EXIT_FAILURE); }