/* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 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. * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * */ /*********************************************************** * * Test program: tcoords * * Test the element coordinates for dataspace selection. For * chunked dataset, when the hyperslab selection of some * dimensions is full, the library optimize it by "flattening" * the fully selected dimensions. This program tests if the * coordinates of selected elements are correctly calculated. * *************************************************************/ #include "testhdf5.h" #define FILENAME "coord.h5" #define SINGLE_END_DSET "single_end_dset" #define MULTI_ENDS_SEL_HYPER_DSET "multiple_ends_dset" #define NAME_LEN 128 /* Data written to the dataset for single block test. Global variable * for convenience. */ static int da_buffer[2][3][6][2]; /*********************************************************** ** ** test_singleEnd_selElements(): Test element selection of only ** one block. ** *************************************************************/ static void test_singleEnd_selElements(hid_t file, bool is_chunked) { hid_t sid, plid, did, msid; char dset_name[NAME_LEN]; /* Dataset name */ size_t elmts_numb; herr_t ret; /* Generic error return */ int i, j, k; hsize_t da_dims[4] = {2, 3, 6, 2}; hsize_t da_chunksize[4] = {1, 3, 3, 2}; /* For testing the full selection in the fastest-growing end */ int mem1_buffer[1][1][6][2]; hsize_t mem1_dims[4] = {1, 1, 6, 2}; hsize_t da_elmts1[12][4] = {{0, 0, 0, 0}, {0, 0, 0, 1}, {0, 0, 1, 0}, {0, 0, 1, 1}, {0, 0, 2, 0}, {0, 0, 2, 1}, {0, 0, 3, 0}, {0, 0, 3, 1}, {0, 0, 4, 0}, {0, 0, 4, 1}, {0, 0, 5, 0}, {0, 0, 5, 1}}; /* For testing the full selection in the slowest-growing end */ int mem2_buffer[2][3][1][1]; hsize_t mem2_dims[4] = {2, 3, 1, 1}; hsize_t da_elmts2[6][4] = {{0, 0, 0, 0}, {0, 1, 0, 0}, {0, 2, 0, 0}, {1, 0, 0, 0}, {1, 1, 0, 0}, {1, 2, 0, 0}}; /* For testing the full selection in the middle dimensions */ int mem3_buffer[1][3][6][1]; hsize_t mem3_dims[4] = {1, 3, 6, 1}; hsize_t da_elmts3[18][4] = {{0, 0, 0, 0}, {0, 0, 1, 0}, {0, 0, 2, 0}, {0, 0, 3, 0}, {0, 0, 4, 0}, {0, 0, 5, 0}, {0, 1, 0, 0}, {0, 1, 1, 0}, {0, 1, 2, 0}, {0, 1, 3, 0}, {0, 1, 4, 0}, {0, 1, 5, 0}, {0, 2, 0, 0}, {0, 2, 1, 0}, {0, 2, 2, 0}, {0, 2, 3, 0}, {0, 2, 4, 0}, {0, 2, 5, 0}}; /* Create and write the dataset */ sid = H5Screate_simple(4, da_dims, da_dims); CHECK(sid, FAIL, "H5Screate_simple"); plid = H5Pcreate(H5P_DATASET_CREATE); CHECK(plid, FAIL, "H5Pcreate"); if (is_chunked) { ret = H5Pset_chunk(plid, 4, da_chunksize); CHECK(ret, FAIL, "H5Pset_chunk"); } /* Construct dataset's name */ memset(dset_name, 0, (size_t)NAME_LEN); strcat(dset_name, SINGLE_END_DSET); if (is_chunked) strcat(dset_name, "_chunked"); did = H5Dcreate2(file, dset_name, H5T_NATIVE_INT, sid, H5P_DEFAULT, plid, H5P_DEFAULT); CHECK(did, FAIL, "H5Dcreate2"); /* Initialize the data to be written to file */ for (i = 0; i < 2; i++) { for (j = 0; j < 3; j++) { for (k = 0; k < 6; k++) { da_buffer[i][j][k][0] = i * 100 + j * 10 + k; da_buffer[i][j][k][1] = i * 100 + j * 10 + k + 1; } } } ret = H5Dwrite(did, H5T_NATIVE_INT, sid, sid, H5P_DEFAULT, da_buffer); CHECK(ret, FAIL, "H5Dwrite"); ret = H5Dclose(did); CHECK(ret, FAIL, "H5Dclose"); /* ****** Case 1: ****** * Testing the full selection in the fastest-growing end */ did = H5Dopen2(file, dset_name, H5P_DEFAULT); CHECK(did, FAIL, "H5Dopen2"); /* Select the elements in the dataset */ elmts_numb = 12; ret = H5Sselect_elements(sid, H5S_SELECT_SET, elmts_numb, (const hsize_t *)da_elmts1); CHECK(ret, FAIL, "H5Sselect_elements"); /* Dataspace for memory buffer */ msid = H5Screate_simple(4, mem1_dims, mem1_dims); CHECK(msid, FAIL, "H5Screate_simple"); ret = H5Sselect_all(msid); CHECK(ret, FAIL, "H5Sselect_all"); ret = H5Dread(did, H5T_NATIVE_INT, msid, sid, H5P_DEFAULT, mem1_buffer); CHECK(ret, FAIL, "H5Dread"); ret = H5Dclose(did); CHECK(ret, FAIL, "H5Dclose"); ret = H5Sclose(msid); CHECK(ret, FAIL, "H5Sclose"); for (i = 0; i < 6; i++) for (j = 0; j < 2; j++) if (da_buffer[0][0][i][j] != mem1_buffer[0][0][i][j]) { TestErrPrintf("%u: Read different values than written at index 0,0,%d,%d\n", __LINE__, i, j); } /* ****** Case 2: ****** * Testing the full selection in the slowest-growing end */ did = H5Dopen2(file, dset_name, H5P_DEFAULT); CHECK(did, FAIL, "H5Dopen2"); /* Select the elements in the dataset */ elmts_numb = 6; ret = H5Sselect_elements(sid, H5S_SELECT_SET, elmts_numb, (const hsize_t *)da_elmts2); CHECK(ret, FAIL, "H5Sselect_elements"); /* Dataspace for memory buffer */ msid = H5Screate_simple(4, mem2_dims, mem2_dims); CHECK(msid, FAIL, "H5Screate_simple"); ret = H5Sselect_all(msid); CHECK(ret, FAIL, "H5Sselect_all"); ret = H5Dread(did, H5T_NATIVE_INT, msid, sid, H5P_DEFAULT, mem2_buffer); CHECK(ret, FAIL, "H5Dread"); ret = H5Dclose(did); CHECK(ret, FAIL, "H5Dclose"); ret = H5Sclose(msid); CHECK(ret, FAIL, "H5Sclose"); for (i = 0; i < 2; i++) for (j = 0; j < 3; j++) if (da_buffer[i][j][0][0] != mem2_buffer[i][j][0][0]) { TestErrPrintf("%u: Read different values than written at index %d,%d,0,0, da_buffer = %d, " "mem2_buffer = %d\n", __LINE__, i, j, da_buffer[i][j][0][0], mem2_buffer[i][j][0][0]); } /* ****** Case 3: ****** * Testing the full selection in the middle dimensions */ did = H5Dopen2(file, dset_name, H5P_DEFAULT); CHECK(did, FAIL, "H5Dopen2"); /* Select the elements in the dataset */ elmts_numb = 18; ret = H5Sselect_elements(sid, H5S_SELECT_SET, elmts_numb, (const hsize_t *)da_elmts3); CHECK(ret, FAIL, "H5Sselect_elements"); /* Dataspace for memory buffer */ msid = H5Screate_simple(4, mem3_dims, mem3_dims); CHECK(msid, FAIL, "H5Screate_simple"); ret = H5Sselect_all(msid); CHECK(ret, FAIL, "H5Sselect_all"); ret = H5Dread(did, H5T_NATIVE_INT, msid, sid, H5P_DEFAULT, mem3_buffer); CHECK(ret, FAIL, "H5Dread"); ret = H5Dclose(did); CHECK(ret, FAIL, "H5Dclose"); ret = H5Sclose(msid); CHECK(ret, FAIL, "H5Sclose"); for (i = 0; i < 3; i++) for (j = 0; j < 6; j++) if (da_buffer[0][i][j][0] != mem3_buffer[0][i][j][0]) { TestErrPrintf("%u: Read different values than written at index 0,%d,%d,0\n", __LINE__, i, j); } ret = H5Sclose(sid); CHECK(ret, FAIL, "H5Sclose"); ret = H5Pclose(plid); CHECK(ret, FAIL, "H5Pclose"); } /*********************************************************** ** ** test_singleEnd_selHyperslab(): Test full hyperslab selection ** of only one block. ** *************************************************************/ static void test_singleEnd_selHyperslab(hid_t file, bool is_chunked) { hid_t sid, did, msid; char dset_name[NAME_LEN]; /* Dataset name */ herr_t ret; /* Generic error return */ int i, j; hsize_t da_dims[4] = {2, 3, 6, 2}; /* For testing the full selection in the fastest-growing end */ int mem1_buffer[1][1][6][2]; hsize_t mem1_dims[4] = {1, 1, 6, 2}; hsize_t mem1_start[4] = {0, 0, 0, 0}; hsize_t mem1_count[4] = {1, 1, 1, 1}; hsize_t mem1_stride[4] = {1, 1, 1, 1}; hsize_t mem1_block[4] = {1, 1, 6, 2}; /* For testing the full selection in the slowest-growing end */ int mem2_buffer[2][3][1][1]; hsize_t mem2_dims[4] = {2, 3, 1, 1}; hsize_t mem2_start[4] = {0, 0, 0, 0}; hsize_t mem2_count[4] = {1, 1, 1, 1}; hsize_t mem2_stride[4] = {1, 1, 1, 1}; hsize_t mem2_block[4] = {2, 3, 1, 1}; /* For testing the full selection in the middle dimensions */ int mem3_buffer[1][3][6][1]; hsize_t mem3_dims[4] = {1, 3, 6, 1}; hsize_t mem3_start[4] = {0, 0, 0, 0}; hsize_t mem3_count[4] = {1, 1, 1, 1}; hsize_t mem3_stride[4] = {1, 1, 1, 1}; hsize_t mem3_block[4] = {1, 3, 6, 1}; /* Construct dataset's name */ memset(dset_name, 0, NAME_LEN); strcat(dset_name, SINGLE_END_DSET); if (is_chunked) strcat(dset_name, "_chunked"); /* Dataspace for the dataset in file */ sid = H5Screate_simple(4, da_dims, da_dims); CHECK(sid, FAIL, "H5Screate_simple"); /* ****** Case 1: ****** * Testing the full selection in the fastest-growing end */ did = H5Dopen2(file, dset_name, H5P_DEFAULT); CHECK(did, FAIL, "H5Dopen2"); /* Select the elements in the dataset */ ret = H5Sselect_hyperslab(sid, H5S_SELECT_SET, mem1_start, mem1_stride, mem1_count, mem1_block); CHECK(ret, FAIL, "H5Sselect_hyperslab"); /* Dataspace for memory buffer */ msid = H5Screate_simple(4, mem1_dims, mem1_dims); CHECK(msid, FAIL, "H5Screate_simple"); ret = H5Sselect_all(msid); CHECK(ret, FAIL, "H5Sselect_all"); ret = H5Dread(did, H5T_NATIVE_INT, msid, sid, H5P_DEFAULT, mem1_buffer); CHECK(ret, FAIL, "H5Dread"); ret = H5Dclose(did); CHECK(ret, FAIL, "H5Dclose"); ret = H5Sclose(msid); CHECK(ret, FAIL, "H5Sclose"); for (i = 0; i < 6; i++) for (j = 0; j < 2; j++) if (da_buffer[0][0][i][j] != mem1_buffer[0][0][i][j]) { TestErrPrintf("%u: Read different values than written at index 0,0,%d,%d\n", __LINE__, i, j); } /* ****** Case 2: ****** * Testing the full selection in the slowest-growing end */ did = H5Dopen2(file, dset_name, H5P_DEFAULT); CHECK(did, FAIL, "H5Dopen2"); /* Select the elements in the dataset */ ret = H5Sselect_hyperslab(sid, H5S_SELECT_SET, mem2_start, mem2_stride, mem2_count, mem2_block); CHECK(ret, FAIL, "H5Sselect_hyperslab"); /* Dataspace for memory buffer */ msid = H5Screate_simple(4, mem2_dims, mem2_dims); CHECK(msid, FAIL, "H5Screate_simple"); ret = H5Sselect_all(msid); CHECK(ret, FAIL, "H5Sselect_all"); ret = H5Dread(did, H5T_NATIVE_INT, msid, sid, H5P_DEFAULT, mem2_buffer); CHECK(ret, FAIL, "H5Dread"); ret = H5Dclose(did); CHECK(ret, FAIL, "H5Dclose"); ret = H5Sclose(msid); CHECK(ret, FAIL, "H5Sclose"); for (i = 0; i < 2; i++) for (j = 0; j < 3; j++) if (da_buffer[i][j][0][0] != mem2_buffer[i][j][0][0]) { TestErrPrintf("%u: Read different values than written at index %d,%d,0,0\n", __LINE__, i, j); } /* ****** Case 3: ****** * Testing the full selection in the middle dimensions */ did = H5Dopen2(file, dset_name, H5P_DEFAULT); CHECK(did, FAIL, "H5Dopen2"); /* Select the elements in the dataset */ ret = H5Sselect_hyperslab(sid, H5S_SELECT_SET, mem3_start, mem3_stride, mem3_count, mem3_block); CHECK(ret, FAIL, "H5Sselect_hyperslab"); /* Dataspace for memory buffer */ msid = H5Screate_simple(4, mem3_dims, mem3_dims); CHECK(msid, FAIL, "H5Screate_simple"); ret = H5Sselect_all(msid); CHECK(ret, FAIL, "H5Sselect_all"); ret = H5Dread(did, H5T_NATIVE_INT, msid, sid, H5P_DEFAULT, mem3_buffer); CHECK(ret, FAIL, "H5Dread"); ret = H5Dclose(did); CHECK(ret, FAIL, "H5Dclose"); ret = H5Sclose(msid); CHECK(ret, FAIL, "H5Sclose"); for (i = 0; i < 3; i++) for (j = 0; j < 6; j++) if (da_buffer[0][i][j][0] != mem3_buffer[0][i][j][0]) { TestErrPrintf("%u: Read different values than written at index 0,%d,%d,0\n", __LINE__, i, j); } ret = H5Sclose(sid); CHECK(ret, FAIL, "H5Sclose"); } /*********************************************************** ** ** test_multiple_end(): Test full hyperslab selection of ** multiple blocks. ** *************************************************************/ static void test_multiple_ends(hid_t file, bool is_chunked) { hid_t sid, plid, did, msid; char dset_name[NAME_LEN]; /* Dataset name */ herr_t ret; /* Generic error return */ int i, j, k, l, m, n, p; hsize_t da_dims[8] = {4, 5, 3, 4, 2, 3, 6, 2}; hsize_t da_chunksize[8] = {1, 5, 3, 2, 2, 3, 3, 2}; struct { int arr[4][5][3][4][2][3][6][2]; } *data_buf = NULL; /* For testing the full selections in the fastest-growing end and in the middle dimensions */ struct { int arr[1][1][1][4][2][1][6][2]; } *mem1_buffer = NULL; hsize_t mem1_dims[8] = {1, 1, 1, 4, 2, 1, 6, 2}; hsize_t mem1_start[8] = {0, 0, 0, 0, 0, 0, 0, 0}; hsize_t mem1_count[8] = {1, 1, 1, 1, 1, 1, 1, 1}; hsize_t mem1_stride[8] = {1, 1, 1, 1, 1, 1, 1, 1}; hsize_t mem1_block[8] = {1, 1, 1, 4, 2, 1, 6, 2}; /* For testing the full selections in the slowest-growing end and in the middle dimensions */ struct { int arr[4][5][1][4][2][1][1][1]; } *mem2_buffer = NULL; hsize_t mem2_dims[8] = {4, 5, 1, 4, 2, 1, 1, 1}; hsize_t mem2_start[8] = {0, 0, 0, 0, 0, 0, 0, 0}; hsize_t mem2_count[8] = {1, 1, 1, 1, 1, 1, 1, 1}; hsize_t mem2_stride[8] = {1, 1, 1, 1, 1, 1, 1, 1}; hsize_t mem2_block[8] = {4, 5, 1, 4, 2, 1, 1, 1}; /* For testing two unadjacent full selections in the middle dimensions */ struct { int arr[1][5][3][1][1][3][6][1]; } *mem3_buffer = NULL; hsize_t mem3_dims[8] = {1, 5, 3, 1, 1, 3, 6, 1}; hsize_t mem3_start[8] = {0, 0, 0, 0, 0, 0, 0, 0}; hsize_t mem3_count[8] = {1, 1, 1, 1, 1, 1, 1, 1}; hsize_t mem3_stride[8] = {1, 1, 1, 1, 1, 1, 1, 1}; hsize_t mem3_block[8] = {1, 5, 3, 1, 1, 3, 6, 1}; /* For testing the full selections in the fastest-growing end and the slowest-growing end */ struct { int arr[4][5][1][1][1][1][6][2]; } *mem4_buffer = NULL; hsize_t mem4_dims[8] = {4, 5, 1, 1, 1, 1, 6, 2}; hsize_t mem4_start[8] = {0, 0, 0, 0, 0, 0, 0, 0}; hsize_t mem4_count[8] = {1, 1, 1, 1, 1, 1, 1, 1}; hsize_t mem4_stride[8] = {1, 1, 1, 1, 1, 1, 1, 1}; hsize_t mem4_block[8] = {4, 5, 1, 1, 1, 1, 6, 2}; /* For testing the full selections in the fastest-growing end and slowest-growing end, * also in the middle dimensions */ struct { int arr[4][5][1][4][2][1][6][2]; } *mem5_buffer = NULL; hsize_t mem5_dims[8] = {4, 5, 1, 4, 2, 1, 6, 2}; hsize_t mem5_start[8] = {0, 0, 0, 0, 0, 0, 0, 0}; hsize_t mem5_count[8] = {1, 1, 1, 1, 1, 1, 1, 1}; hsize_t mem5_stride[8] = {1, 1, 1, 1, 1, 1, 1, 1}; hsize_t mem5_block[8] = {4, 5, 1, 4, 2, 1, 6, 2}; /* Initialize dynamic arrays */ data_buf = calloc(1, sizeof(*data_buf)); CHECK_PTR(data_buf, "calloc"); mem1_buffer = calloc(1, sizeof(*mem1_buffer)); CHECK_PTR(data_buf, "calloc"); mem2_buffer = calloc(1, sizeof(*mem2_buffer)); CHECK_PTR(data_buf, "calloc"); mem3_buffer = calloc(1, sizeof(*mem3_buffer)); CHECK_PTR(data_buf, "calloc"); mem4_buffer = calloc(1, sizeof(*mem4_buffer)); CHECK_PTR(data_buf, "calloc"); mem5_buffer = calloc(1, sizeof(*mem5_buffer)); CHECK_PTR(data_buf, "calloc"); /* Create and write the dataset */ sid = H5Screate_simple(8, da_dims, da_dims); CHECK(sid, FAIL, "H5Screate_simple"); plid = H5Pcreate(H5P_DATASET_CREATE); CHECK(plid, FAIL, "H5Pcreate"); if (is_chunked) { ret = H5Pset_chunk(plid, 8, da_chunksize); CHECK(ret, FAIL, "H5Pset_chunk"); } /* Construct dataset's name */ memset(dset_name, 0, NAME_LEN); strcat(dset_name, MULTI_ENDS_SEL_HYPER_DSET); if (is_chunked) strcat(dset_name, "_chunked"); did = H5Dcreate2(file, dset_name, H5T_NATIVE_INT, sid, H5P_DEFAULT, plid, H5P_DEFAULT); CHECK(did, FAIL, "H5Dcreate2"); for (i = 0; i < 4; i++) for (j = 0; j < 5; j++) for (k = 0; k < 3; k++) for (l = 0; l < 4; l++) for (m = 0; m < 2; m++) for (n = 0; n < 3; n++) for (p = 0; p < 6; p++) { data_buf->arr[i][j][k][l][m][n][p][0] = i * 1000000 + j * 100000 + k * 10000 + l * 1000 + m * 100 + n * 10 + p; data_buf->arr[i][j][k][l][m][n][p][1] = i * 1000000 + j * 100000 + k * 10000 + l * 1000 + m * 100 + n * 10 + p + 1; } ret = H5Dwrite(did, H5T_NATIVE_INT, sid, sid, H5P_DEFAULT, data_buf); CHECK(ret, FAIL, "H5Dwrite"); ret = H5Dclose(did); CHECK(ret, FAIL, "H5Dclose"); /* ****** Case 1: ****** * Testing the full selections in the fastest-growing end and in the middle dimensions*/ did = H5Dopen2(file, dset_name, H5P_DEFAULT); CHECK(did, FAIL, "H5Dopen2"); /* Select the elements in the dataset */ ret = H5Sselect_hyperslab(sid, H5S_SELECT_SET, mem1_start, mem1_stride, mem1_count, mem1_block); CHECK(ret, FAIL, "H5Sselect_hyperslab"); msid = H5Screate_simple(8, mem1_dims, mem1_dims); CHECK(msid, FAIL, "H5Screate_simple"); ret = H5Sselect_all(msid); CHECK(ret, FAIL, "H5Sselect_all"); ret = H5Dread(did, H5T_NATIVE_INT, msid, sid, H5P_DEFAULT, mem1_buffer); CHECK(ret, FAIL, "H5Dread"); ret = H5Dclose(did); CHECK(ret, FAIL, "H5Dclose"); ret = H5Sclose(msid); CHECK(ret, FAIL, "H5Sclose"); for (i = 0; i < 4; i++) for (j = 0; j < 2; j++) for (k = 0; k < 6; k++) for (l = 0; l < 2; l++) if (data_buf->arr[0][0][0][i][j][0][k][l] != mem1_buffer->arr[0][0][0][i][j][0][k][l]) { TestErrPrintf("%u: Read different values than written at index 0,0,0,%d,%d,0,%d,%d\n", __LINE__, i, j, k, l); } /* ****** Case 2: ****** * Testing the full selections in the slowest-growing end and in the middle dimensions*/ did = H5Dopen2(file, dset_name, H5P_DEFAULT); CHECK(did, FAIL, "H5Dopen2"); /* Select the elements in the dataset */ ret = H5Sselect_hyperslab(sid, H5S_SELECT_SET, mem2_start, mem2_stride, mem2_count, mem2_block); CHECK(ret, FAIL, "H5Sselect_hyperslab"); msid = H5Screate_simple(8, mem2_dims, mem2_dims); CHECK(msid, FAIL, "H5Screate_simple"); ret = H5Sselect_all(msid); CHECK(ret, FAIL, "H5Sselect_all"); ret = H5Dread(did, H5T_NATIVE_INT, msid, sid, H5P_DEFAULT, mem2_buffer); CHECK(ret, FAIL, "H5Dread"); ret = H5Dclose(did); CHECK(ret, FAIL, "H5Dclose"); ret = H5Sclose(msid); CHECK(ret, FAIL, "H5Sclose"); for (i = 0; i < 4; i++) for (j = 0; j < 5; j++) for (k = 0; k < 4; k++) for (l = 0; l < 2; l++) if (data_buf->arr[i][j][0][k][l][0][0][0] != mem2_buffer->arr[i][j][0][k][l][0][0][0]) { TestErrPrintf("%u: Read different values than written at index %d,%d,0,%d,%d,0,0,0\n", __LINE__, i, j, k, l); } /* ****** Case 3: ****** * Testing two unadjacent full selections in the middle dimensions */ did = H5Dopen2(file, dset_name, H5P_DEFAULT); CHECK(did, FAIL, "H5Dopen2"); /* Select the elements in the dataset */ ret = H5Sselect_hyperslab(sid, H5S_SELECT_SET, mem3_start, mem3_stride, mem3_count, mem3_block); CHECK(ret, FAIL, "H5Sselect_hyperslab"); msid = H5Screate_simple(8, mem3_dims, mem3_dims); CHECK(msid, FAIL, "H5Screate_simple"); ret = H5Sselect_all(msid); CHECK(ret, FAIL, "H5Sselect_all"); ret = H5Dread(did, H5T_NATIVE_INT, msid, sid, H5P_DEFAULT, mem3_buffer); CHECK(ret, FAIL, "H5Dread"); ret = H5Dclose(did); CHECK(ret, FAIL, "H5Dclose"); ret = H5Sclose(msid); CHECK(ret, FAIL, "H5Sclose"); for (i = 0; i < 5; i++) for (j = 0; j < 3; j++) for (k = 0; k < 3; k++) for (l = 0; l < 6; l++) if (data_buf->arr[0][i][j][0][0][k][l][0] != mem3_buffer->arr[0][i][j][0][0][k][l][0]) { TestErrPrintf("%u: Read different values than written at index 0,%d,%d,0,0,%d,%d,0\n", __LINE__, i, j, k, l); } /* ****** Case 4: ****** * Testing the full selections in the fastest-growing end and the slowest-growing end */ did = H5Dopen2(file, dset_name, H5P_DEFAULT); CHECK(did, FAIL, "H5Dopen2"); /* Select the elements in the dataset */ ret = H5Sselect_hyperslab(sid, H5S_SELECT_SET, mem4_start, mem4_stride, mem4_count, mem4_block); CHECK(ret, FAIL, "H5Sselect_hyperslab"); msid = H5Screate_simple(8, mem4_dims, mem4_dims); CHECK(msid, FAIL, "H5Screate_simple"); ret = H5Sselect_all(msid); CHECK(ret, FAIL, "H5Sselect_all"); ret = H5Dread(did, H5T_NATIVE_INT, msid, sid, H5P_DEFAULT, mem4_buffer); CHECK(ret, FAIL, "H5Dread"); ret = H5Dclose(did); CHECK(ret, FAIL, "H5Dclose"); ret = H5Sclose(msid); CHECK(ret, FAIL, "H5Sclose"); for (i = 0; i < 4; i++) for (j = 0; j < 5; j++) for (k = 0; k < 6; k++) for (l = 0; l < 2; l++) if (data_buf->arr[i][j][0][0][0][0][k][l] != mem4_buffer->arr[i][j][0][0][0][0][k][l]) { TestErrPrintf("%u: Read different values than written at index %d,%d,0,0,0,0,%d,%d\n", __LINE__, i, j, k, l); } /* ****** Case 5: ****** * Testing the full selections in the fastest-growing end and the slowest-growing end, * and also in the middle dimensions */ did = H5Dopen2(file, dset_name, H5P_DEFAULT); CHECK(did, FAIL, "H5Dopen2"); /* Select the elements in the dataset */ ret = H5Sselect_hyperslab(sid, H5S_SELECT_SET, mem5_start, mem5_stride, mem5_count, mem5_block); CHECK(ret, FAIL, "H5Sselect_hyperslab"); msid = H5Screate_simple(8, mem5_dims, mem5_dims); CHECK(msid, FAIL, "H5Screate_simple"); ret = H5Sselect_all(msid); CHECK(ret, FAIL, "H5Sselect_all"); ret = H5Dread(did, H5T_NATIVE_INT, msid, sid, H5P_DEFAULT, mem5_buffer); CHECK(ret, FAIL, "H5Dread"); ret = H5Dclose(did); CHECK(ret, FAIL, "H5Dclose"); ret = H5Sclose(msid); CHECK(ret, FAIL, "H5Sclose"); for (i = 0; i < 4; i++) for (j = 0; j < 5; j++) for (k = 0; k < 4; k++) for (l = 0; l < 2; l++) for (m = 0; m < 6; m++) for (n = 0; n < 2; n++) if (data_buf->arr[i][j][0][k][l][0][m][n] != mem5_buffer->arr[i][j][0][k][l][0][m][n]) { TestErrPrintf( "%u: Read different values than written at index %d,%d,0,%d,%d,0,%d,%d\n", __LINE__, i, j, k, l, m, n); } ret = H5Sclose(sid); CHECK(ret, FAIL, "H5Sclose"); ret = H5Pclose(plid); CHECK(ret, FAIL, "H5Pclose"); free(data_buf); free(mem1_buffer); free(mem2_buffer); free(mem3_buffer); free(mem4_buffer); free(mem5_buffer); } /**************************************************************** ** ** test_coords(): Main testing routine. ** ****************************************************************/ void test_coords(void) { hid_t fid; bool is_chunk[2] = {true, false}; int i; herr_t ret; /* Generic error return */ fid = H5Fcreate(FILENAME, H5F_ACC_TRUNC, H5P_DEFAULT, H5P_DEFAULT); CHECK(fid, FAIL, "H5Fcreate"); for (i = 0; i < 2; i++) { /* Test H5Sselect_elements with selection of one block of data */ test_singleEnd_selElements(fid, is_chunk[i]); /* Test H5Sselect_hyperslab with selection of one block of data */ test_singleEnd_selHyperslab(fid, is_chunk[i]); /* Test H5Sselect_hyperslab with selection of multiple blocks of data */ test_multiple_ends(fid, is_chunk[i]); } ret = H5Fclose(fid); CHECK(ret, FAIL, "H5Fclose"); } /*------------------------------------------------------------------------- * Function: cleanup_coords * * Purpose: Cleanup temporary test files * * Return: none * *------------------------------------------------------------------------- */ void cleanup_coords(void) { H5E_BEGIN_TRY { H5Fdelete(FILENAME, H5P_DEFAULT); } H5E_END_TRY }