mirror/hdf5
2
0
Fork 0
mirror of https://github.com/HDFGroup/hdf5.git synced 2024-11-27 02:10:55 +08:00
Code Issues Packages Projects Releases Wiki Activity
c72ac3ec13
hdf5/test/dtransform.c
Dana Robinson 7f1e49206d
Renamed COPYING to LICENSE (#4978) 
This is where most people will expect to find license information. The
COPYING_LBNL_HDF5 file has also been renamed to LICENSE_LBNL_HDF5.
The licenses are unchanged.
2024-10-18 21:13:04 -07:00

986 lines
52 KiB
C
Raw Blame History

/* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
* Copyright by The HDF Group. *
* All rights reserved. *
* *
* This file is part of HDF5. The full HDF5 copyright notice, including *
* terms governing use, modification, and redistribution, is contained in *
* the LICENSE file, which can be found at the root of the source code *
* distribution tree, or in https://www.hdfgroup.org/licenses. *
* If you do not have access to either file, you may request a copy from *
* help@hdfgroup.org. *
* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * */
#include "h5test.h"
#define ROWS 12
#define COLS 18
#define FLOAT_TOL 0.0001F
static int init_test(hid_t file_id);
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(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 = H5I_INVALID_HID;
hid_t dset_id_float = H5I_INVALID_HID;
hid_t dset_id_int_chunk = H5I_INVALID_HID;
hid_t dset_id_float_chunk = H5I_INVALID_HID;
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) \
do { \
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(); \
} while (0)
#define COMPARE(TYPE, VAR1, VAR2, TOL) \
do { \
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(); \
} while (0)
#define COMPARE_INT(VAR1, VAR2) \
do { \
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; \
} \
} \
} while (0)
#define TEST_TYPE_CONTIG(XFORM, TYPE, HDF_TYPE, TEST_STR, COMPARE_DATA, SIGNED) \
do { \
struct { \
TYPE arr[ROWS][COLS]; \
} *array = NULL; \
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}; \
\
/* NOTE: If this macro encounters errors, this memory will leak */ \
if (NULL == (array = calloc(1, sizeof(*array)))) \
TEST_ERROR; \
\
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->arr, COMPARE_DATA, 2); \
else \
UCOMPARE(TYPE, array->arr, 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->arr, COMPARE_DATA, 2); \
else \
UCOMPARE(TYPE, array->arr, 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->arr, COMPARE_DATA, 2); \
} \
\
if (H5Dclose(dset_nn) < 0) \
TEST_ERROR; \
if (H5Dclose(dset) < 0) \
TEST_ERROR; \
if (H5Sclose(dataspace) < 0) \
TEST_ERROR; \
\
free(array); \
} while (0)
#define TEST_TYPE_CHUNK(XFORM, TYPE, HDF_TYPE, TEST_STR, COMPARE_DATA, SIGNED) \
do { \
struct { \
TYPE arr[ROWS][COLS]; \
} *array = NULL; \
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}; \
\
/* NOTE: If this macro encounters errors, this memory will leak */ \
if (NULL == (array = calloc(1, sizeof(*array)))) \
TEST_ERROR; \
\
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->arr, COMPARE_DATA, 2); \
else \
UCOMPARE(TYPE, array->arr, 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->arr, 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; \
\
free(array); \
} while (0)
#define INVALID_SET_TEST(TRANSFORM) \
do { \
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; \
} \
} while (0)
int
main(void)
{
hid_t dxpl_id_c_to_f = H5I_INVALID_HID;
hid_t dxpl_id_c_to_f_copy = 1;
hid_t dxpl_id_simple = H5I_INVALID_HID;
hid_t dxpl_id_polynomial = H5I_INVALID_HID;
hid_t dxpl_id_polynomial_copy = H5I_INVALID_HID;
hid_t dxpl_id_utrans_inv = H5I_INVALID_HID;
hid_t file_id = H5I_INVALID_HID;
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 - 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;
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);
TEST_TYPE_CONTIG(dxpl_id_c_to_f, short, H5T_NATIVE_SHORT, "short", windchillFfloat, 1);
TEST_TYPE_CONTIG(dxpl_id_utrans_inv, unsigned short, H5T_NATIVE_USHORT, "ushort", transformData, 0);
TEST_TYPE_CONTIG(dxpl_id_c_to_f, int, H5T_NATIVE_INT, "int", windchillFfloat, 1);
TEST_TYPE_CONTIG(dxpl_id_utrans_inv, unsigned int, H5T_NATIVE_UINT, "uint", transformData, 0);
TEST_TYPE_CONTIG(dxpl_id_c_to_f, long, H5T_NATIVE_LONG, "long", windchillFfloat, 1);
TEST_TYPE_CONTIG(dxpl_id_utrans_inv, unsigned long, H5T_NATIVE_ULONG, "ulong", transformData, 0);
TEST_TYPE_CONTIG(dxpl_id_c_to_f, long long, H5T_NATIVE_LLONG, "llong", windchillFfloat, 1);
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);
TEST_TYPE_CONTIG(dxpl_id_c_to_f, long double, H5T_NATIVE_LDOUBLE, "ldouble", windchillFfloat, 1);
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);
TEST_TYPE_CHUNK(dxpl_id_c_to_f, signed char, H5T_NATIVE_SCHAR, "schar", windchillFfloat, 1);
TEST_TYPE_CHUNK(dxpl_id_c_to_f, short, H5T_NATIVE_SHORT, "short", windchillFfloat, 1);
TEST_TYPE_CHUNK(dxpl_id_utrans_inv, unsigned short, H5T_NATIVE_USHORT, "ushort", transformData, 0);
TEST_TYPE_CHUNK(dxpl_id_c_to_f, int, H5T_NATIVE_INT, "int", windchillFfloat, 1);
TEST_TYPE_CHUNK(dxpl_id_utrans_inv, unsigned int, H5T_NATIVE_UINT, "uint", transformData, 0);
TEST_TYPE_CHUNK(dxpl_id_c_to_f, long, H5T_NATIVE_LONG, "long", windchillFfloat, 1);
TEST_TYPE_CHUNK(dxpl_id_utrans_inv, unsigned long, H5T_NATIVE_ULONG, "ulong", transformData, 0);
TEST_TYPE_CHUNK(dxpl_id_c_to_f, long long, H5T_NATIVE_LLONG, "llong", windchillFfloat, 1);
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);
TEST_TYPE_CHUNK(dxpl_id_c_to_f, long double, H5T_NATIVE_LDOUBLE, "ldouble", windchillFfloat, 1);
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;
return 0;
error:
H5E_BEGIN_TRY
{
H5Dclose(dset_id_int);
H5Dclose(dset_id_int_chunk);
H5Dclose(dset_id_float);
H5Dclose(dset_id_float_chunk);
H5Fclose(file_id);
H5Pclose(dxpl_id_c_to_f);
H5Pclose(dxpl_id_c_to_f_copy);
H5Pclose(dxpl_id_polynomial);
H5Pclose(dxpl_id_polynomial_copy);
H5Pclose(dxpl_id_simple);
H5Pclose(dxpl_id_utrans_inv);
}
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 char types: numbers are restricted from -128 to 127, fits into char */
const char *utrans = "(x/4+25)*3";
hid_t dataspace = H5I_INVALID_HID;
hid_t dxpl_id_f_to_c = H5I_INVALID_HID;
hid_t dxpl_id_utrans = H5I_INVALID_HID;
hid_t cparms = H5I_INVALID_HID;
hid_t filespace = H5I_INVALID_HID;
hsize_t dim[2] = {ROWS, COLS};
hsize_t offset[2] = {0, 0};
if ((dxpl_id_f_to_c = H5Pcreate(H5P_DATASET_XFER)) < 0)
TEST_ERROR;
if ((dxpl_id_utrans = H5Pcreate(H5P_DATASET_XFER)) < 0)
TEST_ERROR;
if (H5Pset_data_transform(dxpl_id_f_to_c, f_to_c) < 0)
TEST_ERROR;
if (H5Pset_data_transform(dxpl_id_utrans, utrans) < 0)
TEST_ERROR;
cparms = H5Pcreate(H5P_DATASET_CREATE);
if (H5Pset_chunk(cparms, 2, dim) < 0)
TEST_ERROR;
if ((dataspace = H5Screate_simple(2, dim, NULL)) < 0)
TEST_ERROR;
TESTING("Initializing test...");
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)
TEST_ERROR;
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)
TEST_ERROR;
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)
TEST_ERROR;
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)
TEST_ERROR;
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)
TEST_ERROR;
if (H5Pclose(cparms) < 0)
TEST_ERROR;
if (H5Pclose(dxpl_id_f_to_c) < 0)
TEST_ERROR;
if (H5Pclose(dxpl_id_utrans) < 0)
TEST_ERROR;
if (H5Sclose(dataspace) < 0)
TEST_ERROR;
if (H5Sclose(filespace) < 0)
TEST_ERROR;
PASSED();
return 0;
error:
H5E_BEGIN_TRY
{
H5Pclose(cparms);
H5Pclose(dxpl_id_f_to_c);
H5Pclose(dxpl_id_utrans);
H5Sclose(dataspace);
H5Sclose(filespace);
}
H5E_END_TRY
return -1;
}
static int
test_poly(const hid_t dxpl_id_polynomial)
{
float polyflres[ROWS][COLS];
int polyintread[ROWS][COLS];
float polyflread[ROWS][COLS];
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));
}
TESTING("data transform, polynomial transform (int->float)");
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));
}
TESTING("data transform, polynomial transform (float->int)");
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;
}
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";
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)
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;
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)
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;
/*-----------------------------
* Operation 2: 100-x
*----------------------------*/
if (H5Pset_data_transform(dxpl_id, special2) < 0)
TEST_ERROR;
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)
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;
/*-----------------------------
* Operation 3: 1000/x
*----------------------------*/
if (H5Pset_data_transform(dxpl_id, special3) < 0)
TEST_ERROR;
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)
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;
/*-----------------------------
* Operation 4: -x
*----------------------------*/
if (H5Pset_data_transform(dxpl_id, special4) < 0)
TEST_ERROR;
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)
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;
/*-----------------------------
* Operation 5: +x
*----------------------------*/
if (H5Pset_data_transform(dxpl_id, special5) < 0)
TEST_ERROR;
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)
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;
/*-----------------------------
* 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)
TEST_ERROR;
COMPARE_INT(read_buf, data_res);
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 ((dset_id = H5Dcreate2(file, "/special7", 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 (H5Pclose(dxpl_id) < 0)
TEST_ERROR;
if (H5Sclose(dataspace) < 0)
TEST_ERROR;
PASSED();
return 0;
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;
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));
}
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)
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)
TEST_ERROR;
COMPARE(int, polyintread, polyflres, 2);
return 0;
error:
return -1;
}
static int
test_trivial(const hid_t dxpl_id_simple)
{
float windchillFfloatread[ROWS][COLS];
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)
TEST_ERROR;
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)
TEST_ERROR;
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();
return 0;
error:
return -1;
}
static int
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;
TESTING("H5Pget_data_transform");
if (NULL == (ptrgetTest = (char *)malloc(strlen(simple) + 1)))
TEST_ERROR;
if (H5Pget_data_transform(dxpl_id_c_to_f, ptrgetTest, strlen(c_to_f) + 1) < 0)
TEST_ERROR;
if (strcmp(c_to_f, ptrgetTest) != 0)
FAIL_PUTS_ERROR(" ERROR: Data transform failed to match what was set\n");
PASSED();
free(ptrgetTest);
ptrgetTest = NULL;
TESTING("data transform, read after resetting of transform property");
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)
TEST_ERROR;
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("H5Pget_data_transform, after resetting transform property");
if (NULL == (ptrgetTest = (char *)calloc((size_t)1, strlen(simple) + 1)))
TEST_ERROR;
if (H5Pget_data_transform(dxpl_id_c_to_f, ptrgetTest, strlen(simple) + 1) < 0)
TEST_ERROR;
if (strcmp(simple, ptrgetTest) != 0)
FAIL_PUTS_ERROR(" ERROR: Data transform failed to match what was set\n");
PASSED();
free(ptrgetTest);
ptrgetTest = NULL;
return 0;
error:
if (ptrgetTest)
free(ptrgetTest);
return -1;
}
static int
test_set(void)
{
hid_t dxpl_id = H5I_INVALID_HID;
H5E_auto2_t func;
const char *str = "(9/5.0)*x + 32";
char *ptrgetTest = NULL;
TESTING("H5Pget_data_transform (get before set)");
if (NULL == (ptrgetTest = (char *)malloc(strlen(str) + 1)))
TEST_ERROR;
if ((dxpl_id = H5Pcreate(H5P_DATASET_XFER)) < 0)
TEST_ERROR;
/* Test get before set */
H5Eget_auto2(H5E_DEFAULT, &func, NULL);
H5Eset_auto2(H5E_DEFAULT, NULL, NULL);
if (H5Pget_data_transform(dxpl_id, ptrgetTest, strlen(str) + 1) < 0)
PASSED();
else
FAIL_PUTS_ERROR(" ERROR: Data transform get before set succeeded (it shouldn't have)\n");
free(ptrgetTest);
ptrgetTest = NULL;
TESTING("H5Pset_data_transform (set with NULL transform)");
INVALID_SET_TEST(NULL);
TESTING("H5Pset_data_transform (set with invalid transform 1)");
INVALID_SET_TEST("\0");
TESTING("H5Pset_data_transform (set with invalid transform 2)");
INVALID_SET_TEST(" ");
TESTING("H5Pset_data_transform (set with invalid transform 3)");
INVALID_SET_TEST("x+");
TESTING("H5Pset_data_transform (set with invalid transform 4)");
INVALID_SET_TEST("(x+5");
TESTING("H5Pset_data_transform (set with invalid transform 5)");
INVALID_SET_TEST("+");
TESTING("H5Pset_data_transform (set with invalid transform 6)");
INVALID_SET_TEST("(9/5)*x + x**2");
TESTING("H5Pset_data_transform (set with invalid transform 7)");
INVALID_SET_TEST("(9/5)x");
TESTING("H5Pset_data_transform (set with invalid transform 8)");
INVALID_SET_TEST("(9/5)*x + x^2");
H5Eset_auto2(H5E_DEFAULT, func, NULL);
if (H5Pclose(dxpl_id) < 0)
TEST_ERROR;
return 0;
error:
if (ptrgetTest)
free(ptrgetTest);
H5E_BEGIN_TRY
{
H5Pclose(dxpl_id);
}
H5E_END_TRY
return -1;
}
Reference in New Issue View Git Blame Copy Permalink