diff options
Diffstat (limited to 'test/dtransform.c')
-rw-r--r-- | test/dtransform.c | 1002 |
1 files changed, 532 insertions, 470 deletions
diff --git a/test/dtransform.c b/test/dtransform.c index b2dcd6f..d2f53fe 100644 --- a/test/dtransform.c +++ b/test/dtransform.c @@ -6,15 +6,15 @@ * 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); @@ -26,266 +26,324 @@ static int test_set(void); static int test_getset(const 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(); \ - fprintf(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(); \ - fprintf(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(); \ - fprintf(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(); \ - fprintf(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 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) { - 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)*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; /* 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,7 +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 +#if H5_SIZEOF_LONG_DOUBLE != 0 TEST_TYPE_CONTIG(dxpl_id_c_to_f, long double, H5T_NATIVE_LDOUBLE, "ldouble", windchillFfloat, 1); #endif @@ -316,34 +374,50 @@ 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 +#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,100 +429,102 @@ 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 unsigned types: no negative numbers involved and results are < 255 to fit + * into uchar */ + const char *utrans = "((x+100)/4)*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...") - 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(); return 0; - -error: - H5E_BEGIN_TRY { + +error: + 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,229 +533,220 @@ 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] = (float) ((2.0f + windchillC) * ((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) - 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) return 0; -error: - return -1; +error: + 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"; 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) + 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 - - if(H5Pclose(dxpl_id) < 0) + if (H5Pclose(dxpl_id) < 0) TEST_ERROR - if(H5Sclose(dataspace) < 0) + if (H5Sclose(dataspace) < 0) TEST_ERROR PASSED(); return 0; -error: - return -1; +error: + 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) return 0; -error: +error: return -1; } @@ -687,28 +754,26 @@ 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() 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") } @@ -725,37 +790,36 @@ 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() - + HDfree(ptrgetTest); ptrgetTest = NULL; - TESTING("data transform, read after reseting of transform property") + 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") } @@ -763,11 +827,11 @@ test_getset(const hid_t dxpl_id_c_to_f) 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() @@ -778,7 +842,7 @@ test_getset(const hid_t dxpl_id_c_to_f) return 0; error: - if(ptrgetTest) + if (ptrgetTest) HDfree(ptrgetTest); return -1; @@ -787,17 +851,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,9 +869,9 @@ test_set(void) H5Eset_auto2(H5E_DEFAULT, NULL, NULL); - if(H5Pget_data_transform(dxpl_id, ptrgetTest, HDstrlen(str) + 1) < 0) + if (H5Pget_data_transform(dxpl_id, ptrgetTest, HDstrlen(str) + 1) < 0) PASSED() - else + else FAIL_PUTS_ERROR(" ERROR: Data transform get before set succeeded (it shouldn't have)\n"); HDfree(ptrgetTest); @@ -842,18 +906,16 @@ 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 { - H5Pclose(dxpl_id); - } H5E_END_TRY + H5E_BEGIN_TRY { H5Pclose(dxpl_id); } + H5E_END_TRY return -1; } - |