[svn-r16614] 3. #1501 (B1) tools bug if dataset is larger than H5TOOLS_BUFSIZE limit.

ISSUE : the tools use the following formula to read by hyperslabs: hyperslab_size[i] = MIN( dim_size[i], H5TOOLS_BUFSIZE / datum_size) where H5TOOLS_BUFSIZE is a constant defined of 1024K. This is OK as long as the datum_size does not exceed 1024K, otherwise we have a hyperslab size of 0 (since 1024K/(greater than 1024K) = 0). This affects h5dump. h5repack, h5diff
SOLUTION: add a check for a 0 size and define as 1 if so. 
TEST FOR H5DUMP: Defined a case in the h5dump test generator program of such a type (an array type of doubles with a large array dimension, that was the case the user reported).  Since the written file commited in svn would be around 1024K, opted for not writing the data (the part of the code where the hyperslab is defined is executed, since h5dump always reads the files). Defined a macro WRITE_ARRAY to enable such writing if needed. Added a run on the h5dump shell script. Added 2 new files to svn: tools/testfiles/tarray8.ddl, tools/testfiles/tarray8.h5. NOTE: while doing this I thought of adding this dataset case to an existing file, but that would add the large array output to those files (the ddls). The issue is that the file list is increasing.
TEST FOR H5DIFF: for h5diff the check for reading by hyperslabs is H5TOOLS_MALLOCSIZE (128 * H5TOOLS_BUFSIZE) or 128 Mb. This makes it not possible to add such a file to svn, so used the same method as h5dump (only write the dataset if WRITE_ARRAY is defined). As opposed to h5dump, the hyperslab code is NOT executed when the dataset is empty (dataset is not read). Added the new dataset to existing files and shell run (tools/h5diff/testfiles/h5diff_dset1.h5 and tools/h5diff/testfiles/h5diff_dset2.h5 and output in tools/h5diff/testfiles/h5diff_80.txt).
TEST FOR H5REPACK: similar issue as h5diff with the difference that the hyperslab code is run. Added a run to the shell script (with a filter, otherwise the code uses H5Ocopy). 

tested: linux (h5commitest failed , apparently it did not detect the code changes in /tools/lib that fix the bug: the error in an assertion in the hyperslab of 0. I am sure that making h5ccomitest --distclean will detect the new code , but don't want to wait more 3 hours :-) )

1 files changed, 55 insertions, 13 deletions

diff --git a/tools/h5diff/h5diffgentest.c b/tools/h5diff/h5diffgentest.c
index 181e242..92e5493 100644
--- a/tools/h5diff/h5diffgentest.c
+++ b/tools/h5diff/h5diffgentest.c
@@ -17,6 +17,8 @@
 #include <stdlib.h>
 #include "hdf5.h"
 #include "H5private.h"
+#include "h5tools.h"
+
 
 /*-------------------------------------------------------------------------
  * Program: h5diffgentest
@@ -153,12 +155,12 @@ int test_basic(const char *fname1, const char *fname2, const char *fname3)
         goto out;
     if (( fid2 = H5Fcreate (fname2, H5F_ACC_TRUNC, H5P_DEFAULT, H5P_DEFAULT)) < 0 )
         goto out;
-    
-        /*-------------------------------------------------------------------------
-        * create groups
-        *-------------------------------------------------------------------------
+
+    /*-------------------------------------------------------------------------
+    * create groups
+    *-------------------------------------------------------------------------
     */
-    
+
     gid1 = H5Gcreate2(fid1, "g1", H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
     gid2 = H5Gcreate2(fid2, "g1", H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
     gid3 = H5Gcreate2(fid2, "g2", H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
@@ -600,42 +602,42 @@ int test_datatypes(const char *fname)
     write_dset(fid1,2,dims,"dset0b",H5T_STD_I32LE,buf3b);
     
     /*-------------------------------------------------------------------------
-    * Check H5T_NATIVE_CHAR
+    * H5T_NATIVE_CHAR
     *-------------------------------------------------------------------------
     */
     write_dset(fid1,2,dims,"dset1a",H5T_NATIVE_CHAR,buf1a);
     write_dset(fid1,2,dims,"dset1b",H5T_NATIVE_CHAR,buf1b);
     
     /*-------------------------------------------------------------------------
-    * Check H5T_NATIVE_SHORT
+    * H5T_NATIVE_SHORT
     *-------------------------------------------------------------------------
     */
     write_dset(fid1,2,dims,"dset2a",H5T_NATIVE_SHORT,buf2a);
     write_dset(fid1,2,dims,"dset2b",H5T_NATIVE_SHORT,buf2b);
     
     /*-------------------------------------------------------------------------
-    * Check H5T_NATIVE_INT
+    * H5T_NATIVE_INT
     *-------------------------------------------------------------------------
     */
     write_dset(fid1,2,dims,"dset3a",H5T_NATIVE_INT,buf3a);
     write_dset(fid1,2,dims,"dset3b",H5T_NATIVE_INT,buf3b);
     
     /*-------------------------------------------------------------------------
-    * Check H5T_NATIVE_LONG
+    * H5T_NATIVE_LONG
     *-------------------------------------------------------------------------
     */
     write_dset(fid1,2,dims,"dset4a",H5T_NATIVE_LONG,buf4a);
     write_dset(fid1,2,dims,"dset4b",H5T_NATIVE_LONG,buf4b);
     
     /*-------------------------------------------------------------------------
-    * Check H5T_NATIVE_FLOAT
+    * H5T_NATIVE_FLOAT
     *-------------------------------------------------------------------------
     */
     write_dset(fid1,2,dims,"dset5a",H5T_NATIVE_FLOAT,buf5a);
     write_dset(fid1,2,dims,"dset5b",H5T_NATIVE_FLOAT,buf5b);
     
     /*-------------------------------------------------------------------------
-    * Check H5T_NATIVE_DOUBLE
+    * H5T_NATIVE_DOUBLE
     *-------------------------------------------------------------------------
     */
     
@@ -1909,7 +1911,7 @@ void write_dset_in(hid_t loc_id,
     typedef enum
     {
         RED,
-            GREEN
+        GREEN
     } e_t;
     
     hid_t   did;
@@ -2155,7 +2157,47 @@ void write_dset_in(hid_t loc_id,
     tid = H5Tarray_create2(H5T_NATIVE_INT, 1, dimarray);
     write_dset(loc_id, 1, dims, "array", tid, buf6);
     status = H5Tclose(tid);
-    
+
+    {
+
+        double  *dbuf;                           /* information to write */
+        hid_t   did;                             /* dataset ID   */
+        hid_t   sid;                             /* dataspace ID   */
+        hid_t   tid;                             /* datatype ID   */
+        size_t  size;
+        hsize_t sdims[] = {1};
+        hsize_t tdims[] = {H5TOOLS_MALLOCSIZE / sizeof(double) + 1};
+        int     j;
+
+        /* allocate and initialize array data to write */
+        size = ( H5TOOLS_MALLOCSIZE / sizeof(double) + 1 ) * sizeof(double);
+        dbuf = malloc( size );
+
+        for( j = 0; j < H5TOOLS_MALLOCSIZE / sizeof(double) + 1; j++)
+            dbuf[j] = j;
+
+        if (make_diffs)
+        {
+            dbuf[5] = 0;
+            dbuf[6] = 0;
+        }
+
+        /* create a type larger than H5TOOLS_MALLOCSIZE */
+        tid = H5Tarray_create2(H5T_NATIVE_DOUBLE, 1, tdims);
+        size = H5Tget_size(tid);
+        sid = H5Screate_simple(1, sdims, NULL);
+        did = H5Dcreate2(loc_id, "arrayd", tid, sid, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
+#if defined(WRITE_ARRAY)
+        H5Dwrite(did, tid, H5S_ALL, H5S_ALL, H5P_DEFAULT, dbuf);
+#endif
+
+        /* close */
+        H5Dclose(did);
+        H5Tclose(tid);
+        H5Sclose(sid);
+        free( dbuf );
+    }
+
     /*-------------------------------------------------------------------------
     * H5T_INTEGER and H5T_FLOAT
     *-------------------------------------------------------------------------