From 01c6249fbb8e777c3da8cb1715a585aed22d5c0d Mon Sep 17 00:00:00 2001
From: Pedro Vicente Nunes <pvn@hdfgroup.org>
Date: Tue, 8 Apr 2008 11:47:09 -0500
Subject: [svn-r14813] modified the test suite for H5Dget_chunk_info so that it
 writes a 4X4 dataset by iterating by 2X2 chunks

tested: windows
---
 test/chunk_info.c | 87 +++++++++++++++++++++++++++++++++++++++++++++++--------
 1 file 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();
-- 
cgit v0.12