diff options
Diffstat (limited to 'test')
-rw-r--r-- | test/trefer.c | 279 |
1 files changed, 279 insertions, 0 deletions
diff --git a/test/trefer.c b/test/trefer.c index db5d558..6bd69ef 100644 --- a/test/trefer.c +++ b/test/trefer.c @@ -39,6 +39,11 @@ #define SPACE2_DIM1 10 #define SPACE2_DIM2 10 +/* Larger 1-D dataset with fixed dimensions */ +#define SPACE3_NAME "Space3" +#define SPACE3_RANK 1 +#define SPACE3_DIM1 100 + /* Element selection information */ #define POINT1_NPOINTS 10 @@ -560,6 +565,279 @@ test_reference_region(void) /**************************************************************** ** +** test_reference_region_1D(): Test H5R (reference) object reference code. +** Tests 1-D references to various kinds of objects +** +****************************************************************/ +static void +test_reference_region_1D(void) +{ + hid_t fid1; /* HDF5 File IDs */ + hid_t dset1, /* Dataset ID */ + dset3; /* Dereferenced dataset ID */ + hid_t sid1, /* Dataspace ID #1 */ + sid3; /* Dataspace ID #3 */ + hsize_t dims1[] = {SPACE1_DIM1}, + dims3[] = {SPACE3_DIM1}; + hssize_t start[SPACE3_RANK]; /* Starting location of hyperslab */ + hsize_t stride[SPACE3_RANK]; /* Stride of hyperslab */ + hsize_t count[SPACE3_RANK]; /* Element count of hyperslab */ + hsize_t block[SPACE3_RANK]; /* Block size of hyperslab */ + hssize_t coord1[POINT1_NPOINTS][SPACE3_RANK]; /* Coordinates for point selection */ + hsize_t * coords; /* Coordinate buffer */ + hsize_t low[SPACE3_RANK]; /* Selection bounds */ + hsize_t high[SPACE3_RANK]; /* Selection bounds */ + hdset_reg_ref_t *wbuf, /* buffer to write to disk */ + *rbuf; /* buffer read from disk */ + uint8_t *dwbuf, /* Buffer for writing numeric data to disk */ + *drbuf; /* Buffer for reading numeric data from disk */ + uint8_t *tu8; /* Temporary pointer to uint8 data */ + int i; /* counting variables */ + herr_t ret; /* Generic return value */ + + /* Output message about test being performed */ + MESSAGE(5, ("Testing 1-D Dataset Region Reference Functions\n")); + + /* Allocate write & read buffers */ + wbuf=calloc(sizeof(hdset_reg_ref_t), SPACE1_DIM1); + rbuf=malloc(sizeof(hdset_reg_ref_t)*SPACE1_DIM1); + dwbuf=malloc(sizeof(uint8_t)*SPACE3_DIM1); + drbuf=calloc(sizeof(uint8_t),SPACE3_DIM1); + + /* Create file */ + fid1 = H5Fcreate(FILE2, H5F_ACC_TRUNC, H5P_DEFAULT, H5P_DEFAULT); + CHECK(fid1, FAIL, "H5Fcreate"); + + /* Create dataspace for datasets */ + sid3 = H5Screate_simple(SPACE3_RANK, dims3, NULL); + CHECK(sid3, FAIL, "H5Screate_simple"); + + /* Create a dataset */ + dset3=H5Dcreate(fid1,"Dataset2",H5T_STD_U8LE,sid3,H5P_DEFAULT); + CHECK(dset3, FAIL, "H5Dcreate"); + + for(tu8=dwbuf,i=0; i<SPACE3_DIM1; i++) + *tu8++=i*3; + + /* Write selection to disk */ + ret=H5Dwrite(dset3,H5T_STD_U8LE,H5S_ALL,H5S_ALL,H5P_DEFAULT,dwbuf); + CHECK(ret, FAIL, "H5Dwrite"); + + /* Close Dataset */ + ret = H5Dclose(dset3); + CHECK(ret, FAIL, "H5Dclose"); + + /* Create dataspace for the reference dataset */ + sid1 = H5Screate_simple(SPACE1_RANK, dims1, NULL); + CHECK(sid1, FAIL, "H5Screate_simple"); + + /* Create a dataset */ + dset1=H5Dcreate(fid1,"Dataset1",H5T_STD_REF_DSETREG,sid1,H5P_DEFAULT); + CHECK(ret, FAIL, "H5Dcreate"); + + /* Create references */ + + /* Select 15 2x1 hyperslabs for first reference */ + start[0]=2; + stride[0]=5; + count[0]=15; + block[0]=2; + ret = H5Sselect_hyperslab(sid3,H5S_SELECT_SET,start,stride,count,block); + CHECK(ret, FAIL, "H5Sselect_hyperslab"); + + ret = (int)H5Sget_select_npoints(sid3); + VERIFY(ret, 30, "H5Sget_select_npoints"); + + /* Store first dataset region */ + ret = H5Rcreate(&wbuf[0],fid1,"/Dataset2",H5R_DATASET_REGION,sid3); + CHECK(ret, FAIL, "H5Rcreate"); + ret = H5Rget_obj_type(dset1,H5R_DATASET_REGION,&wbuf[0]); + VERIFY(ret, H5G_DATASET, "H5Rget_obj_type"); + + /* Select sequence of ten points for second reference */ + coord1[0][0]=16; + coord1[1][0]=22; + coord1[2][0]=38; + coord1[3][0]=41; + coord1[4][0]=52; + coord1[5][0]=63; + coord1[6][0]=70; + coord1[7][0]=89; + coord1[8][0]=97; + coord1[9][0]=03; + ret = H5Sselect_elements(sid3,H5S_SELECT_SET,POINT1_NPOINTS,(const hssize_t **)coord1); + CHECK(ret, FAIL, "H5Sselect_elements"); + + ret = (int)H5Sget_select_npoints(sid3); + VERIFY(ret, 10, "H5Sget_select_npoints"); + + /* Store second dataset region */ + ret = H5Rcreate(&wbuf[1],fid1,"/Dataset2",H5R_DATASET_REGION,sid3); + CHECK(ret, FAIL, "H5Rcreate"); + + /* Write selection to disk */ + ret=H5Dwrite(dset1,H5T_STD_REF_DSETREG,H5S_ALL,H5S_ALL,H5P_DEFAULT,wbuf); + CHECK(ret, FAIL, "H5Dwrite"); + + /* Close disk dataspace */ + ret = H5Sclose(sid1); + CHECK(ret, FAIL, "H5Sclose"); + + /* Close Dataset */ + ret = H5Dclose(dset1); + CHECK(ret, FAIL, "H5Dclose"); + + /* Close uint8 dataset dataspace */ + ret = H5Sclose(sid3); + CHECK(ret, FAIL, "H5Sclose"); + + /* Close file */ + ret = H5Fclose(fid1); + CHECK(ret, FAIL, "H5Fclose"); + + /* Re-open the file */ + fid1 = H5Fopen(FILE2, H5F_ACC_RDWR, H5P_DEFAULT); + CHECK(fid1, FAIL, "H5Fopen"); + + /* Open the dataset */ + dset1=H5Dopen(fid1,"/Dataset1"); + CHECK(dset1, FAIL, "H5Dopen"); + + /* Read selection from disk */ + ret=H5Dread(dset1,H5T_STD_REF_DSETREG,H5S_ALL,H5S_ALL,H5P_DEFAULT,rbuf); + CHECK(ret, FAIL, "H5Dread"); + + /* Try to open objects */ + dset3 = H5Rdereference(dset1,H5R_DATASET_REGION,&rbuf[0]); + CHECK(dset3, FAIL, "H5Rdereference"); + + /* Check what H5Rget_obj_type function returns */ + ret = H5Rget_obj_type(dset1, H5R_DATASET_REGION,&rbuf[0]); + VERIFY(ret, H5G_DATASET, "H5Rget_obj_type"); + + /* Check information in referenced dataset */ + sid1 = H5Dget_space(dset3); + CHECK(sid1, FAIL, "H5Dget_space"); + + ret=(int)H5Sget_simple_extent_npoints(sid1); + VERIFY(ret, 100, "H5Sget_simple_extent_npoints"); + + /* Read from disk */ + ret=H5Dread(dset3,H5T_STD_U8LE,H5S_ALL,H5S_ALL,H5P_DEFAULT,drbuf); + CHECK(ret, FAIL, "H5Dread"); + + for(tu8=(uint8_t *)drbuf,i=0; i<SPACE3_DIM1; i++,tu8++) + VERIFY(*tu8, (uint8_t)(i*3), "Data"); + + /* Get the hyperslab selection */ + sid3=H5Rget_region(dset1,H5R_DATASET_REGION,&rbuf[0]); + CHECK(sid3, FAIL, "H5Rget_region"); + + /* Verify correct hyperslab selected */ + ret = (int)H5Sget_select_npoints(sid3); + VERIFY(ret, 30, "H5Sget_select_npoints"); + ret = (int)H5Sget_select_hyper_nblocks(sid3); + VERIFY(ret, 15, "H5Sget_select_hyper_nblocks"); + coords=HDmalloc(ret*SPACE3_RANK*sizeof(hsize_t)*2); /* allocate space for the hyperslab blocks */ + ret = H5Sget_select_hyper_blocklist(sid3,(hsize_t)0,(hsize_t)ret,coords); + CHECK(ret, FAIL, "H5Sget_select_hyper_blocklist"); + VERIFY(coords[0], 2, "Hyperslab Coordinates"); + VERIFY(coords[1], 3, "Hyperslab Coordinates"); + VERIFY(coords[2], 7, "Hyperslab Coordinates"); + VERIFY(coords[3], 8, "Hyperslab Coordinates"); + VERIFY(coords[4],12, "Hyperslab Coordinates"); + VERIFY(coords[5],13, "Hyperslab Coordinates"); + VERIFY(coords[6],17, "Hyperslab Coordinates"); + VERIFY(coords[7],18, "Hyperslab Coordinates"); + VERIFY(coords[8],22, "Hyperslab Coordinates"); + VERIFY(coords[9],23, "Hyperslab Coordinates"); + VERIFY(coords[10],27, "Hyperslab Coordinates"); + VERIFY(coords[11],28, "Hyperslab Coordinates"); + VERIFY(coords[12],32, "Hyperslab Coordinates"); + VERIFY(coords[13],33, "Hyperslab Coordinates"); + VERIFY(coords[14],37, "Hyperslab Coordinates"); + VERIFY(coords[15],38, "Hyperslab Coordinates"); + VERIFY(coords[16],42, "Hyperslab Coordinates"); + VERIFY(coords[17],43, "Hyperslab Coordinates"); + VERIFY(coords[18],47, "Hyperslab Coordinates"); + VERIFY(coords[19],48, "Hyperslab Coordinates"); + VERIFY(coords[20],52, "Hyperslab Coordinates"); + VERIFY(coords[21],53, "Hyperslab Coordinates"); + VERIFY(coords[22],57, "Hyperslab Coordinates"); + VERIFY(coords[23],58, "Hyperslab Coordinates"); + VERIFY(coords[24],62, "Hyperslab Coordinates"); + VERIFY(coords[25],63, "Hyperslab Coordinates"); + VERIFY(coords[26],67, "Hyperslab Coordinates"); + VERIFY(coords[27],68, "Hyperslab Coordinates"); + VERIFY(coords[28],72, "Hyperslab Coordinates"); + VERIFY(coords[29],73, "Hyperslab Coordinates"); + HDfree(coords); + ret = H5Sget_select_bounds(sid3,low,high); + CHECK(ret, FAIL, "H5Sget_select_bounds"); + VERIFY(low[0], 2, "Selection Bounds"); + VERIFY(high[0], 73, "Selection Bounds"); + + /* Close region space */ + ret = H5Sclose(sid3); + CHECK(ret, FAIL, "H5Sclose"); + + /* Get the element selection */ + sid3=H5Rget_region(dset1,H5R_DATASET_REGION,&rbuf[1]); + CHECK(sid3, FAIL, "H5Rget_region"); + + /* Verify correct elements selected */ + ret = (int)H5Sget_select_npoints(sid3); + VERIFY(ret, 10, "H5Sget_select_npoints"); + ret = (int)H5Sget_select_elem_npoints(sid3); + VERIFY(ret, 10, "H5Sget_select_elem_npoints"); + coords=HDmalloc(ret*SPACE3_RANK*sizeof(hsize_t)); /* allocate space for the element points */ + ret = H5Sget_select_elem_pointlist(sid3,(hsize_t)0,(hsize_t)ret,coords); + CHECK(ret, FAIL, "H5Sget_select_elem_pointlist"); + VERIFY((hssize_t)coords[0], coord1[0][0], "Element Coordinates"); + VERIFY((hssize_t)coords[1], coord1[1][0], "Element Coordinates"); + VERIFY((hssize_t)coords[2], coord1[2][0], "Element Coordinates"); + VERIFY((hssize_t)coords[3], coord1[3][0], "Element Coordinates"); + VERIFY((hssize_t)coords[4], coord1[4][0], "Element Coordinates"); + VERIFY((hssize_t)coords[5], coord1[5][0], "Element Coordinates"); + VERIFY((hssize_t)coords[6], coord1[6][0], "Element Coordinates"); + VERIFY((hssize_t)coords[7], coord1[7][0], "Element Coordinates"); + VERIFY((hssize_t)coords[8], coord1[8][0], "Element Coordinates"); + VERIFY((hssize_t)coords[9], coord1[9][0], "Element Coordinates"); + HDfree(coords); + ret = H5Sget_select_bounds(sid3,low,high); + CHECK(ret, FAIL, "H5Sget_select_bounds"); + VERIFY(low[0], 3, "Selection Bounds"); + VERIFY(high[0], 97, "Selection Bounds"); + + /* Close region space */ + ret = H5Sclose(sid3); + CHECK(ret, FAIL, "H5Sclose"); + + /* Close first space */ + ret = H5Sclose(sid1); + CHECK(ret, FAIL, "H5Sclose"); + + /* Close dereferenced Dataset */ + ret = H5Dclose(dset3); + CHECK(ret, FAIL, "H5Dclose"); + + /* Close Dataset */ + ret = H5Dclose(dset1); + CHECK(ret, FAIL, "H5Dclose"); + + /* Close file */ + ret = H5Fclose(fid1); + CHECK(ret, FAIL, "H5Fclose"); + + /* Free memory buffers */ + free(wbuf); + free(rbuf); + free(dwbuf); + free(drbuf); +} /* test_reference_region_1D() */ + +/**************************************************************** +** ** test_reference_obj_deleted(): Test H5R (reference) object reference code. ** Tests for correct failures for deleted and non-existent objects ** @@ -670,6 +948,7 @@ test_reference(void) test_reference_obj(); /* Test basic H5R object reference code */ test_reference_region(); /* Test basic H5R dataset region reference code */ + test_reference_region_1D(); /* Test H5R dataset region reference code for 1-D datasets */ test_reference_obj_deleted(); /* Test H5R object reference code for deleted objects */ } /* test_reference() */ |