/* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 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: Raymond Lu * October 14, 2001 * * Purpose: Tests the H5Tget_native_type function. */ #include "h5test.h" const char *FILENAME[] = {"ntypes", NULL}; #define DIM0 100 #define DIM1 200 #define DIM3 20 #define DSET_ATOMIC_NAME_1 "atomic_type_1" #define DSET_ATOMIC_NAME_2 "atomic_type_2" #define DSET_ATOMIC_NAME_3 "atomic_type_3" #define DSET_ATOMIC_NAME_4 "atomic_type_4" #define DSET_ATOMIC_NAME_5 "atomic_type_5" #define DSET_COMPOUND_NAME "compound_type" #define DSET_COMPOUND_NAME_2 "compound_type_2" #define DSET_COMPOUND_NAME_3 "compound_type_3" #define DSET_COMPOUND_NAME_4 "compound_type_4" #define DSET_ENUM_NAME "enum_type" #define DSET_ARRAY_NAME "array_type" #define DSET_ARRAY2_NAME "array_type_2" #define DSET_VL_NAME "vl_type" #define DSET_VLSTR_NAME "vlstr_type" #define DSET_STR_NAME "str_type" #define DSET_OPAQUE_NAME "opaque_type" #define DSET1_BITFIELD_NAME "bitfield_type_1" #define DSET2_BITFIELD_NAME "bitfield_type_2" #define SPACE1_DIM1 4 #define SPACE1_RANK 1 #define SPACE2_RANK 2 #define SPACE2_DIM1 10 #define SPACE2_DIM2 10 #define BITFIELD_ENUMB 8 /*------------------------------------------------------------------------- * Function: test_atomic_dtype * * Purpose: Test H5Tget_native_type for atomic datatype * * Return: Success: 0 * * Failure: -1 * * Programmer: Raymond Lu * October 15, 2002 * * Modifications: * *------------------------------------------------------------------------- */ static herr_t test_atomic_dtype(hid_t file) { hid_t dataset = -1, space = -1; hid_t dtype = -1, native_type = -1; int i, j, n; hsize_t dims[2]; void *tmp = NULL; struct { int arr[DIM0][DIM1]; } *ipoints2 = NULL; struct { int arr[DIM0][DIM1]; } *icheck2 = NULL; TESTING("atomic datatype"); if (NULL == (ipoints2 = HDcalloc(1, sizeof(*ipoints2)))) TEST_ERROR; if (NULL == (icheck2 = HDcalloc(1, sizeof(*icheck2)))) TEST_ERROR; /* Initialize the dataset */ for (i = n = 0; i < DIM0; i++) for (j = 0; j < DIM1; j++) ipoints2->arr[i][j] = n++; /* Create the data space */ dims[0] = DIM0; dims[1] = DIM1; if ((space = H5Screate_simple(2, dims, NULL)) < 0) TEST_ERROR; /*------------------- Test data values ------------------------*/ /* Create the dataset */ if ((dataset = H5Dcreate2(file, DSET_ATOMIC_NAME_1, H5T_STD_I32BE, space, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) TEST_ERROR; /* Write the data to the dataset */ if (H5Dwrite(dataset, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, ipoints2) < 0) TEST_ERROR; /* Close dataset */ if (H5Dclose(dataset) < 0) TEST_ERROR; /* Open dataset again to check H5Tget_native_type */ if ((dataset = H5Dopen2(file, DSET_ATOMIC_NAME_1, H5P_DEFAULT)) < 0) TEST_ERROR; if ((dtype = H5Dget_type(dataset)) < 0) TEST_ERROR; if ((native_type = H5Tget_native_type(dtype, H5T_DIR_DEFAULT)) < 0) TEST_ERROR; /* Verify the datatype retrieved and converted */ if (H5Tget_order(native_type) != H5Tget_order(H5T_NATIVE_INT)) TEST_ERROR; if (H5Tget_size(native_type) < H5Tget_size(H5T_STD_I32BE)) TEST_ERROR; if (H5T_INTEGER != H5Tget_class(native_type)) TEST_ERROR; /* Read the dataset back. The temporary buffer is for special platforms * like Cray. */ if (NULL == (tmp = HDmalloc((size_t)(DIM0 * DIM1 * H5Tget_size(native_type))))) TEST_ERROR; if (H5Dread(dataset, native_type, H5S_ALL, H5S_ALL, H5P_DEFAULT, tmp) < 0) TEST_ERROR; /* Copy data from temporary buffer to destination buffer */ HDmemcpy(icheck2, tmp, (size_t)(DIM0 * DIM1 * H5Tget_size(native_type))); HDfree(tmp); tmp = NULL; /* Convert to the integer type */ if (H5Tconvert(native_type, H5T_NATIVE_INT, (DIM0 * DIM1), icheck2, NULL, H5P_DEFAULT) < 0) TEST_ERROR; /* Check that the values read are the same as the values written */ for (i = 0; i < DIM0; i++) for (j = 0; j < DIM1; j++) if (ipoints2->arr[i][j] != icheck2->arr[i][j]) { H5_FAILED(); HDprintf(" Read different values than written.\n"); HDprintf(" At index %d,%d\n", i, j); goto error; } /* end if */ if (H5Dclose(dataset) < 0) TEST_ERROR; if (H5Tclose(native_type) < 0) TEST_ERROR; if (H5Tclose(dtype) < 0) TEST_ERROR; /*------------------ Test different data types ----------------*/ /* Create the dataset of H5T_STD_I64LE */ if ((dataset = H5Dcreate2(file, DSET_ATOMIC_NAME_2, H5T_STD_I64LE, space, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) TEST_ERROR; if ((dtype = H5Dget_type(dataset)) < 0) TEST_ERROR; if ((native_type = H5Tget_native_type(dtype, H5T_DIR_DEFAULT)) < 0) TEST_ERROR; /* Verify the datatype retrieved and converted */ if (H5Tget_order(native_type) != H5Tget_order(H5T_NATIVE_LLONG)) TEST_ERROR; if (H5Tget_size(native_type) < H5Tget_size(H5T_STD_I64LE)) TEST_ERROR; if (H5T_INTEGER != H5Tget_class(native_type)) TEST_ERROR; if (H5Dclose(dataset) < 0) TEST_ERROR; if (H5Tclose(native_type) < 0) TEST_ERROR; if (H5Tclose(dtype) < 0) TEST_ERROR; /* Create the dataset of H5T_STD_I8LE */ if ((dataset = H5Dcreate2(file, DSET_ATOMIC_NAME_3, H5T_STD_I8LE, space, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) TEST_ERROR; if ((dtype = H5Dget_type(dataset)) < 0) TEST_ERROR; if ((native_type = H5Tget_native_type(dtype, H5T_DIR_ASCEND)) < 0) TEST_ERROR; /* Verify the datatype retrieved and converted */ if (H5Tget_order(native_type) != H5Tget_order(H5T_NATIVE_CHAR)) TEST_ERROR; if (H5Tget_size(native_type) < H5Tget_size(H5T_STD_I8LE)) TEST_ERROR; if (H5T_INTEGER != H5Tget_class(native_type)) TEST_ERROR; if (H5Dclose(dataset) < 0) TEST_ERROR; if (H5Tclose(native_type) < 0) TEST_ERROR; if (H5Tclose(dtype) < 0) TEST_ERROR; /* Create the dataset of H5T_IEEE_F32BE */ if ((dataset = H5Dcreate2(file, DSET_ATOMIC_NAME_4, H5T_IEEE_F32BE, space, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) TEST_ERROR; if ((dtype = H5Dget_type(dataset)) < 0) TEST_ERROR; if ((native_type = H5Tget_native_type(dtype, H5T_DIR_DESCEND)) < 0) TEST_ERROR; /* Verify the datatype retrieved and converted */ if (H5Tget_order(native_type) != H5Tget_order(H5T_NATIVE_FLOAT)) TEST_ERROR; if (H5Tget_size(native_type) < H5Tget_size(H5T_IEEE_F32BE)) TEST_ERROR; if (H5T_FLOAT != H5Tget_class(native_type)) TEST_ERROR; if (H5Dclose(dataset) < 0) TEST_ERROR; if (H5Tclose(native_type) < 0) TEST_ERROR; if (H5Tclose(dtype) < 0) TEST_ERROR; /* Create the dataset of H5T_IEEE_F64BE */ if ((dataset = H5Dcreate2(file, DSET_ATOMIC_NAME_5, H5T_IEEE_F64BE, space, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) TEST_ERROR; if ((dtype = H5Dget_type(dataset)) < 0) TEST_ERROR; if ((native_type = H5Tget_native_type(dtype, H5T_DIR_DESCEND)) < 0) TEST_ERROR; /* Verify the datatype retrieved and converted */ if (H5Tget_order(native_type) != H5Tget_order(H5T_NATIVE_DOUBLE)) TEST_ERROR; if (H5Tget_size(native_type) < H5Tget_size(H5T_IEEE_F64BE)) TEST_ERROR; if (H5T_FLOAT != H5Tget_class(native_type)) TEST_ERROR; if (H5Dclose(dataset) < 0) TEST_ERROR; if (H5Tclose(native_type) < 0) TEST_ERROR; if (H5Tclose(dtype) < 0) TEST_ERROR; /* Close dataspace */ if (H5Sclose(space) < 0) TEST_ERROR; HDfree(ipoints2); HDfree(icheck2); PASSED(); return 0; error: if (tmp) HDfree(tmp); H5E_BEGIN_TRY { H5Dclose(dataset); H5Tclose(native_type); H5Tclose(dtype); H5Sclose(space); } H5E_END_TRY; HDfree(ipoints2); HDfree(icheck2); return -1; } /*------------------------------------------------------------------------- * Function: test_compound_dtype2 * * Purpose: Test H5Tget_native_type for compound datatype * * Return: Success: 0 * * Failure: -1 * * Programmer: Raymond Lu * October 15, 2002 * * Modifications: * *------------------------------------------------------------------------- */ static herr_t test_compound_dtype2(hid_t file) { typedef struct s2 { short c2; long l2; long long ll2; } s2; typedef struct s1 { char c; int i; s2 st; unsigned long long l; } s1; hid_t dataset = -1, space = -1; hid_t dtype = -1, native_type = -1, tid = -1, tid2 = -1, tid_m = -1, tid_m2 = -1, mem_id = -1, nest_mem_id = -1; int i, j, n; hsize_t dims[2]; s1 *temp_point = NULL, *temp_check = NULL; s1 *points = NULL, *check = NULL; void *tmp = NULL, *bkg = NULL; TESTING("nested compound datatype"); /* Allocate space for the points & check arrays */ if (NULL == (points = (s1 *)HDmalloc(sizeof(s1) * DIM0 * DIM1))) TEST_ERROR; if (NULL == (check = (s1 *)HDcalloc(sizeof(s1), DIM0 * DIM1))) TEST_ERROR; /* Initialize the dataset */ for (i = n = 0, temp_point = points; i < DIM0; i++) { for (j = 0; j < DIM1; j++, temp_point++) { temp_point->c = 't'; temp_point->i = n++; temp_point->st.c2 = (short)(i + j); temp_point->st.l2 = (i * 5 + j * 50) * n; temp_point->st.ll2 = (i * 10 + j * 100) * n; temp_point->l = (unsigned long long)((i * 40 + j * 400) * n); } /* end for */ } /* end for */ /* Create the data space */ dims[0] = DIM0; dims[1] = DIM1; if ((space = H5Screate_simple(2, dims, NULL)) < 0) TEST_ERROR; /* Create compound datatype for disk storage */ if ((tid2 = H5Tcreate(H5T_COMPOUND, sizeof(s2))) < 0) TEST_ERROR; if ((tid = H5Tcreate(H5T_COMPOUND, sizeof(s1))) < 0) TEST_ERROR; /* Insert and pack members */ if (H5Tinsert(tid2, "c2", HOFFSET(s2, c2), H5T_STD_I16BE) < 0) TEST_ERROR; #if H5_SIZEOF_LONG == 4 if (H5Tinsert(tid2, "l2", HOFFSET(s2, l2), H5T_STD_I32LE) < 0) TEST_ERROR; #elif H5_SIZEOF_LONG == 8 if (H5Tinsert(tid2, "l2", HOFFSET(s2, l2), H5T_STD_I64LE) < 0) TEST_ERROR; #else #error "Unknown 'long' size" #endif #if H5_SIZEOF_LONG_LONG == 4 if (H5Tinsert(tid2, "ll2", HOFFSET(s2, ll2), H5T_STD_I32BE) < 0) TEST_ERROR; #elif H5_SIZEOF_LONG_LONG == 8 if (H5Tinsert(tid2, "ll2", HOFFSET(s2, ll2), H5T_STD_I64BE) < 0) TEST_ERROR; #else #error "Unknown 'long long' size" #endif if (H5Tinsert(tid, "c", HOFFSET(s1, c), H5T_STD_U8LE) < 0) TEST_ERROR; if (H5Tinsert(tid, "i", HOFFSET(s1, i), H5T_STD_I32LE) < 0) TEST_ERROR; if (H5Tinsert(tid, "st", HOFFSET(s1, st), tid2) < 0) TEST_ERROR; #if H5_SIZEOF_LONG_LONG == 4 if (H5Tinsert(tid, "l", HOFFSET(s1, l), H5T_STD_U32BE) < 0) TEST_ERROR; #elif H5_SIZEOF_LONG_LONG == 8 if (H5Tinsert(tid, "l", HOFFSET(s1, l), H5T_STD_U64BE) < 0) TEST_ERROR; #else #error "Unknown 'long long' size" #endif /* Take away the paddings */ if (H5Tpack(tid) < 0) TEST_ERROR; /* Create the dataset */ if ((dataset = H5Dcreate2(file, DSET_COMPOUND_NAME_2, tid, space, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) TEST_ERROR; /* Create compound datatype for memory */ if ((tid_m2 = H5Tcreate(H5T_COMPOUND, sizeof(s2))) < 0) TEST_ERROR; if ((tid_m = H5Tcreate(H5T_COMPOUND, sizeof(s1))) < 0) TEST_ERROR; /* Insert members */ if (H5Tinsert(tid_m2, "c2", HOFFSET(s2, c2), H5T_NATIVE_SHORT) < 0) TEST_ERROR; if (H5Tinsert(tid_m2, "l2", HOFFSET(s2, l2), H5T_NATIVE_LONG) < 0) TEST_ERROR; if (H5Tinsert(tid_m2, "ll2", HOFFSET(s2, ll2), H5T_NATIVE_LLONG) < 0) TEST_ERROR; if (H5Tinsert(tid_m, "c", HOFFSET(s1, c), H5T_NATIVE_UCHAR) < 0) TEST_ERROR; if (H5Tinsert(tid_m, "i", HOFFSET(s1, i), H5T_NATIVE_INT) < 0) TEST_ERROR; if (H5Tinsert(tid_m, "st", HOFFSET(s1, st), tid_m2) < 0) TEST_ERROR; if (H5Tinsert(tid_m, "l", HOFFSET(s1, l), H5T_NATIVE_ULLONG) < 0) TEST_ERROR; /* Write the data to the dataset */ if (H5Dwrite(dataset, tid_m, H5S_ALL, H5S_ALL, H5P_DEFAULT, points) < 0) TEST_ERROR; /* Close dataset */ if (H5Dclose(dataset) < 0) TEST_ERROR; /* Close dataspace */ if (H5Sclose(space) < 0) TEST_ERROR; /* Open dataset again to check H5Tget_native_type */ if ((dataset = H5Dopen2(file, DSET_COMPOUND_NAME_2, H5P_DEFAULT)) < 0) TEST_ERROR; if ((dtype = H5Dget_type(dataset)) < 0) TEST_ERROR; if ((native_type = H5Tget_native_type(dtype, H5T_DIR_DEFAULT)) < 0) TEST_ERROR; if (H5Tequal(native_type, tid_m) != TRUE) TEST_ERROR; /* Verify the datatype of each field retrieved and converted */ /* check the char member */ if ((mem_id = H5Tget_member_type(native_type, 0)) < 0) TEST_ERROR; if (H5Tget_order(mem_id) != H5Tget_order(H5T_NATIVE_SCHAR)) TEST_ERROR; if (H5Tget_size(mem_id) < H5Tget_size(H5T_STD_I8LE)) TEST_ERROR; if (H5T_INTEGER != H5Tget_class(mem_id)) TEST_ERROR; H5Tclose(mem_id); /* check the integer member */ if ((mem_id = H5Tget_member_type(native_type, 1)) < 0) TEST_ERROR; if (H5Tget_order(mem_id) != H5Tget_order(H5T_NATIVE_INT)) TEST_ERROR; if (H5Tget_size(mem_id) < H5Tget_size(H5T_STD_I32LE)) TEST_ERROR; if (H5T_INTEGER != H5Tget_class(mem_id)) TEST_ERROR; H5Tclose(mem_id); /* check the unsigned long long member */ if ((mem_id = H5Tget_member_type(native_type, 3)) < 0) TEST_ERROR; if (H5Tget_order(mem_id) != H5Tget_order(H5T_NATIVE_ULLONG)) TEST_ERROR; if (H5Tget_size(mem_id) < H5Tget_size(H5T_STD_U64BE)) TEST_ERROR; if (H5T_INTEGER != H5Tget_class(mem_id)) TEST_ERROR; H5Tclose(mem_id); /* check the nested compound member */ if ((nest_mem_id = H5Tget_member_type(native_type, 2)) < 0) TEST_ERROR; if ((mem_id = H5Tget_member_type(nest_mem_id, 0)) < 0) TEST_ERROR; if (H5Tget_order(mem_id) != H5Tget_order(H5T_NATIVE_SHORT)) TEST_ERROR; if (H5Tget_size(mem_id) < H5Tget_size(H5T_STD_I16BE)) TEST_ERROR; if (H5T_INTEGER != H5Tget_class(mem_id)) TEST_ERROR; H5Tclose(mem_id); if ((mem_id = H5Tget_member_type(nest_mem_id, 1)) < 0) TEST_ERROR; if (H5Tget_order(mem_id) != H5Tget_order(H5T_NATIVE_LONG)) TEST_ERROR; #if H5_SIZEOF_LONG == 4 if (H5Tget_size(mem_id) < H5Tget_size(H5T_STD_I32LE)) TEST_ERROR; #elif H5_SIZEOF_LONG == 8 if (H5Tget_size(mem_id) < H5Tget_size(H5T_STD_I64LE)) TEST_ERROR; #else #error "Unknown 'long' size" #endif if (H5T_INTEGER != H5Tget_class(mem_id)) TEST_ERROR; H5Tclose(mem_id); if ((mem_id = H5Tget_member_type(nest_mem_id, 2)) < 0) TEST_ERROR; if (H5Tget_order(mem_id) != H5Tget_order(H5T_NATIVE_LLONG)) TEST_ERROR; #if H5_SIZEOF_LONG_LONG == 4 if (H5Tget_size(mem_id) < H5Tget_size(H5T_STD_I32LE)) TEST_ERROR; #elif H5_SIZEOF_LONG_LONG == 8 if (H5Tget_size(mem_id) < H5Tget_size(H5T_STD_I64LE)) TEST_ERROR; #else #error "Unknown 'long long' size" #endif if (H5T_INTEGER != H5Tget_class(mem_id)) TEST_ERROR; H5Tclose(mem_id); /* Read the dataset back. Temporary buffer is for special platforms like * Cray */ if (NULL == (tmp = HDmalloc(DIM0 * DIM1 * H5Tget_size(native_type)))) TEST_ERROR; if (NULL == (bkg = HDcalloc(sizeof(s1), DIM0 * DIM1))) TEST_ERROR; if (H5Dread(dataset, native_type, H5S_ALL, H5S_ALL, H5P_DEFAULT, tmp) < 0) TEST_ERROR; HDmemcpy(check, tmp, DIM0 * DIM1 * H5Tget_size(native_type)); HDfree(tmp); tmp = NULL; if (H5Tconvert(native_type, tid_m, (DIM0 * DIM1), check, bkg, H5P_DEFAULT) < 0) TEST_ERROR; HDfree(bkg); bkg = NULL; /* Check that the values read are the same as the values written */ for (i = 0, temp_point = points, temp_check = check; i < DIM0; i++) { for (j = 0; j < DIM1; j++, temp_point++, temp_check++) { if (temp_point->c != temp_check->c || temp_point->i != temp_check->i || temp_point->st.c2 != temp_check->st.c2 || temp_point->st.l2 != temp_check->st.l2 || temp_point->st.ll2 != temp_check->st.ll2 || temp_point->l != temp_check->l) { H5_FAILED(); HDprintf(" Read different values than written.\n"); HDprintf(" At index %d,%d\n", i, j); goto error; } /* end if */ } /* end for */ } /* end for */ /* Close temporary datatypes */ if (H5Tclose(tid2) < 0) TEST_ERROR; if (H5Tclose(tid) < 0) TEST_ERROR; if (H5Tclose(tid_m2) < 0) TEST_ERROR; /* Close HDF5 objects */ H5Dclose(dataset); H5Tclose(dtype); H5Tclose(native_type); H5Tclose(tid_m); /* Free memory for test data */ HDfree(points); HDfree(check); PASSED(); return 0; error: if (tmp) HDfree(tmp); if (bkg) HDfree(bkg); if (points) HDfree(points); if (check) HDfree(check); H5E_BEGIN_TRY { H5Tclose(tid); H5Tclose(tid2); H5Tclose(tid_m); H5Tclose(tid_m2); H5Sclose(space); H5Tclose(mem_id); H5Tclose(nest_mem_id); H5Dclose(dataset); H5Tclose(dtype); H5Tclose(native_type); } H5E_END_TRY; return -1; } /*------------------------------------------------------------------------- * Function: test_compound_dtype * * Purpose: Test H5Tget_native_type for compound datatype * * Return: Success: 0 * * Failure: -1 * * Programmer: Raymond Lu * October 15, 2002 * * Modifications: * *------------------------------------------------------------------------- */ static herr_t test_compound_dtype(hid_t file) { typedef struct { char c; unsigned int i; long long l; } s1; hid_t dataset = -1, space = -1; hid_t dtype = -1, native_type = -1, tid = -1, tid2 = -1, mem_id = -1; int i, j, n; hsize_t dims[2]; s1 *temp_point = NULL; s1 *temp_check = NULL; s1 *points = NULL; s1 *check = NULL; void *tmp = NULL; void *bkg = NULL; TESTING("compound datatype"); /* Allocate space for the points & check arrays */ if (NULL == (points = (s1 *)HDmalloc(sizeof(s1) * DIM0 * DIM1))) TEST_ERROR; if (NULL == (check = (s1 *)HDcalloc(sizeof(s1), DIM0 * DIM1))) TEST_ERROR; /* Initialize the dataset */ for (i = n = 0, temp_point = points; i < DIM0; i++) for (j = 0; j < DIM1; j++, temp_point++) { temp_point->c = 't'; temp_point->i = (unsigned int)(n++); temp_point->l = (i * 10 + j * 100) * n; } /* end for */ /* Create the data space */ dims[0] = DIM0; dims[1] = DIM1; if ((space = H5Screate_simple(2, dims, NULL)) < 0) TEST_ERROR; /* Create compound datatype for disk storage */ if ((tid = H5Tcreate(H5T_COMPOUND, sizeof(s1))) < 0) TEST_ERROR; /* Insert members */ if (H5Tinsert(tid, "c", 0, H5T_STD_U8LE) < 0) TEST_ERROR; if (H5Tinsert(tid, "i", 1, H5T_STD_U32LE) < 0) TEST_ERROR; if (H5Tinsert(tid, "l", 5, H5T_STD_I64BE) < 0) TEST_ERROR; /* Create the dataset */ if ((dataset = H5Dcreate2(file, DSET_COMPOUND_NAME, tid, space, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) TEST_ERROR; /* Create compound datatype for datatype in memory */ if ((tid2 = H5Tcreate(H5T_COMPOUND, sizeof(s1))) < 0) TEST_ERROR; if (H5Tinsert(tid2, "c", HOFFSET(s1, c), H5T_NATIVE_UCHAR) < 0) TEST_ERROR; if (H5Tinsert(tid2, "i", HOFFSET(s1, i), H5T_NATIVE_UINT) < 0) TEST_ERROR; if (H5Tinsert(tid2, "l", HOFFSET(s1, l), H5T_NATIVE_LLONG) < 0) TEST_ERROR; /* Write the data to the dataset */ if (H5Dwrite(dataset, tid2, H5S_ALL, H5S_ALL, H5P_DEFAULT, points) < 0) TEST_ERROR; /* Close dataset */ if (H5Dclose(dataset) < 0) TEST_ERROR; /* Close dataspace */ if (H5Sclose(space) < 0) TEST_ERROR; /* Open dataset again to check H5Tget_native_type */ if ((dataset = H5Dopen2(file, DSET_COMPOUND_NAME, H5P_DEFAULT)) < 0) TEST_ERROR; if ((dtype = H5Dget_type(dataset)) < 0) TEST_ERROR; if ((native_type = H5Tget_native_type(dtype, H5T_DIR_DEFAULT)) < 0) TEST_ERROR; if (H5Tequal(native_type, tid2) != TRUE) TEST_ERROR; /* Verify the datatype of each field retrieved and converted */ if ((mem_id = H5Tget_member_type(native_type, 0)) < 0) TEST_ERROR; if (H5Tget_order(mem_id) != H5Tget_order(H5T_NATIVE_UCHAR)) TEST_ERROR; if (H5Tget_size(mem_id) < H5Tget_size(H5T_STD_U8LE)) TEST_ERROR; if (H5T_INTEGER != H5Tget_class(mem_id)) TEST_ERROR; H5Tclose(mem_id); if ((mem_id = H5Tget_member_type(native_type, 1)) < 0) TEST_ERROR; if (H5Tget_order(mem_id) != H5Tget_order(H5T_NATIVE_UINT)) TEST_ERROR; if (H5Tget_size(mem_id) < H5Tget_size(H5T_STD_U32LE)) TEST_ERROR; if (H5T_INTEGER != H5Tget_class(mem_id)) TEST_ERROR; H5Tclose(mem_id); if ((mem_id = H5Tget_member_type(native_type, 2)) < 0) TEST_ERROR; if (H5Tget_order(mem_id) != H5Tget_order(H5T_NATIVE_LLONG)) TEST_ERROR; if (H5Tget_size(mem_id) < H5Tget_size(H5T_STD_I64BE)) TEST_ERROR; if (H5T_INTEGER != H5Tget_class(mem_id)) TEST_ERROR; H5Tclose(mem_id); /* Read the dataset back. Temporary buffer is for special platforms like * Cray */ if (NULL == (tmp = HDmalloc(DIM0 * DIM1 * H5Tget_size(native_type)))) TEST_ERROR; if (NULL == (bkg = HDcalloc(sizeof(s1), DIM0 * DIM1))) TEST_ERROR; if (H5Dread(dataset, native_type, H5S_ALL, H5S_ALL, H5P_DEFAULT, tmp) < 0) TEST_ERROR; HDmemcpy(check, tmp, DIM0 * DIM1 * H5Tget_size(native_type)); HDfree(tmp); tmp = NULL; if (H5Tconvert(native_type, tid2, (DIM0 * DIM1), check, bkg, H5P_DEFAULT) < 0) TEST_ERROR; HDfree(bkg); bkg = NULL; /* Check that the values read are the same as the values written */ for (i = 0, temp_point = points, temp_check = check; i < DIM0; i++) for (j = 0; j < DIM1; j++, temp_point++, temp_check++) if (temp_point->c != temp_check->c || temp_point->i != temp_check->i || temp_point->l != temp_check->l) { H5_FAILED(); HDprintf(" Read different values than written.\n"); HDprintf(" At index %d,%d\n", i, j); goto error; } /* end if */ /* Close datatype */ if (H5Tclose(tid) < 0) TEST_ERROR; H5Dclose(dataset); H5Tclose(dtype); H5Tclose(native_type); H5Tclose(tid2); /* Free memory for test data */ HDfree(points); HDfree(check); PASSED(); return 0; error: /* Free memory for test data */ if (tmp) HDfree(tmp); if (bkg) HDfree(bkg); if (points) HDfree(points); if (check) HDfree(check); H5E_BEGIN_TRY { H5Tclose(tid); H5Sclose(space); H5Tclose(mem_id); H5Dclose(dataset); H5Tclose(dtype); H5Tclose(native_type); H5Tclose(tid2); } H5E_END_TRY; return -1; } /*------------------------------------------------------------------------- * Function: test_compound_dtype3 * * Purpose: Test H5Tget_native_type for compound datatype * * Return: Success: 0 * * Failure: -1 * * Programmer: Raymond Lu * October 15, 2002 * * Modifications: * *------------------------------------------------------------------------- */ static herr_t test_compound_dtype3(hid_t file) { typedef struct { char c; int a[5]; long long l; } s1; hid_t dataset = -1, space = -1; hid_t dtype = -1, native_type = -1, tid = -1, tid2 = -1, tid_m = -1, tid_m2 = -1, mem_id = -1, nest_mem_id = -1; hsize_t array_dims[1] = {5}; int i, j, k, n; hsize_t dims[2]; s1 *temp_point = NULL, *temp_check = NULL; s1 *points = NULL, *check = NULL; void *tmp = NULL, *bkg = NULL; TESTING("compound datatype with array as field"); /* Allocate space for the points & check arrays */ if (NULL == (points = (s1 *)HDmalloc(sizeof(s1) * DIM0 * DIM1))) TEST_ERROR; if (NULL == (check = (s1 *)HDcalloc(sizeof(s1), DIM0 * DIM1))) TEST_ERROR; /* Initialize the dataset */ for (i = n = 0, temp_point = points; i < DIM0; i++) for (j = 0; j < DIM1; j++, temp_point++) { temp_point->c = 't'; temp_point->l = (i * 10 + j * 100) * n; for (k = 0; k < 5; k++) (temp_point->a)[k] = n++; } /* end for */ /* Create the data space */ dims[0] = DIM0; dims[1] = DIM1; if ((space = H5Screate_simple(2, dims, NULL)) < 0) TEST_ERROR; /* Create array datatype */ if ((tid2 = H5Tarray_create2(H5T_STD_I32LE, 1, array_dims)) < 0) TEST_ERROR; /* Create compound datatype for disk storage */ if ((tid = H5Tcreate(H5T_COMPOUND, 29)) < 0) TEST_ERROR; /* Insert members */ if (H5Tinsert(tid, "c", 0, H5T_STD_U8LE) < 0) TEST_ERROR; if (H5Tinsert(tid, "a", 1, tid2) < 0) TEST_ERROR; if (H5Tinsert(tid, "l", 21, H5T_STD_I64BE) < 0) TEST_ERROR; /* Create the dataset */ if ((dataset = H5Dcreate2(file, DSET_COMPOUND_NAME_3, tid, space, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) TEST_ERROR; /* Create array datatype */ if ((tid_m2 = H5Tarray_create2(H5T_NATIVE_INT, 1, array_dims)) < 0) TEST_ERROR; /* Create compound datatype for datatype in memory */ if ((tid_m = H5Tcreate(H5T_COMPOUND, sizeof(s1))) < 0) TEST_ERROR; if (H5Tinsert(tid_m, "c", HOFFSET(s1, c), H5T_NATIVE_UCHAR) < 0) TEST_ERROR; if (H5Tinsert(tid_m, "a", HOFFSET(s1, a), tid_m2) < 0) TEST_ERROR; if (H5Tinsert(tid_m, "l", HOFFSET(s1, l), H5T_NATIVE_LLONG) < 0) TEST_ERROR; /* Write the data to the dataset */ if (H5Dwrite(dataset, tid_m, H5S_ALL, H5S_ALL, H5P_DEFAULT, points) < 0) TEST_ERROR; /* Close dataset */ if (H5Dclose(dataset) < 0) TEST_ERROR; /* Close datatype */ if (H5Tclose(tid) < 0) TEST_ERROR; if (H5Tclose(tid2) < 0) TEST_ERROR; /* Close dataspace */ if (H5Sclose(space) < 0) TEST_ERROR; /* Open dataset again to check H5Tget_native_type */ if ((dataset = H5Dopen2(file, DSET_COMPOUND_NAME_3, H5P_DEFAULT)) < 0) TEST_ERROR; if ((dtype = H5Dget_type(dataset)) < 0) TEST_ERROR; if ((native_type = H5Tget_native_type(dtype, H5T_DIR_DEFAULT)) < 0) TEST_ERROR; /* Verify the datatype of each field retrieved and converted */ /* check the char member */ if ((mem_id = H5Tget_member_type(native_type, 0)) < 0) TEST_ERROR; if (H5Tget_order(mem_id) != H5Tget_order(H5T_NATIVE_UCHAR)) TEST_ERROR; if (H5Tget_size(mem_id) < H5Tget_size(H5T_STD_U8LE)) TEST_ERROR; if (H5T_INTEGER != H5Tget_class(mem_id)) TEST_ERROR; H5Tclose(mem_id); /* check the array member */ if ((mem_id = H5Tget_member_type(native_type, 1)) < 0) TEST_ERROR; if (H5T_ARRAY != H5Tget_class(mem_id)) TEST_ERROR; if ((nest_mem_id = H5Tget_super(mem_id)) < 0) TEST_ERROR; if (H5Tget_order(nest_mem_id) != H5Tget_order(H5T_NATIVE_INT)) TEST_ERROR; if (H5Tget_size(nest_mem_id) < H5Tget_size(H5T_STD_I32LE)) TEST_ERROR; if (H5T_INTEGER != H5Tget_class(nest_mem_id)) TEST_ERROR; H5Tclose(nest_mem_id); H5Tclose(mem_id); /* check the long long member */ if ((mem_id = H5Tget_member_type(native_type, 2)) < 0) TEST_ERROR; if (H5Tget_order(mem_id) != H5Tget_order(H5T_NATIVE_LLONG)) TEST_ERROR; if (H5Tget_size(mem_id) < H5Tget_size(H5T_STD_I64BE)) TEST_ERROR; if (H5T_INTEGER != H5Tget_class(mem_id)) TEST_ERROR; H5Tclose(mem_id); /* Read the dataset back. Temporary buffer is for special platforms like * Cray */ if (NULL == (tmp = HDmalloc(DIM0 * DIM1 * H5Tget_size(native_type)))) TEST_ERROR; if (NULL == (bkg = HDcalloc(sizeof(s1), DIM0 * DIM1))) TEST_ERROR; if (H5Dread(dataset, native_type, H5S_ALL, H5S_ALL, H5P_DEFAULT, tmp) < 0) TEST_ERROR; HDmemcpy(check, tmp, DIM0 * DIM1 * H5Tget_size(native_type)); HDfree(tmp); tmp = NULL; if (H5Tconvert(native_type, tid_m, (DIM0 * DIM1), check, bkg, H5P_DEFAULT) < 0) TEST_ERROR; HDfree(bkg); bkg = NULL; /* Check that the values read are the same as the values written */ for (i = 0, temp_point = points, temp_check = check; i < DIM0; i++) for (j = 0; j < DIM1; j++, temp_point++, temp_check++) { if (temp_point->c != temp_check->c || temp_point->l != temp_check->l) { H5_FAILED(); HDprintf(" Read different values than written.\n"); HDprintf(" At index %d,%d\n", i, j); goto error; } /* end if */ for (k = 0; k < 5; k++) { if (temp_point->a[k] != temp_check->a[k]) { H5_FAILED(); HDprintf(" Read different values than written.\n"); HDprintf(" At index %d,%d,%d\n", i, j, k); goto error; } /* end if */ } /* end for */ } /* end for */ H5Dclose(dataset); H5Tclose(dtype); H5Tclose(native_type); H5Tclose(tid_m); H5Tclose(tid_m2); /* Free memory for test data */ HDfree(points); HDfree(check); PASSED(); return 0; error: /* Free memory for test data */ if (tmp) HDfree(tmp); if (bkg) HDfree(bkg); if (points) HDfree(points); if (check) HDfree(check); H5E_BEGIN_TRY { H5Sclose(space); H5Tclose(mem_id); H5Tclose(nest_mem_id); H5Dclose(dataset); H5Tclose(dtype); H5Tclose(native_type); H5Tclose(tid); H5Tclose(tid2); H5Tclose(tid_m); H5Tclose(tid_m2); } H5E_END_TRY; return -1; } /*------------------------------------------------------------------------- * Function: test_compound_opaque * * Purpose: Test H5Tget_native_type for compound datatype with opaque field * * Return: Success: 0 * * Failure: -1 * * Programmer: Quincey Koziol * January 31, 2004 * * Modifications: * *------------------------------------------------------------------------- */ static herr_t test_compound_opaque(hid_t file) { typedef struct { char c; unsigned char o[5]; long long l; } s1; hid_t dataset = -1, space = -1; hid_t dtype = -1, native_type = -1, tid = -1, tid2 = -1, tid_m = -1, mem_id = -1; int i, j, k, n; hsize_t dims[2]; s1 *temp_point = NULL, *temp_check = NULL; s1 *points = NULL, *check = NULL; void *tmp = NULL, *bkg = NULL; TESTING("compound datatype with opaque field"); /* Allocate space for the points & check arrays */ if (NULL == (points = (s1 *)HDmalloc(sizeof(s1) * DIM0 * DIM1))) TEST_ERROR; if (NULL == (check = (s1 *)HDcalloc(sizeof(s1), DIM0 * DIM1))) TEST_ERROR; /* Initialize the dataset */ for (i = n = 0, temp_point = points; i < DIM0; i++) for (j = 0; j < DIM1; j++, temp_point++) { temp_point->c = 't'; temp_point->l = (i * 10 + j * 100) * n; for (k = 0; k < 5; k++) (temp_point->o)[k] = (unsigned char)(n++); } /* end for */ /* Create the data space */ dims[0] = DIM0; dims[1] = DIM1; if ((space = H5Screate_simple(2, dims, NULL)) < 0) TEST_ERROR; /* Create opaque datatype */ if ((tid2 = H5Tcreate(H5T_OPAQUE, sizeof(temp_point->o))) < 0) TEST_ERROR; if (H5Tset_tag(tid2, "testing opaque field") < 0) TEST_ERROR; /* Create compound datatype for disk storage */ if ((tid = H5Tcreate(H5T_COMPOUND, 14)) < 0) TEST_ERROR; /* Insert members */ if (H5Tinsert(tid, "c", 0, H5T_STD_U8LE) < 0) TEST_ERROR; if (H5Tinsert(tid, "o", 1, tid2) < 0) TEST_ERROR; if (H5Tinsert(tid, "l", 6, H5T_STD_I64BE) < 0) TEST_ERROR; /* Create the dataset */ if ((dataset = H5Dcreate2(file, DSET_COMPOUND_NAME_4, tid, space, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) TEST_ERROR; /* Create compound datatype for datatype in memory */ if ((tid_m = H5Tcreate(H5T_COMPOUND, sizeof(s1))) < 0) TEST_ERROR; if (H5Tinsert(tid_m, "c", HOFFSET(s1, c), H5T_NATIVE_UCHAR) < 0) TEST_ERROR; if (H5Tinsert(tid_m, "o", HOFFSET(s1, o), tid2) < 0) TEST_ERROR; if (H5Tinsert(tid_m, "l", HOFFSET(s1, l), H5T_NATIVE_LLONG) < 0) TEST_ERROR; /* Write the data to the dataset */ if (H5Dwrite(dataset, tid_m, H5S_ALL, H5S_ALL, H5P_DEFAULT, points) < 0) TEST_ERROR; /* Close dataset */ if (H5Dclose(dataset) < 0) TEST_ERROR; /* Close datatype */ if (H5Tclose(tid) < 0) TEST_ERROR; if (H5Tclose(tid2) < 0) TEST_ERROR; /* Close dataspace */ if (H5Sclose(space) < 0) TEST_ERROR; /* Open dataset again to check H5Tget_native_type */ if ((dataset = H5Dopen2(file, DSET_COMPOUND_NAME_4, H5P_DEFAULT)) < 0) TEST_ERROR; if ((dtype = H5Dget_type(dataset)) < 0) TEST_ERROR; if ((native_type = H5Tget_native_type(dtype, H5T_DIR_DEFAULT)) < 0) TEST_ERROR; /* Verify the datatype of each field retrieved and converted */ /* check the char member */ if ((mem_id = H5Tget_member_type(native_type, 0)) < 0) TEST_ERROR; if (H5Tget_order(mem_id) != H5Tget_order(H5T_NATIVE_UCHAR)) TEST_ERROR; if (H5Tget_size(mem_id) < H5Tget_size(H5T_STD_U8LE)) TEST_ERROR; if (H5T_INTEGER != H5Tget_class(mem_id)) TEST_ERROR; H5Tclose(mem_id); /* check the array member */ if ((mem_id = H5Tget_member_type(native_type, 1)) < 0) TEST_ERROR; if (H5T_OPAQUE != H5Tget_class(mem_id)) TEST_ERROR; if (H5Tget_size(mem_id) != sizeof(temp_point->o)) TEST_ERROR; H5Tclose(mem_id); /* check the long long member */ if ((mem_id = H5Tget_member_type(native_type, 2)) < 0) TEST_ERROR; if (H5Tget_order(mem_id) != H5Tget_order(H5T_NATIVE_LLONG)) TEST_ERROR; if (H5Tget_size(mem_id) < H5Tget_size(H5T_STD_I64BE)) TEST_ERROR; if (H5T_INTEGER != H5Tget_class(mem_id)) TEST_ERROR; H5Tclose(mem_id); /* Read the dataset back. Temporary buffer is for special platforms like * Cray */ if (NULL == (tmp = HDmalloc(DIM0 * DIM1 * H5Tget_size(native_type)))) TEST_ERROR; if (NULL == (bkg = HDcalloc(sizeof(s1), DIM0 * DIM1))) TEST_ERROR; if (H5Dread(dataset, native_type, H5S_ALL, H5S_ALL, H5P_DEFAULT, tmp) < 0) TEST_ERROR; HDmemcpy(check, tmp, DIM0 * DIM1 * H5Tget_size(native_type)); HDfree(tmp); tmp = NULL; if (H5Tconvert(native_type, tid_m, (DIM0 * DIM1), check, bkg, H5P_DEFAULT) < 0) TEST_ERROR; HDfree(bkg); bkg = NULL; /* Check that the values read are the same as the values written */ for (i = 0, temp_point = points, temp_check = check; i < DIM0; i++) for (j = 0; j < DIM1; j++, temp_point++, temp_check++) { if (temp_point->c != temp_check->c || temp_point->l != temp_check->l) { H5_FAILED(); HDprintf(" Read different values than written.\n"); HDprintf(" At index %d,%d\n", i, j); goto error; } /* end if */ for (k = 0; k < 5; k++) { if (temp_point->o[k] != temp_check->o[k]) { H5_FAILED(); HDprintf(" Read different values than written.\n"); HDprintf(" At index %d,%d,%d\n", i, j, k); goto error; } /* end if */ } /* end for */ } /* end for */ H5Dclose(dataset); H5Tclose(dtype); H5Tclose(native_type); H5Tclose(tid_m); /* Free memory for test data */ HDfree(points); HDfree(check); PASSED(); return 0; error: /* Free memory for test data */ if (tmp) HDfree(tmp); if (bkg) HDfree(bkg); if (points) HDfree(points); if (check) HDfree(check); H5E_BEGIN_TRY { H5Sclose(space); H5Tclose(mem_id); H5Dclose(dataset); H5Tclose(dtype); H5Tclose(native_type); H5Tclose(tid); H5Tclose(tid2); H5Tclose(tid_m); } H5E_END_TRY; return -1; } /*------------------------------------------------------------------------- * Function: test_enum_dtype * * Purpose: Test H5Tget_native_type for enumerate datatype * * Return: Success: 0 * * Failure: -1 * * Programmer: Raymond Lu * October 15, 2002 * * Modifications: * *------------------------------------------------------------------------- */ static herr_t test_enum_dtype(hid_t file) { hid_t dataset = -1, space = -1; hid_t tid = -1, tid_m = -1, dtype = -1, native_type = -1; int i, j, n; hsize_t dims[2]; void *tmp = NULL; short colors[8]; unsigned char sub_colors[16]; const char *mname[] = {"RED", "GREEN", "BLUE", "YELLOW", "PINK", "PURPLE", "ORANGE", "WHITE"}; struct { short arr[DIM0][DIM1]; } *spoints2 = NULL; struct { short arr[DIM0][DIM1]; } *scheck2 = NULL; TESTING("enum datatype"); if (NULL == (spoints2 = HDcalloc(1, sizeof(*spoints2)))) TEST_ERROR; if (NULL == (scheck2 = HDcalloc(1, sizeof(*scheck2)))) TEST_ERROR; /* Initialize the dataset */ for (i = 0; i < DIM0; i++) for (j = 0, n = 0; j < DIM1; j++, n++) spoints2->arr[i][j] = (short)((i * 10 + j * 100 + n) % 8); /* Create the data space */ dims[0] = DIM0; dims[1] = DIM1; if ((space = H5Screate_simple(2, dims, NULL)) < 0) TEST_ERROR; /* Construct enum type based on native type */ if ((tid = H5Tenum_create(H5T_STD_I16LE)) < 0) TEST_ERROR; for (i = 0; i < 8; i++) { sub_colors[i * 2] = (unsigned char)i; sub_colors[i * 2 + 1] = 0; if (H5Tenum_insert(tid, mname[i], &(sub_colors[i * 2])) < 0) TEST_ERROR; } /* end for */ /* Create the dataset */ if ((dataset = H5Dcreate2(file, DSET_ENUM_NAME, tid, space, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) TEST_ERROR; /* Construct enum type based on native type in memory */ if ((tid_m = H5Tenum_create(H5T_NATIVE_SHORT)) < 0) TEST_ERROR; for (i = 0; i < 8; i++) { colors[i] = (short)i; if (H5Tenum_insert(tid_m, mname[i], &(colors[i])) < 0) TEST_ERROR; } /* end for */ /* Write the data to the dataset */ if (H5Dwrite(dataset, tid_m, H5S_ALL, H5S_ALL, H5P_DEFAULT, spoints2) < 0) TEST_ERROR; /* Close dataset */ if (H5Dclose(dataset) < 0) TEST_ERROR; /* Close datatype */ if (H5Tclose(tid) < 0) TEST_ERROR; /* Close dataspace */ if (H5Sclose(space) < 0) TEST_ERROR; /* Open dataset again to check H5Tget_native_type */ if ((dataset = H5Dopen2(file, DSET_ENUM_NAME, H5P_DEFAULT)) < 0) TEST_ERROR; if ((dtype = H5Dget_type(dataset)) < 0) TEST_ERROR; if ((native_type = H5Tget_native_type(dtype, H5T_DIR_DEFAULT)) < 0) TEST_ERROR; /* Read the dataset back. Temporary buffer is for special platforms like * Cray */ if (NULL == (tmp = HDmalloc(DIM0 * DIM1 * H5Tget_size(native_type)))) TEST_ERROR; if (H5Dread(dataset, native_type, H5S_ALL, H5S_ALL, H5P_DEFAULT, tmp) < 0) TEST_ERROR; HDmemcpy(scheck2, tmp, DIM0 * DIM1 * H5Tget_size(native_type)); if (H5Tconvert(native_type, tid_m, (DIM0 * DIM1), scheck2, NULL, H5P_DEFAULT) < 0) TEST_ERROR; /* Check that the values read are the same as the values written */ for (i = 0; i < DIM0; i++) for (j = 0; j < DIM1; j++) if (spoints2->arr[i][j] != scheck2->arr[i][j]) { H5_FAILED(); HDprintf(" Read different values than written.\n"); HDprintf(" At index %d,%d\n", i, j); HDprintf(" spoints2[i][j]=%hd, scheck2[i][j]=%hd\n", spoints2->arr[i][j], scheck2->arr[i][j]); goto error; } /* end if */ H5Dclose(dataset); H5Tclose(dtype); H5Tclose(native_type); H5Tclose(tid_m); HDfree(tmp); HDfree(spoints2); HDfree(scheck2); PASSED(); return 0; error: /* Free memory for test data */ HDfree(tmp); HDfree(spoints2); HDfree(scheck2); H5E_BEGIN_TRY { H5Sclose(space); H5Dclose(dataset); H5Tclose(dtype); H5Tclose(native_type); H5Tclose(tid); H5Tclose(tid_m); } H5E_END_TRY; return -1; } /*------------------------------------------------------------------------- * Function: test_array_dtype * * Purpose: Test H5Tget_native_type for array datatype * * Return: Success: 0 * * Failure: -1 * * Programmer: Raymond Lu * October 15, 2002 * * Modifications: * *------------------------------------------------------------------------- */ static herr_t test_array_dtype(hid_t file) { typedef struct { char c; int i; long long l; } s1; hid_t dataset = -1, space = -1; hid_t dtype = -1, native_type = -1, tid = -1, tid2 = -1, tid3 = -1, tid_m = -1; int i, j, k, n; hsize_t space_dims[2], array_dims[1] = {5}; s1 *temp_point = NULL, *temp_check = NULL; s1 *points = NULL, *check = NULL; void *tmp = NULL; TESTING("array of compound datatype"); /* Allocate space for the points & check arrays */ if (NULL == (points = (s1 *)HDmalloc(sizeof(s1) * DIM0 * DIM1 * 5))) TEST_ERROR; if (NULL == (check = (s1 *)HDcalloc(sizeof(s1), DIM0 * DIM1 * 5))) TEST_ERROR; /* Initialize the dataset */ for (i = n = 0, temp_point = points; i < DIM0; i++) for (j = 0; j < DIM1; j++) for (k = 0; k < 5; k++, temp_point++) { temp_point->c = 't'; temp_point->i = n++; temp_point->l = (i * 10 + j * 100) * n; } /* end for */ /* Create the data space */ space_dims[0] = DIM0; space_dims[1] = DIM1; if ((space = H5Screate_simple(2, space_dims, NULL)) < 0) TEST_ERROR; /* Create compound datatype for disk storage */ if ((tid2 = H5Tcreate(H5T_COMPOUND, 13)) < 0) TEST_ERROR; /* Insert members */ if (H5Tinsert(tid2, "c", 0, H5T_STD_U8BE) < 0) TEST_ERROR; if (H5Tinsert(tid2, "i", 1, H5T_STD_U32LE) < 0) TEST_ERROR; if (H5Tinsert(tid2, "l", 5, H5T_STD_I64BE) < 0) TEST_ERROR; /* Create array datatype for disk storage */ if ((tid = H5Tarray_create2(tid2, 1, array_dims)) < 0) TEST_ERROR; /* Create the dataset */ if ((dataset = H5Dcreate2(file, DSET_ARRAY_NAME, tid, space, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) TEST_ERROR; /* Create compound datatype for datatype in memory */ if ((tid3 = H5Tcreate(H5T_COMPOUND, sizeof(s1))) < 0) TEST_ERROR; if (H5Tinsert(tid3, "c", HOFFSET(s1, c), H5T_NATIVE_UCHAR) < 0) TEST_ERROR; if (H5Tinsert(tid3, "i", HOFFSET(s1, i), H5T_NATIVE_UINT) < 0) TEST_ERROR; if (H5Tinsert(tid3, "l", HOFFSET(s1, l), H5T_NATIVE_LLONG) < 0) TEST_ERROR; /* Create array datatype for memory */ if ((tid_m = H5Tarray_create2(tid3, 1, array_dims)) < 0) TEST_ERROR; /* Write the data to the dataset */ if (H5Dwrite(dataset, tid_m, H5S_ALL, H5S_ALL, H5P_DEFAULT, points) < 0) TEST_ERROR; /* Close dataset */ if (H5Dclose(dataset) < 0) TEST_ERROR; /* Close datatype */ if (H5Tclose(tid) < 0) TEST_ERROR; if (H5Tclose(tid2) < 0) TEST_ERROR; /* Close dataspace */ if (H5Sclose(space) < 0) TEST_ERROR; /* Open dataset again to check H5Tget_native_type */ if ((dataset = H5Dopen2(file, DSET_ARRAY_NAME, H5P_DEFAULT)) < 0) TEST_ERROR; if ((dtype = H5Dget_type(dataset)) < 0) TEST_ERROR; if ((native_type = H5Tget_native_type(dtype, H5T_DIR_DEFAULT)) < 0) TEST_ERROR; /* Read the dataset back. Temporary buffer is for special platforms like * Cray */ if (NULL == (tmp = HDmalloc(DIM0 * DIM1 * H5Tget_size(native_type)))) TEST_ERROR; if (H5Dread(dataset, native_type, H5S_ALL, H5S_ALL, H5P_DEFAULT, tmp) < 0) TEST_ERROR; HDmemcpy(check, tmp, DIM0 * DIM1 * H5Tget_size(native_type)); HDfree(tmp); tmp = NULL; if (H5Tconvert(native_type, tid_m, (DIM0 * DIM1), check, NULL, H5P_DEFAULT) < 0) TEST_ERROR; /* Check that the values read are the same as the values written */ for (i = 0, temp_point = points, temp_check = check; i < DIM0; i++) for (j = 0; j < DIM1; j++) for (k = 0; k < 5; k++, temp_point++, temp_check++) if (temp_point->c != temp_check->c || temp_point->i != temp_check->i || temp_point->l != temp_check->l) { H5_FAILED(); HDprintf(" Read different values than written.\n"); HDprintf(" At index %d,%d\n", i, j); goto error; } /* end if */ /* Close HDF5 objects */ if (H5Dclose(dataset)) TEST_ERROR; if (H5Tclose(native_type)) TEST_ERROR; if (H5Tclose(dtype)) TEST_ERROR; if (H5Tclose(tid_m) < 0) TEST_ERROR; if (H5Tclose(tid3) < 0) TEST_ERROR; /* Free memory for test data */ HDfree(points); HDfree(check); PASSED(); return 0; error: /* Free memory for test data */ if (tmp) HDfree(tmp); if (points) HDfree(points); if (check) HDfree(check); H5E_BEGIN_TRY { H5Sclose(space); H5Dclose(dataset); H5Tclose(dtype); H5Tclose(native_type); H5Tclose(tid); H5Tclose(tid2); H5Tclose(tid3); H5Tclose(tid_m); } H5E_END_TRY; return -1; } /*------------------------------------------------------------------------- * Function: test_array_dtype2 * * Purpose: Test H5Tget_native_type for array datatype * * Return: Success: 0 * * Failure: -1 * * Programmer: Raymond Lu * October 15, 2002 * * Modifications: * *------------------------------------------------------------------------- */ static herr_t test_array_dtype2(hid_t file) { hid_t dataset = -1, space = -1; hid_t dtype = -1, native_type = -1, tid = -1, tid_m = -1; int i, j, k, n; hsize_t space_dims[2], array_dims[1] = {5}; void *tmp = NULL; struct { int arr[DIM0][DIM1][5]; } *ipoints3 = NULL; struct { int arr[DIM0][DIM1][5]; } *icheck3 = NULL; TESTING("array of atomic datatype"); if (NULL == (ipoints3 = HDcalloc(1, sizeof(*ipoints3)))) goto error; if (NULL == (icheck3 = HDcalloc(1, sizeof(*icheck3)))) goto error; /* Initialize the dataset */ for (i = n = 0; i < DIM0; i++) for (j = 0; j < DIM1; j++) for (k = 0; k < 5; k++) ipoints3->arr[i][j][k] = n++; /* Create the data space */ space_dims[0] = DIM0; space_dims[1] = DIM1; if ((space = H5Screate_simple(2, space_dims, NULL)) < 0) TEST_ERROR; /* Create array datatype for disk storage */ if ((tid = H5Tarray_create2(H5T_STD_I32LE, 1, array_dims)) < 0) TEST_ERROR; /* Create the dataset */ if ((dataset = H5Dcreate2(file, DSET_ARRAY2_NAME, tid, space, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) TEST_ERROR; /* Create array datatype for memory */ if ((tid_m = H5Tarray_create2(H5T_NATIVE_INT, 1, array_dims)) < 0) TEST_ERROR; /* Write the data to the dataset */ if (H5Dwrite(dataset, tid_m, H5S_ALL, H5S_ALL, H5P_DEFAULT, ipoints3) < 0) TEST_ERROR; /* Close dataset */ if (H5Dclose(dataset) < 0) TEST_ERROR; /* Close datatype */ if (H5Tclose(tid) < 0) TEST_ERROR; /* Close dataspace */ if (H5Sclose(space) < 0) TEST_ERROR; /* Open dataset again to check H5Tget_native_type */ if ((dataset = H5Dopen2(file, DSET_ARRAY2_NAME, H5P_DEFAULT)) < 0) TEST_ERROR; if ((dtype = H5Dget_type(dataset)) < 0) TEST_ERROR; if ((native_type = H5Tget_native_type(dtype, H5T_DIR_DEFAULT)) < 0) TEST_ERROR; /* Read the dataset back. Temporary buffer is for special platforms like * Cray */ if (NULL == (tmp = HDmalloc(DIM0 * DIM1 * H5Tget_size(native_type)))) TEST_ERROR; if (H5Dread(dataset, native_type, H5S_ALL, H5S_ALL, H5P_DEFAULT, tmp) < 0) TEST_ERROR; HDmemcpy(icheck3, tmp, DIM0 * DIM1 * H5Tget_size(native_type)); HDfree(tmp); tmp = NULL; if (H5Tconvert(native_type, tid_m, (DIM0 * DIM1), icheck3, NULL, H5P_DEFAULT) < 0) TEST_ERROR; /* Check that the values read are the same as the values written */ for (i = 0; i < DIM0; i++) for (j = 0; j < DIM1; j++) for (k = 0; k < 5; k++) if (icheck3->arr[i][j][k] != ipoints3->arr[i][j][k]) { H5_FAILED(); HDprintf(" Read different values than written.\n"); HDprintf(" At index %d,%d\n", i, j); goto error; } /* end if */ /* Close HDF5 objects */ if (H5Dclose(dataset)) TEST_ERROR; if (H5Tclose(native_type)) TEST_ERROR; if (H5Tclose(dtype)) TEST_ERROR; if (H5Tclose(tid_m) < 0) TEST_ERROR; HDfree(ipoints3); HDfree(icheck3); PASSED(); return 0; error: /* Free memory for test data */ if (tmp) HDfree(tmp); H5E_BEGIN_TRY { H5Sclose(space); H5Dclose(dataset); H5Tclose(dtype); H5Tclose(native_type); H5Tclose(tid); H5Tclose(tid_m); } H5E_END_TRY; HDfree(ipoints3); HDfree(icheck3); return -1; } /*------------------------------------------------------------------------- * Function: test_vl_dtype * * Purpose: Test H5Tget_native_type for variable length datatype * * Return: Success: 0 * * Failure: -1 * * Programmer: Raymond Lu * October 15, 2002 * * Modifications: * *------------------------------------------------------------------------- */ static herr_t test_vl_dtype(hid_t file) { hvl_t wdata[SPACE1_DIM1]; /* Information to write */ hvl_t rdata[SPACE1_DIM1]; /* Information read in */ hvl_t *t1, *t2; /* Temporary pointer to VL information */ hsize_t dims1[] = {SPACE1_DIM1}; hid_t dataset = -1, space = -1; hid_t dtype = -1, native_type = -1, nat_super_type = -1, tid = -1, tid2 = -1, tid_m = -1, tid_m2 = -1; size_t i, j, k; void **tmp = NULL; TESTING("variable length datatype"); /* Allocate and initialize VL data to write */ for (i = 0; i < SPACE1_DIM1; i++) { wdata[i].p = HDmalloc((i + 1) * sizeof(hvl_t)); if (NULL == wdata[i].p) { H5_FAILED(); HDprintf(" Cannot allocate memory for VL data! i=%u\n", (unsigned)i); goto error; } /* end if */ wdata[i].len = i + 1; for (t1 = (hvl_t *)wdata[i].p, j = 0; j < (i + 1); j++, t1++) { t1->p = HDmalloc((j + 1) * sizeof(unsigned int)); if (NULL == t1->p) { H5_FAILED(); HDprintf(" Cannot allocate memory for VL data! i=%u, j=%u\n", (unsigned)i, (unsigned)j); goto error; } /* end if */ t1->len = j + 1; for (k = 0; k < (j + 1); k++) ((unsigned int *)t1->p)[k] = (unsigned int)(i * 100 + j * 10 + k); } /* end for */ } /* end for */ /* Create dataspace for datasets */ if ((space = H5Screate_simple(SPACE1_RANK, dims1, NULL)) < 0) TEST_ERROR; /* Create the base VL type */ if ((tid2 = H5Tvlen_create(H5T_STD_U32LE)) < 0) TEST_ERROR; /* Create a VL datatype for disk storage */ if ((tid = H5Tvlen_create(tid2)) < 0) TEST_ERROR; /* Create a dataset */ if ((dataset = H5Dcreate2(file, DSET_VL_NAME, tid, space, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) TEST_ERROR; /* Create a base VL datatype for memory */ if ((tid_m2 = H5Tvlen_create(H5T_NATIVE_UINT)) < 0) TEST_ERROR; /* Create a VL datatype for memory */ if ((tid_m = H5Tvlen_create(tid_m2)) < 0) TEST_ERROR; /* Write dataset to disk */ if (H5Dwrite(dataset, tid_m, H5S_ALL, H5S_ALL, H5P_DEFAULT, wdata) < 0) TEST_ERROR; /* Close Dataset */ if (H5Dclose(dataset) < 0) TEST_ERROR; /* Close datatype */ if (H5Tclose(tid2) < 0) TEST_ERROR; if (H5Tclose(tid) < 0) TEST_ERROR; /* Open a dataset */ if ((dataset = H5Dopen2(file, DSET_VL_NAME, H5P_DEFAULT)) < 0) TEST_ERROR; /* Get native datatype for dataset */ if ((dtype = H5Dget_type(dataset)) < 0) TEST_ERROR; if ((native_type = H5Tget_native_type(dtype, H5T_DIR_DEFAULT)) < 0) TEST_ERROR; /* Also get native base type for this nested VL type. Should be an integer type. */ if ((nat_super_type = H5Tget_super(native_type)) < 0) TEST_ERROR; if ((nat_super_type = H5Tget_super(nat_super_type)) < 0) TEST_ERROR; /* Read dataset from disk */ if (H5Dread(dataset, native_type, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata) < 0) TEST_ERROR; /* Compare data read in */ for (i = 0; i < SPACE1_DIM1; i++) { if (wdata[i].len != rdata[i].len) { H5_FAILED(); HDprintf(" VL data length don't match!, wdata[%d].len=%d, rdata[%d].len=%d\n", (int)i, (int)wdata[i].len, (int)i, (int)rdata[i].len); goto error; } /* end if */ for (t1 = (hvl_t *)wdata[i].p, t2 = (hvl_t *)rdata[i].p, j = 0; j < rdata[i].len; j++, t1++, t2++) { if (t1->len != t2->len) { H5_FAILED(); HDprintf(" VL data length don't match!, wdata[%d].len=%d, rdata[%d].len=%d\n", (int)i, (int)wdata[i].len, (int)i, (int)rdata[i].len); goto error; } /* end if */ /* use temporary buffer to convert datatype. This is for special * platforms like Cray */ if (NULL == (tmp = (void **)HDmalloc(t2->len * sizeof(unsigned int)))) TEST_ERROR; HDmemcpy(tmp, t2->p, t2->len * H5Tget_size(nat_super_type)); if (H5Tconvert(nat_super_type, H5T_NATIVE_UINT, t2->len, tmp, NULL, H5P_DEFAULT) < 0) TEST_ERROR; for (k = 0; k < t2->len; k++) { if (((unsigned int *)t1->p)[k] != ((unsigned int *)tmp)[k]) { H5_FAILED(); HDprintf(" VL data don't match!, wdata[%u].p=%d, rdata[%u].p=%u\n", (unsigned)i, ((unsigned int *)t1->p)[k], (unsigned)i, ((unsigned int *)tmp)[k]); goto error; } /* end if */ } /* end for */ HDfree(tmp); tmp = NULL; } /* end for */ } /* end for */ /* Reclaim the read VL data */ if (H5Treclaim(native_type, space, H5P_DEFAULT, rdata) < 0) TEST_ERROR; /* Reclaim the write VL data */ if (H5Treclaim(native_type, space, H5P_DEFAULT, wdata) < 0) TEST_ERROR; /* Close Dataset */ if (H5Dclose(dataset) < 0) TEST_ERROR; /* Close datatype */ if (H5Tclose(native_type) < 0) TEST_ERROR; native_type = -1; /* reset so that error handling can check for VL reclaim */ if (H5Tclose(dtype) < 0) TEST_ERROR; if (H5Tclose(tid_m) < 0) TEST_ERROR; if (H5Tclose(tid_m2) < 0) TEST_ERROR; /* Close disk dataspace */ if (H5Sclose(space) < 0) TEST_ERROR; PASSED(); return 0; error: /* Free memory for test data */ if (tmp) HDfree(tmp); H5E_BEGIN_TRY { if (native_type > 0) { H5Treclaim(native_type, space, H5P_DEFAULT, rdata); H5Treclaim(native_type, space, H5P_DEFAULT, wdata); } /* end if */ H5Sclose(space); H5Dclose(dataset); H5Tclose(dtype); H5Tclose(native_type); H5Tclose(nat_super_type); H5Tclose(tid); H5Tclose(tid2); H5Tclose(tid_m); H5Tclose(tid_m2); } H5E_END_TRY; return -1; } /* end test_vl_type() */ /*------------------------------------------------------------------------- * Function: test_vlstr_dtype * * Purpose: Test H5Tget_native_type for variable length string datatype * * Return: Success: 0 * * Failure: -1 * * Programmer: Raymond Lu * October 15, 2002 * * Modifications: * *------------------------------------------------------------------------- */ static herr_t test_vlstr_dtype(hid_t file) { /* Information to write */ const char *wdata[SPACE1_DIM1] = { "Four score and seven years ago our forefathers brought forth on this continent a new nation,", "conceived in liberty and dedicated to the proposition that all men are created equal.", "Now we are engaged in a great civil war,", "testing whether that nation or any nation so conceived and so dedicated can long endure."}; char *rdata[SPACE1_DIM1]; /* Information read in */ hbool_t rdata_alloc = FALSE; /* Whether the read data is allocated */ hid_t dataset = -1; /* Dataset ID */ hid_t sid1 = -1; /* Dataspace ID */ hid_t tid1 = -1, dtype = -1, native_type = -1; /* Datatype ID */ hsize_t dims1[] = {SPACE1_DIM1}; unsigned i; /* counting variable */ /* Output message about test being performed */ TESTING("variable length string datatype"); /* Create dataspace for datasets */ if ((sid1 = H5Screate_simple(SPACE1_RANK, dims1, NULL)) < 0) TEST_ERROR; /* Create a datatype to refer to */ if ((tid1 = H5Tcopy(H5T_C_S1)) < 0) TEST_ERROR; if (H5Tset_size(tid1, H5T_VARIABLE) < 0) TEST_ERROR; if (H5T_STRING != H5Tget_class(tid1) || !H5Tis_variable_str(tid1)) TEST_ERROR; /* Create a dataset */ if ((dataset = H5Dcreate2(file, DSET_VLSTR_NAME, tid1, sid1, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) TEST_ERROR; /* Write dataset to disk */ if (H5Dwrite(dataset, tid1, H5S_ALL, H5S_ALL, H5P_DEFAULT, wdata) < 0) TEST_ERROR; /* Close Dataset */ if (H5Dclose(dataset) < 0) TEST_ERROR; /* Open a dataset */ if ((dataset = H5Dopen2(file, DSET_VLSTR_NAME, H5P_DEFAULT)) < 0) TEST_ERROR; /* Get datatype for dataset */ if ((dtype = H5Dget_type(dataset)) < 0) TEST_ERROR; /* Construct native type */ if ((native_type = H5Tget_native_type(dtype, H5T_DIR_DEFAULT)) < 0) TEST_ERROR; /* Check if the data type is equal */ if (!H5Tequal(native_type, tid1)) TEST_ERROR; /* Read dataset from disk */ if (H5Dread(dataset, native_type, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata) < 0) TEST_ERROR; rdata_alloc = TRUE; /* Compare data read in */ for (i = 0; i < SPACE1_DIM1; i++) { if (HDstrlen(wdata[i]) != HDstrlen(rdata[i])) { H5_FAILED(); HDprintf(" VL data length don't match!, strlen(wdata[%d])=%d, strlen(rdata[%d])=%d\n", (int)i, (int)HDstrlen(wdata[i]), (int)i, (int)HDstrlen(rdata[i])); goto error; } /* end if */ if (HDstrcmp(wdata[i], rdata[i]) != 0) { H5_FAILED(); HDprintf(" VL data values don't match!, wdata[%d]=%s, rdata[%d]=%s\n", (int)i, wdata[i], (int)i, rdata[i]); goto error; } /* end if */ } /* end for */ /* Close Dataset */ if (H5Dclose(dataset) < 0) TEST_ERROR; /* Close datatype */ if (H5Tclose(tid1) < 0) TEST_ERROR; if (H5Tclose(native_type) < 0) TEST_ERROR; /* Close disk dataspace */ if (H5Sclose(sid1) < 0) TEST_ERROR; /* Free memory for rdata */ for (i = 0; i < SPACE1_DIM1; i++) HDfree(rdata[i]); rdata_alloc = FALSE; PASSED(); return 0; error: if (rdata_alloc) { /* Free memory for rdata */ for (i = 0; i < SPACE1_DIM1; i++) HDfree(rdata[i]); } /* end if */ H5E_BEGIN_TRY { H5Dclose(dataset); H5Tclose(dtype); H5Tclose(native_type); H5Tclose(sid1); H5Tclose(tid1); } H5E_END_TRY; return -1; } /* end test_vlstr_dtype() */ /*------------------------------------------------------------------------- * Function: test_str_dtype * * Purpose: Test H5Tget_native_type for fixed-length string datatype * * Return: Success: 0 * * Failure: -1 * * Programmer: Raymond Lu * October 15, 2002 * * Modifications: * *------------------------------------------------------------------------- */ static herr_t test_str_dtype(hid_t file) { const char wdata[SPACE1_DIM1][4] = {"one", "two", "3rd", "4th"}; /* Information to write */ char rdata[SPACE1_DIM1][4]; /* Information read in */ hid_t dataset = -1; /* Dataset ID */ hid_t sid1 = -1; /* Dataspace ID */ hid_t tid1 = -1, dtype = -1, native_type = -1; /* Datatype ID */ hsize_t dims1[] = {SPACE1_DIM1}; unsigned i; /* counting variable */ /* Output message about test being performed */ TESTING("fixed-length string datatype"); /* Create dataspace for datasets */ if ((sid1 = H5Screate_simple(SPACE1_RANK, dims1, NULL)) < 0) TEST_ERROR; /* Create a datatype to refer to */ if ((tid1 = H5Tcopy(H5T_C_S1)) < 0) TEST_ERROR; if (H5Tset_size(tid1, 4) < 0) TEST_ERROR; if (H5T_STRING != H5Tget_class(tid1)) TEST_ERROR; /* Create a dataset */ if ((dataset = H5Dcreate2(file, DSET_STR_NAME, tid1, sid1, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) TEST_ERROR; /* Write dataset to disk */ if (H5Dwrite(dataset, tid1, H5S_ALL, H5S_ALL, H5P_DEFAULT, wdata) < 0) TEST_ERROR; /* Close Dataset */ if (H5Dclose(dataset) < 0) TEST_ERROR; /* Open a dataset */ if ((dataset = H5Dopen2(file, DSET_STR_NAME, H5P_DEFAULT)) < 0) TEST_ERROR; /* Get datatype for dataset */ if ((dtype = H5Dget_type(dataset)) < 0) TEST_ERROR; /* Construct native type */ if ((native_type = H5Tget_native_type(dtype, H5T_DIR_DEFAULT)) < 0) TEST_ERROR; /* Check if the data type is equal */ if (!H5Tequal(native_type, tid1) || H5T_STRING != H5Tget_class(native_type)) TEST_ERROR; /* Read dataset from disk */ if (H5Dread(dataset, native_type, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata) < 0) TEST_ERROR; /* Compare data read in */ for (i = 0; i < SPACE1_DIM1; i++) { if (HDstrlen(wdata[i]) != HDstrlen(rdata[i])) { H5_FAILED(); HDprintf(" data length don't match!, strlen(wdata[%d])=%d, strlen(rdata[%d])=%d\n", (int)i, (int)HDstrlen(wdata[i]), (int)i, (int)HDstrlen(rdata[i])); goto error; } /* end if */ if (HDstrcmp(wdata[i], rdata[i]) != 0) { H5_FAILED(); HDprintf(" data values don't match!, wdata[%d]=%s, rdata[%d]=%s\n", (int)i, wdata[i], (int)i, rdata[i]); goto error; } /* end if */ } /* end for */ /* Close Dataset */ if (H5Dclose(dataset) < 0) TEST_ERROR; /* Close datatype */ if (H5Tclose(tid1) < 0) TEST_ERROR; if (H5Tclose(native_type) < 0) TEST_ERROR; /* Close disk dataspace */ if (H5Sclose(sid1) < 0) TEST_ERROR; PASSED(); return 0; error: H5E_BEGIN_TRY { H5Dclose(dataset); H5Tclose(dtype); H5Tclose(native_type); H5Tclose(tid1); H5Tclose(sid1); } H5E_END_TRY; return -1; } /* end test_str_dtype() */ /*------------------------------------------------------------------------- * Function: test_refer_dtype * * Purpose: Test H5Tget_native_type for reference datatype * * Return: Success: 0 * * Failure: -1 * * Programmer: Raymond Lu * October 15, 2002 * * Modifications: * *------------------------------------------------------------------------- */ static herr_t test_refer_dtype(hid_t file) { /* Compound datatype */ typedef struct s1_t { unsigned int a; unsigned int b; float c; } s1_t; hid_t dataset = -1; /* Dataset ID */ hid_t group = -1; /* Group ID */ hid_t sid1 = -1; /* Dataspace ID */ hid_t tid1 = -1, dtype = -1, native_type = -1; /* Datatype ID */ hsize_t dims1[] = {1}; H5O_type_t obj_type; /* Object type */ hobj_ref_t *wbuf = NULL, /* buffer to write to disk */ *rbuf = NULL; /* buffer read from disk */ /* Output message about test being performed */ TESTING("reference datatype"); /* Allocate write & read buffers */ if (NULL == (wbuf = (hobj_ref_t *)HDmalloc(MAX(sizeof(unsigned), sizeof(hobj_ref_t))))) TEST_ERROR; if (NULL == (rbuf = (hobj_ref_t *)HDmalloc(MAX(sizeof(unsigned), sizeof(hobj_ref_t))))) TEST_ERROR; /* Create dataspace for datasets */ if ((sid1 = H5Screate_simple(SPACE1_RANK, dims1, NULL)) < 0) TEST_ERROR; /* Create a group */ if ((group = H5Gcreate2(file, "Group1", H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) FAIL_STACK_ERROR; /* Create a datatype to refer to */ if ((tid1 = H5Tcreate(H5T_COMPOUND, sizeof(s1_t))) < 0) TEST_ERROR; /* Insert fields */ if (H5Tinsert(tid1, "a", HOFFSET(s1_t, a), H5T_NATIVE_INT) < 0) TEST_ERROR; if (H5Tinsert(tid1, "b", HOFFSET(s1_t, b), H5T_NATIVE_INT) < 0) TEST_ERROR; if (H5Tinsert(tid1, "c", HOFFSET(s1_t, c), H5T_NATIVE_FLOAT) < 0) TEST_ERROR; /* Save datatype for later */ if (H5Tcommit2(group, "Datatype1", tid1, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT) < 0) TEST_ERROR; /* Close datatype */ if (H5Tclose(tid1) < 0) TEST_ERROR; /* Close group */ if (H5Gclose(group) < 0) TEST_ERROR; /* Create a dataset */ if ((dataset = H5Dcreate2(file, "Dataset3", H5T_STD_REF_OBJ, sid1, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) TEST_ERROR; /* Create reference to named datatype */ if (H5Rcreate(wbuf, file, "/Group1/Datatype1", H5R_OBJECT, (hid_t)-1) < 0) TEST_ERROR; if (H5Rget_obj_type2(dataset, H5R_OBJECT, wbuf, &obj_type) < 0) TEST_ERROR; if (obj_type != H5O_TYPE_NAMED_DATATYPE) TEST_ERROR; /* Write selection to disk */ if (H5Dwrite(dataset, H5T_STD_REF_OBJ, H5S_ALL, H5S_ALL, H5P_DEFAULT, wbuf) < 0) TEST_ERROR; /* Close disk dataspace */ if (H5Sclose(sid1) < 0) TEST_ERROR; /* Close Dataset */ if (H5Dclose(dataset) < 0) TEST_ERROR; /* Open the dataset */ if ((dataset = H5Dopen2(file, "/Dataset3", H5P_DEFAULT)) < 0) TEST_ERROR; /* Get datatype for dataset */ if ((dtype = H5Dget_type(dataset)) < 0) TEST_ERROR; /* Construct native type */ if ((native_type = H5Tget_native_type(dtype, H5T_DIR_DEFAULT)) < 0) TEST_ERROR; /* Check if the data type is equal */ if (!H5Tequal(native_type, H5T_STD_REF_OBJ)) TEST_ERROR; /* Read selection from disk */ if (H5Dread(dataset, native_type, H5S_ALL, H5S_ALL, H5P_DEFAULT, rbuf) < 0) TEST_ERROR; /* Open datatype object */ if ((tid1 = H5Rdereference2(dataset, H5P_DEFAULT, H5R_OBJECT, rbuf)) < 0) TEST_ERROR; /* Verify correct datatype */ if (H5Tget_class(tid1) != H5T_COMPOUND) TEST_ERROR; if (H5Tget_nmembers(tid1) != 3) TEST_ERROR; /* Close datatype */ if (H5Tclose(tid1) < 0) TEST_ERROR; if (H5Tclose(native_type) < 0) TEST_ERROR; /* Close Dataset */ if (H5Dclose(dataset) < 0) TEST_ERROR; /* Free memory buffers */ HDfree(wbuf); HDfree(rbuf); PASSED(); return 0; error: if (wbuf) HDfree(wbuf); if (rbuf) HDfree(rbuf); H5E_BEGIN_TRY { H5Sclose(sid1); H5Gclose(group); H5Tclose(tid1); H5Tclose(dtype); H5Tclose(native_type); H5Dclose(dataset); } H5E_END_TRY; return -1; } /* test_refer_dtype() */ /*------------------------------------------------------------------------- * Function: test_refer_dtype2 * * Purpose: Test H5Tget_native_type for reference * * Return: Success: 0 * * Failure: -1 * * Programmer: Raymond Lu * October 15, 2002 * * Modifications: * *------------------------------------------------------------------------- */ static herr_t test_refer_dtype2(hid_t file) { hid_t dset1 = -1, /* Dataset ID */ dset2 = -1; /* Dereferenced dataset ID */ hid_t sid1 = -1, /* Dataspace ID #1 */ sid2 = -1; /* Dataspace ID #2 */ hid_t dtype = -1, native_type = -1; hsize_t dims1[] = {1}, dims2[] = {SPACE2_DIM1, SPACE2_DIM2}; hsize_t start[SPACE2_RANK]; /* Starting location of hyperslab */ hsize_t stride[SPACE2_RANK]; /* Stride of hyperslab */ hsize_t count[SPACE2_RANK]; /* Element count of hyperslab */ hsize_t block[SPACE2_RANK]; /* Block size of hyperslab */ hdset_reg_ref_t wbuf, /* buffer to write to disk */ rbuf; /* buffer read from disk */ uint8_t *dwbuf = NULL, /* Buffer for writing numeric data to disk */ *drbuf = NULL; /* Buffer for reading numeric data from disk */ uint8_t *tu8 = NULL; /* Temporary pointer to uint8 data */ H5O_type_t obj_type; /* Object type */ int i; /* counting variables */ /* Output message about test being performed */ TESTING("dataset region reference"); /* Allocate write & read buffers */ if (NULL == (dwbuf = (uint8_t *)HDmalloc(sizeof(uint8_t) * SPACE2_DIM1 * SPACE2_DIM2))) TEST_ERROR; if (NULL == (drbuf = (uint8_t *)HDcalloc(sizeof(uint8_t), SPACE2_DIM1 * SPACE2_DIM2))) TEST_ERROR; /* Create dataspace for datasets */ if ((sid2 = H5Screate_simple(SPACE2_RANK, dims2, NULL)) < 0) TEST_ERROR; /* Create a dataset */ if ((dset2 = H5Dcreate2(file, "Dataset2", H5T_STD_U8LE, sid2, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) TEST_ERROR; for (tu8 = dwbuf, i = 0; i < SPACE2_DIM1 * SPACE2_DIM2; i++) *tu8++ = (uint8_t)(i * 3); /* Write selection to disk */ if (H5Dwrite(dset2, H5T_STD_U8LE, H5S_ALL, H5S_ALL, H5P_DEFAULT, dwbuf) < 0) TEST_ERROR; /* Close Dataset */ if (H5Dclose(dset2) < 0) TEST_ERROR; /* Create dataspace for the reference dataset */ if ((sid1 = H5Screate_simple(SPACE1_RANK, dims1, NULL)) < 0) TEST_ERROR; /* Create a reference dataset */ if ((dset1 = H5Dcreate2(file, "Dataset1", H5T_STD_REF_DSETREG, sid1, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) TEST_ERROR; /* Create references */ /* Select 6x6 hyperslab for first reference */ start[0] = 2; start[1] = 2; stride[0] = 1; stride[1] = 1; count[0] = 1; count[1] = 1; block[0] = 6; block[1] = 6; if (H5Sselect_hyperslab(sid2, H5S_SELECT_SET, start, stride, count, block) < 0) TEST_ERROR; if ((int)H5Sget_select_npoints(sid2) != 36) TEST_ERROR; /* Store first dataset region */ if (H5Rcreate(&wbuf, file, "/Dataset2", H5R_DATASET_REGION, sid2) < 0) TEST_ERROR; if (H5Rget_obj_type2(dset1, H5R_DATASET_REGION, &wbuf, &obj_type) < 0) TEST_ERROR; if (obj_type != H5O_TYPE_DATASET) TEST_ERROR; /* Write selection to disk */ if (H5Dwrite(dset1, H5T_STD_REF_DSETREG, H5S_ALL, H5S_ALL, H5P_DEFAULT, &wbuf) < 0) TEST_ERROR; /* Close disk dataspace */ if (H5Sclose(sid1) < 0) TEST_ERROR; /* Close Dataset */ if (H5Dclose(dset1) < 0) TEST_ERROR; /* Close uint8 dataset dataspace */ if (H5Sclose(sid2) < 0) TEST_ERROR; /* Open the dataset */ if ((dset1 = H5Dopen2(file, "/Dataset1", H5P_DEFAULT)) < 0) TEST_ERROR; /* Get datatype for dataset */ if ((dtype = H5Dget_type(dset1)) < 0) TEST_ERROR; /* Construct native type */ if ((native_type = H5Tget_native_type(dtype, H5T_DIR_DEFAULT)) < 0) TEST_ERROR; /* Check if the data type is equal */ if (!H5Tequal(native_type, H5T_STD_REF_DSETREG)) TEST_ERROR; /* Read selection from disk */ if (H5Dread(dset1, H5T_STD_REF_DSETREG, H5S_ALL, H5S_ALL, H5P_DEFAULT, &rbuf) < 0) TEST_ERROR; /* Try to open objects */ if ((dset2 = H5Rdereference2(dset1, H5P_DEFAULT, H5R_DATASET_REGION, &rbuf)) < 0) TEST_ERROR; /* Check what H5Rget_obj_type2 function returns */ if (H5Rget_obj_type2(dset1, H5R_DATASET_REGION, &rbuf, &obj_type) < 0) TEST_ERROR; if (obj_type != H5O_TYPE_DATASET) TEST_ERROR; /* Check information in referenced dataset */ if ((sid1 = H5Dget_space(dset2)) < 0) TEST_ERROR; if ((int)H5Sget_simple_extent_npoints(sid1) != 100) TEST_ERROR; /* Read from disk */ if (H5Dread(dset2, H5T_STD_U8LE, H5S_ALL, H5S_ALL, H5P_DEFAULT, drbuf) < 0) TEST_ERROR; for (tu8 = (uint8_t *)drbuf, i = 0; i < (SPACE2_DIM1 * SPACE2_DIM2); i++, tu8++) if (*tu8 != (uint8_t)(i * 3)) TEST_ERROR; /* Get the hyperslab selection */ if ((sid2 = H5Rget_region(dset1, H5R_DATASET_REGION, &rbuf)) < 0) TEST_ERROR; /* Verify correct hyperslab selected */ if ((int)H5Sget_select_npoints(sid2) != 36) TEST_ERROR; if ((int)H5Sget_select_hyper_nblocks(sid2) != 1) TEST_ERROR; /* Close region space */ if (H5Sclose(sid2) < 0) TEST_ERROR; /* Close first space */ if (H5Sclose(sid1) < 0) TEST_ERROR; /* Close dereferenced Dataset */ if (H5Dclose(dset2) < 0) TEST_ERROR; /* Close Dataset */ if (H5Dclose(dset1) < 0) TEST_ERROR; /* Free memory buffers */ HDfree(dwbuf); HDfree(drbuf); PASSED(); return 0; error: /* Free memory buffers */ if (dwbuf) HDfree(dwbuf); if (drbuf) HDfree(drbuf); H5E_BEGIN_TRY { H5Sclose(sid2); H5Sclose(sid1); H5Tclose(dtype); H5Tclose(native_type); H5Dclose(dset2); H5Dclose(dset1); } H5E_END_TRY; return -1; } /* test_refer_dtype2() */ /*------------------------------------------------------------------------- * Function: test_opaque_dtype * * Purpose: Test H5Tget_native_type for opaque datatype * * Return: Success: 0 * * Failure: -1 * * Programmer: Raymond Lu * October 15, 2002 * * Modifications: * *------------------------------------------------------------------------- */ static herr_t test_opaque_dtype(hid_t file) { hid_t type = -1, space = -1, dset = -1; hid_t dataset = -1, dtype = -1, native_type = -1; size_t i; unsigned char wbuf[32], rbuf[32]; hsize_t nelmts; TESTING("opaque datatype"); /* opaque_1 */ nelmts = sizeof(wbuf); if ((type = H5Tcreate(H5T_OPAQUE, 1)) < 0) TEST_ERROR; if (H5Tset_tag(type, "testing 1-byte opaque type") < 0) TEST_ERROR; if ((space = H5Screate_simple(1, &nelmts, NULL)) < 0) TEST_ERROR; if ((dset = H5Dcreate2(file, DSET_OPAQUE_NAME, type, space, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) TEST_ERROR; for (i = 0; i < sizeof(wbuf); i++) wbuf[i] = (unsigned char)(0xff ^ i); if (H5Dwrite(dset, type, H5S_ALL, H5S_ALL, H5P_DEFAULT, wbuf) < 0) TEST_ERROR; if (H5Sclose(space) < 0) TEST_ERROR; if (H5Dclose(dset) < 0) TEST_ERROR; /* Open dataset again to check H5Tget_native_type */ if ((dataset = H5Dopen2(file, DSET_OPAQUE_NAME, H5P_DEFAULT)) < 0) TEST_ERROR; if ((dtype = H5Dget_type(dataset)) < 0) TEST_ERROR; if ((native_type = H5Tget_native_type(dtype, H5T_DIR_DEFAULT)) < 0) TEST_ERROR; if (!H5Tequal(native_type, type)) TEST_ERROR; if (H5Dread(dataset, native_type, H5S_ALL, H5S_ALL, H5P_DEFAULT, rbuf) < 0) TEST_ERROR; for (i = 0; i < sizeof(rbuf); i++) if (rbuf[i] != wbuf[i]) { H5_FAILED(); HDprintf(" Read different values than written.\n"); HDprintf(" At index %u\n", (unsigned)i); goto error; } /* end if */ if (H5Tclose(type) < 0) TEST_ERROR; if (H5Tclose(dtype) < 0) TEST_ERROR; if (H5Tclose(native_type) < 0) TEST_ERROR; if (H5Dclose(dataset) < 0) TEST_ERROR; PASSED(); return 0; error: H5E_BEGIN_TRY { H5Sclose(space); H5Tclose(type); H5Tclose(dtype); H5Tclose(native_type); H5Dclose(dset); H5Dclose(dataset); } H5E_END_TRY; return -1; } /* test_opaque_dtype */ /*------------------------------------------------------------------------- * Function: test_bitfield_dtype * * Purpose: Test H5Tget_native_type for bitfield datatype * * Return: Success: 0 * * Failure: -1 * * Programmer: Raymond Lu * October 15, 2002 * * Modifications: * Raymond Lu * 1 December 2009 * I added the support for bitfield and changed the test to * compare the data being read back. *------------------------------------------------------------------------- */ static herr_t test_bitfield_dtype(hid_t file) { hid_t type = -1, space = -1, dset1 = -1, dset2 = -1; hid_t dataset1 = -1, dataset2 = -1, dtype = -1, native_type = -1; size_t ntype_size, i; unsigned char wbuf[BITFIELD_ENUMB * sizeof(int)]; unsigned char *p = NULL; void *rbuf = NULL; unsigned int intw[BITFIELD_ENUMB], intr[BITFIELD_ENUMB]; hsize_t nelmts; TESTING("bitfield datatype"); nelmts = BITFIELD_ENUMB; if ((type = H5Tcopy(H5T_STD_B32BE)) < 0) TEST_ERROR; if ((space = H5Screate_simple(1, &nelmts, NULL)) < 0) TEST_ERROR; /* Create and write to dataset1 with a unsigned char buffer */ if ((dset1 = H5Dcreate2(file, DSET1_BITFIELD_NAME, type, space, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) TEST_ERROR; for (i = 0; i < BITFIELD_ENUMB * sizeof(int); i++) wbuf[i] = (unsigned char)((unsigned int)0xff ^ (unsigned int)i); if (H5Dwrite(dset1, H5T_NATIVE_B32, H5S_ALL, H5S_ALL, H5P_DEFAULT, wbuf) < 0) TEST_ERROR; if (H5Dclose(dset1) < 0) TEST_ERROR; /* Create and write to dataset2 with a unsigned int buffer */ if ((dset2 = H5Dcreate2(file, DSET2_BITFIELD_NAME, type, space, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) TEST_ERROR; for (i = 0; i < BITFIELD_ENUMB; i++) intw[i] = (unsigned int)0xff << (unsigned int)((i * 8) % 32); if (H5Dwrite(dset2, H5T_NATIVE_B32, H5S_ALL, H5S_ALL, H5P_DEFAULT, intw) < 0) TEST_ERROR; if (H5Dclose(dset2) < 0) TEST_ERROR; if (H5Sclose(space) < 0) TEST_ERROR; if (H5Tclose(type) < 0) TEST_ERROR; /* Open dataset1 again to check H5Tget_native_type */ if ((dataset1 = H5Dopen2(file, DSET1_BITFIELD_NAME, H5P_DEFAULT)) < 0) TEST_ERROR; if ((dtype = H5Dget_type(dataset1)) < 0) TEST_ERROR; if ((native_type = H5Tget_native_type(dtype, H5T_DIR_DEFAULT)) < 0) TEST_ERROR; if ((ntype_size = H5Tget_size(native_type)) == 0) TEST_ERROR; rbuf = HDmalloc((size_t)nelmts * ntype_size); /* Read the data and compare them */ if (H5Dread(dataset1, native_type, H5S_ALL, H5S_ALL, H5P_DEFAULT, rbuf) < 0) TEST_ERROR; p = (unsigned char *)rbuf; for (i = 0; i < BITFIELD_ENUMB * 4; i++) { if (*p != wbuf[i]) { H5_FAILED(); HDprintf(" Read different values than written.\n"); HDprintf(" At index %zu\n", i); TEST_ERROR; } p++; } if (H5Tclose(dtype) < 0) TEST_ERROR; if (H5Tclose(native_type) < 0) TEST_ERROR; if (H5Dclose(dataset1) < 0) TEST_ERROR; if (rbuf) HDfree(rbuf); /* Open dataset2 again to check H5Tget_native_type */ if ((dataset2 = H5Dopen2(file, DSET2_BITFIELD_NAME, H5P_DEFAULT)) < 0) TEST_ERROR; if ((dtype = H5Dget_type(dataset2)) < 0) TEST_ERROR; if ((native_type = H5Tget_native_type(dtype, H5T_DIR_DEFAULT)) < 0) TEST_ERROR; /* Read the data and compare them */ if (H5Dread(dataset2, native_type, H5S_ALL, H5S_ALL, H5P_DEFAULT, intr) < 0) TEST_ERROR; for (i = 0; i < BITFIELD_ENUMB; i++) { if (intr[i] != intw[i]) { H5_FAILED(); HDprintf(" Read different values than written.\n"); HDprintf(" At index %zu\n", i); TEST_ERROR; } } if (H5Tclose(dtype) < 0) TEST_ERROR; if (H5Tclose(native_type) < 0) TEST_ERROR; if (H5Dclose(dataset2) < 0) TEST_ERROR; PASSED(); return 0; error: H5E_BEGIN_TRY { H5Sclose(space); H5Tclose(type); H5Tclose(dtype); H5Tclose(native_type); H5Dclose(dset1); H5Dclose(dset2); H5Dclose(dataset1); H5Dclose(dataset2); } H5E_END_TRY; return -1; } /* test_bitfield_dtype */ /*------------------------------------------------------------------------- * Function: test_ninteger * * Purpose: Test the native integer function; made to check the case * like the Cray SV1, where the size of short is 8 but precision is 32 * * Return: Success: 0 * Failure: -1 * * Programmer: Pedro Vicente * September 3, 2004 * * Modifications: * *------------------------------------------------------------------------- */ static herr_t test_ninteger(void) { hid_t fid1 = -1; /* file ID */ hid_t fid2 = -1; /* file ID */ hid_t did1 = -1; /* dataset ID */ hid_t did2 = -1; /* dataset ID */ hid_t sid1 = -1; /* dataspace ID */ hid_t dcpl1 = -1; /* dataset creation property list ID */ hid_t dcpl2 = -1; /* dataset creation property list ID */ hid_t tid1 = -1; /* file datatype */ hid_t tid2 = -1; /* file datatype */ hid_t nid1 = -1; /* native datatype */ hid_t nid2 = -1; /* native datatype */ hsize_t dims[1] = {DIM3}; /* dataspace dimensions */ int rank = 1; /* rank of dataset */ int buf[DIM3]; int chk[DIM3]; int i; TESTING("native integer "); for (i = 0; i < DIM3; i++) buf[i] = i; /*------------------------------------------------------------------------- * step1: create a file *------------------------------------------------------------------------- */ /* create a file using default properties */ if ((fid1 = H5Fcreate("tstint1.h5", H5F_ACC_TRUNC, H5P_DEFAULT, H5P_DEFAULT)) < 0) FAIL_STACK_ERROR; /* create a data space */ if ((sid1 = H5Screate_simple(rank, dims, NULL)) < 0) FAIL_STACK_ERROR; /* create dcpl */ if ((dcpl1 = H5Pcreate(H5P_DATASET_CREATE)) < 0) FAIL_STACK_ERROR; /* create a dataset */ if ((did1 = H5Dcreate2(fid1, "dset", H5T_NATIVE_INT, sid1, H5P_DEFAULT, dcpl1, H5P_DEFAULT)) < 0) FAIL_STACK_ERROR; /* write */ if (H5Dwrite(did1, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, buf) < 0) FAIL_STACK_ERROR; /* close */ if (H5Sclose(sid1) < 0) FAIL_STACK_ERROR; if (H5Pclose(dcpl1) < 0) FAIL_STACK_ERROR; if (H5Dclose(did1) < 0) FAIL_STACK_ERROR; if (H5Fclose(fid1) < 0) FAIL_STACK_ERROR; /*------------------------------------------------------------------------- * step 2: open and create another file copying the data from file1 *------------------------------------------------------------------------- */ /* open */ if ((fid1 = H5Fopen("tstint1.h5", H5F_ACC_RDONLY, H5P_DEFAULT)) < 0) FAIL_STACK_ERROR; /* open dataset */ if ((did1 = H5Dopen2(fid1, "dset", H5P_DEFAULT)) < 0) FAIL_STACK_ERROR; if ((sid1 = H5Dget_space(did1)) < 0) FAIL_STACK_ERROR; /* get dcpl */ if ((dcpl1 = H5Dget_create_plist(did1)) < 0) FAIL_STACK_ERROR; /* get file datatype */ if ((tid1 = H5Dget_type(did1)) < 0) FAIL_STACK_ERROR; /* get native datatype */ if ((nid1 = H5Tget_native_type(tid1, H5T_DIR_DEFAULT)) < 0) FAIL_STACK_ERROR; /* get size */ if (H5Tget_size(nid1) == 0) FAIL_STACK_ERROR; /* get rank */ if ((rank = H5Sget_simple_extent_ndims(sid1)) < 0) FAIL_STACK_ERROR; HDmemset(dims, 0, sizeof dims); /* get dimension */ if (H5Sget_simple_extent_dims(sid1, dims, NULL) < 0) FAIL_STACK_ERROR; /* read */ if (H5Dread(did1, nid1, H5S_ALL, H5S_ALL, H5P_DEFAULT, chk) < 0) FAIL_STACK_ERROR; /* create a file using default properties */ if ((fid2 = H5Fcreate("tstint2.h5", H5F_ACC_TRUNC, H5P_DEFAULT, H5P_DEFAULT)) < 0) FAIL_STACK_ERROR; /* create a dataset using the native type */ if ((did2 = H5Dcreate2(fid2, "dset", nid1, sid1, H5P_DEFAULT, dcpl1, H5P_DEFAULT)) < 0) FAIL_STACK_ERROR; /* write */ if (H5Dwrite(did2, nid1, H5S_ALL, H5S_ALL, H5P_DEFAULT, chk) < 0) FAIL_STACK_ERROR; /* get dcpl */ if ((dcpl2 = H5Dget_create_plist(did2)) < 0) FAIL_STACK_ERROR; /* get file datatype */ if ((tid2 = H5Dget_type(did2)) < 0) FAIL_STACK_ERROR; /* get native datatype */ if ((nid2 = H5Tget_native_type(tid2, H5T_DIR_DEFAULT)) < 0) FAIL_STACK_ERROR; /* check */ if (H5Tget_precision(nid1) != H5Tget_precision(nid2)) { HDprintf(" Precision differ.\n"); TEST_ERROR; } /* end if */ /* compare dataset creation property lists */ if (H5Pequal(dcpl1, dcpl2) <= 0) { HDprintf(" Property lists differ.\n"); TEST_ERROR; } /* end if */ /* check */ for (i = 0; i < DIM3; i++) if (buf[i] != chk[i]) { H5_FAILED(); HDprintf(" Read different values than written.\n"); HDprintf(" At index %d\n", i); TEST_ERROR; } /* end if */ /* close */ if (H5Sclose(sid1) < 0) FAIL_STACK_ERROR; if (H5Pclose(dcpl1) < 0) FAIL_STACK_ERROR; if (H5Pclose(dcpl2) < 0) FAIL_STACK_ERROR; if (H5Tclose(tid1) < 0) FAIL_STACK_ERROR; if (H5Tclose(tid2) < 0) FAIL_STACK_ERROR; if (H5Tclose(nid1) < 0) FAIL_STACK_ERROR; if (H5Tclose(nid2) < 0) FAIL_STACK_ERROR; if (H5Dclose(did1) < 0) FAIL_STACK_ERROR; if (H5Dclose(did2) < 0) FAIL_STACK_ERROR; if (H5Fclose(fid1) < 0) FAIL_STACK_ERROR; if (H5Fclose(fid2) < 0) FAIL_STACK_ERROR; PASSED(); return 0; error: H5E_BEGIN_TRY { H5Pclose(dcpl1); H5Pclose(dcpl2); H5Tclose(tid1); H5Tclose(tid2); H5Tclose(nid1); H5Tclose(nid2); H5Dclose(did1); H5Dclose(did2); H5Sclose(sid1); H5Fclose(fid1); H5Fclose(fid2); } H5E_END_TRY; return -1; } /* end test_ninteger() */ /*------------------------------------------------------------------------- * Function: main * * Purpose: Test H5Tget_native_type for different datatype * * Programmer: Raymond Lu * October 15, 2002 * * Modifications: * *------------------------------------------------------------------------- */ int main(void) { hid_t file, fapl; int nerrors = 0; char filename[1024]; hbool_t driver_is_parallel; h5_reset(); fapl = h5_fileaccess(); h5_fixname(FILENAME[0], fapl, filename, sizeof filename); if ((file = H5Fcreate(filename, H5F_ACC_TRUNC, H5P_DEFAULT, fapl)) < 0) goto error; if (h5_using_parallel_driver(fapl, &driver_is_parallel) < 0) goto error; nerrors += test_atomic_dtype(file) < 0 ? 1 : 0; nerrors += test_compound_dtype(file) < 0 ? 1 : 0; nerrors += test_compound_dtype2(file) < 0 ? 1 : 0; nerrors += test_compound_dtype3(file) < 0 ? 1 : 0; nerrors += test_compound_opaque(file) < 0 ? 1 : 0; nerrors += test_enum_dtype(file) < 0 ? 1 : 0; nerrors += test_array_dtype(file) < 0 ? 1 : 0; nerrors += test_array_dtype2(file) < 0 ? 1 : 0; if (!driver_is_parallel) { nerrors += test_vl_dtype(file) < 0 ? 1 : 0; nerrors += test_vlstr_dtype(file) < 0 ? 1 : 0; } nerrors += test_str_dtype(file) < 0 ? 1 : 0; if (!driver_is_parallel) { nerrors += test_refer_dtype(file) < 0 ? 1 : 0; nerrors += test_refer_dtype2(file) < 0 ? 1 : 0; } nerrors += test_opaque_dtype(file) < 0 ? 1 : 0; nerrors += test_bitfield_dtype(file) < 0 ? 1 : 0; nerrors += test_ninteger() < 0 ? 1 : 0; if (H5Fclose(file) < 0) goto error; /* Verify symbol table messages are cached */ nerrors += (h5_verify_cached_stabs(FILENAME, fapl) < 0 ? 1 : 0); if (nerrors) goto error; HDprintf("All native datatype tests passed.\n"); h5_cleanup(FILENAME, fapl); return 0; error: H5E_BEGIN_TRY { H5Fclose(file); h5_cleanup(FILENAME, fapl); } H5E_END_TRY; nerrors = MAX(1, nerrors); HDprintf("***** %d DATASET TEST%s FAILED! *****\n", nerrors, 1 == nerrors ? "" : "S"); return 1; }