summaryrefslogtreecommitdiffstats
path: root/test
diff options
context:
space:
mode:
Diffstat (limited to 'test')
-rw-r--r--test/chunk_info.c87
1 files changed, 75 insertions, 12 deletions
diff --git a/test/chunk_info.c b/test/chunk_info.c
index d56cf88..fe585e0 100644
--- a/test/chunk_info.c
+++ b/test/chunk_info.c
@@ -18,12 +18,15 @@
* April 7, 2008
*
* Purpose: Tests the H5Dget_chunk_info API function
+ * This program writes a 4x4 did by iterating on 2x2 chunks
+ * at a time
*/
#include "hdf5.h"
#include "h5test.h"
+#define PRINT_DATA
#define H5FILE_NAME "chunk_info.h5"
#define DATASETNAME "2d"
#define RANK 2
@@ -32,33 +35,91 @@
int main( void )
{
- hid_t fid;
- hid_t did;
- hid_t sid;
- hid_t pid;
+ hid_t fid; /* file ID */
+ hid_t did; /* did ID */
+ hid_t f_sid; /* file space ID */
+ hid_t m_sid; /* memory space ID */
+ hid_t pid; /* property list ID */
+ hsize_t start[2]; /* chunk location to start writing */
hsize_t dims[2] = { 4, 4};
hsize_t chunk_dims[2] = { 2, 2 };
+ int chunk_data[2][2] = { {1, 1}, {1, 1} };
+ int buf[4][4];
int fillvalue = 0;
+ int i, j, ii, jj;
/* create a new file using default properties. */
if ((fid = H5Fcreate(H5FILE_NAME, H5F_ACC_TRUNC, H5P_DEFAULT, H5P_DEFAULT)) < 0) TEST_ERROR;
- /* create the data space */
- if ((sid = H5Screate_simple(RANK, dims, dims)) < 0) TEST_ERROR;
+ /* create the file space */
+ if ((f_sid = H5Screate_simple(RANK, dims, dims)) < 0) TEST_ERROR;
+ /* create the memory space with chunk dimensions */
+ if ((m_sid = H5Screate_simple(RANK, chunk_dims, chunk_dims)) < 0) TEST_ERROR;
+ start[0] = 0;
+ start[1] = 0;
+ if (H5Sselect_hyperslab(m_sid, H5S_SELECT_SET, start, NULL, chunk_dims, NULL) < 0) TEST_ERROR;
TESTING("chunk info");
/*-------------------------------------------------------------------------
- * create a dataset
+ * create a did
*-------------------------------------------------------------------------
*/
- /* modify dataset creation properties, i.e. enable chunking. */
+ /* modify did creation properties, i.e. enable chunking. */
if ((pid = H5Pcreate (H5P_DATASET_CREATE)) < 0) TEST_ERROR;
if (H5Pset_chunk(pid, RANK, chunk_dims) < 0) TEST_ERROR;
if (H5Pset_fill_value(pid, H5T_NATIVE_INT, &fillvalue) < 0) TEST_ERROR;
- /* create a new dataset */
- if((did = H5Dcreate2(fid , DATASETNAME, H5T_NATIVE_INT, sid, H5P_DEFAULT, pid, H5P_DEFAULT)) < 0) TEST_ERROR;
+ /* create a new did */
+ if((did = H5Dcreate2(fid , DATASETNAME, H5T_NATIVE_INT, f_sid, H5P_DEFAULT, pid, H5P_DEFAULT)) < 0) TEST_ERROR;
+
+
+ /*-------------------------------------------------------------------------
+ * write the did in 2x2 chunks
+ *-------------------------------------------------------------------------
+ */
+
+ /* iterate in dim 0 */
+ for (j = 0; j < chunk_dims[0]; j++)
+ {
+
+ /* reset start in dim 1 */
+ start[1] = 0;
+
+ /* iterate in dim 1 */
+ for (i = 0; i < chunk_dims[1]; i++)
+ {
+
+ /* select file hyperslab to save a 2x2 chunk */
+ if (H5Sselect_hyperslab(f_sid, H5S_SELECT_SET, start, NULL, chunk_dims, NULL) < 0) TEST_ERROR;
+
+ /* write the data to the hyperslab. */
+ if (H5Dwrite(did, H5T_NATIVE_INT, m_sid, f_sid, H5P_DEFAULT, chunk_data) < 0) TEST_ERROR;
+
+ /* read back and display complete dataset 4x4 */
+ if (H5Dread(did, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, buf) < 0) TEST_ERROR;
+
+
+#if defined (PRINT_DATA)
+ printf("\n");
+ printf("dataset: \n");
+ for (jj = 0; jj < dims[0]; jj++) {
+ for (ii = 0; ii < dims[1]; ii++) printf("%d ", buf[jj][ii]);
+ printf("\n");
+ }
+#endif
+
+
+ /* increment start in dim 1 */
+ start[1] += 2;
+
+
+ }
+
+ /* increment start in dim 0 */
+ start[0] += 2;
+ }
+
/*-------------------------------------------------------------------------
@@ -67,7 +128,8 @@ int main( void )
*/
if (H5Dclose(did) < 0) TEST_ERROR
- if (H5Sclose(sid) < 0) TEST_ERROR
+ if (H5Sclose(f_sid) < 0) TEST_ERROR
+ if (H5Sclose(m_sid) < 0) TEST_ERROR
if (H5Pclose(pid) < 0) TEST_ERROR
if (H5Fclose(fid) < 0) TEST_ERROR
@@ -79,7 +141,8 @@ int main( void )
error:
H5Dclose( did );
- H5Sclose( sid );
+ H5Sclose( f_sid );
+ H5Sclose( m_sid );
H5Pclose( pid );
H5Fclose( fid );
H5_FAILED();