diff options
| -rw-r--r-- | test/Makefile.in | 27 | ||||
| -rw-r--r-- | test/srb_append.c | 133 | ||||
| -rw-r--r-- | test/srb_read.c | 164 | ||||
| -rw-r--r-- | test/srb_write.c | 131 |
4 files changed, 447 insertions, 8 deletions
diff --git a/test/Makefile.in b/test/Makefile.in index 0039159..5d11b48 100644 --- a/test/Makefile.in +++ b/test/Makefile.in @@ -16,9 +16,10 @@ CPPFLAGS=-I. -I$(srcdir) -I../src -I$(top_srcdir)/src @CPPFLAGS@ ## These are our main targets. They should be listed in the order to be ## executed, generally most specific tests to least specific tests. RUNTEST=$(LT_RUN) -TEST_PROGS=testhdf5 lheap ohdr stab gheap hyperslab istore bittests dtypes \ - dsets cmpd_dset extend external links unlink big mtime fillval mount \ - flush1 flush2 enum gass_write gass_read gass_append dpss_read dpss_write +TEST_PROGS=testhdf5 lheap ohdr stab gheap hyperslab istore bittests dtypes \ + dsets cmpd_dset extend external links unlink big mtime fillval mount \ + flush1 flush2 enum gass_write gass_read gass_append dpss_write \ + dpss_read srb_write srb_append srb_read TIMINGS=iopipe chunk ragged overhead ## The libh5test.a library provides common support code for the tests. We link @@ -53,11 +54,12 @@ CLEAN=$(TIMINGS) ## other source lists are for the individual tests, the files of which may ## overlap with other tests. TEST_SRC=big.c bittests.c chunk.c cmpd_dset.c dsets.c dtypes.c extend.c \ - external.c fillval.c flush1.c flush2.c gheap.c h5test.c hyperslab.c \ - iopipe.c istore.c lheap.c links.c mount.c mtime.c ohdr.c overhead.c \ - ragged.c stab.c tattr.c testhdf5.c tfile.c th5s.c titerate.c tmeta.c \ - trefer.c tselect.c tvltypes.c tvlstr.c unlink.c enum.c gass_write.c \ - gass_read.c gass_append.c dpss_read.c dpss_write.c + external.c fillval.c flush1.c flush2.c gheap.c h5test.c hyperslab.c \ + iopipe.c istore.c lheap.c links.c mount.c mtime.c ohdr.c overhead.c \ + ragged.c stab.c tattr.c testhdf5.c tfile.c th5s.c tmeta.c trefer.c \ + tselect.c tvltypes.c tvlstr.c unlink.c enum.c gass_write.c \ + gass_read.c gass_append.c dpss_read.c dpss_write.c srb_read.c \ + srb_write.c srb_append.c TEST_OBJ=$(TEST_SRC:.c=.lo) @@ -170,4 +172,13 @@ dpss_read: dpss_read.lo dpss_write: dpss_write.lo @$(LT_LINK_EXE) $(CFLAGS) -o $@ dpss_write.lo $(LIB) $(LIBHDF5) $(LDFLAGS) $(LIBS) +srb_read: srb_read.lo + @$(LT_LINK_EXE) $(CFLAGS) -o $@ srb_read.lo $(LIB) $(LIBHDF5) $(LDFLAGS) $(LIBS) + +srb_write: srb_write.lo + @$(LT_LINK_EXE) $(CFLAGS) -o $@ srb_write.lo $(LIB) $(LIBHDF5) $(LDFLAGS) $(LIBS) + +srb_append: srb_append.lo + @$(LT_LINK_EXE) $(CFLAGS) -o $@ srb_append.lo $(LIB) $(LIBHDF5) $(LDFLAGS) $(LIBS) + @CONCLUDE@ diff --git a/test/srb_append.c b/test/srb_append.c new file mode 100644 index 0000000..fff7ac4 --- /dev/null +++ b/test/srb_append.c @@ -0,0 +1,133 @@ +#include <h5test.h> + +#ifndef H5_HAVE_SRB +int main(void) +{ + printf("Test skipped because SRB driver not available\n"); + return 0; +} +#else + +#define fileName "/projects/mdas/srb/SRBVault/slu.ncsa/a.h5" +#define DATASETNAME "Int1Array" +#define NX 5 /* dataset dimensions */ +#define NY 6 +#define RANK 2 + +int main(void) +{ + SRB_Info srb_info={ NULL, + NULL, + NULL, + 0, + 0600, + -1 + }; + hid_t fapl =-1, file; + hid_t dataspace, datatype, dataset; + hsize_t dimsf[2]; + + herr_t status = 0; + int data[NX][NY]; /* data to write */ + int i, j; + + /* + * Data and output buffer initialization. + */ + for (j = 0; j < NX; j++) { + for (i = 0; i < NY; i++) + data[j][i] = i*i + j*j; + } + /* + * 0 1 4 9 16 25 + * 1 2 5 10 17 26 + * 4 5 8 13 20 29 + * 9 10 13 18 25 34 + * 16 17 20 25 32 41 + */ + + /* Create access property list and set the driver to GASS */ + fapl = H5Pcreate (H5P_FILE_ACCESS); + if (fapl < 0) { + printf (" H5Pcreate failed. \n"); + return -1; + } + status = H5Pset_fapl_srb (fapl, srb_info); + if (status < 0) { + printf ("H5Pset_fapl_gass failed. \n"); + return -1; + } + + /* + * Open an existing file using H5F_ACC_RDWR access, + * and srb file access properties. + */ + file = H5Fopen(fileName, H5F_ACC_RDWR, fapl); + if (file < 0) { + printf ("H5Fopen failed. \n"); + return -1; + } + + /* + * Describe the size of the array and create the data space for fixed + * size dataset. + */ + dimsf[0] = NX; + dimsf[1] = NY; + dataspace = H5Screate_simple(RANK, dimsf, NULL); + if (dataspace < 0) { + printf ("H5Screate failed. \n"); + return -1; + } + + /* + * Define datatype for the data in the file. + * We will store little endian INT numbers. + */ + datatype = H5Tcopy(H5T_NATIVE_INT); + if (datatype < 0) { + printf ("H5Tcopy failed. \n"); + return -1; + } + + status = H5Tset_order(datatype, H5T_ORDER_LE); + if (status < 0) { + printf ("H5Tset_order failed. \n"); + return -1; + } + + /* + * Create a new dataset within the file using defined dataspace and + * datatype and default dataset creation properties. + */ + dataset = H5Dcreate(file, DATASETNAME, datatype, dataspace, + H5P_DEFAULT); + if (dataset < 0) { + printf ("H5Dcreate failed. \n"); + return -1; + } + + /* + * Write the data to the dataset using default transfer properties. + */ + status = H5Dwrite(dataset, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, + H5P_DEFAULT, data); + if (status < 0) { + printf ("H5Dwrite failed. \n"); + return -1; + } + + /* + * Close/release resources. + */ + H5Sclose(dataspace); + H5Tclose(datatype); + H5Dclose(dataset); + H5Fclose(file); + H5Pclose(fapl); + + printf("Test finished!\n"); + return 0; +} + +#endif diff --git a/test/srb_read.c b/test/srb_read.c new file mode 100644 index 0000000..6269243 --- /dev/null +++ b/test/srb_read.c @@ -0,0 +1,164 @@ +#include <h5test.h> + +#ifndef H5_HAVE_SRB +int main(void) +{ + printf("Test skipped because SRB driver not available\n"); + return 0; +} +#else + +#define fileName "/projects/mdas/srb/SRBVault/slu.ncsa/a.h5" +#define DATASETNAME "IntArray" +#define NX_SUB 3 /* hyperslab dimensions */ +#define NY_SUB 4 +#define NX 7 /* output buffer dimensions */ +#define NY 7 +#define NZ 3 +#define RANK 2 +#define RANK_OUT 3 + +int main(void) +{ + hid_t fapl=-1, fid = -1, dataset; + hid_t datatype, dataspace; + hid_t memspace; + H5T_class_t class; /* data type class */ + H5T_order_t order; /* data order */ + size_t size; /* + * size of the data element + * stored in file + */ + hsize_t dimsm[3]; /* memory space dimensions */ + hsize_t dims_out[2]; /* dataset dimensions */ + herr_t status; + + int data_out[NX][NY][NZ ]; /* output buffer */ + + hsize_t count[2]; /* size of the hyperslab in the file */ + hssize_t offset[2]; /* hyperslab offset in the file */ + hsize_t count_out[3]; /* size of the hyperslab in memory */ + hssize_t offset_out[3]; /* hyperslab offset in memory */ + int i, j, k, status_n, rank; + SRB_Info srb_info={ NULL, + NULL, + NULL, + 0, + 0600, + -1 + }; + + for (j = 0; j < NX; j++) { + for (i = 0; i < NY; i++) { + for (k = 0; k < NZ ; k++) + data_out[j][i][k] = 0; + } + } + + fapl = H5Pcreate(H5P_FILE_ACCESS); + if (fapl < 0) { + printf (" H5Pcreate failed. \n"); + return -1; + } + + status = H5Pset_fapl_srb(fapl, srb_info); + if (status < 0) { + printf ("H5Pset_fapl_gass failed. \n"); + return -1; + } + + fid = H5Fopen(fileName, H5F_ACC_RDONLY, fapl); + /*fid = H5Fopen(fileName, H5F_ACC_RDWR, H5P_DEFAULT);*/ + if (fid < 0) { + printf ("H5Fopen failed. \n"); + return -1; + } + dataset = H5Dopen(fid, DATASETNAME); + if(dataset<0) { + printf ("H5Dopen failed. \n"); + return -1; + } + /* + * Get datatype and dataspace handles and then query + * dataset class, order, size, rank and dimensions. + */ + datatype = H5Dget_type(dataset); /* datatype handle */ + class = H5Tget_class(datatype); + if (class == H5T_INTEGER) printf("Data set has INTEGER type \n"); + order = H5Tget_order(datatype); + if (order == H5T_ORDER_LE) printf("Little endian order \n"); + + size = H5Tget_size(datatype); + printf(" Data size is %d \n", size); + + dataspace = H5Dget_space(dataset); /* dataspace handle */ + rank = H5Sget_simple_extent_ndims(dataspace); + status_n = H5Sget_simple_extent_dims(dataspace, dims_out, NULL); + printf("rank %d, dimensions %lu x %lu \n", rank, + (unsigned long)(dims_out[0]), (unsigned long)(dims_out[1])); + + /* + * Define hyperslab in the dataset. + */ + offset[0] = 1; + offset[1] = 2; + count[0] = NX_SUB; + count[1] = NY_SUB; + status = H5Sselect_hyperslab(dataspace, H5S_SELECT_SET, offset, NULL, + count, NULL); + + /* + * Define the memory dataspace. + */ + dimsm[0] = NX; + dimsm[1] = NY; + dimsm[2] = NZ ; + memspace = H5Screate_simple(RANK_OUT,dimsm,NULL); + + /* + * Define memory hyperslab. + */ + offset_out[0] = 3; + offset_out[1] = 0; + offset_out[2] = 0; + count_out[0] = NX_SUB; + count_out[1] = NY_SUB; + count_out[2] = 1; + status = H5Sselect_hyperslab(memspace, H5S_SELECT_SET, offset_out, NULL, + count_out, NULL); + + /* + * Read data from hyperslab in the file into the hyperslab in + * memory and display. + */ + status = H5Dread(dataset, H5T_NATIVE_INT, memspace, dataspace, + H5P_DEFAULT, data_out); + for (j = 0; j < NX; j++) { + for (i = 0; i < NY; i++) printf("%d ", data_out[j][i][0]); + printf("\n"); + } + /* + * 0 0 0 0 0 0 0 + * 0 0 0 0 0 0 0 + * 0 0 0 0 0 0 0 + * 3 4 5 6 0 0 0 + * 4 5 6 7 0 0 0 + * 5 6 7 8 0 0 0 + * 0 0 0 0 0 0 0 + */ + + /* + * Close/release resources. + */ + H5Tclose(datatype); + H5Dclose(dataset); + H5Sclose(dataspace); + H5Sclose(memspace); + H5Fclose(fid); + H5Pclose(fapl); + + printf("Test finished!\n"); + return 0; +} + +#endif diff --git a/test/srb_write.c b/test/srb_write.c new file mode 100644 index 0000000..a3e58a5 --- /dev/null +++ b/test/srb_write.c @@ -0,0 +1,131 @@ +#include <h5test.h> + +#ifndef H5_HAVE_SRB +int main(void) +{ + printf("Test skipped because SRB driver not available\n"); + return 0; +} +#else + +#define fileName "/projects/mdas/srb/SRBVault/slu.ncsa/a.h5" +#define DATASETNAME "IntArray" +#define NX 5 /* dataset dimensions */ +#define NY 6 +#define RANK 2 + +int main(void) +{ + hid_t fapl=-1, fid = -1; + hid_t dataspace, datatype, dataset; + hsize_t dimsf[2]; + herr_t status = 0; + int data[NX][NY]; /* data to write */ + int i, j; + SRB_Info srb_info={ "ghidorah.sdsc.edu", /* SRB host address of server */ + "5557", /* SRB host port number */ + NULL, /* SRB Authentication-password */ + 0, /* Storage Type: 0=Unix, + * 1=UniTree, + * 2=HPSS, + * 3=FTP, + * 4=HTTP */ + 0600, /* File mode-Unix access mode */ + -1 /* File Size-Only valid for HPSS, + * -1 is default. */ + }; + + /* + * Data and output buffer initialization. + */ + for (j = 0; j < NX; j++) { + for (i = 0; i < NY; i++) + data[j][i] = i + j; + } + /* + * 0 1 2 3 4 5 + * 1 2 3 4 5 6 + * 2 3 4 5 6 7 + * 3 4 5 6 7 8 + * 4 5 6 7 8 9 + */ + + fapl = H5Pcreate(H5P_FILE_ACCESS); + if (fapl < 0) { + printf (" H5Pcreate failed. \n"); + return -1; + } + + status = H5Pset_fapl_srb(fapl, srb_info); + if (status < 0) { + printf ("H5Pset_fapl_srb failed. \n"); + return -1; + } + + fid = H5Fcreate(fileName, H5F_ACC_TRUNC, H5P_DEFAULT, fapl); + /*fid = H5Fcreate(fileName, H5F_ACC_TRUNC, H5P_DEFAULT, H5P_DEFAULT);*/ + if (fid < 0) { + printf ("H5Fcreate failed. \n"); + return -1; + } + + /* + * Describe the size of the array and create the data space for fixed + * size dataset. + */ + dimsf[0] = NX; + dimsf[1] = NY; + dataspace = H5Screate_simple(RANK, dimsf, NULL); + if (dataspace < 0) { + printf ("H5Screate failed. \n"); + return -1; + } + + /* + * Define datatype for the data in the file. + * We will store little endian INT numbers. + */ + datatype = H5Tcopy(H5T_NATIVE_INT); + if (datatype < 0) { + printf ("H5Tcopy failed. \n"); + return -1; + } + + status = H5Tset_order(datatype, H5T_ORDER_LE); + if (status < 0) { + printf ("H5Tset_order failed. \n"); + return -1; + } + /* + * Create a new dataset within the file using defined dataspace and + * datatype and default dataset creation properties. + */ + dataset = H5Dcreate(fid, DATASETNAME, datatype, dataspace, + H5P_DEFAULT); + if (dataset < 0) { + printf ("H5Dcreate failed. \n"); + return -1; + } + + + /* + * Write the data to the dataset using default transfer properties. + */ + status = H5Dwrite(dataset, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, + H5P_DEFAULT, data); + if (status < 0) { + printf ("H5Dwrite failed. \n"); + return -1; + } + + H5Sclose(dataspace); + H5Tclose(datatype); + H5Dclose(dataset); + H5Fclose(fid); + H5Pclose(fapl); + + printf("Test finished!\n"); + return 0; +} + +#endif |