summaryrefslogtreecommitdiffstats
path: root/test/big.c
diff options
context:
space:
mode:
authorRobb Matzke <matzke@llnl.gov>1998-04-14 16:44:46 (GMT)
committerRobb Matzke <matzke@llnl.gov>1998-04-14 16:44:46 (GMT)
commitb6fc6ede10fee24e9ef0a21df81cc181f6d7fec0 (patch)
tree66321165e4bd58a9fe3722a1a3baee2f835d0544 /test/big.c
parent3f571b6b4499e09aa6db29cb389e1121624a081d (diff)
downloadhdf5-b6fc6ede10fee24e9ef0a21df81cc181f6d7fec0.zip
hdf5-b6fc6ede10fee24e9ef0a21df81cc181f6d7fec0.tar.gz
hdf5-b6fc6ede10fee24e9ef0a21df81cc181f6d7fec0.tar.bz2
[svn-r346] Changes since 19980410
---------------------- ./html/H5.format.html ./src/H5E.c ./src/H5Epublic.h ./src/H5F.c ./src/H5G.c ./src/H5Gent.c ./src/H5Gnode.c ./src/H5Gprivate.h ./src/h5ls.c ./test/Makefile.in Symbolic links are now supported. The ./test/links program will create a `links.h5' file that demonstrates hard links, soft links, dangling links, and recursive links. A symbolic link is a symbol table entity and doesn't have an object header. It's format is described in H5.format.hml. ./src/H5G.c ./src/H5Gpublic.h ./src/h5ls.c Implemented H5Gstat() and H5Gget_linkval() as documented by Quincey. The `H5G_type_t type' field of `H5G_stat_t' was changed to `int type' because H5G_type_t was already used and the `type' data type should be open-ended. ./src/H5Ffamily.c Removed an incorrect diagnostic message. ./test/big.c Added read/write calls to test partial I/O to big contiguous datasets. With no arguments it writes to a dataset and prints the list of points written which should be redirected to a file. With an argument (the name of a file containing the stdout of a run with no arguments) values are read from the dataset. One would typically say `big >x && big x'.
Diffstat (limited to 'test/big.c')
-rw-r--r--test/big.c160
1 files changed, 150 insertions, 10 deletions
diff --git a/test/big.c b/test/big.c
index 44d550f..a75eb5d 100644
--- a/test/big.c
+++ b/test/big.c
@@ -6,7 +6,26 @@
* Wednesday, April 8, 1998
*/
#include <assert.h>
+#include <ctype.h>
#include <hdf5.h>
+#include <math.h>
+#include <stdlib.h>
+
+#include <H5private.h> /*needed for HDfprintf() */
+
+#define FNAME "big%05d.h5"
+#define WRT_N 50
+#define WRT_SIZE 4*1024
+
+static hsize_t
+randll (hsize_t limit)
+{
+ hsize_t acc = rand ();
+ acc *= rand ();
+
+ return acc % limit;
+}
+
/*-------------------------------------------------------------------------
@@ -25,12 +44,16 @@
*
*-------------------------------------------------------------------------
*/
-int
-main (void)
+static void
+writer (int wrt_n)
{
hsize_t size1[4] = {8, 1024, 1024, 1024};
hsize_t size2[1] = {8589934592LL};
- hid_t plist, file, space1, space2, dset;
+ hssize_t hs_start[1];
+ hsize_t hs_size[1];
+ hid_t plist, file, space1, space2, mem_space, d1, d2;
+ int *buf = malloc (sizeof(int) * WRT_SIZE);
+ int i, j;
/*
* Make sure that `hsize_t' is large enough to represent the entire data
@@ -44,20 +67,137 @@ main (void)
*/
plist = H5Pcreate (H5P_FILE_ACCESS);
H5Pset_family (plist, 30, H5P_DEFAULT);
- file = H5Fcreate ("big%05d.h5", H5F_ACC_TRUNC, H5P_DEFAULT, plist);
+ file = H5Fcreate (FNAME, H5F_ACC_TRUNC, H5P_DEFAULT, plist);
/* Create simple data spaces according to the size specified above. */
space1 = H5Screate_simple (4, size1, size1);
space2 = H5Screate_simple (1, size2, size2);
/* Create the datasets */
- dset = H5Dcreate (file, "d1", H5T_NATIVE_INT, space1, H5P_DEFAULT);
- H5Dclose (dset);
- dset = H5Dcreate (file, "d2", H5T_NATIVE_INT, space2, H5P_DEFAULT);
- H5Dclose (dset);
-
+ d1 = H5Dcreate (file, "d1", H5T_NATIVE_INT, space1, H5P_DEFAULT);
+ d2 = H5Dcreate (file, "d2", H5T_NATIVE_INT, space2, H5P_DEFAULT);
+
+ /* Write some things to them randomly */
+ hs_size[0] = WRT_SIZE;
+ mem_space = H5Screate_simple (1, hs_size, hs_size);
+ for (i=0; i<wrt_n; i++) {
+ hs_start[0] = randll (size2[0]);
+ HDfprintf (stdout, "#%03d 0x%016Hx\n", i, hs_start[0]);
+ H5Sset_hyperslab (space2, hs_start, hs_size, NULL);
+ for (j=0; j<WRT_SIZE; j++) {
+ buf[j] = i+1;
+ }
+ H5Dwrite (d2, H5T_NATIVE_INT, mem_space, space2, H5P_DEFAULT, buf);
+ }
+
+ H5Dclose (d1);
+ H5Dclose (d2);
+ H5Sclose (mem_space);
H5Sclose (space1);
H5Sclose (space2);
H5Fclose (file);
- exit (0);
+ free (buf);
+}
+
+
+/*-------------------------------------------------------------------------
+ * Function: reader
+ *
+ * Purpose: Reads some data from random locations in the dataset.
+ *
+ * Return: void
+ *
+ * Programmer: Robb Matzke
+ * Friday, April 10, 1998
+ *
+ * Modifications:
+ *
+ *-------------------------------------------------------------------------
+ */
+static void
+reader (const char *script_name)
+{
+ FILE *script;
+ hid_t plist, file, mspace, fspace, d2;
+ char ln[64], *s;
+ hssize_t hs_offset[1];
+ hsize_t hs_size[1] = {WRT_SIZE};
+ int *buf = malloc (sizeof(int) * WRT_SIZE);
+ int i, j, zero, wrong;
+
+ /* Open script file */
+ script = fopen (script_name, "r");
+
+ /* Open HDF5 file */
+ plist = H5Pcreate (H5P_FILE_ACCESS);
+ H5Pset_family (plist, 30, H5P_DEFAULT);
+ file = H5Fopen (FNAME, H5F_ACC_RDONLY, plist);
+
+ /* Open the dataset */
+ d2 = H5Dopen (file, "d2");
+ fspace = H5Dget_space (d2);
+
+ /* Describe `buf' */
+ mspace = H5Screate_simple (1, hs_size, hs_size);
+
+ /* Read each region */
+ while (fgets (ln, sizeof(ln), script)) {
+ if ('#'!=ln[0]) break;
+ i = strtol (ln+1, &s, 10);
+ hs_offset[0] = HDstrtoll (s, NULL, 0);
+ HDfprintf (stdout, "#%03d 0x%016Hx", i, hs_offset[0]);
+ fflush (stdout);
+
+ H5Sset_hyperslab (fspace, hs_offset, hs_size, NULL);
+ H5Dread (d2, H5T_NATIVE_INT, mspace, fspace, H5P_DEFAULT, buf);
+
+ /* Check */
+ for (j=zero=wrong=0; j<WRT_SIZE; j++) {
+ if (0==buf[j]) zero++;
+ else if (buf[j]!=i+1) wrong++;
+ }
+ if (zero) {
+ printf (" *FAILED* (%d zeros)\n", zero);
+ } else if (wrong) {
+ printf (" *SKIPPED* (possible overlap with another region)\n");
+ } else {
+ printf (" PASSED\n");
+ }
+ }
+
+ H5Dclose (d2);
+ H5Sclose (mspace);
+ H5Sclose (fspace);
+ H5Fclose (file);
+ fclose (script);
+}
+
+
+/*-------------------------------------------------------------------------
+ * Function: main
+ *
+ * Purpose:
+ *
+ * Return: Success:
+ *
+ * Failure:
+ *
+ * Programmer: Robb Matzke
+ * Friday, April 10, 1998
+ *
+ * Modifications:
+ *
+ *-------------------------------------------------------------------------
+ */
+int
+main (int argc, char *argv[])
+{
+ if (1==argc) {
+ writer (WRT_N);
+ } else if (isdigit (argv[1][0])) {
+ writer (strtol (argv[1], NULL, 0));
+ } else {
+ reader (argv[1]);
+ }
+ return 0;
}