mirror of
https://github.com/HDFGroup/hdf5.git
synced 2024-11-27 02:10:55 +08:00
427ff7da28
Bug Fix/Code Cleanup/Doc Cleanup/Optimization/Branch Sync :-)
Description:
Generally speaking, this is the "signed->unsigned" change to selections.
However, in the process of merging code back, things got stickier and stickier
until I ended up doing a big "sync the two branches up" operation. So... I
brought back all the "infrastructure" fixes from the development branch to the
release branch (which I think were actually making some improvement in
performance) as well as fixed several bugs which had been fixed in one branch,
but not the other.
I've also tagged the repository before making this checkin with the label
"before_signed_unsigned_changes".
Platforms tested:
FreeBSD 4.10 (sleipnir) w/parallel & fphdf5
FreeBSD 4.10 (sleipnir) w/threadsafe
FreeBSD 4.10 (sleipnir) w/backward compatibility
Solaris 2.7 (arabica) w/"purify options"
Solaris 2.8 (sol) w/FORTRAN & C++
AIX 5.x (copper) w/parallel & FORTRAN
IRIX64 6.5 (modi4) w/FORTRAN
Linux 2.4 (heping) w/FORTRAN & C++
Misc. update:
590 lines
22 KiB
C
590 lines
22 KiB
C
#include "h5test.h"
|
|
|
|
#define ROWS 12
|
|
#define COLS 18
|
|
#define FLOAT_TOL 0.0001
|
|
|
|
int init_test(hid_t* file_id);
|
|
int test_copy(const hid_t dxpl_id_c_to_f_copy, const hid_t dxpl_id_polynomial_copy);
|
|
int test_trivial(const hid_t dxpl_id_simple);
|
|
int test_poly(const hid_t dxpl_id_polynomial);
|
|
int test_set(void);
|
|
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;
|
|
hid_t dset_id_float;
|
|
hid_t dset_id_int_chunk;
|
|
hid_t dset_id_float_chunk;
|
|
|
|
|
|
const float windchillFfloat[ROWS][COLS] =
|
|
{ {36.0, 31.0, 25.0, 19.0, 13.0, 7.0, 1.0, -5.0, -11.0, -16.0, -22.0, -28.0, -34.0, -40.0, -46.0, -52.0, -57.0, -63.0 },
|
|
{34.0, 27.0, 21.0, 15.0, 9.0, 3.0, -4.0, -10.0, -16.0, -22.0, -28.0, -35.0, -41.0, -47.0, -53.0, -59.0, -66.0, -72.0 } ,
|
|
{32.0, 25.0, 19.0, 13.0, 6.0, 0.0, -7.0, -13.0, -19.0, -26.0, -32.0, -39.0, -45.0, -51.0, -58.0, -64.0, -71.0, -77.0 },
|
|
{30.0, 24.0, 17.0, 11.0, 4.0, -2.0, -9.0, -15.0, -22.0, -29.0, -35.0, -42.0, -48.0, -55.0, -61.0, -68.0, -74.0, -81.0 },
|
|
{29.0, 23.0, 16.0, 9.0, 3.0, -4.0, -11.0, -17.0, -24.0, -31.0, -37.0, -44.0, -51.0, -58.0, -64.0, -71.0, -78.0, -84.0 },
|
|
{28.0, 22.0, 15.0, 8.0, 1.0, -5.0, -12.0, -19.0, -26.0, -33.0, -39.0, -46.0, -53.0, -60.0, -67.0, -73.0, -80.0, -87.0 },
|
|
{28.0, 21.0, 14.0, 7.0, 0.0, -7.0, -14.0, -21.0, -27.0, -34.0, -41.0, -48.0, -55.0, -62.0, -69.0, -76.0, -82.0, -89.0 },
|
|
{27.0, 20.0, 13.0, 6.0, -1.0, -8.0, -15.0, -22.0, -29.0, -36.0, -43.0, -50.0, -57.0, -64.0, -71.0, -78.0, -84.0, -91.0 },
|
|
{26.0, 19.0, 12.0, 5.0, -2.0, -9.0, -16.0, -23.0, -30.0, -37.0, -44.0, -51.0, -58.0, -65.0, -72.0, -79.0, -86.0, -93.0 },
|
|
{26.0, 19.0, 12.0, 4.0, -3.0, -10.0, -17.0, -24.0, -31.0, -38.0, -45.0, -52.0, -60.0, -67.0, -74.0, -81.0, -88.0, -95.0},
|
|
{25.0, 18.0, 11.0, 4.0, -3.0, -11.0, -18.0, -25.0, -32.0, -39.0, -46.0, -54.0, -61.0, -68.0, -75.0, -82.0, -89.0, -97.0},
|
|
{25.0, 17.0, 10.0, 3.0, -4.0, -11.0, -19.0, -26.0, -33.0, -40.0, -48.0, -55.0, -62.0, -69.0, -76.0, -84.0, -91.0, -98.0}
|
|
};
|
|
|
|
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, 72 } ,
|
|
{32, 25, 19, 13, 6, 0, 7, 13, 19, 26, 32, 39, 45, 51, 58, 64, 71, 77 },
|
|
{30, 24, 17, 11, 4, 2, 9, 15, 22, 29, 35, 42, 48, 55, 61, 68, 74, 81 },
|
|
{29, 23, 16, 9, 3, 4, 11, 17, 24, 31, 37, 44, 51, 58, 64, 71, 78, 84 },
|
|
{28, 22, 15, 8, 1, 5, 12, 19, 26, 33, 39, 46, 53, 60, 67, 73, 80, 87 },
|
|
{28, 21, 14, 7, 0, 7, 14, 21, 27, 34, 41, 48, 55, 62, 69, 76, 82, 89 },
|
|
{27, 20, 13, 6, 1, 8, 15, 22, 29, 36, 43, 50, 57, 64, 71, 78, 84, 91 },
|
|
{26, 19, 12, 5, 2, 9, 16, 23, 30, 37, 44, 51, 58, 65, 72, 79, 86, 93 },
|
|
{26, 19, 12, 4, 3, 10, 17, 24, 31, 38, 45, 52, 60, 67, 74, 81, 88, 95},
|
|
{25, 18, 11, 4, 3, 11, 18, 25, 32, 39, 46, 54, 61, 68, 75, 82, 89, 97},
|
|
{25, 17, 10, 3, 4, 11, 19, 26, 33, 40, 48, 55, 62, 69, 76, 84, 91, 98}
|
|
};
|
|
|
|
#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 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; \
|
|
hsize_t dim[2] = {ROWS, COLS}; \
|
|
\
|
|
if((dataspace = H5Screate_simple(2, dim, NULL))<0) TEST_ERROR; \
|
|
if((dset = H5Dcreate(file_id, "/transformtest_"TEST_STR, HDF_TYPE, dataspace, H5P_DEFAULT))<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(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(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) \
|
|
\
|
|
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; \
|
|
if(SIGNED) \
|
|
COMPARE(TYPE, array, COMPARE_DATA, 2) \
|
|
else \
|
|
UCOMPARE(TYPE, array, COMPARE_DATA, 4) \
|
|
} \
|
|
\
|
|
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; \
|
|
\
|
|
cparms = H5Pcreate(H5P_DATASET_CREATE); \
|
|
if(H5Pset_chunk(cparms, 2, dim) < 0) TEST_ERROR; \
|
|
\
|
|
if((dset_chunk = H5Dcreate(file_id, "/transformtest_chunk_"TEST_STR, HDF_TYPE, dataspace, cparms))<0) TEST_ERROR; \
|
|
filespace = H5Dget_space (dset_chunk); \
|
|
memspace = H5Screate_simple(2, dim, NULL); \
|
|
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; \
|
|
if(SIGNED) \
|
|
COMPARE(TYPE, array, COMPARE_DATA, 2) \
|
|
else \
|
|
UCOMPARE(TYPE, array, COMPARE_DATA, 4) \
|
|
} \
|
|
\
|
|
\
|
|
\
|
|
if(H5Dclose(dset_chunk)<0) TEST_ERROR; \
|
|
if(H5Sclose(dataspace)<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)
|
|
{
|
|
hid_t dxpl_id_c_to_f, dxpl_id_c_to_f_copy, dxpl_id_simple, dxpl_id_polynomial, dxpl_id_polynomial_copy, dxpl_id_utrans_inv, file_id;
|
|
|
|
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((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);
|
|
#ifdef H5_ULLONG_TO_FP_CAST_WORKS
|
|
TEST_TYPE_CONTIG(dxpl_id_utrans_inv, unsigned long_long, H5T_NATIVE_ULLONG, "ullong", transformData, 0);
|
|
#endif
|
|
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);
|
|
#ifdef H5_ULLONG_TO_FP_CAST_WORKS
|
|
TEST_TYPE_CHUNK(dxpl_id_utrans_inv, unsigned long_long, H5T_NATIVE_ULLONG, "ullong", transformData, 0);
|
|
#endif
|
|
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;
|
|
|
|
/* 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:
|
|
return -1;
|
|
}
|
|
|
|
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, dxpl_id_f_to_c, dxpl_id_utrans, cparms, filespace;
|
|
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;
|
|
|
|
if((*file_id = H5Fcreate("dtransform.h5", H5F_ACC_TRUNC, H5P_DEFAULT, H5P_DEFAULT))<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 = H5Dcreate(*file_id, "/default_int", H5T_NATIVE_INT, dataspace, 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 = H5Dcreate(*file_id, "/default_float", H5T_NATIVE_FLOAT, dataspace, 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 = H5Dcreate(*file_id, "/default_chunk_int", H5T_NATIVE_INT, dataspace, cparms))<0) TEST_ERROR;
|
|
|
|
filespace = H5Dget_space (dset_id_int_chunk);
|
|
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 = H5Dcreate(*file_id, "/default_chunk_float", H5T_NATIVE_FLOAT, dataspace, cparms))<0) TEST_ERROR;
|
|
if(H5Dwrite(dset_id_float_chunk, H5T_NATIVE_FLOAT, dataspace, filespace, dxpl_id_f_to_c, windchillFfloat)<0) TEST_ERROR;
|
|
|
|
|
|
PASSED();
|
|
|
|
|
|
if(H5Pclose(dxpl_id_f_to_c)<0) TEST_ERROR;
|
|
if(H5Pclose(dxpl_id_utrans)<0) TEST_ERROR;
|
|
if(H5Sclose(dataspace)<0) TEST_ERROR;
|
|
|
|
return 0;
|
|
error:
|
|
return -1;
|
|
}
|
|
|
|
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 = (5/9.0)*(windchillFfloat[row][col] - 32);
|
|
polyflres[row][col] = (2.0+windchillC)*((windchillC-8.0)/2.0);
|
|
}
|
|
}
|
|
|
|
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.0)
|
|
|
|
|
|
for(row = 0; row<ROWS; row++)
|
|
{
|
|
for(col = 0; col<COLS; col++)
|
|
{
|
|
windchillC = (5/9.0)*(windchillFfloat[row][col] - 32);
|
|
polyflres[row][col] = (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;
|
|
}
|
|
|
|
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 = (5/9.0)*(windchillFfloat[row][col] - 32);
|
|
polyflres[row][col] = (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;
|
|
}
|
|
|
|
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.8) > FLOAT_TOL)
|
|
{
|
|
H5_FAILED();
|
|
fprintf(stderr, " ERROR: Conversion failed to match computed data\n");
|
|
goto error;
|
|
}
|
|
}
|
|
}
|
|
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)
|
|
{
|
|
H5_FAILED();
|
|
fprintf(stderr, " ERROR: Conversion failed to match computed data\n");
|
|
goto error;
|
|
}
|
|
}
|
|
}
|
|
PASSED();
|
|
|
|
return 0;
|
|
error:
|
|
return -1;
|
|
}
|
|
|
|
int test_getset(const hid_t dxpl_id_c_to_f)
|
|
{
|
|
int row, 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 = HDmalloc(HDstrlen(simple)+1);
|
|
|
|
TESTING("H5Pget_data_transform")
|
|
H5Pget_data_transform(dxpl_id_c_to_f, ptrgetTest, HDstrlen(c_to_f)+1);
|
|
if(HDstrcmp(c_to_f, ptrgetTest) != 0)
|
|
{
|
|
H5_FAILED();
|
|
fprintf(stderr, " ERROR: Data transform failed to match what was set\n");
|
|
goto error;
|
|
}
|
|
else
|
|
PASSED();
|
|
|
|
if(H5Pset_data_transform(dxpl_id_c_to_f, simple)<0) TEST_ERROR;
|
|
|
|
TESTING("data transform, read after reseting of transform property")
|
|
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.8) > FLOAT_TOL)
|
|
{
|
|
H5_FAILED();
|
|
fprintf(stderr, " ERROR: Conversion failed to match computed data\n");
|
|
goto error;
|
|
}
|
|
}
|
|
}
|
|
PASSED();
|
|
|
|
HDmemset(ptrgetTest, 0, strlen(c_to_f)+1);
|
|
free(ptrgetTest);
|
|
|
|
ptrgetTest = malloc(strlen(simple)+1);
|
|
|
|
HDmemset(ptrgetTest, 0, strlen(simple)+1);
|
|
TESTING("H5Pget_data_transform, after resetting transform property")
|
|
H5Pget_data_transform(dxpl_id_c_to_f, ptrgetTest, strlen(simple)+1);
|
|
if(strcmp(simple, ptrgetTest) != 0)
|
|
{
|
|
H5_FAILED();
|
|
fprintf(stderr, " ERROR: Data transform failed to match what was set\n");
|
|
goto error;
|
|
}
|
|
else
|
|
PASSED();
|
|
|
|
free(ptrgetTest);
|
|
|
|
return 0;
|
|
|
|
error:
|
|
return -1;
|
|
}
|
|
|
|
int test_set(void)
|
|
{
|
|
hid_t dxpl_id;
|
|
H5E_auto_stack_t func;
|
|
const char* str = "(9/5.0)*x + 32";
|
|
char* ptrgetTest = malloc(strlen(str)+1);
|
|
|
|
if((dxpl_id = H5Pcreate(H5P_DATASET_XFER))<0) TEST_ERROR;
|
|
|
|
/* Test get before set */
|
|
H5Eget_auto_stack(H5E_DEFAULT,&func,NULL);
|
|
|
|
H5Eset_auto_stack(H5E_DEFAULT, NULL, NULL);
|
|
TESTING("H5Pget_data_transform (get before set)")
|
|
if(H5Pget_data_transform(dxpl_id, ptrgetTest, strlen(str)+1) < 0)
|
|
{
|
|
PASSED();
|
|
}
|
|
else
|
|
{
|
|
H5_FAILED();
|
|
fprintf(stderr, " ERROR: Data transform get before set succeeded (it shouldn't have)\n");
|
|
goto error;
|
|
}
|
|
|
|
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_auto_stack(H5E_DEFAULT, func, NULL);
|
|
|
|
return 0;
|
|
|
|
error:
|
|
return -1;
|
|
|
|
}
|