diff options
author | Quincey Koziol <koziol@hdfgroup.org> | 1998-01-16 22:23:43 (GMT) |
---|---|---|
committer | Quincey Koziol <koziol@hdfgroup.org> | 1998-01-16 22:23:43 (GMT) |
commit | c2c94c31878dc42926661c9cb7e71be620196fc1 (patch) | |
tree | ee9ce2ae309882348dc1eafc1d105b5e6646605e /test/istore.c | |
parent | 903e677366a86ea385d5cfe1241f3f53132941de (diff) | |
download | hdf5-c2c94c31878dc42926661c9cb7e71be620196fc1.zip hdf5-c2c94c31878dc42926661c9cb7e71be620196fc1.tar.gz hdf5-c2c94c31878dc42926661c9cb7e71be620196fc1.tar.bz2 |
[svn-r157] Reformatted code with indent...
Diffstat (limited to 'test/istore.c')
-rw-r--r-- | test/istore.c | 1079 |
1 files changed, 538 insertions, 541 deletions
diff --git a/test/istore.c b/test/istore.c index 8f246c8..6115a5f 100644 --- a/test/istore.c +++ b/test/istore.c @@ -2,10 +2,10 @@ * Copyright (C) 1997 NCSA * All rights reserved. * - * Programmer: Robb Matzke <matzke@llnl.gov> - * Wednesday, October 15, 1997 + * Programmer: Robb Matzke <matzke@llnl.gov> + * Wednesday, October 15, 1997 * - * Purpose: Tests various aspects of indexed raw data storage. + * Purpose: Tests various aspects of indexed raw data storage. */ #include <H5private.h> #include <H5Cprivate.h> @@ -16,38 +16,38 @@ #include <H5Vprivate.h> #if 1 -# define FILETYPE H5F_LOW_DFLT -# define FILENAME "istore.h5" +# define FILETYPE H5F_LOW_DFLT +# define FILENAME "istore.h5" #elif 0 -# define FILETYPE H5F_LOW_FAM -# define FILENAME "istore-%05d.h5" -# define TEST_FAMILY 1 +# define FILETYPE H5F_LOW_FAM +# define FILENAME "istore-%05d.h5" +# define TEST_FAMILY 1 #else -# define FILETYPE H5F_LOW_SPLIT -# define FILENAME "istore-split" +# define FILETYPE H5F_LOW_SPLIT +# define FILENAME "istore-split" #endif -#define TEST_SMALL 0x0001 -#define TEST_MEDIUM 0x0002 -#define TEST_LARGE 0x0004 +#define TEST_SMALL 0x0001 +#define TEST_MEDIUM 0x0002 +#define TEST_LARGE 0x0004 #ifndef HAVE_FUNCTION #define __FUNCTION__ "" #endif -#define AT() printf (" at %s:%d in %s()...\n", \ - __FILE__, __LINE__, __FUNCTION__); - -size_t align_g[3] = {50, 50, 50}; +#define AT() printf (" at %s:%d in %s()...\n", \ + __FILE__, __LINE__, __FUNCTION__); +size_t align_g[3] = +{50, 50, 50}; /*------------------------------------------------------------------------- - * Function: print_array + * Function: print_array * - * Purpose: Prints the values in an array + * Purpose: Prints the values in an array * - * Return: void + * Return: void * - * Programmer: Robb Matzke + * Programmer: Robb Matzke * Friday, October 10, 1997 * * Modifications: @@ -55,41 +55,41 @@ size_t align_g[3] = {50, 50, 50}; *------------------------------------------------------------------------- */ static void -print_array (uint8 *array, size_t nx, size_t ny, size_t nz) +print_array(uint8 *array, size_t nx, size_t ny, size_t nz) { - int i, j, k; - - for (i=0; i<nx; i++) { - if (nz>1) { - printf ("i=%d:\n", i); - } else { - printf ("%03d:", i); - } - - for (j=0; j<ny; j++) { - if (nz>1) printf ("%03d:", j); - for (k=0; k<nz; k++) { - printf (" %3d", *array++); - } - if (nz>1) printf ("\n"); - } - printf ("\n"); - } + int i, j, k; + + for (i = 0; i < nx; i++) { + if (nz > 1) { + printf("i=%d:\n", i); + } else { + printf("%03d:", i); + } + + for (j = 0; j < ny; j++) { + if (nz > 1) + printf("%03d:", j); + for (k = 0; k < nz; k++) { + printf(" %3d", *array++); + } + if (nz > 1) + printf("\n"); + } + printf("\n"); + } } - - /*------------------------------------------------------------------------- - * Function: new_object + * Function: new_object * - * Purpose: Creates a new object that refers to a indexed storage of raw - * data. No raw data is stored. + * Purpose: Creates a new object that refers to a indexed storage of raw + * data. No raw data is stored. * - * Return: Success: Handle to a new open object. + * Return: Success: Handle to a new open object. * - * Failure: NULL, error message printed. + * Failure: NULL, error message printed. * - * Programmer: Robb Matzke + * Programmer: Robb Matzke * Wednesday, October 15, 1997 * * Modifications: @@ -97,69 +97,65 @@ print_array (uint8 *array, size_t nx, size_t ny, size_t nz) *------------------------------------------------------------------------- */ static int -new_object (H5F_t *f, const char *name, size_t ndims, H5G_entry_t *ent/*out*/) +new_object(H5F_t *f, const char *name, size_t ndims, H5G_entry_t *ent /*out */ ) { - H5O_layout_t layout; - intn i; - - /* Create the object header */ - if (H5O_create (f, 64, ent)) { - puts ("*FAILED*"); - if (!isatty (1)) { - AT (); - printf (" H5O_create() = NULL\n"); - } - return -1; - } - - /* create chunked storage */ - layout.type = H5D_CHUNKED; - layout.ndims = ndims; - for (i=0; i<ndims; i++) { - if (i<NELMTS (align_g)) { - layout.dim[i] = align_g[i]; - } else { - layout.dim[i] = 2; - } - } - H5F_arr_create (f, &layout/*in,out*/); - if (H5O_modify (ent, H5O_LAYOUT, H5O_NEW_MESG, 0, &layout)<0) { - printf ("*FAILED*\n"); - if (!isatty (1)) { - AT(); - printf (" H5O_modify istore message failure\n"); - } - return -1; - } - - /* Give the object header a name */ - if (H5G_insert (name, ent)<0) { - printf ("*FAILED*\n"); - if (!isatty (1)) { - AT (); - printf (" H5G_insert(f, name=\"%s\", ent) failed\n", name); - } - return -1; - } - - /* Close the header */ - H5O_close (ent); - - return 0; + H5O_layout_t layout; + intn i; + + /* Create the object header */ + if (H5O_create(f, 64, ent)) { + puts("*FAILED*"); + if (!isatty(1)) { + AT(); + printf(" H5O_create() = NULL\n"); + } + return -1; + } + /* create chunked storage */ + layout.type = H5D_CHUNKED; + layout.ndims = ndims; + for (i = 0; i < ndims; i++) { + if (i < NELMTS(align_g)) { + layout.dim[i] = align_g[i]; + } else { + layout.dim[i] = 2; + } + } + H5F_arr_create(f, &layout /*in,out */ ); + if (H5O_modify(ent, H5O_LAYOUT, H5O_NEW_MESG, 0, &layout) < 0) { + printf("*FAILED*\n"); + if (!isatty(1)) { + AT(); + printf(" H5O_modify istore message failure\n"); + } + return -1; + } + /* Give the object header a name */ + if (H5G_insert(name, ent) < 0) { + printf("*FAILED*\n"); + if (!isatty(1)) { + AT(); + printf(" H5G_insert(f, name=\"%s\", ent) failed\n", name); + } + return -1; + } + /* Close the header */ + H5O_close(ent); + + return 0; } - /*------------------------------------------------------------------------- - * Function: test_create + * Function: test_create * - * Purpose: Creates a named object that refers to indexed storage of raw - * data. No raw data is stored. + * Purpose: Creates a named object that refers to indexed storage of raw + * data. No raw data is stored. * - * Return: Success: SUCCEED + * Return: Success: SUCCEED * - * Failure: FAIL + * Failure: FAIL * - * Programmer: Robb Matzke + * Programmer: Robb Matzke * Wednesday, October 15, 1997 * * Modifications: @@ -167,37 +163,36 @@ new_object (H5F_t *f, const char *name, size_t ndims, H5G_entry_t *ent/*out*/) *------------------------------------------------------------------------- */ static herr_t -test_create (H5F_t *f, const char *prefix) +test_create(H5F_t *f, const char *prefix) { - H5G_entry_t handle; - intn i; - char name[256]; - - printf ("%-70s", "Testing istore create"); - - for (i=1; i<=H5O_LAYOUT_NDIMS; i++) { - sprintf (name, "%s_%02d", prefix, i); - if (new_object (f, name, i, &handle)<0) return FAIL; - } - - puts (" PASSED"); - return SUCCEED; -} + H5G_entry_t handle; + intn i; + char name[256]; + printf("%-70s", "Testing istore create"); + for (i = 1; i <= H5O_LAYOUT_NDIMS; i++) { + sprintf(name, "%s_%02d", prefix, i); + if (new_object(f, name, i, &handle) < 0) + return FAIL; + } + + puts(" PASSED"); + return SUCCEED; +} /*------------------------------------------------------------------------- - * Function: test_extend + * Function: test_extend * - * Purpose: Creates an empty object and then writes to it in such a way - * as to always extend the object's domain without creating - * holes and without causing the object to become concave. + * Purpose: Creates an empty object and then writes to it in such a way + * as to always extend the object's domain without creating + * holes and without causing the object to become concave. * - * Return: Success: SUCCEED + * Return: Success: SUCCEED * - * Failure: FAIL + * Failure: FAIL * - * Programmer: Robb Matzke + * Programmer: Robb Matzke * Wednesday, October 15, 1997 * * Modifications: @@ -205,214 +200,217 @@ test_create (H5F_t *f, const char *prefix) *------------------------------------------------------------------------- */ static herr_t -test_extend (H5F_t *f, const char *prefix, - size_t nx, size_t ny, size_t nz) +test_extend(H5F_t *f, const char *prefix, + size_t nx, size_t ny, size_t nz) { - H5G_entry_t handle; - int i, j, k, ndims, ctr; - uint8 *buf=NULL, *check=NULL, *whole=NULL; - char dims[64], s[256], name[256]; - size_t offset[3]; - size_t max_corner[3]; - size_t size[3]; - size_t whole_size[3]; - size_t nelmts; - H5O_layout_t layout; - - if (!nz) { - if (!ny) { - ndims = 1; - ny = nz = 1; - sprintf (dims, "%lu", (unsigned long)nx); - } else { - ndims = 2; - nz = 1; - sprintf (dims, "%lux%lu", (unsigned long)nx, (unsigned long)ny); - } - } else { - ndims = 3; - sprintf (dims, "%lux%lux%lu", - (unsigned long)nx, (unsigned long)ny, (unsigned long)nz); - } - - - sprintf (s, "Testing istore extend: %s", dims); - printf ("%-70s", s); - buf = H5MM_xmalloc (nx*ny*nz); - check = H5MM_xmalloc (nx*ny*nz); - whole = H5MM_xcalloc (nx*ny*nz, 1); - - /* Build the new empty object */ - sprintf (name, "%s_%s", prefix, dims); - if (new_object (f, name, ndims, &handle)<0) { - if (!isatty (1)) { - AT (); - printf (" Cannot create %d-d object `%s'\n", ndims, name); - } - goto error; - } - if (NULL==H5O_read (&handle, H5O_LAYOUT, 0, &layout)) { - puts ("*FAILED*"); - if (!isatty (1)) { - AT (); - printf (" Unable to read istore message\n"); - } - goto error; - } - if (ndims!=layout.ndims) { - puts ("*FAILED*"); - if (!isatty (1)) { - AT (); - printf (" Header read error: istore.ndims != %d\n", ndims); - } - goto error; - } - - whole_size[0] = nx; - whole_size[1] = ny; - whole_size[2] = nz; - max_corner[0] = 0; - max_corner[1] = 0; - max_corner[2] = 0; - - for (ctr=0; H5V_vector_lt (ndims, max_corner, whole_size); ctr++) { - - /* Size and location */ - if (0==ctr) { - offset[0] = offset[1] = offset[2] = 0; - size[0] = size[1] = size[2] = 1; - nelmts = 1; - } else { - for (i=0, nelmts=1; i<ndims; i++) { - if (ctr % ndims == i) { - offset[i] = max_corner[i]; - size[i] = MIN (1, whole_size[i]-offset[i]); - } else { - offset[i] = 0; - size[i] = max_corner[i]; - } - nelmts *= size[i]; - } - } + H5G_entry_t handle; + int i, j, k, ndims, ctr; + uint8 *buf = NULL, *check = NULL, *whole = NULL; + char dims[64], s[256], name[256]; + size_t offset[3]; + size_t max_corner[3]; + size_t size[3]; + size_t whole_size[3]; + size_t nelmts; + H5O_layout_t layout; + + if (!nz) { + if (!ny) { + ndims = 1; + ny = nz = 1; + sprintf(dims, "%lu", (unsigned long) nx); + } else { + ndims = 2; + nz = 1; + sprintf(dims, "%lux%lu", (unsigned long) nx, (unsigned long) ny); + } + } else { + ndims = 3; + sprintf(dims, "%lux%lux%lu", + (unsigned long) nx, (unsigned long) ny, (unsigned long) nz); + } + + sprintf(s, "Testing istore extend: %s", dims); + printf("%-70s", s); + buf = H5MM_xmalloc(nx * ny * nz); + check = H5MM_xmalloc(nx * ny * nz); + whole = H5MM_xcalloc(nx * ny * nz, 1); + + /* Build the new empty object */ + sprintf(name, "%s_%s", prefix, dims); + if (new_object(f, name, ndims, &handle) < 0) { + if (!isatty(1)) { + AT(); + printf(" Cannot create %d-d object `%s'\n", ndims, name); + } + goto error; + } + if (NULL == H5O_read(&handle, H5O_LAYOUT, 0, &layout)) { + puts("*FAILED*"); + if (!isatty(1)) { + AT(); + printf(" Unable to read istore message\n"); + } + goto error; + } + if (ndims != layout.ndims) { + puts("*FAILED*"); + if (!isatty(1)) { + AT(); + printf(" Header read error: istore.ndims != %d\n", ndims); + } + goto error; + } + whole_size[0] = nx; + whole_size[1] = ny; + whole_size[2] = nz; + max_corner[0] = 0; + max_corner[1] = 0; + max_corner[2] = 0; + + for (ctr = 0; H5V_vector_lt(ndims, max_corner, whole_size); ctr++) { + + /* Size and location */ + if (0 == ctr) { + offset[0] = offset[1] = offset[2] = 0; + size[0] = size[1] = size[2] = 1; + nelmts = 1; + } else { + for (i = 0, nelmts = 1; i < ndims; i++) { + if (ctr % ndims == i) { + offset[i] = max_corner[i]; + size[i] = MIN(1, whole_size[i] - offset[i]); + } else { + offset[i] = 0; + size[i] = max_corner[i]; + } + nelmts *= size[i]; + } + } #if 0 - if (0==ctr) printf ("\n"); - printf (" Insert: ctr=%d, corner=(%d", ctr, offset[0]); - if (ndims>1) printf (",%d", offset[1]); - if (ndims>2) printf (",%d", offset[2]); - printf ("), size=(%d", size[0]); - if (ndims>1) printf (",%d", size[1]); - if (ndims>2) printf (",%d", size[2]); - printf ("), %d element%s", nelmts, 1==nelmts?"":"s"); - if (0==nelmts) printf (" *SKIPPED*"); - printf ("\n"); + if (0 == ctr) + printf("\n"); + printf(" Insert: ctr=%d, corner=(%d", ctr, offset[0]); + if (ndims > 1) + printf(",%d", offset[1]); + if (ndims > 2) + printf(",%d", offset[2]); + printf("), size=(%d", size[0]); + if (ndims > 1) + printf(",%d", size[1]); + if (ndims > 2) + printf(",%d", size[2]); + printf("), %d element%s", nelmts, 1 == nelmts ? "" : "s"); + if (0 == nelmts) + printf(" *SKIPPED*"); + printf("\n"); #endif - - /* Fill the source array */ - if (0==nelmts) continue; - memset (buf, 128+ctr, nelmts); - - /* Write to disk */ - if (H5F_arr_write (f, &layout, offset, size, buf)<0) { - puts ("*FAILED*"); - if (!isatty (1)) { - AT (); - printf (" Write failed: ctr=%d\n", ctr); - } - goto error; - } - - /* Read from disk */ - memset (check, 0xff, nelmts); - if (H5F_arr_read (f, &layout, offset, size, check)<0) { - puts ("*FAILED*"); - if (!isatty (1)) { - AT (); - printf (" Read failed: ctr=%d\n", ctr); - } - goto error; - } - if (memcmp (buf, check, nelmts)) { - puts ("*FAILED*"); - if (!isatty (1)) { - AT (); - printf (" Read check failed: ctr=%d\n", ctr); - printf (" Wrote:\n"); - print_array (buf, size[0], size[1], size[2]); - printf (" Read:\n"); - print_array (check, size[0], size[1], size[2]); - } - goto error; - } - - /* Write to `whole' buffer for later checking */ - H5V_hyper_copy (ndims, size, - whole_size, offset, whole, /*dst*/ - size, H5V_ZERO, buf); /*src*/ - - /* Update max corner */ - for (i=0; i<ndims; i++) { - max_corner[i] = MAX (max_corner[i], offset[i]+size[i]); - } - } - - /* Now read the entire array back out and check it */ - memset (buf, 0xff, nx*ny*nz); - if (H5F_arr_read (f, &layout, H5V_ZERO, whole_size, buf)<0) { - puts ("*FAILED*"); - if (!isatty (1)) { - AT (); - printf (" Read failed for whole array\n"); - } - goto error; - } - for (i=0; i<nx; i++) { - for (j=0; j<ny; j++) { - for (k=0; k<nz; k++) { - if (whole[i*ny*nz + j*nz + k] != buf[i*ny*nz + j*nz + k]) { - puts ("*FAILED*"); - if (!isatty (1)) { - AT (); - printf (" Check failed at i=%d", i); - if (ndims>1) printf (", j=%d", j); - if (ndims>2) printf (", k=%d", k); - printf ("\n Check array is:\n"); - print_array (whole, nx, ny, nz); - printf (" Value read is:\n"); - print_array (buf, nx, ny, nz); - } - goto error; - } - } - } - } - - puts (" PASSED"); - H5MM_xfree (buf); - H5MM_xfree (check); - H5MM_xfree (whole); - return SUCCEED; - - error: - H5MM_xfree (buf); - H5MM_xfree (check); - H5MM_xfree (whole); - return FAIL; -} - + /* Fill the source array */ + if (0 == nelmts) + continue; + memset(buf, 128 + ctr, nelmts); + + /* Write to disk */ + if (H5F_arr_write(f, &layout, offset, size, buf) < 0) { + puts("*FAILED*"); + if (!isatty(1)) { + AT(); + printf(" Write failed: ctr=%d\n", ctr); + } + goto error; + } + /* Read from disk */ + memset(check, 0xff, nelmts); + if (H5F_arr_read(f, &layout, offset, size, check) < 0) { + puts("*FAILED*"); + if (!isatty(1)) { + AT(); + printf(" Read failed: ctr=%d\n", ctr); + } + goto error; + } + if (memcmp(buf, check, nelmts)) { + puts("*FAILED*"); + if (!isatty(1)) { + AT(); + printf(" Read check failed: ctr=%d\n", ctr); + printf(" Wrote:\n"); + print_array(buf, size[0], size[1], size[2]); + printf(" Read:\n"); + print_array(check, size[0], size[1], size[2]); + } + goto error; + } + /* Write to `whole' buffer for later checking */ + H5V_hyper_copy(ndims, size, + whole_size, offset, whole, /*dst */ + size, H5V_ZERO, buf); /*src */ + + /* Update max corner */ + for (i = 0; i < ndims; i++) { + max_corner[i] = MAX(max_corner[i], offset[i] + size[i]); + } + } + + /* Now read the entire array back out and check it */ + memset(buf, 0xff, nx * ny * nz); + if (H5F_arr_read(f, &layout, H5V_ZERO, whole_size, buf) < 0) { + puts("*FAILED*"); + if (!isatty(1)) { + AT(); + printf(" Read failed for whole array\n"); + } + goto error; + } + for (i = 0; i < nx; i++) { + for (j = 0; j < ny; j++) { + for (k = 0; k < nz; k++) { + if (whole[i * ny * nz + j * nz + k] != buf[i * ny * nz + j * nz + k]) { + puts("*FAILED*"); + if (!isatty(1)) { + AT(); + printf(" Check failed at i=%d", i); + if (ndims > 1) + printf(", j=%d", j); + if (ndims > 2) + printf(", k=%d", k); + printf("\n Check array is:\n"); + print_array(whole, nx, ny, nz); + printf(" Value read is:\n"); + print_array(buf, nx, ny, nz); + } + goto error; + } + } + } + } + + puts(" PASSED"); + H5MM_xfree(buf); + H5MM_xfree(check); + H5MM_xfree(whole); + return SUCCEED; + + error: + H5MM_xfree(buf); + H5MM_xfree(check); + H5MM_xfree(whole); + return FAIL; +} /*------------------------------------------------------------------------- - * Function: test_sparse + * Function: test_sparse * - * Purpose: Creates a sparse matrix consisting of NBLOCKS randomly placed - * blocks each of size NX,NY,NZ. + * Purpose: Creates a sparse matrix consisting of NBLOCKS randomly placed + * blocks each of size NX,NY,NZ. * - * Return: Success: SUCCEED + * Return: Success: SUCCEED * - * Failure: FAIL + * Failure: FAIL * - * Programmer: Robb Matzke + * Programmer: Robb Matzke * Wednesday, October 22, 1997 * * Modifications: @@ -420,108 +418,109 @@ test_extend (H5F_t *f, const char *prefix, *------------------------------------------------------------------------- */ static herr_t -test_sparse (H5F_t *f, const char *prefix, size_t nblocks, - size_t nx, size_t ny, size_t nz) +test_sparse(H5F_t *f, const char *prefix, size_t nblocks, + size_t nx, size_t ny, size_t nz) { - intn ndims, ctr; - char dims[64], s[256], name[256]; - size_t offset[3], size[3], total=0; - H5G_entry_t handle; - H5O_layout_t layout; - uint8 *buf = NULL; - - if (!nz) { - if (!ny) { - ndims = 1; - ny = nz = 1; - sprintf (dims, "%lu", (unsigned long)nx); - } else { - ndims = 2; - nz = 1; - sprintf (dims, "%lux%lu", (unsigned long)nx, (unsigned long)ny); - } - } else { - ndims = 3; - sprintf (dims, "%lux%lux%lu", - (unsigned long)nx, (unsigned long)ny, (unsigned long)nz); - } - - sprintf (s, "Testing istore sparse: %s", dims); - printf ("%-70s", s); - buf = H5MM_xmalloc (nx*ny*nz); - - /* Build the new empty object */ - sprintf (name, "%s_%s", prefix, dims); - if (new_object (f, name, ndims, &handle)<0) { - if (!isatty (1)) { - AT (); - printf (" Cannot create %d-d object `%s'\n", ndims, name); - } - goto error; - } - if (NULL==H5O_read (&handle, H5O_LAYOUT, 0, &layout)) { - puts ("*FAILED*"); - if (!isatty (1)) { - AT (); - printf (" Unable to read istore message\n"); - } - goto error; - } - - for (ctr=0; ctr<nblocks; ctr++) { - offset[0] = rand () % 1000000; - offset[1] = rand () % 1000000; - offset[2] = rand () % 1000000; - size[0] = nx; - size[1] = ny; - size[2] = nz; - memset (buf, 128+ctr, nx*ny*nz); - - /* write to disk */ - if (H5F_arr_write (f, &layout, offset, size, buf)<0) { - puts ("*FAILED*"); - if (!isatty (1)) { - AT (); - printf (" Write failed: ctr=%d\n", ctr); - printf (" offset=(%lu", (unsigned long)(offset[0])); - if (ndims>1) printf (",%lu", (unsigned long)(offset[1])); - if (ndims>2) printf (",%lu", (unsigned long)(offset[2])); - printf ("), size=(%lu", (unsigned long)(size[0])); - if (ndims>1) printf (",%lu", (unsigned long)(size[1])); - if (ndims>2) printf (",%lu", (unsigned long)(size[2])); - printf (")\n"); - } - goto error; - } - total += nx*ny*nz; + intn ndims, ctr; + char dims[64], s[256], name[256]; + size_t offset[3], size[3], total = 0; + H5G_entry_t handle; + H5O_layout_t layout; + uint8 *buf = NULL; + + if (!nz) { + if (!ny) { + ndims = 1; + ny = nz = 1; + sprintf(dims, "%lu", (unsigned long) nx); + } else { + ndims = 2; + nz = 1; + sprintf(dims, "%lux%lu", (unsigned long) nx, (unsigned long) ny); + } + } else { + ndims = 3; + sprintf(dims, "%lux%lux%lu", + (unsigned long) nx, (unsigned long) ny, (unsigned long) nz); + } + + sprintf(s, "Testing istore sparse: %s", dims); + printf("%-70s", s); + buf = H5MM_xmalloc(nx * ny * nz); + + /* Build the new empty object */ + sprintf(name, "%s_%s", prefix, dims); + if (new_object(f, name, ndims, &handle) < 0) { + if (!isatty(1)) { + AT(); + printf(" Cannot create %d-d object `%s'\n", ndims, name); + } + goto error; + } + if (NULL == H5O_read(&handle, H5O_LAYOUT, 0, &layout)) { + puts("*FAILED*"); + if (!isatty(1)) { + AT(); + printf(" Unable to read istore message\n"); + } + goto error; + } + for (ctr = 0; ctr < nblocks; ctr++) { + offset[0] = rand() % 1000000; + offset[1] = rand() % 1000000; + offset[2] = rand() % 1000000; + size[0] = nx; + size[1] = ny; + size[2] = nz; + memset(buf, 128 + ctr, nx * ny * nz); + + /* write to disk */ + if (H5F_arr_write(f, &layout, offset, size, buf) < 0) { + puts("*FAILED*"); + if (!isatty(1)) { + AT(); + printf(" Write failed: ctr=%d\n", ctr); + printf(" offset=(%lu", (unsigned long) (offset[0])); + if (ndims > 1) + printf(",%lu", (unsigned long) (offset[1])); + if (ndims > 2) + printf(",%lu", (unsigned long) (offset[2])); + printf("), size=(%lu", (unsigned long) (size[0])); + if (ndims > 1) + printf(",%lu", (unsigned long) (size[1])); + if (ndims > 2) + printf(",%lu", (unsigned long) (size[2])); + printf(")\n"); + } + goto error; + } + total += nx * ny * nz; #if 0 - printf ("ctr: ctr=%d, total=%lu\n", ctr, (unsigned long)total); + printf("ctr: ctr=%d, total=%lu\n", ctr, (unsigned long) total); #endif - /* We don't test reading yet.... */ - } - - - puts (" PASSED"); - H5MM_xfree (buf); - return SUCCEED; - - error: - H5MM_xfree (buf); - return FAIL; -} + /* We don't test reading yet.... */ + } + puts(" PASSED"); + H5MM_xfree(buf); + return SUCCEED; + + error: + H5MM_xfree(buf); + return FAIL; +} /*------------------------------------------------------------------------- - * Function: main + * Function: main * - * Purpose: Tests indexed storage stuff. + * Purpose: Tests indexed storage stuff. * - * Return: Success: exit(0) + * Return: Success: exit(0) * - * Failure: exit(non-zero) + * Failure: exit(non-zero) * - * Programmer: Robb Matzke + * Programmer: Robb Matzke * Wednesday, October 15, 1997 * * Modifications: @@ -529,136 +528,134 @@ test_sparse (H5F_t *f, const char *prefix, size_t nblocks, *------------------------------------------------------------------------- */ int -main (int argc, char *argv[]) +main(int argc, char *argv[]) { - H5F_t *f; - herr_t status; - int nerrors = 0; - uintn size_of_test; - hid_t template_id; - H5F_create_t *creation_template = NULL; - H5G_t *dir = NULL; - - setbuf (stdout, NULL); - - /* Parse arguments or assume `small' */ - if (1==argc) { - size_of_test = TEST_SMALL; - } else { - intn i; - for (i=1,size_of_test=0; i<argc; i++) { - if (!strcmp (argv[i], "small")) { - size_of_test |= TEST_SMALL; - } else if (!strcmp (argv[i], "medium")) { - size_of_test |= TEST_MEDIUM; - } else if (!strcmp (argv[i], "large")) { - size_of_test |= TEST_LARGE; - } else { - printf ("unrecognized argument: %s\n", argv[i]); - exit (1); - } - } - } - printf ("Test sizes: "); - if (size_of_test & TEST_SMALL) printf (" SMALL"); - if (size_of_test & TEST_MEDIUM) printf (" MEDIUM"); - if (size_of_test & TEST_LARGE) printf (" LARGE"); - printf ("\n"); - - /* - * Use larger file addresses... - */ - template_id = H5Ccreate (H5C_FILE_CREATE); - H5Cset_sizes (template_id, 8, 0); - creation_template = H5Aatom_object (template_id); - - /* Create the test file */ - if (NULL==(f=H5F_open (FILETYPE, FILENAME, - (H5F_ACC_CREAT|H5F_ACC_WRITE|H5F_ACC_TRUNC| - H5F_ACC_DEBUG), - creation_template))) { - printf ("Cannot create file %s; test aborted\n", FILENAME); - exit (1); - } - + H5F_t *f; + herr_t status; + int nerrors = 0; + uintn size_of_test; + hid_t template_id; + H5F_create_t *creation_template = NULL; + H5G_t *dir = NULL; + + setbuf(stdout, NULL); + + /* Parse arguments or assume `small' */ + if (1 == argc) { + size_of_test = TEST_SMALL; + } else { + intn i; + for (i = 1, size_of_test = 0; i < argc; i++) { + if (!strcmp(argv[i], "small")) { + size_of_test |= TEST_SMALL; + } else if (!strcmp(argv[i], "medium")) { + size_of_test |= TEST_MEDIUM; + } else if (!strcmp(argv[i], "large")) { + size_of_test |= TEST_LARGE; + } else { + printf("unrecognized argument: %s\n", argv[i]); + exit(1); + } + } + } + printf("Test sizes: "); + if (size_of_test & TEST_SMALL) + printf(" SMALL"); + if (size_of_test & TEST_MEDIUM) + printf(" MEDIUM"); + if (size_of_test & TEST_LARGE) + printf(" LARGE"); + printf("\n"); + + /* + * Use larger file addresses... + */ + template_id = H5Ccreate(H5C_FILE_CREATE); + H5Cset_sizes(template_id, 8, 0); + creation_template = H5Aatom_object(template_id); + + /* Create the test file */ + if (NULL == (f = H5F_open(FILETYPE, FILENAME, + (H5F_ACC_CREAT | H5F_ACC_WRITE | H5F_ACC_TRUNC | + H5F_ACC_DEBUG), + creation_template))) { + printf("Cannot create file %s; test aborted\n", FILENAME); + exit(1); + } #ifdef TEST_FAMILY - { - /* - * For testing file families, fool the library into thinking it already - * allocated a whole bunch of data. - */ - haddr_t addr; - addr.offset = 8 * ((uint64)1<<30); /*8 GB*/ - H5F_low_seteof (f->shared->lf, &addr); - } + { + /* + * For testing file families, fool the library into thinking it already + * allocated a whole bunch of data. + */ + haddr_t addr; + addr.offset = 8 * ((uint64) 1 << 30); /*8 GB */ + H5F_low_seteof(f->shared->lf, &addr); + } #endif - /* - * By creating a group we cause the library to emit it's debugging - * diagnostic messages before we begin testing... - */ - dir = H5G_create (f, "flushing_diagnostics", 0); - H5G_close (dir); - dir = NULL; - - - /* - * Creation test: Creates empty objects with various raw data sizes - * and alignments. - */ - status = test_create (f, "create"); - nerrors += status<0 ? 1 : 0; - - if (size_of_test & TEST_SMALL) { - status = test_extend (f, "extend", 10, 0, 0); - nerrors += status<0 ? 1 : 0; - status = test_extend (f, "extend", 10, 10, 0); - nerrors += status<0 ? 1 : 0; - status = test_extend (f, "extend", 10, 10, 10); - nerrors += status<0 ? 1 : 0; - } - if (size_of_test & TEST_MEDIUM) { - status = test_extend (f, "extend", 10000, 0, 0); - nerrors += status<0 ? 1 : 0; - status = test_extend (f, "extend", 2500, 10, 0); - nerrors += status<0 ? 1 : 0; - status = test_extend (f, "extend", 10, 400, 10); - nerrors += status<0 ? 1 : 0; - } - if (size_of_test & TEST_SMALL) { - status = test_sparse (f, "sparse", 100, 5, 0, 0); - nerrors += status<0 ? 1 : 0; - status = test_sparse (f, "sparse", 100, 3, 4, 0); - nerrors += status<0 ? 1 : 0; - status = test_sparse (f, "sparse", 100, 2, 3, 4); - nerrors += status<0 ? 1 : 0; - } - if (size_of_test & TEST_MEDIUM) { - status = test_sparse (f, "sparse", 1000, 30, 0, 0); - nerrors += status<0 ? 1 : 0; - status = test_sparse (f, "sparse", 2000, 7, 3, 0); - nerrors += status<0 ? 1 : 0; - status = test_sparse (f, "sparse", 2000, 4, 2, 3); - nerrors += status<0 ? 1 : 0; - } - if (size_of_test & TEST_LARGE) { - status = test_sparse (f, "sparse", 800, 50, 50, 50); - nerrors += status<0 ? 1 : 0; - } - - - /* Close the test file and exit */ - H5F_close (f); - if (nerrors) { - printf ("***** %d I-STORE TEST%s FAILED! *****\n", - nerrors, 1==nerrors?"":"S"); - if (isatty (1)) { - printf ("(Redirect output to a pager or a file to see " - "debug output)\n"); - } - exit (1); - } - - printf ("All i-store tests passed.\n"); - exit (0); + /* + * By creating a group we cause the library to emit it's debugging + * diagnostic messages before we begin testing... + */ + dir = H5G_create(f, "flushing_diagnostics", 0); + H5G_close(dir); + dir = NULL; + + /* + * Creation test: Creates empty objects with various raw data sizes + * and alignments. + */ + status = test_create(f, "create"); + nerrors += status < 0 ? 1 : 0; + + if (size_of_test & TEST_SMALL) { + status = test_extend(f, "extend", 10, 0, 0); + nerrors += status < 0 ? 1 : 0; + status = test_extend(f, "extend", 10, 10, 0); + nerrors += status < 0 ? 1 : 0; + status = test_extend(f, "extend", 10, 10, 10); + nerrors += status < 0 ? 1 : 0; + } + if (size_of_test & TEST_MEDIUM) { + status = test_extend(f, "extend", 10000, 0, 0); + nerrors += status < 0 ? 1 : 0; + status = test_extend(f, "extend", 2500, 10, 0); + nerrors += status < 0 ? 1 : 0; + status = test_extend(f, "extend", 10, 400, 10); + nerrors += status < 0 ? 1 : 0; + } + if (size_of_test & TEST_SMALL) { + status = test_sparse(f, "sparse", 100, 5, 0, 0); + nerrors += status < 0 ? 1 : 0; + status = test_sparse(f, "sparse", 100, 3, 4, 0); + nerrors += status < 0 ? 1 : 0; + status = test_sparse(f, "sparse", 100, 2, 3, 4); + nerrors += status < 0 ? 1 : 0; + } + if (size_of_test & TEST_MEDIUM) { + status = test_sparse(f, "sparse", 1000, 30, 0, 0); + nerrors += status < 0 ? 1 : 0; + status = test_sparse(f, "sparse", 2000, 7, 3, 0); + nerrors += status < 0 ? 1 : 0; + status = test_sparse(f, "sparse", 2000, 4, 2, 3); + nerrors += status < 0 ? 1 : 0; + } + if (size_of_test & TEST_LARGE) { + status = test_sparse(f, "sparse", 800, 50, 50, 50); + nerrors += status < 0 ? 1 : 0; + } + /* Close the test file and exit */ + H5F_close(f); + if (nerrors) { + printf("***** %d I-STORE TEST%s FAILED! *****\n", + nerrors, 1 == nerrors ? "" : "S"); + if (isatty(1)) { + printf("(Redirect output to a pager or a file to see " + "debug output)\n"); + } + exit(1); + } + printf("All i-store tests passed.\n"); + exit(0); } |