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/gen_cross.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/gen_cross.c')
-rwxr-xr-x | test/gen_cross.c | 251 |
1 files changed, 227 insertions, 24 deletions
diff --git a/test/gen_cross.c b/test/gen_cross.c index 46afc35..9b9b425 100755 --- a/test/gen_cross.c +++ b/test/gen_cross.c @@ -27,17 +27,39 @@ #include "h5test.h" #define H5FILE_NAME "data.h5" -#define DATASETNAME "Array" -#define NX 5 /* dataset 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 +#define RANK 2 +#define CHUNK0 3 +#define CHUNK1 3 +int create_normal_dset(hid_t fid, hid_t sid); +int create_scale_offset_dset_float(hid_t fid, hid_t sid); +int create_scale_offset_dset_int(hid_t fid, hid_t sid); + + +/*------------------------------------------------------------------------- + * Function: create_normal_dset + * + * Purpose: Create a regular dataset of FLOAT datatype. + * + * Return: Success: 0 + * Failure: -1 + * + * Programmer: Raymond Lu + * Some time ago + * + * Modifications: + * + *------------------------------------------------------------------------- + */ int -main (void) +create_normal_dset(hid_t fid, hid_t sid) { - hid_t file, dataset; /* file and dataset handles */ - hid_t datatype, dataspace; /* handles */ - hsize_t dimsf[2]; /* dataset dimensions */ + hid_t dataset; /* file and dataset handles */ herr_t status; float data[NX][NY]; /* data to write */ int i, j; @@ -55,35 +77,79 @@ main (void) * 2 3 4 5 6 7 * 3 4 5 6 7 8 * 4 5 6 7 8 9 + * 5 6 7 8 9 10 */ /* - * Create a new file using H5F_ACC_TRUNC access, - * default file creation properties, and default file - * access properties. + * Create a new dataset within the file using defined dataspace and + * datatype and default dataset creation properties. */ - file = H5Fcreate(H5FILE_NAME, H5F_ACC_TRUNC, H5P_DEFAULT, H5P_DEFAULT); + dataset = H5Dcreate2(fid, DATASETNAME, H5T_NATIVE_FLOAT, sid, + H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT); /* - * Describe the size of the array and create the data space for fixed - * size dataset. + * Write the data to the dataset using default transfer properties. */ - dimsf[0] = NX; - dimsf[1] = NY; - dataspace = H5Screate_simple(RANK, dimsf, NULL); + status = H5Dwrite(dataset, H5T_NATIVE_FLOAT, H5S_ALL, H5S_ALL, + H5P_DEFAULT, data); /* - * Define datatype for the data in the file. - * We will store little endian INT numbers. + * Close/release resources. */ - datatype = H5Tcopy(H5T_NATIVE_DOUBLE); + H5Dclose(dataset); + + return 0; +} + + +/*------------------------------------------------------------------------- + * Function: create_scale_offset_dset_float + * + * Purpose: Create a dataset of FLOAT datatype with scale-offset filter + * + * Return: Success: 0 + * Failure: -1 + * + * Programmer: Raymond Lu + * 21 January 2011 + * + * Modifications: + * + *------------------------------------------------------------------------- + */ +int +create_scale_offset_dset_float(hid_t fid, hid_t sid) +{ +#ifdef H5_HAVE_FILTER_SCALEOFFSET + hid_t dataset; /* dataset handles */ + hid_t dcpl; + herr_t status; + float data[NX][NY]; /* data to write */ + hsize_t chunk[RANK] = {CHUNK0, CHUNK1}; + int i, j; + + /* + * Data and output buffer initialization. + */ + for (j = 0; j < NX; j++) { + for (i = 0; i < NY; i++) + data[j][i] = ((float)(i + j + 1))/3; + } + + /* + * Create the dataset creation property list, add the Scale-Offset + * filter and set the chunk size. + */ + dcpl = H5Pcreate (H5P_DATASET_CREATE); + status = H5Pset_scaleoffset (dcpl, H5Z_SO_FLOAT_DSCALE, 3); + status = H5Pset_chunk (dcpl, RANK, chunk); /* * Create a new dataset within the file using defined dataspace and * datatype and default dataset creation properties. */ - dataset = H5Dcreate2(file, DATASETNAME, datatype, dataspace, - H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT); + dataset = H5Dcreate2(fid, DATASETNAME2, H5T_NATIVE_FLOAT, sid, + H5P_DEFAULT, dcpl, H5P_DEFAULT); /* * Write the data to the dataset using default transfer properties. @@ -94,9 +160,146 @@ main (void) /* * Close/release resources. */ - H5Sclose(dataspace); - H5Tclose(datatype); + H5Pclose(dcpl); + H5Dclose(dataset); + +#else /* H5_HAVE_FILTER_SCALEOFFSET */ + const char *not_supported= "Scaleoffset filter is not enabled. Can't create the dataset."; + + puts(not_supported); +#endif /* H5_HAVE_FILTER_SCALEOFFSET */ + + return 0; +} + + +/*------------------------------------------------------------------------- + * Function: create_scale_offset_dset_int + * + * Purpose: Create a dataset of INT datatype with scale-offset filter + * + * Return: Success: 0 + * Failure: -1 + * + * Programmer: Raymond Lu + * 21 January 2011 + * + * Modifications: + * + *------------------------------------------------------------------------- + */ +int +create_scale_offset_dset_int(hid_t fid, hid_t sid) +{ +#ifdef H5_HAVE_FILTER_SCALEOFFSET + hid_t dataset; /* dataset handles */ + hid_t dcpl; + herr_t status; + int data[NX][NY]; /* data to write */ + hsize_t chunk[RANK] = {CHUNK0, CHUNK1}; + int i, j; + + /* + * Data and output buffer initialization. + */ + for (j = 0; j < NX; j++) { + for (i = 0; i < NY; i++) + data[j][i] = i + j; + } + /* + * 0 1 2 3 4 5 + * 1 2 3 4 5 6 + * 2 3 4 5 6 7 + * 3 4 5 6 7 8 + * 4 5 6 7 8 9 + * 5 6 7 8 9 10 + */ + + /* + * Create the dataset creation property list, add the Scale-Offset + * filter and set the chunk size. + */ + dcpl = H5Pcreate (H5P_DATASET_CREATE); + status = H5Pset_scaleoffset (dcpl, H5Z_SO_INT, H5Z_SO_INT_MINBITS_DEFAULT); + status = H5Pset_chunk (dcpl, RANK, chunk); + + /* + * Create a new dataset within the file using defined dataspace and + * datatype and default dataset creation properties. + */ + dataset = H5Dcreate2(fid, DATASETNAME3, H5T_NATIVE_INT, sid, + H5P_DEFAULT, dcpl, H5P_DEFAULT); + + /* + * Write the data to the dataset using default transfer properties. + */ + status = H5Dwrite(dataset, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, + H5P_DEFAULT, data); + + /* + * Close/release resources. + */ + H5Pclose(dcpl); H5Dclose(dataset); + +#else /* H5_HAVE_FILTER_SCALEOFFSET */ + const char *not_supported= "Scaleoffset filter is not enabled. Can't create the dataset."; + + puts(not_supported); +#endif /* H5_HAVE_FILTER_SCALEOFFSET */ + + return 0; +} + + +/*------------------------------------------------------------------------- + * Function: main + * + * Purpose: Create a file for cross_read.c test. + * + * Return: Success: exit(0) + * Failure: exit(1) + * + * Programmer: Raymond Lu + * Some time ago + * + *------------------------------------------------------------------------- + */ +int +main (void) +{ + hid_t file; /* file and dataset handles */ + hid_t dataspace; + hsize_t dimsf[RANK]; + + /* + * Create a new file using H5F_ACC_TRUNC access, + * default file creation properties, and default file + * access properties. + */ + file = H5Fcreate(H5FILE_NAME, H5F_ACC_TRUNC, H5P_DEFAULT, H5P_DEFAULT); + + /* + * Describe the size of the array and create the data space for fixed + * size dataset. + */ + dimsf[0] = NX; + dimsf[1] = NY; + dataspace = H5Screate_simple(RANK, dimsf, NULL); + + /* Create a regular dataset */ + create_normal_dset(file, dataspace); + + /* Create a dataset of FLOAT with scale-offset filter */ + create_scale_offset_dset_float(file, dataspace); + + /* Create a dataset of INT with scale-offset filter */ + create_scale_offset_dset_int(file, dataspace); + + /* + * Close/release resources. + */ + H5Sclose(dataspace); H5Fclose(file); return 0; |