path: root/test/dtransform.c

/* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
 * Copyright by The HDF Group.                                               *
 * Copyright by the Board of Trustees of the University of Illinois.         *
 * All rights reserved.                                                      *
 *                                                                           *
 * This file is part of HDF5.  The full HDF5 copyright notice, including     *
 * terms governing use, modification, and redistribution, is contained in    *
 * the COPYING file, which can be found at the root of the source code       *
 * distribution tree, or in https://www.hdfgroup.org/licenses.               *
 * If you do not have access to either file, you may request a copy from     *
 * help@hdfgroup.org.                                                        *
 * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * */

#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(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_specials(hid_t file);
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_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)
{
    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;

    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;

    /* 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);
#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);
    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);
#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;

    /* 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 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   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};

    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("Intializing 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";

    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

    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 *)HDmalloc(HDstrlen(simple) + 1)))
        TEST_ERROR

    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)
        FAIL_PUTS_ERROR("    ERROR: Data transform failed to match what was set\n")

    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)
        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 *)HDcalloc((size_t)1, HDstrlen(simple) + 1)))
        TEST_ERROR
    if (H5Pget_data_transform(dxpl_id_c_to_f, ptrgetTest, HDstrlen(simple) + 1) < 0)
        TEST_ERROR
    if (HDstrcmp(simple, ptrgetTest) != 0)
        FAIL_PUTS_ERROR("    ERROR: Data transform failed to match what was set\n")

    PASSED()

    HDfree(ptrgetTest);
    ptrgetTest = NULL;

    return 0;

error:
    if (ptrgetTest)
        HDfree(ptrgetTest);

    return -1;
}

static int
test_set(void)
{
    hid_t       dxpl_id = -1;
    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 *)HDmalloc(HDstrlen(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, HDstrlen(str) + 1) < 0)
        PASSED()
    else
        FAIL_PUTS_ERROR("    ERROR: Data transform get before set succeeded (it shouldn't have)\n");

    HDfree(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)
        HDfree(ptrgetTest);
    H5E_BEGIN_TRY { H5Pclose(dxpl_id); }
    H5E_END_TRY

    return -1;
}