diff options
author | Raymond Lu <songyulu@hdfgroup.org> | 2011-01-24 20:39:17 (GMT) |
---|---|---|
committer | Raymond Lu <songyulu@hdfgroup.org> | 2011-01-24 20:39:17 (GMT) |
commit | c7782121eedc332660a85ce1a4fa97ed3289cf49 (patch) | |
tree | d10614114b716815ca84291afe7e08a8360370d1 /test/cross_read.c | |
parent | f4117c14361c96f99ccff22de1939bb2454e8c99 (diff) | |
download | hdf5-c7782121eedc332660a85ce1a4fa97ed3289cf49.zip hdf5-c7782121eedc332660a85ce1a4fa97ed3289cf49.tar.gz hdf5-c7782121eedc332660a85ce1a4fa97ed3289cf49.tar.bz2 |
[svn-r19983] I added a test case for dataset with scale-offset filter into cross_read.c and updated the data files from BE, LE, and VMS. This is the test for bug 2131 that dataset of FLOAT with scale-offset filter can't be read across BE and LE.
Tested on jam and linew.
Diffstat (limited to 'test/cross_read.c')
-rwxr-xr-x | test/cross_read.c | 215 |
1 files changed, 199 insertions, 16 deletions
diff --git a/test/cross_read.c b/test/cross_read.c index 3c64f26..92802af 100755 --- a/test/cross_read.c +++ b/test/cross_read.c @@ -17,8 +17,8 @@ * Programmer: Raymond Lu <slu@ncsa.uiuc.edu> * Thursday, March 23, 2006 * - * Purpose: Check if floating-point data created on OpenVMS (VAX type), Solaris, - * and Linux machines can be read on the machine running this test. + * Purpose: Check if floating-point data created on OpenVMS, big-endian, and + * little-endian machines can be read on the machine running this test. */ #include "h5test.h" @@ -31,11 +31,28 @@ const char *FILENAME[] = { NULL }; -#define DATASETNAME "Array" -#define NX 5 /* output buffer dimensions */ -#define NY 6 -#define RANK 2 +#define DATASETNAME "Array" +#define DATASETNAME2 "Scale_offset_float_data" +#define DATASETNAME3 "Scale_offset_int_data" +#define NX 6 +#define NY 6 + +/*------------------------------------------------------------------------- + * Function: read_data + * + * Purpose: Read data from a data file. + * + * Return: Success: 0 + * Failure: -1 + * + * Programmer: Raymond Lu + * 21 January 2011 + * + * Modifications: + * + *------------------------------------------------------------------------- + */ static int read_data(char *fname) { char pathname[1024]; @@ -43,10 +60,14 @@ static int read_data(char *fname) hid_t file, dataset; /* handles */ hid_t datatype; hid_t dt; - double data_in[NX][NY]; /* input buffer */ - double data_out[NX][NY]; /* output buffer */ + float data_in[NX][NY]; /* input buffer */ + float data_out[NX][NY]; /* output buffer */ + int int_data_in[NX][NY]; /* input buffer */ + int int_data_out[NX][NY]; /* output buffer */ int i, j; unsigned nerrors = 0; + const char *not_supported= " Scaleoffset filter is not enabled."; + const char *not_fixed= " Scaleoffset filter bug (2131) is not fixed yet."; pathname[0] = '\0'; /* Generate correct name for test file by prepending the source path */ @@ -57,6 +78,20 @@ static int read_data(char *fname) strcat(pathname, fname); /* + * Open the file. + */ + if((file = H5Fopen(pathname, H5F_ACC_RDONLY, H5P_DEFAULT)) < 0) + TEST_ERROR; + + TESTING(" regular dataset"); + + /* + * Open the regular dataset. + */ + if((dataset = H5Dopen2(file, DATASETNAME, H5P_DEFAULT)) < 0) + TEST_ERROR; + + /* * Data and output buffer initialization. */ for (j = 0; j < NX; j++) { @@ -71,17 +106,74 @@ static int read_data(char *fname) * 2 3 4 5 6 7 * 3 4 5 6 7 8 * 4 5 6 7 8 9 + * 5 6 7 8 9 10 + */ + + /* + * Get datatype and dataspace handles and then query + * dataset class, order, size, rank and dimensions. */ + if((dt = H5Dget_type(dataset)) < 0) /* datatype handle */ + TEST_ERROR; + if((datatype = H5Tget_native_type(dt, H5T_DIR_DEFAULT)) < 0) + TEST_ERROR; /* - * Open the file and the dataset. + * Read data from hyperslab in the file into the hyperslab in + * memory and display. */ - if((file = H5Fopen(pathname, H5F_ACC_RDONLY, H5P_DEFAULT)) < 0) + if(H5Dread(dataset, datatype, H5S_ALL, H5S_ALL, H5P_DEFAULT, data_out) < 0) TEST_ERROR; - if((dataset = H5Dopen2(file, DATASETNAME, H5P_DEFAULT)) < 0) + + /* Check results */ + for (j=0; j<NX; j++) { + for (i=0; i<NY; i++) { + /* if (data_out[j][i] != data_in[j][i]) { */ + if (!DBL_ABS_EQUAL(data_out[j][i], data_in[j][i])) { + if (!nerrors++) { + H5_FAILED(); + printf("element [%d][%d] is %g but should have been %g\n", + j, i, data_out[j][i], data_in[j][i]); + } + } + } + } + + /* + * Close/release resources. + */ + H5Tclose(dt); + H5Tclose(datatype); + H5Dclose(dataset); + + /* Failure */ + if (nerrors) { + printf("total of %d errors out of %d elements\n", nerrors, NX*NY); + return 1; + } + + PASSED(); + + TESTING(" dataset of FLOAT with scale-offset filter"); +#ifdef TMP +#ifdef H5_HAVE_FILTER_SCALEOFFSET + /* + * Open the dataset with scale-offset filter. + */ + if((dataset = H5Dopen2(file, DATASETNAME2, H5P_DEFAULT)) < 0) TEST_ERROR; /* + * Data and output buffer initialization. + */ + for (j = 0; j < NX; j++) { + for (i = 0; i < NY; i++) { + data_in[j][i] = ((float)(i + j + 1))/3; + data_out[j][i] = 0; + } + } + + /* * Get datatype and dataspace handles and then query * dataset class, order, size, rank and dimensions. */ @@ -100,7 +192,7 @@ static int read_data(char *fname) /* Check results */ for (j=0; j<NX; j++) { for (i=0; i<NY; i++) { - if (data_out[j][i] != data_in[j][i]) { + if (!DBL_REL_EQUAL(data_out[j][i], data_in[j][i], 0.001)) { if (!nerrors++) { H5_FAILED(); printf("element [%d][%d] is %g but should have been %g\n", @@ -116,7 +208,6 @@ static int read_data(char *fname) H5Tclose(dt); H5Tclose(datatype); H5Dclose(dataset); - H5Fclose(file); /* Failure */ if (nerrors) { @@ -125,6 +216,84 @@ static int read_data(char *fname) } PASSED(); +#else /*H5_HAVE_FILTER_SCALEOFFSET*/ + SKIPPED(); + puts(not_supported); +#endif /*H5_HAVE_FILTER_SCALEOFFSET*/ +#else /*TMP*/ + SKIPPED(); + puts(not_fixed); +#endif /*TMP*/ + + TESTING(" dataset of INT with scale-offset filter"); + +#ifdef H5_HAVE_FILTER_SCALEOFFSET + /* + * Open the dataset with scale-offset filter. + */ + if((dataset = H5Dopen2(file, DATASETNAME3, H5P_DEFAULT)) < 0) + TEST_ERROR; + + /* + * Data and output buffer initialization. + */ + for (j = 0; j < NX; j++) { + for (i = 0; i < NY; i++) { + int_data_in[j][i] = i + j; + int_data_out[j][i] = 0; + } + } + + /* + * Get datatype and dataspace handles and then query + * dataset class, order, size, rank and dimensions. + */ + if((dt = H5Dget_type(dataset)) < 0) /* datatype handle */ + TEST_ERROR; + if((datatype = H5Tget_native_type(dt, H5T_DIR_DEFAULT)) < 0) + TEST_ERROR; + + /* + * Read data from hyperslab in the file into the hyperslab in + * memory and display. + */ + if(H5Dread(dataset, datatype, H5S_ALL, H5S_ALL, H5P_DEFAULT, int_data_out) < 0) + TEST_ERROR; + + /* Check results */ + for (j=0; j<NX; j++) { + for (i=0; i<NY; i++) { + if (int_data_out[j][i] != int_data_in[j][i]) { + if (!nerrors++) { + H5_FAILED(); + printf("element [%d][%d] is %d but should have been %d\n", + j, i, int_data_out[j][i], int_data_in[j][i]); + } + } + } + } + + /* + * Close/release resources. + */ + H5Tclose(dt); + H5Tclose(datatype); + H5Dclose(dataset); + + /* Failure */ + if (nerrors) { + printf("total of %d errors out of %d elements\n", nerrors, NX*NY); + return 1; + } + + PASSED(); + +#else /*H5_HAVE_FILTER_SCALEOFFSET*/ + SKIPPED(); + puts(not_supported); +#endif /*H5_HAVE_FILTER_SCALEOFFSET*/ + + H5Fclose(file); return 0; error: @@ -134,6 +303,20 @@ error: return 1; } + +/*------------------------------------------------------------------------- + * Function: main + * + * Purpose: Tests the basic features of Virtual File Drivers + * + * Return: Success: exit(0) + * Failure: exit(1) + * + * Programmer: Raymond Lu + * Tuesday, Sept 24, 2002 + * + *------------------------------------------------------------------------- + */ int main(void) { char filename[1024]; @@ -141,15 +324,15 @@ int main(void) h5_reset(); - TESTING("reading data created on OpenVMS"); + puts("Testing reading data created on OpenVMS"); h5_fixname(FILENAME[0], H5P_DEFAULT, filename, sizeof filename); nerrors += read_data(filename); - TESTING("reading data created on Linux"); + puts("Testing reading data created on Linux"); h5_fixname(FILENAME[1], H5P_DEFAULT, filename, sizeof filename); nerrors += read_data(filename); - TESTING("reading data created on Solaris"); + puts("Testing reading data created on Solaris"); h5_fixname(FILENAME[2], H5P_DEFAULT, filename, sizeof filename); nerrors += read_data(filename); |