diff options
Diffstat (limited to 'test/dtransform.c')
-rw-r--r-- | test/dtransform.c | 1037 |
1 files changed, 574 insertions, 463 deletions
diff --git a/test/dtransform.c b/test/dtransform.c index 6f7e8a4..1f6d48b 100644 --- a/test/dtransform.c +++ b/test/dtransform.c @@ -6,286 +6,344 @@ * 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://support.hdfgroup.org/ftp/HDF5/releases. * + * 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. * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * */ #include "h5test.h" -#define ROWS 12 -#define COLS 18 +#define ROWS 12 +#define COLS 18 #define FLOAT_TOL 0.0001F static int init_test(hid_t file_id); -static int test_copy(const hid_t dxpl_id_c_to_f_copy, const hid_t dxpl_id_polynomial_copy); -static int test_trivial(const hid_t dxpl_id_simple); -static int test_poly(const hid_t dxpl_id_polynomial); +static int test_copy(hid_t dxpl_id_c_to_f_copy, hid_t dxpl_id_polynomial_copy); +static int test_trivial(hid_t dxpl_id_simple); +static int test_poly(hid_t dxpl_id_polynomial); static int test_specials(hid_t file); static int test_set(void); -static int test_getset(const hid_t dxpl_id_simple); +static int test_getset(hid_t dxpl_id_simple); /* These are needed for multiple tests, so are declared here globally and are init'ed in init_test */ -hid_t dset_id_int = -1; -hid_t dset_id_float = -1; -hid_t dset_id_int_chunk = -1; +hid_t dset_id_int = -1; +hid_t dset_id_float = -1; +hid_t dset_id_int_chunk = -1; hid_t dset_id_float_chunk = -1; - -const float windchillFfloat[ROWS][COLS] = - { {36.0f, 31.0f, 25.0f, 19.0f, 13.0f, 7.0f, 1.0f, -5.0f, -11.0f, -16.0f, -22.0f, -28.0f, -34.0f, -40.0f, -46.0f, -52.0f, -57.0f, -63.0f}, - {34.0f, 27.0f, 21.0f, 15.0f, 9.0f, 3.0f, -4.0f, -10.0f, -16.0f, -22.0f, -28.0f, -35.0f, -41.0f, -47.0f, -53.0f, -59.0f, -66.0f, -72.0f} , - {32.0f, 25.0f, 19.0f, 13.0f, 6.0f, 0.0f, -7.0f, -13.0f, -19.0f, -26.0f, -32.0f, -39.0f, -45.0f, -51.0f, -58.0f, -64.0f, -71.0f, -77.0f}, - {30.0f, 24.0f, 17.0f, 11.0f, 4.0f, -2.0f, -9.0f, -15.0f, -22.0f, -29.0f, -35.0f, -42.0f, -48.0f, -55.0f, -61.0f, -68.0f, -74.0f, -81.0f}, - {29.0f, 23.0f, 16.0f, 9.0f, 3.0f, -4.0f, -11.0f, -17.0f, -24.0f, -31.0f, -37.0f, -44.0f, -51.0f, -58.0f, -64.0f, -71.0f, -78.0f, -84.0f}, - {28.0f, 22.0f, 15.0f, 8.0f, 1.0f, -5.0f, -12.0f, -19.0f, -26.0f, -33.0f, -39.0f, -46.0f, -53.0f, -60.0f, -67.0f, -73.0f, -80.0f, -87.0f}, - {28.0f, 21.0f, 14.0f, 7.0f, 0.0f, -7.0f, -14.0f, -21.0f, -27.0f, -34.0f, -41.0f, -48.0f, -55.0f, -62.0f, -69.0f, -76.0f, -82.0f, -89.0f}, - {27.0f, 20.0f, 13.0f, 6.0f, -1.0f, -8.0f, -15.0f, -22.0f, -29.0f, -36.0f, -43.0f, -50.0f, -57.0f, -64.0f, -71.0f, -78.0f, -84.0f, -91.0f}, - {26.0f, 19.0f, 12.0f, 5.0f, -2.0f, -9.0f, -16.0f, -23.0f, -30.0f, -37.0f, -44.0f, -51.0f, -58.0f, -65.0f, -72.0f, -79.0f, -86.0f, -93.0f}, - {26.0f, 19.0f, 12.0f, 4.0f, -3.0f, -10.0f, -17.0f, -24.0f, -31.0f, -38.0f, -45.0f, -52.0f, -60.0f, -67.0f, -74.0f, -81.0f, -88.0f, -95.0f}, - {25.0f, 18.0f, 11.0f, 4.0f, -3.0f, -11.0f, -18.0f, -25.0f, -32.0f, -39.0f, -46.0f, -54.0f, -61.0f, -68.0f, -75.0f, -82.0f, -89.0f, -97.0f}, - {25.0f, 17.0f, 10.0f, 3.0f, -4.0f, -11.0f, -19.0f, -26.0f, -33.0f, -40.0f, -48.0f, -55.0f, -62.0f, -69.0f, -76.0f, -84.0f, -91.0f, -98.0f} - }; - -const int transformData[ROWS][COLS] = - { {36, 31, 25, 19, 13, 7, 1, 5, 11, 16, 22, 28, 34, 40, 46, 52, 57, 63 }, - {34, 27, 21, 15, 9, 3, 4, 10, 16, 22, 28, 35, 41, 47, 53, 59, 66, 1 } , - {32, 25, 19, 13, 6, 2, 7, 13, 19, 26, 32, 39, 45, 51, 58, 64, 71, 5 }, - {30, 24, 17, 11, 4, 2, 9, 15, 22, 29, 35, 42, 48, 55, 61, 68, 2, 9 }, - {29, 23, 16, 9, 3, 4, 11, 17, 24, 31, 37, 44, 51, 58, 64, 71, 6, 12 }, - {28, 22, 15, 8, 1, 5, 12, 19, 26, 33, 39, 46, 53, 60, 67, 1, 8, 15 }, - {28, 21, 14, 7, 6, 7, 14, 21, 27, 34, 41, 48, 55, 62, 69, 4, 10, 17 }, - {27, 20, 13, 6, 1, 8, 15, 22, 29, 36, 43, 50, 57, 64, 71, 6, 12, 19 }, - {26, 19, 12, 5, 2, 9, 16, 23, 30, 37, 44, 51, 58, 65, 5, 7, 14, 21 }, - {26, 19, 12, 4, 3, 10, 17, 24, 31, 38, 45, 52, 60, 67, 2, 9, 16, 23}, - {25, 18, 11, 4, 3, 11, 18, 25, 32, 39, 46, 54, 61, 68, 3, 10, 17, 25}, - {25, 17, 10, 3, 4, 11, 19, 26, 33, 40, 48, 55, 62, 69, 4, 12, 19, 26} - }; - -#define UCOMPARE(TYPE,VAR1,VAR2,TOL) \ -{ \ - size_t i,j; \ - \ - for(i=0; i<ROWS; i++) \ - for(j=0; j<COLS; j++) \ - { \ - if(!( (((VAR1)[i][j] >= (TYPE)((VAR2)[i][j])) && ( ((VAR1)[i][j] - TOL) < (TYPE)((VAR2)[i][j]))) || ( ((VAR1)[i][j] <= (TYPE)((VAR2)[i][j])) && ( ((VAR1)[i][j] + TOL) > (TYPE)((VAR2)[i][j]))))) \ - { \ - H5_FAILED(); \ - HDfprintf(stderr, " ERROR: Conversion failed to match computed data\n"); \ - goto error; \ - } \ - } \ - PASSED(); \ -} - -#define COMPARE(TYPE,VAR1,VAR2,TOL) \ -{ \ - size_t i,j; \ - \ - for(i=0; i<ROWS; i++) \ - for(j=0; j<COLS; j++) \ - { \ - if( !(((VAR1)[i][j] <= ((TYPE)(VAR2)[i][j] + TOL)) && ((VAR1)[i][j] >= ((TYPE)(VAR2)[i][j] - TOL))) ) \ - { \ - H5_FAILED(); \ - HDfprintf(stderr, " ERROR: Conversion failed to match computed data\n"); \ - goto error; \ - } \ - } \ - PASSED(); \ -} - -#define COMPARE_INT(VAR1,VAR2) \ -{ \ - size_t i,j; \ - \ - for(i=0; i<ROWS; i++) \ - for(j=0; j<COLS; j++) \ - { \ - if( (VAR1)[i][j] != (VAR2)[i][j] ) \ - { \ - H5_FAILED(); \ - HDfprintf(stderr, " ERROR: data failed to match computed data\n"); \ - goto error; \ - } \ - } \ -} - -#define TEST_TYPE_CONTIG(XFORM, TYPE, HDF_TYPE, TEST_STR, COMPARE_DATA, SIGNED) \ -{ \ - TYPE array[ROWS][COLS]; \ - const char* f_to_c = "(5/9.0)*(x-32)"; \ - /* utrans is a transform for unsigned types: no negative numbers involved and results are < 255 to fit into uchar */ \ - const char* utrans = "((x+100)/4)*3"; \ - \ - hid_t dataspace, dxpl_id_f_to_c, dxpl_id_utrans, dset, dset_nn, dt_nn; \ - H5T_order_t order; \ - hsize_t dim[2] = {ROWS, COLS}; \ - \ - if((dataspace = H5Screate_simple(2, dim, NULL)) < 0) TEST_ERROR; \ - if((dset = H5Dcreate2(file_id, "/transformtest_"TEST_STR, HDF_TYPE, dataspace, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) TEST_ERROR; \ - \ - if((dt_nn = H5Tcopy(HDF_TYPE)) < 0) TEST_ERROR \ - if((order = H5Tget_order(dt_nn)) == H5T_ORDER_ERROR) TEST_ERROR \ - if(H5Tset_order(dt_nn, order == H5T_ORDER_LE ? H5T_ORDER_BE : H5T_ORDER_LE) < 0) TEST_ERROR \ - if((dset_nn = H5Dcreate2(file_id, "/nonnative_transformtest_"TEST_STR, dt_nn, dataspace, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) TEST_ERROR \ - if(H5Tclose(dt_nn) < 0) TEST_ERROR \ - \ - if(SIGNED) \ - { \ - if((dxpl_id_f_to_c = H5Pcreate(H5P_DATASET_XFER)) < 0) TEST_ERROR; \ - if(H5Pset_data_transform(dxpl_id_f_to_c, f_to_c) < 0) TEST_ERROR; \ - if(H5Dwrite(dset, H5T_NATIVE_FLOAT, H5S_ALL, H5S_ALL, dxpl_id_f_to_c, windchillFfloat) < 0) TEST_ERROR; \ - if(H5Dwrite(dset_nn, H5T_NATIVE_FLOAT, H5S_ALL, H5S_ALL, dxpl_id_f_to_c, windchillFfloat) < 0) TEST_ERROR; \ - if(H5Pclose(dxpl_id_f_to_c) < 0) TEST_ERROR; \ - } \ - else \ - { \ - if((dxpl_id_utrans = H5Pcreate(H5P_DATASET_XFER)) < 0) TEST_ERROR; \ - if(H5Pset_data_transform(dxpl_id_utrans, utrans) < 0) TEST_ERROR; \ - if(H5Dwrite(dset, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, dxpl_id_utrans, transformData) < 0) TEST_ERROR; \ - if(H5Dwrite(dset_nn, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, dxpl_id_utrans, transformData) < 0) TEST_ERROR; \ - if(H5Pclose(dxpl_id_utrans) < 0) TEST_ERROR; \ - } \ - \ - \ - TESTING("contiguous, no data type conversion ("TEST_STR"->"TEST_STR")") \ - \ - if(H5Dread(dset, HDF_TYPE, H5S_ALL, H5S_ALL, XFORM, array) < 0) TEST_ERROR; \ - if(SIGNED) \ - COMPARE(TYPE, array, COMPARE_DATA, 2) \ - else \ - UCOMPARE(TYPE, array, COMPARE_DATA, 4) \ - \ - TESTING("contiguous, byte order conversion ("TEST_STR"->"TEST_STR")") \ - \ - if(H5Dread(dset_nn, HDF_TYPE, H5S_ALL, H5S_ALL, XFORM, array) < 0) TEST_ERROR; \ - if(SIGNED) \ - COMPARE(TYPE, array, COMPARE_DATA, 2) \ - else \ - UCOMPARE(TYPE, array, COMPARE_DATA, 4) \ - \ - if(SIGNED) \ - { \ - TESTING("contiguous, with type conversion (float->"TEST_STR")") \ - \ - if(H5Dread(dset_id_float, HDF_TYPE, H5S_ALL, H5S_ALL, XFORM, array) < 0) TEST_ERROR; \ - COMPARE(TYPE, array, COMPARE_DATA, 2) \ - } \ - \ - if(H5Dclose(dset) < 0) TEST_ERROR; \ - if(H5Sclose(dataspace) < 0) TEST_ERROR; \ -} - -#define TEST_TYPE_CHUNK(XFORM, TYPE, HDF_TYPE, TEST_STR, COMPARE_DATA, SIGNED) \ -{ \ - TYPE array[ROWS][COLS]; \ - const char* f_to_c = "(5/9.0)*(x-32)"; \ - /* utrans is a transform for unsigned types: no negative numbers involved and results are < 255 to fit into uchar */ \ - const char* utrans = "((x+100)/4)*3"; \ - \ - hid_t dataspace, dxpl_id_f_to_c, dxpl_id_utrans, cparms, memspace, dset_chunk, filespace; \ - hsize_t dim[2] = {ROWS, COLS}; \ - hsize_t offset[2] = {0, 0}; \ - \ - \ - if((dataspace = H5Screate_simple(2, dim, NULL)) < 0) TEST_ERROR; \ - \ - if((cparms = H5Pcreate(H5P_DATASET_CREATE)) < 0) TEST_ERROR; \ - if(H5Pset_chunk(cparms, 2, dim) < 0) TEST_ERROR; \ - \ - if((dset_chunk = H5Dcreate2(file_id, "/transformtest_chunk_"TEST_STR, HDF_TYPE, dataspace, H5P_DEFAULT, cparms, H5P_DEFAULT)) < 0) TEST_ERROR; \ - if((filespace = H5Dget_space(dset_chunk)) < 0) TEST_ERROR \ - if((memspace = H5Screate_simple(2, dim, NULL)) < 0) TEST_ERROR \ - if(H5Sselect_hyperslab(filespace, H5S_SELECT_SET, offset, NULL, dim, NULL) < 0) TEST_ERROR; \ - \ - if(SIGNED) \ - { \ - if((dxpl_id_f_to_c = H5Pcreate(H5P_DATASET_XFER)) < 0) TEST_ERROR; \ - if(H5Pset_data_transform(dxpl_id_f_to_c, f_to_c) < 0) TEST_ERROR; \ - if(H5Dwrite(dset_chunk, H5T_NATIVE_FLOAT, dataspace, filespace, dxpl_id_f_to_c, windchillFfloat) < 0) TEST_ERROR; \ - if(H5Pclose(dxpl_id_f_to_c) < 0) TEST_ERROR; \ - } \ - else \ - { \ - if((dxpl_id_utrans = H5Pcreate(H5P_DATASET_XFER)) < 0) TEST_ERROR; \ - if(H5Pset_data_transform(dxpl_id_utrans, utrans) < 0) TEST_ERROR; \ - if(H5Dwrite(dset_chunk, H5T_NATIVE_INT, dataspace, filespace, dxpl_id_utrans, transformData) < 0) TEST_ERROR; \ - if(H5Pclose(dxpl_id_utrans) < 0) TEST_ERROR; \ - } \ - \ - \ - TESTING("chunked, no data type conversion ("TEST_STR"->"TEST_STR")") \ - \ - if(H5Dread(dset_chunk, HDF_TYPE, memspace, filespace, XFORM, array) < 0) TEST_ERROR; \ - if(SIGNED) \ - COMPARE(TYPE, array, COMPARE_DATA, 2) \ - else \ - UCOMPARE(TYPE, array, COMPARE_DATA, 4) \ - \ - if(SIGNED) \ - { \ - TESTING("chunked, with type conversion (float->"TEST_STR")") \ - \ - if(H5Dread(dset_id_float_chunk, HDF_TYPE, memspace, filespace, XFORM, array) < 0) TEST_ERROR; \ - COMPARE(TYPE, array, COMPARE_DATA, 2) \ - } \ - \ - \ - if(H5Pclose(cparms) < 0) TEST_ERROR; \ - if(H5Dclose(dset_chunk) < 0) TEST_ERROR; \ - if(H5Sclose(dataspace) < 0) TEST_ERROR; \ - if(H5Sclose(memspace) < 0) TEST_ERROR; \ -} - -#define INVALID_SET_TEST(TRANSFORM) \ -{ \ - if(H5Pset_data_transform(dxpl_id, TRANSFORM) < 0) \ - { \ - PASSED(); \ - } \ - else \ - { \ - H5_FAILED(); \ - HDfprintf(stderr, " ERROR: Data transform allowed invalid TRANSFORM transform to be set\n"); \ - goto error; \ - } \ -} - -int main(void) +const float windchillFfloat[ROWS][COLS] = { + {36.0F, 31.0F, 25.0F, 19.0F, 13.0F, 7.0F, 1.0F, -5.0F, -11.0F, -16.0F, -22.0F, -28.0F, -34.0F, -40.0F, + -46.0F, -52.0F, -57.0F, -63.0F}, + {34.0F, 27.0F, 21.0F, 15.0F, 9.0F, 3.0F, -4.0F, -10.0F, -16.0F, -22.0F, -28.0F, -35.0F, -41.0F, -47.0F, + -53.0F, -59.0F, -66.0F, -72.0F}, + {32.0F, 25.0F, 19.0F, 13.0F, 6.0F, 0.0F, -7.0F, -13.0F, -19.0F, -26.0F, -32.0F, -39.0F, -45.0F, -51.0F, + -58.0F, -64.0F, -71.0F, -77.0F}, + {30.0F, 24.0F, 17.0F, 11.0F, 4.0F, -2.0F, -9.0F, -15.0F, -22.0F, -29.0F, -35.0F, -42.0F, -48.0F, -55.0F, + -61.0F, -68.0F, -74.0F, -81.0F}, + {29.0F, 23.0F, 16.0F, 9.0F, 3.0F, -4.0F, -11.0F, -17.0F, -24.0F, -31.0F, -37.0F, -44.0F, -51.0F, -58.0F, + -64.0F, -71.0F, -78.0F, -84.0F}, + {28.0F, 22.0F, 15.0F, 8.0F, 1.0F, -5.0F, -12.0F, -19.0F, -26.0F, -33.0F, -39.0F, -46.0F, -53.0F, -60.0F, + -67.0F, -73.0F, -80.0F, -87.0F}, + {28.0F, 21.0F, 14.0F, 7.0F, 0.0F, -7.0F, -14.0F, -21.0F, -27.0F, -34.0F, -41.0F, -48.0F, -55.0F, -62.0F, + -69.0F, -76.0F, -82.0F, -89.0F}, + {27.0F, 20.0F, 13.0F, 6.0F, -1.0F, -8.0F, -15.0F, -22.0F, -29.0F, -36.0F, -43.0F, -50.0F, -57.0F, -64.0F, + -71.0F, -78.0F, -84.0F, -91.0F}, + {26.0F, 19.0F, 12.0F, 5.0F, -2.0F, -9.0F, -16.0F, -23.0F, -30.0F, -37.0F, -44.0F, -51.0F, -58.0F, -65.0F, + -72.0F, -79.0F, -86.0F, -93.0F}, + {26.0F, 19.0F, 12.0F, 4.0F, -3.0F, -10.0F, -17.0F, -24.0F, -31.0F, -38.0F, -45.0F, -52.0F, -60.0F, -67.0F, + -74.0F, -81.0F, -88.0F, -95.0F}, + {25.0F, 18.0F, 11.0F, 4.0F, -3.0F, -11.0F, -18.0F, -25.0F, -32.0F, -39.0F, -46.0F, -54.0F, -61.0F, -68.0F, + -75.0F, -82.0F, -89.0F, -97.0F}, + {25.0F, 17.0F, 10.0F, 3.0F, -4.0F, -11.0F, -19.0F, -26.0F, -33.0F, -40.0F, -48.0F, -55.0F, -62.0F, -69.0F, + -76.0F, -84.0F, -91.0F, -98.0F}}; + +const int transformData[ROWS][COLS] = {{36, 31, 25, 19, 13, 7, 1, 5, 11, 16, 22, 28, 34, 40, 46, 52, 57, 63}, + {34, 27, 21, 15, 9, 3, 4, 10, 16, 22, 28, 35, 41, 47, 53, 59, 66, 1}, + {32, 25, 19, 13, 6, 2, 7, 13, 19, 26, 32, 39, 45, 51, 58, 64, 71, 5}, + {30, 24, 17, 11, 4, 2, 9, 15, 22, 29, 35, 42, 48, 55, 61, 68, 2, 9}, + {29, 23, 16, 9, 3, 4, 11, 17, 24, 31, 37, 44, 51, 58, 64, 71, 6, 12}, + {28, 22, 15, 8, 1, 5, 12, 19, 26, 33, 39, 46, 53, 60, 67, 1, 8, 15}, + {28, 21, 14, 7, 6, 7, 14, 21, 27, 34, 41, 48, 55, 62, 69, 4, 10, 17}, + {27, 20, 13, 6, 1, 8, 15, 22, 29, 36, 43, 50, 57, 64, 71, 6, 12, 19}, + {26, 19, 12, 5, 2, 9, 16, 23, 30, 37, 44, 51, 58, 65, 5, 7, 14, 21}, + {26, 19, 12, 4, 3, 10, 17, 24, 31, 38, 45, 52, 60, 67, 2, 9, 16, 23}, + {25, 18, 11, 4, 3, 11, 18, 25, 32, 39, 46, 54, 61, 68, 3, 10, 17, 25}, + {25, 17, 10, 3, 4, 11, 19, 26, 33, 40, 48, 55, 62, 69, 4, 12, 19, 26}}; + +#define UCOMPARE(TYPE, VAR1, VAR2, TOL) \ + { \ + size_t i, j; \ + \ + for (i = 0; i < ROWS; i++) \ + for (j = 0; j < COLS; j++) { \ + if (!((((VAR1)[i][j] >= (TYPE)((VAR2)[i][j])) && \ + (((VAR1)[i][j] - (TOL)) < (TYPE)((VAR2)[i][j]))) || \ + (((VAR1)[i][j] <= (TYPE)((VAR2)[i][j])) && \ + (((VAR1)[i][j] + (TOL)) > (TYPE)((VAR2)[i][j]))))) { \ + H5_FAILED(); \ + HDfprintf(stderr, " ERROR: Conversion failed to match computed data\n"); \ + goto error; \ + } \ + } \ + PASSED(); \ + } + +#define COMPARE(TYPE, VAR1, VAR2, TOL) \ + { \ + size_t i, j; \ + \ + for (i = 0; i < ROWS; i++) \ + for (j = 0; j < COLS; j++) { \ + if (!(((VAR1)[i][j] <= ((TYPE)(VAR2)[i][j] + (TOL))) && \ + ((VAR1)[i][j] >= ((TYPE)(VAR2)[i][j] - (TOL))))) { \ + H5_FAILED(); \ + HDfprintf(stderr, " ERROR: Conversion failed to match computed data\n"); \ + goto error; \ + } \ + } \ + PASSED(); \ + } + +#define COMPARE_INT(VAR1, VAR2) \ + { \ + size_t i, j; \ + \ + for (i = 0; i < ROWS; i++) \ + for (j = 0; j < COLS; j++) { \ + if ((VAR1)[i][j] != (VAR2)[i][j]) { \ + H5_FAILED(); \ + HDfprintf(stderr, " ERROR: data failed to match computed data\n"); \ + goto error; \ + } \ + } \ + } + +#define TEST_TYPE_CONTIG(XFORM, TYPE, HDF_TYPE, TEST_STR, COMPARE_DATA, SIGNED) \ + { \ + TYPE array[ROWS][COLS]; \ + const char *f_to_c = "(5/9.0)*(x-32)"; \ + /* utrans is a transform for char types: numbers are restricted from -128 to 127, fits into char */ \ + const char *utrans = "(x/4+25)*3"; \ + \ + hid_t dataspace, dxpl_id_f_to_c, dxpl_id_utrans, dset, dset_nn, dt_nn; \ + H5T_order_t order; \ + hsize_t dim[2] = {ROWS, COLS}; \ + \ + if ((dataspace = H5Screate_simple(2, dim, NULL)) < 0) \ + TEST_ERROR; \ + if ((dset = H5Dcreate2(file_id, "/transformtest_" TEST_STR, HDF_TYPE, dataspace, H5P_DEFAULT, \ + H5P_DEFAULT, H5P_DEFAULT)) < 0) \ + TEST_ERROR; \ + \ + if ((dt_nn = H5Tcopy(HDF_TYPE)) < 0) \ + TEST_ERROR \ + if ((order = H5Tget_order(dt_nn)) == H5T_ORDER_ERROR) \ + TEST_ERROR \ + if (H5Tset_order(dt_nn, order == H5T_ORDER_LE ? H5T_ORDER_BE : H5T_ORDER_LE) < 0) \ + TEST_ERROR \ + if ((dset_nn = H5Dcreate2(file_id, "/nonnative_transformtest_" TEST_STR, dt_nn, dataspace, \ + H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) \ + TEST_ERROR \ + if (H5Tclose(dt_nn) < 0) \ + TEST_ERROR \ + \ + if (SIGNED) { \ + if ((dxpl_id_f_to_c = H5Pcreate(H5P_DATASET_XFER)) < 0) \ + TEST_ERROR; \ + if (H5Pset_data_transform(dxpl_id_f_to_c, f_to_c) < 0) \ + TEST_ERROR; \ + if (H5Dwrite(dset, H5T_NATIVE_FLOAT, H5S_ALL, H5S_ALL, dxpl_id_f_to_c, windchillFfloat) < 0) \ + TEST_ERROR; \ + if (H5Dwrite(dset_nn, H5T_NATIVE_FLOAT, H5S_ALL, H5S_ALL, dxpl_id_f_to_c, windchillFfloat) < 0) \ + TEST_ERROR; \ + if (H5Pclose(dxpl_id_f_to_c) < 0) \ + TEST_ERROR; \ + } \ + else { \ + if ((dxpl_id_utrans = H5Pcreate(H5P_DATASET_XFER)) < 0) \ + TEST_ERROR; \ + if (H5Pset_data_transform(dxpl_id_utrans, utrans) < 0) \ + TEST_ERROR; \ + if (H5Dwrite(dset, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, dxpl_id_utrans, transformData) < 0) \ + TEST_ERROR; \ + if (H5Dwrite(dset_nn, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, dxpl_id_utrans, transformData) < 0) \ + TEST_ERROR; \ + if (H5Pclose(dxpl_id_utrans) < 0) \ + TEST_ERROR; \ + } \ + \ + TESTING("contiguous, no data type conversion (" TEST_STR "->" TEST_STR ")") \ + \ + if (H5Dread(dset, HDF_TYPE, H5S_ALL, H5S_ALL, XFORM, array) < 0) \ + TEST_ERROR; \ + if (SIGNED) \ + COMPARE(TYPE, array, COMPARE_DATA, 2) \ + else \ + UCOMPARE(TYPE, array, COMPARE_DATA, 4) \ + \ + TESTING("contiguous, byte order conversion (" TEST_STR "->" TEST_STR ")") \ + \ + if (H5Dread(dset_nn, HDF_TYPE, H5S_ALL, H5S_ALL, XFORM, array) < 0) \ + TEST_ERROR; \ + if (SIGNED) \ + COMPARE(TYPE, array, COMPARE_DATA, 2) \ + else \ + UCOMPARE(TYPE, array, COMPARE_DATA, 4) \ + \ + if (SIGNED) { \ + TESTING("contiguous, with type conversion (float->" TEST_STR ")") \ + \ + if (H5Dread(dset_id_float, HDF_TYPE, H5S_ALL, H5S_ALL, XFORM, array) < 0) \ + TEST_ERROR; \ + COMPARE(TYPE, array, COMPARE_DATA, 2) \ + } \ + \ + if (H5Dclose(dset_nn) < 0) \ + TEST_ERROR; \ + if (H5Dclose(dset) < 0) \ + TEST_ERROR; \ + if (H5Sclose(dataspace) < 0) \ + TEST_ERROR; \ + } + +#define TEST_TYPE_CHUNK(XFORM, TYPE, HDF_TYPE, TEST_STR, COMPARE_DATA, SIGNED) \ + { \ + TYPE array[ROWS][COLS]; \ + const char *f_to_c = "(5/9.0)*(x-32)"; \ + /* utrans is a transform for char types: numbers are restricted from -128 to 127, fits into char */ \ + const char *utrans = "(x/4+25)*3"; \ + \ + hid_t dataspace, dxpl_id_f_to_c, dxpl_id_utrans, cparms, memspace, dset_chunk, filespace; \ + hsize_t dim[2] = {ROWS, COLS}; \ + hsize_t offset[2] = {0, 0}; \ + \ + if ((dataspace = H5Screate_simple(2, dim, NULL)) < 0) \ + TEST_ERROR; \ + \ + if ((cparms = H5Pcreate(H5P_DATASET_CREATE)) < 0) \ + TEST_ERROR; \ + if (H5Pset_chunk(cparms, 2, dim) < 0) \ + TEST_ERROR; \ + \ + if ((dset_chunk = H5Dcreate2(file_id, "/transformtest_chunk_" TEST_STR, HDF_TYPE, dataspace, \ + H5P_DEFAULT, cparms, H5P_DEFAULT)) < 0) \ + TEST_ERROR; \ + if ((filespace = H5Dget_space(dset_chunk)) < 0) \ + TEST_ERROR \ + if ((memspace = H5Screate_simple(2, dim, NULL)) < 0) \ + TEST_ERROR \ + if (H5Sselect_hyperslab(filespace, H5S_SELECT_SET, offset, NULL, dim, NULL) < 0) \ + TEST_ERROR; \ + \ + if (SIGNED) { \ + if ((dxpl_id_f_to_c = H5Pcreate(H5P_DATASET_XFER)) < 0) \ + TEST_ERROR; \ + if (H5Pset_data_transform(dxpl_id_f_to_c, f_to_c) < 0) \ + TEST_ERROR; \ + if (H5Dwrite(dset_chunk, H5T_NATIVE_FLOAT, dataspace, filespace, dxpl_id_f_to_c, \ + windchillFfloat) < 0) \ + TEST_ERROR; \ + if (H5Pclose(dxpl_id_f_to_c) < 0) \ + TEST_ERROR; \ + } \ + else { \ + if ((dxpl_id_utrans = H5Pcreate(H5P_DATASET_XFER)) < 0) \ + TEST_ERROR; \ + if (H5Pset_data_transform(dxpl_id_utrans, utrans) < 0) \ + TEST_ERROR; \ + if (H5Dwrite(dset_chunk, H5T_NATIVE_INT, dataspace, filespace, dxpl_id_utrans, transformData) < \ + 0) \ + TEST_ERROR; \ + if (H5Pclose(dxpl_id_utrans) < 0) \ + TEST_ERROR; \ + } \ + \ + TESTING("chunked, no data type conversion (" TEST_STR "->" TEST_STR ")") \ + \ + if (H5Dread(dset_chunk, HDF_TYPE, memspace, filespace, XFORM, array) < 0) \ + TEST_ERROR; \ + if (SIGNED) \ + COMPARE(TYPE, array, COMPARE_DATA, 2) \ + else \ + UCOMPARE(TYPE, array, COMPARE_DATA, 4) \ + \ + if (SIGNED) { \ + TESTING("chunked, with type conversion (float->" TEST_STR ")") \ + \ + if (H5Dread(dset_id_float_chunk, HDF_TYPE, memspace, filespace, XFORM, array) < 0) \ + TEST_ERROR; \ + COMPARE(TYPE, array, COMPARE_DATA, 2) \ + } \ + \ + if (H5Pclose(cparms) < 0) \ + TEST_ERROR; \ + if (H5Dclose(dset_chunk) < 0) \ + TEST_ERROR; \ + if (H5Sclose(dataspace) < 0) \ + TEST_ERROR; \ + if (H5Sclose(memspace) < 0) \ + TEST_ERROR; \ + } + +#define INVALID_SET_TEST(TRANSFORM) \ + { \ + if (H5Pset_data_transform(dxpl_id, TRANSFORM) < 0) { \ + PASSED(); \ + } \ + else { \ + H5_FAILED(); \ + HDfprintf(stderr, " ERROR: Data transform allowed invalid TRANSFORM transform to be set\n"); \ + goto error; \ + } \ + } + +int +main(void) { - hid_t dxpl_id_c_to_f = -1; - hid_t dxpl_id_c_to_f_copy = 1; - hid_t dxpl_id_simple = -1; - hid_t dxpl_id_polynomial = -1; + hid_t dxpl_id_c_to_f = -1; + hid_t dxpl_id_c_to_f_copy = 1; + hid_t dxpl_id_simple = -1; + hid_t dxpl_id_polynomial = -1; hid_t dxpl_id_polynomial_copy = -1; - hid_t dxpl_id_utrans_inv = -1; - hid_t file_id = -1; + hid_t dxpl_id_utrans_inv = -1; + hid_t file_id = -1; - const char* c_to_f = "(9/5.0)*x + 32"; - const char* simple = "(4/2) * ( (2 + 4)/(5 - 2.5))"; /* this equals 4.8 */ - const char* polynomial = "(2+x)* ((x-8)/2)"; + const char *c_to_f = "(9/5.0)*x + 32"; + const char *simple = "(4/2) * ( (2 + 4)/(5 - 2.5))"; /* this equals 4.8 */ + const char *polynomial = "(2+x)* ((x-8)/2)"; /* inverses the utrans transform in init_test to get back original array */ - const char* utrans_inv = "(x/3)*4 - 100"; - - if((file_id = H5Fcreate("dtransform.h5", H5F_ACC_TRUNC, H5P_DEFAULT, H5P_DEFAULT)) < 0) TEST_ERROR; - - if((dxpl_id_c_to_f = H5Pcreate(H5P_DATASET_XFER)) < 0) TEST_ERROR; - if((dxpl_id_simple = H5Pcreate(H5P_DATASET_XFER)) < 0) TEST_ERROR; - if((dxpl_id_utrans_inv = H5Pcreate(H5P_DATASET_XFER)) < 0) TEST_ERROR; - if((dxpl_id_polynomial = H5Pcreate(H5P_DATASET_XFER)) < 0) TEST_ERROR; - if(H5Pset_data_transform(dxpl_id_c_to_f, c_to_f) < 0) TEST_ERROR; - if(H5Pset_data_transform(dxpl_id_polynomial, polynomial) < 0) TEST_ERROR; - if(H5Pset_data_transform(dxpl_id_simple, simple) < 0) TEST_ERROR; - if(H5Pset_data_transform(dxpl_id_utrans_inv, utrans_inv) < 0) TEST_ERROR; - if((dxpl_id_polynomial_copy = H5Pcopy(dxpl_id_polynomial)) < 0) TEST_ERROR; - if((dxpl_id_c_to_f_copy = H5Pcopy(dxpl_id_c_to_f)) < 0) TEST_ERROR; + const char *utrans_inv = "(x/3 - 25)*4"; + + if ((file_id = H5Fcreate("dtransform.h5", H5F_ACC_TRUNC, H5P_DEFAULT, H5P_DEFAULT)) < 0) + TEST_ERROR; + + if ((dxpl_id_c_to_f = H5Pcreate(H5P_DATASET_XFER)) < 0) + TEST_ERROR; + if ((dxpl_id_simple = H5Pcreate(H5P_DATASET_XFER)) < 0) + TEST_ERROR; + if ((dxpl_id_utrans_inv = H5Pcreate(H5P_DATASET_XFER)) < 0) + TEST_ERROR; + if ((dxpl_id_polynomial = H5Pcreate(H5P_DATASET_XFER)) < 0) + TEST_ERROR; + if (H5Pset_data_transform(dxpl_id_c_to_f, c_to_f) < 0) + TEST_ERROR; + if (H5Pset_data_transform(dxpl_id_polynomial, polynomial) < 0) + TEST_ERROR; + if (H5Pset_data_transform(dxpl_id_simple, simple) < 0) + TEST_ERROR; + if (H5Pset_data_transform(dxpl_id_utrans_inv, utrans_inv) < 0) + TEST_ERROR; + if ((dxpl_id_polynomial_copy = H5Pcopy(dxpl_id_polynomial)) < 0) + TEST_ERROR; + if ((dxpl_id_c_to_f_copy = H5Pcopy(dxpl_id_c_to_f)) < 0) + TEST_ERROR; /* Run all the tests */ - if(init_test(file_id) < 0) TEST_ERROR; - if(test_set() < 0) TEST_ERROR; + if (init_test(file_id) < 0) + TEST_ERROR; + if (test_set() < 0) + TEST_ERROR; TEST_TYPE_CONTIG(dxpl_id_utrans_inv, char, H5T_NATIVE_CHAR, "char", transformData, 0); TEST_TYPE_CONTIG(dxpl_id_utrans_inv, unsigned char, H5T_NATIVE_UCHAR, "uchar", transformData, 0); TEST_TYPE_CONTIG(dxpl_id_c_to_f, signed char, H5T_NATIVE_SCHAR, "schar", windchillFfloat, 1); @@ -299,9 +357,7 @@ int main(void) TEST_TYPE_CONTIG(dxpl_id_utrans_inv, unsigned long long, H5T_NATIVE_ULLONG, "ullong", transformData, 0); TEST_TYPE_CONTIG(dxpl_id_c_to_f, float, H5T_NATIVE_FLOAT, "float", windchillFfloat, 1); TEST_TYPE_CONTIG(dxpl_id_c_to_f, double, H5T_NATIVE_DOUBLE, "double", windchillFfloat, 1); -#if H5_SIZEOF_LONG_DOUBLE!=0 TEST_TYPE_CONTIG(dxpl_id_c_to_f, long double, H5T_NATIVE_LDOUBLE, "ldouble", windchillFfloat, 1); -#endif TEST_TYPE_CHUNK(dxpl_id_utrans_inv, char, H5T_NATIVE_CHAR, "char", transformData, 0); TEST_TYPE_CHUNK(dxpl_id_utrans_inv, unsigned char, H5T_NATIVE_UCHAR, "uchar", transformData, 0); @@ -316,34 +372,48 @@ int main(void) TEST_TYPE_CHUNK(dxpl_id_utrans_inv, unsigned long long, H5T_NATIVE_ULLONG, "ullong", transformData, 0); TEST_TYPE_CHUNK(dxpl_id_c_to_f, float, H5T_NATIVE_FLOAT, "float", windchillFfloat, 1); TEST_TYPE_CHUNK(dxpl_id_c_to_f, double, H5T_NATIVE_DOUBLE, "double", windchillFfloat, 1); -#if H5_SIZEOF_LONG_DOUBLE!=0 TEST_TYPE_CHUNK(dxpl_id_c_to_f, long double, H5T_NATIVE_LDOUBLE, "ldouble", windchillFfloat, 1); -#endif - if(test_copy(dxpl_id_c_to_f_copy, dxpl_id_polynomial_copy) < 0) TEST_ERROR; - if(test_trivial(dxpl_id_simple) < 0) TEST_ERROR; - if(test_poly(dxpl_id_polynomial) < 0) TEST_ERROR; - if(test_getset(dxpl_id_c_to_f) < 0) TEST_ERROR; - if(test_specials(file_id) < 0) TEST_ERROR; + if (test_copy(dxpl_id_c_to_f_copy, dxpl_id_polynomial_copy) < 0) + TEST_ERROR; + if (test_trivial(dxpl_id_simple) < 0) + TEST_ERROR; + if (test_poly(dxpl_id_polynomial) < 0) + TEST_ERROR; + if (test_getset(dxpl_id_c_to_f) < 0) + TEST_ERROR; + if (test_specials(file_id) < 0) + TEST_ERROR; /* Close the objects we opened/created */ - if(H5Dclose(dset_id_int) < 0) TEST_ERROR; - if(H5Dclose(dset_id_int_chunk) < 0) TEST_ERROR; - if(H5Dclose(dset_id_float) < 0) TEST_ERROR; - if(H5Dclose(dset_id_float_chunk) < 0) TEST_ERROR; - if(H5Fclose(file_id) < 0) TEST_ERROR; - if(H5Pclose(dxpl_id_c_to_f) < 0) TEST_ERROR; - if(H5Pclose(dxpl_id_c_to_f_copy) < 0) TEST_ERROR; - if(H5Pclose(dxpl_id_polynomial) < 0) TEST_ERROR; - if(H5Pclose(dxpl_id_polynomial_copy) < 0) TEST_ERROR; - if(H5Pclose(dxpl_id_simple) < 0) TEST_ERROR; - if(H5Pclose(dxpl_id_utrans_inv) < 0) TEST_ERROR; - + if (H5Dclose(dset_id_int) < 0) + TEST_ERROR; + if (H5Dclose(dset_id_int_chunk) < 0) + TEST_ERROR; + if (H5Dclose(dset_id_float) < 0) + TEST_ERROR; + if (H5Dclose(dset_id_float_chunk) < 0) + TEST_ERROR; + if (H5Fclose(file_id) < 0) + TEST_ERROR; + if (H5Pclose(dxpl_id_c_to_f) < 0) + TEST_ERROR; + if (H5Pclose(dxpl_id_c_to_f_copy) < 0) + TEST_ERROR; + if (H5Pclose(dxpl_id_polynomial) < 0) + TEST_ERROR; + if (H5Pclose(dxpl_id_polynomial_copy) < 0) + TEST_ERROR; + if (H5Pclose(dxpl_id_simple) < 0) + TEST_ERROR; + if (H5Pclose(dxpl_id_utrans_inv) < 0) + TEST_ERROR; - return 0; + return 0; error: - H5E_BEGIN_TRY { + H5E_BEGIN_TRY + { H5Dclose(dset_id_int); H5Dclose(dset_id_int_chunk); H5Dclose(dset_id_float); @@ -355,86 +425,85 @@ error: H5Pclose(dxpl_id_polynomial_copy); H5Pclose(dxpl_id_simple); H5Pclose(dxpl_id_utrans_inv); - } H5E_END_TRY - return -1; + } + H5E_END_TRY + return -1; } static int init_test(hid_t file_id) { - const char* f_to_c = "(5/9.0)*(x-32)"; - /* utrans is a transform for unsigned types: no negative numbers involved and results are < 255 to fit into uchar */ - const char* utrans = "((x+100)/4)*3"; + const char *f_to_c = "(5/9.0)*(x-32)"; + /* utrans is a transform for char types: numbers are restricted from -128 to 127, fits into char */ + const char *utrans = "(x/4+25)*3"; - hid_t dataspace = -1; + hid_t dataspace = -1; hid_t dxpl_id_f_to_c = -1; hid_t dxpl_id_utrans = -1; - hid_t cparms = -1; - hid_t filespace = -1; - hsize_t dim[2] = { ROWS, COLS }; - hsize_t offset[2] = { 0, 0 }; + hid_t cparms = -1; + hid_t filespace = -1; + hsize_t dim[2] = {ROWS, COLS}; + hsize_t offset[2] = {0, 0}; - if((dxpl_id_f_to_c = H5Pcreate(H5P_DATASET_XFER)) < 0) + if ((dxpl_id_f_to_c = H5Pcreate(H5P_DATASET_XFER)) < 0) TEST_ERROR - if((dxpl_id_utrans = H5Pcreate(H5P_DATASET_XFER)) < 0) + if ((dxpl_id_utrans = H5Pcreate(H5P_DATASET_XFER)) < 0) TEST_ERROR - if(H5Pset_data_transform(dxpl_id_f_to_c, f_to_c) < 0) + if (H5Pset_data_transform(dxpl_id_f_to_c, f_to_c) < 0) TEST_ERROR - if(H5Pset_data_transform(dxpl_id_utrans, utrans) < 0) + if (H5Pset_data_transform(dxpl_id_utrans, utrans) < 0) TEST_ERROR cparms = H5Pcreate(H5P_DATASET_CREATE); - if(H5Pset_chunk(cparms, 2, dim) < 0) + if (H5Pset_chunk(cparms, 2, dim) < 0) TEST_ERROR - if((dataspace = H5Screate_simple(2, dim, NULL)) < 0) + if ((dataspace = H5Screate_simple(2, dim, NULL)) < 0) TEST_ERROR - TESTING("Intializing test...") + TESTING("Initializing test...") - if((dset_id_int = H5Dcreate2(file_id, "/default_int", H5T_NATIVE_INT, - dataspace, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) + if ((dset_id_int = H5Dcreate2(file_id, "/default_int", H5T_NATIVE_INT, dataspace, H5P_DEFAULT, + H5P_DEFAULT, H5P_DEFAULT)) < 0) TEST_ERROR - if(H5Dwrite(dset_id_int, H5T_NATIVE_FLOAT, H5S_ALL, H5S_ALL, - dxpl_id_f_to_c, windchillFfloat) < 0) + if (H5Dwrite(dset_id_int, H5T_NATIVE_FLOAT, H5S_ALL, H5S_ALL, dxpl_id_f_to_c, windchillFfloat) < 0) TEST_ERROR - if((dset_id_float = H5Dcreate2(file_id, "/default_float", - H5T_NATIVE_FLOAT, dataspace, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) + if ((dset_id_float = H5Dcreate2(file_id, "/default_float", H5T_NATIVE_FLOAT, dataspace, H5P_DEFAULT, + H5P_DEFAULT, H5P_DEFAULT)) < 0) TEST_ERROR - if(H5Dwrite(dset_id_float, H5T_NATIVE_FLOAT, H5S_ALL, H5S_ALL, - dxpl_id_f_to_c, windchillFfloat) < 0) + if (H5Dwrite(dset_id_float, H5T_NATIVE_FLOAT, H5S_ALL, H5S_ALL, dxpl_id_f_to_c, windchillFfloat) < 0) TEST_ERROR - if((dset_id_int_chunk = H5Dcreate2(file_id, "/default_chunk_int", - H5T_NATIVE_INT, dataspace, H5P_DEFAULT, cparms, H5P_DEFAULT)) < 0) + if ((dset_id_int_chunk = H5Dcreate2(file_id, "/default_chunk_int", H5T_NATIVE_INT, dataspace, H5P_DEFAULT, + cparms, H5P_DEFAULT)) < 0) TEST_ERROR - if((filespace = H5Dget_space(dset_id_int_chunk)) < 0) + if ((filespace = H5Dget_space(dset_id_int_chunk)) < 0) TEST_ERROR - if(H5Sselect_hyperslab(filespace, H5S_SELECT_SET, offset, NULL, dim, NULL) < 0) + if (H5Sselect_hyperslab(filespace, H5S_SELECT_SET, offset, NULL, dim, NULL) < 0) TEST_ERROR - if(H5Dwrite(dset_id_int_chunk, H5T_NATIVE_FLOAT, dataspace, filespace, - dxpl_id_f_to_c, windchillFfloat) < 0) + if (H5Dwrite(dset_id_int_chunk, H5T_NATIVE_FLOAT, dataspace, filespace, dxpl_id_f_to_c, windchillFfloat) < + 0) TEST_ERROR - if((dset_id_float_chunk = H5Dcreate2(file_id, "/default_chunk_float", - H5T_NATIVE_FLOAT, dataspace, H5P_DEFAULT, cparms, H5P_DEFAULT)) < 0) + if ((dset_id_float_chunk = H5Dcreate2(file_id, "/default_chunk_float", H5T_NATIVE_FLOAT, dataspace, + H5P_DEFAULT, cparms, H5P_DEFAULT)) < 0) TEST_ERROR - if(H5Dwrite(dset_id_float_chunk, H5T_NATIVE_FLOAT, dataspace, filespace, - dxpl_id_f_to_c, windchillFfloat) < 0) + if (H5Dwrite(dset_id_float_chunk, H5T_NATIVE_FLOAT, dataspace, filespace, dxpl_id_f_to_c, + windchillFfloat) < 0) TEST_ERROR - if(H5Pclose(cparms) < 0) + if (H5Pclose(cparms) < 0) TEST_ERROR - if(H5Pclose(dxpl_id_f_to_c) < 0) + if (H5Pclose(dxpl_id_f_to_c) < 0) TEST_ERROR - if(H5Pclose(dxpl_id_utrans) < 0) + if (H5Pclose(dxpl_id_utrans) < 0) TEST_ERROR - if(H5Sclose(dataspace) < 0) + if (H5Sclose(dataspace) < 0) TEST_ERROR - if(H5Sclose(filespace) < 0) + if (H5Sclose(filespace) < 0) TEST_ERROR PASSED(); @@ -442,13 +511,15 @@ init_test(hid_t file_id) return 0; error: - H5E_BEGIN_TRY { + H5E_BEGIN_TRY + { H5Pclose(cparms); H5Pclose(dxpl_id_f_to_c); H5Pclose(dxpl_id_utrans); H5Sclose(dataspace); H5Sclose(filespace); - } H5E_END_TRY + } + H5E_END_TRY return -1; } @@ -457,33 +528,31 @@ static int test_poly(const hid_t dxpl_id_polynomial) { float polyflres[ROWS][COLS]; - int polyintread[ROWS][COLS]; + int polyintread[ROWS][COLS]; float polyflread[ROWS][COLS]; - int windchillC; - int row, col; + int windchillC; + int row, col; - for(row = 0; row < ROWS; row++) - for(col = 0; col < COLS; col++) { - windchillC = (int) ((5.0f / 9.0f) * (windchillFfloat[row][col] - 32)); - polyflres[row][col] = ((2.0f + (float)windchillC) * (((float)windchillC - 8.0f) / 2.0f)); + for (row = 0; row < ROWS; row++) + for (col = 0; col < COLS; col++) { + windchillC = (int)((5.0F / 9.0F) * (windchillFfloat[row][col] - 32)); + polyflres[row][col] = ((2.0F + (float)windchillC) * (((float)windchillC - 8.0F) / 2.0F)); } TESTING("data transform, polynomial transform (int->float)") - if(H5Dread(dset_id_int, H5T_NATIVE_FLOAT, H5S_ALL, H5S_ALL, - dxpl_id_polynomial, polyflread) < 0) + if (H5Dread(dset_id_int, H5T_NATIVE_FLOAT, H5S_ALL, H5S_ALL, dxpl_id_polynomial, polyflread) < 0) TEST_ERROR - COMPARE(float, polyflread, polyflres, 2.0f) + COMPARE(float, polyflread, polyflres, 2.0F) - for(row = 0; row < ROWS; row++) - for(col = 0; col < COLS; col++) { - windchillC = (int) ((5.0f / 9.0f) * (windchillFfloat[row][col] - 32)); - polyflres[row][col] = (float) ((2 + windchillC) * ((windchillC - 8) / 2)); + for (row = 0; row < ROWS; row++) + for (col = 0; col < COLS; col++) { + windchillC = (int)((5.0F / 9.0F) * (windchillFfloat[row][col] - 32)); + polyflres[row][col] = (float)((2 + windchillC) * ((windchillC - 8) / 2)); } TESTING("data transform, polynomial transform (float->int)") - if(H5Dread(dset_id_float, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, - dxpl_id_polynomial, polyintread) < 0) + if (H5Dread(dset_id_float, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, dxpl_id_polynomial, polyintread) < 0) TEST_ERROR COMPARE(int, polyintread, polyflres, 4) @@ -491,188 +560,232 @@ test_poly(const hid_t dxpl_id_polynomial) return 0; error: - return -1; + return -1; } static int test_specials(hid_t file) { - hid_t dxpl_id, dset_id, dataspace; - hsize_t dim[2] = { ROWS, COLS }; - int read_buf[ROWS][COLS]; - int data_res[ROWS][COLS]; - int row, col; - const char* special1 = "x*-100"; - const char* special2 = "100-x"; - const char* special3 = "1000/x"; - const char* special4 = "-x"; - const char* special5 = "+x"; + hid_t dxpl_id, dset_id, dataspace; + hsize_t dim[2] = {ROWS, COLS}; + int read_buf[ROWS][COLS]; + int data_res[ROWS][COLS]; + int row, col; + const char *special1 = "x*-100"; + const char *special2 = "100-x"; + const char *special3 = "1000/x"; + const char *special4 = "-x"; + const char *special5 = "+x"; + const char *special6 = "2e+1*x"; + const char *special7 = "x"; TESTING("data transform of some special cases") - if((dataspace = H5Screate_simple(2, dim, NULL)) < 0) + if ((dataspace = H5Screate_simple(2, dim, NULL)) < 0) TEST_ERROR - if((dxpl_id = H5Pcreate(H5P_DATASET_XFER)) < 0) TEST_ERROR; + if ((dxpl_id = H5Pcreate(H5P_DATASET_XFER)) < 0) + TEST_ERROR; /*----------------------------- * Operation 1: x*-100 *----------------------------*/ - if(H5Pset_data_transform(dxpl_id, special1) < 0) TEST_ERROR; + if (H5Pset_data_transform(dxpl_id, special1) < 0) + TEST_ERROR; - for(row = 0; row < ROWS; row++) - for(col = 0; col < COLS; col++) + for (row = 0; row < ROWS; row++) + for (col = 0; col < COLS; col++) data_res[row][col] = transformData[row][col] * -100; - if((dset_id = H5Dcreate2(file, "/special1", H5T_NATIVE_INT, - dataspace, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) + if ((dset_id = H5Dcreate2(file, "/special1", H5T_NATIVE_INT, dataspace, H5P_DEFAULT, H5P_DEFAULT, + H5P_DEFAULT)) < 0) TEST_ERROR - if(H5Dwrite(dset_id, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, - dxpl_id, transformData) < 0) + if (H5Dwrite(dset_id, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, dxpl_id, transformData) < 0) TEST_ERROR - if(H5Dread(dset_id, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, - H5P_DEFAULT, read_buf) < 0) + if (H5Dread(dset_id, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, read_buf) < 0) TEST_ERROR COMPARE_INT(read_buf, data_res) - if(H5Dclose(dset_id) < 0) + if (H5Dclose(dset_id) < 0) TEST_ERROR /*----------------------------- * Operation 2: 100-x *----------------------------*/ - if(H5Pset_data_transform(dxpl_id, special2) < 0) TEST_ERROR; + if (H5Pset_data_transform(dxpl_id, special2) < 0) + TEST_ERROR; - for(row = 0; row < ROWS; row++) - for(col = 0; col < COLS; col++) + for (row = 0; row < ROWS; row++) + for (col = 0; col < COLS; col++) data_res[row][col] = 100 - transformData[row][col]; - if((dset_id = H5Dcreate2(file, "/special2", H5T_NATIVE_INT, - dataspace, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) + if ((dset_id = H5Dcreate2(file, "/special2", H5T_NATIVE_INT, dataspace, H5P_DEFAULT, H5P_DEFAULT, + H5P_DEFAULT)) < 0) TEST_ERROR - if(H5Dwrite(dset_id, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, - dxpl_id, transformData) < 0) + if (H5Dwrite(dset_id, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, dxpl_id, transformData) < 0) TEST_ERROR - if(H5Dread(dset_id, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, - H5P_DEFAULT, read_buf) < 0) + if (H5Dread(dset_id, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, read_buf) < 0) TEST_ERROR COMPARE_INT(read_buf, data_res) - if(H5Dclose(dset_id) < 0) + if (H5Dclose(dset_id) < 0) TEST_ERROR /*----------------------------- * Operation 3: 1000/x *----------------------------*/ - if(H5Pset_data_transform(dxpl_id, special3) < 0) TEST_ERROR; + if (H5Pset_data_transform(dxpl_id, special3) < 0) + TEST_ERROR; - for(row = 0; row < ROWS; row++) - for(col = 0; col < COLS; col++) + for (row = 0; row < ROWS; row++) + for (col = 0; col < COLS; col++) data_res[row][col] = 1000 / transformData[row][col]; - if((dset_id = H5Dcreate2(file, "/special3", H5T_NATIVE_INT, - dataspace, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) + if ((dset_id = H5Dcreate2(file, "/special3", H5T_NATIVE_INT, dataspace, H5P_DEFAULT, H5P_DEFAULT, + H5P_DEFAULT)) < 0) TEST_ERROR - if(H5Dwrite(dset_id, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, - dxpl_id, transformData) < 0) + if (H5Dwrite(dset_id, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, dxpl_id, transformData) < 0) TEST_ERROR - if(H5Dread(dset_id, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, - H5P_DEFAULT, read_buf) < 0) + if (H5Dread(dset_id, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, read_buf) < 0) TEST_ERROR COMPARE_INT(read_buf, data_res) - if(H5Dclose(dset_id) < 0) + if (H5Dclose(dset_id) < 0) TEST_ERROR /*----------------------------- * Operation 4: -x *----------------------------*/ - if(H5Pset_data_transform(dxpl_id, special4) < 0) TEST_ERROR; + if (H5Pset_data_transform(dxpl_id, special4) < 0) + TEST_ERROR; - for(row = 0; row < ROWS; row++) - for(col = 0; col < COLS; col++) + for (row = 0; row < ROWS; row++) + for (col = 0; col < COLS; col++) data_res[row][col] = -1 * transformData[row][col]; - if((dset_id = H5Dcreate2(file, "/special4", H5T_NATIVE_INT, - dataspace, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) + if ((dset_id = H5Dcreate2(file, "/special4", H5T_NATIVE_INT, dataspace, H5P_DEFAULT, H5P_DEFAULT, + H5P_DEFAULT)) < 0) TEST_ERROR - if(H5Dwrite(dset_id, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, - dxpl_id, transformData) < 0) + if (H5Dwrite(dset_id, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, dxpl_id, transformData) < 0) TEST_ERROR - if(H5Dread(dset_id, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, - H5P_DEFAULT, read_buf) < 0) + if (H5Dread(dset_id, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, read_buf) < 0) TEST_ERROR COMPARE_INT(read_buf, data_res) - if(H5Dclose(dset_id) < 0) + if (H5Dclose(dset_id) < 0) TEST_ERROR /*----------------------------- * Operation 5: +x *----------------------------*/ - if(H5Pset_data_transform(dxpl_id, special5) < 0) TEST_ERROR; + if (H5Pset_data_transform(dxpl_id, special5) < 0) + TEST_ERROR; - for(row = 0; row < ROWS; row++) - for(col = 0; col < COLS; col++) + for (row = 0; row < ROWS; row++) + for (col = 0; col < COLS; col++) data_res[row][col] = transformData[row][col]; - if((dset_id = H5Dcreate2(file, "/special5", H5T_NATIVE_INT, - dataspace, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) + if ((dset_id = H5Dcreate2(file, "/special5", H5T_NATIVE_INT, dataspace, H5P_DEFAULT, H5P_DEFAULT, + H5P_DEFAULT)) < 0) + TEST_ERROR + if (H5Dwrite(dset_id, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, dxpl_id, transformData) < 0) + TEST_ERROR + if (H5Dread(dset_id, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, read_buf) < 0) + TEST_ERROR + + COMPARE_INT(read_buf, data_res) + + if (H5Dclose(dset_id) < 0) TEST_ERROR - if(H5Dwrite(dset_id, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, - dxpl_id, transformData) < 0) + + /*----------------------------- + * Operation 6: 2e+1*x + *----------------------------*/ + if (H5Pset_data_transform(dxpl_id, special6) < 0) + TEST_ERROR; + + for (row = 0; row < ROWS; row++) + for (col = 0; col < COLS; col++) + data_res[row][col] = transformData[row][col] * 20; + + if ((dset_id = H5Dcreate2(file, "/special6", H5T_NATIVE_INT, dataspace, H5P_DEFAULT, H5P_DEFAULT, + H5P_DEFAULT)) < 0) + TEST_ERROR + if (H5Dwrite(dset_id, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, dxpl_id, transformData) < 0) TEST_ERROR - if(H5Dread(dset_id, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, - H5P_DEFAULT, read_buf) < 0) + if (H5Dread(dset_id, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, read_buf) < 0) TEST_ERROR COMPARE_INT(read_buf, data_res) - if(H5Dclose(dset_id) < 0) + if (H5Dclose(dset_id) < 0) TEST_ERROR + /*----------------------------- + * Operation 7: x + * This operation will be + * treated if no function has + * been specified. + *----------------------------*/ + if (H5Pset_data_transform(dxpl_id, special7) < 0) + TEST_ERROR; + + for (row = 0; row < ROWS; row++) + for (col = 0; col < COLS; col++) + data_res[row][col] = transformData[row][col]; - if(H5Pclose(dxpl_id) < 0) + if ((dset_id = H5Dcreate2(file, "/special7", H5T_NATIVE_INT, dataspace, H5P_DEFAULT, H5P_DEFAULT, + H5P_DEFAULT)) < 0) TEST_ERROR - if(H5Sclose(dataspace) < 0) + if (H5Dwrite(dset_id, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, dxpl_id, transformData) < 0) + TEST_ERROR + if (H5Dread(dset_id, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, read_buf) < 0) + TEST_ERROR + + COMPARE_INT(read_buf, data_res) + + if (H5Dclose(dset_id) < 0) + TEST_ERROR + + if (H5Pclose(dxpl_id) < 0) + TEST_ERROR + if (H5Sclose(dataspace) < 0) TEST_ERROR PASSED(); return 0; error: - return -1; + return -1; } static int test_copy(const hid_t dxpl_id_c_to_f_copy, const hid_t dxpl_id_polynomial_copy) { - int windchillC; + int windchillC; float polyflres[ROWS][COLS]; - int polyintread[ROWS][COLS]; - int windchillFintread[ROWS][COLS]; - int row, col; - - for(row = 0; row < ROWS; row++) - for(col = 0; col < COLS; col++) { - windchillC = (int) ((5.0f / 9.0f) * (windchillFfloat[row][col] - 32)); - polyflres[row][col] = (float) ((2 + windchillC) * ((windchillC - 8) / 2)); + int polyintread[ROWS][COLS]; + int windchillFintread[ROWS][COLS]; + int row, col; + + for (row = 0; row < ROWS; row++) + for (col = 0; col < COLS; col++) { + windchillC = (int)((5.0F / 9.0F) * (windchillFfloat[row][col] - 32)); + polyflres[row][col] = (float)((2 + windchillC) * ((windchillC - 8) / 2)); } TESTING("data transform, linear transform w/ copied property") - if(H5Dread(dset_id_float, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, - dxpl_id_c_to_f_copy, windchillFintread) < 0) + if (H5Dread(dset_id_float, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, dxpl_id_c_to_f_copy, windchillFintread) < 0) TEST_ERROR COMPARE(int, windchillFintread, windchillFfloat, 2) TESTING("data transform, polynomial transform w/ copied property") - if(H5Dread(dset_id_float, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, - dxpl_id_polynomial_copy, polyintread) < 0) + if (H5Dread(dset_id_float, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, dxpl_id_polynomial_copy, polyintread) < 0) TEST_ERROR COMPARE(int, polyintread, polyflres, 2) @@ -687,32 +800,30 @@ static int test_trivial(const hid_t dxpl_id_simple) { float windchillFfloatread[ROWS][COLS]; - int windchillFintread[ROWS][COLS]; - int row, col; + int windchillFintread[ROWS][COLS]; + int row, col; TESTING("data transform, trivial transform, without type conversion") - if(H5Dread(dset_id_float, H5T_NATIVE_FLOAT, H5S_ALL, H5S_ALL, - dxpl_id_simple, windchillFfloatread) < 0) + if (H5Dread(dset_id_float, H5T_NATIVE_FLOAT, H5S_ALL, H5S_ALL, dxpl_id_simple, windchillFfloatread) < 0) TEST_ERROR - for(row = 0; row < ROWS; row++) - for(col = 0; col < COLS; col++) { - if((windchillFfloatread[row][col] - 4.8f) > FLOAT_TOL) + for (row = 0; row < ROWS; row++) + for (col = 0; col < COLS; col++) { + if ((windchillFfloatread[row][col] - 4.8F) > FLOAT_TOL) FAIL_PUTS_ERROR(" ERROR: Conversion failed to match computed data\n"); } - PASSED() + PASSED(); TESTING("data transform, trivial transform, with type conversion") - if(H5Dread(dset_id_float, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, - dxpl_id_simple, windchillFintread) < 0) + if (H5Dread(dset_id_float, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, dxpl_id_simple, windchillFintread) < 0) TEST_ERROR - for(row = 0; row < ROWS; row++) - for(col = 0; col < COLS; col++) { - if(windchillFintread[row][col] != 4) + for (row = 0; row < ROWS; row++) + for (col = 0; col < COLS; col++) { + if (windchillFintread[row][col] != 4) FAIL_PUTS_ERROR(" ERROR: Conversion failed to match computed data\n") } - PASSED() + PASSED(); return 0; error: @@ -725,52 +836,51 @@ test_getset(const hid_t dxpl_id_c_to_f) int row; int col; float windchillFfloatread[ROWS][COLS]; - const char *simple = "(4/2) * ( (2 + 4)/(5 - 2.5))"; /* this equals 4.8 */ - const char *c_to_f = "(9/5.0)*x + 32"; - char *ptrgetTest = NULL; + const char *simple = "(4/2) * ( (2 + 4)/(5 - 2.5))"; /* this equals 4.8 */ + const char *c_to_f = "(9/5.0)*x + 32"; + char * ptrgetTest = NULL; TESTING("H5Pget_data_transform") - if(NULL == (ptrgetTest = (char *)HDmalloc(HDstrlen(simple) + 1))) + if (NULL == (ptrgetTest = (char *)HDmalloc(HDstrlen(simple) + 1))) TEST_ERROR - if(H5Pget_data_transform(dxpl_id_c_to_f, ptrgetTest, HDstrlen(c_to_f) + 1) < 0) + if (H5Pget_data_transform(dxpl_id_c_to_f, ptrgetTest, HDstrlen(c_to_f) + 1) < 0) TEST_ERROR - if(HDstrcmp(c_to_f, ptrgetTest) != 0) + if (HDstrcmp(c_to_f, ptrgetTest) != 0) FAIL_PUTS_ERROR(" ERROR: Data transform failed to match what was set\n") - PASSED() + PASSED(); HDfree(ptrgetTest); ptrgetTest = NULL; TESTING("data transform, read after resetting of transform property") - if(H5Pset_data_transform(dxpl_id_c_to_f, simple) < 0) + if (H5Pset_data_transform(dxpl_id_c_to_f, simple) < 0) TEST_ERROR - if(H5Dread(dset_id_float, H5T_NATIVE_FLOAT, H5S_ALL, H5S_ALL, - dxpl_id_c_to_f, windchillFfloatread) < 0) + if (H5Dread(dset_id_float, H5T_NATIVE_FLOAT, H5S_ALL, H5S_ALL, dxpl_id_c_to_f, windchillFfloatread) < 0) TEST_ERROR - for(row = 0; row < ROWS; row++) - for(col = 0; col < COLS; col++) { - if((windchillFfloatread[row][col] - 4.8f) > FLOAT_TOL) + for (row = 0; row < ROWS; row++) + for (col = 0; col < COLS; col++) { + if ((windchillFfloatread[row][col] - 4.8F) > FLOAT_TOL) FAIL_PUTS_ERROR(" ERROR: Conversion failed to match computed data\n") } - PASSED() + PASSED(); TESTING("H5Pget_data_transform, after resetting transform property") - if(NULL == (ptrgetTest = (char *)HDcalloc((size_t)1, HDstrlen(simple) + 1))) + if (NULL == (ptrgetTest = (char *)HDcalloc((size_t)1, HDstrlen(simple) + 1))) TEST_ERROR - if(H5Pget_data_transform(dxpl_id_c_to_f, ptrgetTest, HDstrlen(simple) + 1) < 0) + if (H5Pget_data_transform(dxpl_id_c_to_f, ptrgetTest, HDstrlen(simple) + 1) < 0) TEST_ERROR - if(HDstrcmp(simple, ptrgetTest) != 0) + if (HDstrcmp(simple, ptrgetTest) != 0) FAIL_PUTS_ERROR(" ERROR: Data transform failed to match what was set\n") - PASSED() + PASSED(); HDfree(ptrgetTest); ptrgetTest = NULL; @@ -778,7 +888,7 @@ test_getset(const hid_t dxpl_id_c_to_f) return 0; error: - if(ptrgetTest) + if (ptrgetTest) HDfree(ptrgetTest); return -1; @@ -787,17 +897,17 @@ error: static int test_set(void) { - hid_t dxpl_id = -1; + hid_t dxpl_id = -1; H5E_auto2_t func; - const char *str = "(9/5.0)*x + 32"; - char *ptrgetTest = NULL; + const char *str = "(9/5.0)*x + 32"; + char * ptrgetTest = NULL; TESTING("H5Pget_data_transform (get before set)") - if(NULL == (ptrgetTest = (char *)HDmalloc(HDstrlen(str) + 1))) + if (NULL == (ptrgetTest = (char *)HDmalloc(HDstrlen(str) + 1))) TEST_ERROR - if((dxpl_id = H5Pcreate(H5P_DATASET_XFER)) < 0) + if ((dxpl_id = H5Pcreate(H5P_DATASET_XFER)) < 0) TEST_ERROR /* Test get before set */ @@ -805,8 +915,8 @@ test_set(void) H5Eset_auto2(H5E_DEFAULT, NULL, NULL); - if(H5Pget_data_transform(dxpl_id, ptrgetTest, HDstrlen(str) + 1) < 0) - PASSED() + if (H5Pget_data_transform(dxpl_id, ptrgetTest, HDstrlen(str) + 1) < 0) + PASSED(); else FAIL_PUTS_ERROR(" ERROR: Data transform get before set succeeded (it shouldn't have)\n"); @@ -842,18 +952,19 @@ test_set(void) H5Eset_auto2(H5E_DEFAULT, func, NULL); - if(H5Pclose(dxpl_id) < 0) + if (H5Pclose(dxpl_id) < 0) TEST_ERROR return 0; error: - if(ptrgetTest) + if (ptrgetTest) HDfree(ptrgetTest); - H5E_BEGIN_TRY { + H5E_BEGIN_TRY + { H5Pclose(dxpl_id); - } H5E_END_TRY + } + H5E_END_TRY return -1; } - |