summaryrefslogtreecommitdiffstats
path: root/test/tselect.c
diff options
context:
space:
mode:
Diffstat (limited to 'test/tselect.c')
-rw-r--r--test/tselect.c201
1 files changed, 201 insertions, 0 deletions
diff --git a/test/tselect.c b/test/tselect.c
index 59150f4..8114c13 100644
--- a/test/tselect.c
+++ b/test/tselect.c
@@ -87,6 +87,9 @@
#define SPACE7_DIM1 10
#define SPACE7_DIM2 10
#define SPACE7_FILL 254
+#define SPACE7_CHUNK_DIM1 5
+#define SPACE7_CHUNK_DIM2 5
+#define SPACE7_NPOINTS 8
/* 4-D dataset with fixed dimensions */
#define SPACE8_NAME "Space8"
@@ -3988,6 +3991,201 @@ test_select_hyper_chunk(hid_t fapl_plist, hid_t xfer_plist)
/****************************************************************
**
+** test_select_point_chunk(): Test basic H5S (dataspace) selection code.
+** Tests combinations of hyperslab and point selections on
+** chunked datasets.
+**
+****************************************************************/
+static void
+test_select_point_chunk(void)
+{
+ hsize_t dimsf[SPACE7_RANK]; /* dataset dimensions */
+ hsize_t chunk_dimsf[SPACE7_RANK] = {SPACE7_CHUNK_DIM1,SPACE7_CHUNK_DIM2}; /* chunk sizes */
+ unsigned *data; /* data to write */
+ unsigned *tmpdata; /* data to write */
+
+ /*
+ * Data and output buffer initialization.
+ */
+ hid_t file, dataset; /* handles */
+ hid_t dataspace;
+ hid_t pnt1_space; /* Dataspace to hold 1st point selection */
+ hid_t pnt2_space; /* Dataspace to hold 2nd point selection */
+ hid_t hyp1_space; /* Dataspace to hold 1st hyperslab selection */
+ hid_t hyp2_space; /* Dataspace to hold 2nd hyperslab selection */
+ hid_t dcpl;
+ herr_t ret; /* Generic return value */
+
+ unsigned *data_out; /* output buffer */
+#ifdef LATER
+ unsigned *tmpdata_out; /* output buffer */
+#endif /* LATER */
+
+ hssize_t start[SPACE7_RANK]; /* hyperslab offset */
+ hsize_t count[SPACE7_RANK]; /* size of the hyperslab */
+ hssize_t points[SPACE7_NPOINTS][SPACE7_RANK]; /* points for selection */
+ unsigned i, j; /* Local index variables */
+
+ /* Output message about test being performed */
+ MESSAGE(5, ("Testing Point Selections on Chunked Datasets\n"));
+
+ /* Allocate the transfer buffers */
+ data = (unsigned*)HDmalloc(sizeof(unsigned)*SPACE7_DIM1*SPACE7_DIM2);
+ data_out = (unsigned*)HDcalloc(SPACE7_DIM1*SPACE7_DIM2,sizeof(unsigned));
+
+ /*
+ * Data buffer initialization.
+ */
+ tmpdata = data;
+ for (i = 0; i < SPACE7_DIM1; i++)
+ for (j = 0; j < SPACE7_DIM1; j++)
+ *tmpdata++ = ((i*SPACE7_DIM2)+j)%256;
+
+ /*
+ * Create a new file using H5F_ACC_TRUNC access,
+ * the default file creation properties and file
+ * access properties.
+ */
+ file = H5Fcreate (FILENAME, H5F_ACC_TRUNC, H5P_DEFAULT, H5P_DEFAULT);
+ CHECK(file, FAIL, "H5Fcreate");
+
+ /* Create file dataspace */
+ dimsf[0] = SPACE7_DIM1;
+ dimsf[1] = SPACE7_DIM2;
+ dataspace = H5Screate_simple (SPACE7_RANK, dimsf, NULL);
+ CHECK(dataspace, FAIL, "H5Screate_simple");
+
+ /*
+ * Create a new dataset within the file using defined dataspace and
+ * chunking properties.
+ */
+ dcpl = H5Pcreate (H5P_DATASET_CREATE);
+ CHECK(dcpl, FAIL, "H5Pcreate");
+ ret = H5Pset_chunk (dcpl, SPACE7_RANK, chunk_dimsf);
+ CHECK(ret, FAIL, "H5Pset_chunk");
+ dataset = H5Dcreate (file, DATASETNAME, H5T_NATIVE_UCHAR, dataspace, dcpl);
+ CHECK(dataset, FAIL, "H5Dcreate");
+
+ /* Create 1st point selection */
+ pnt1_space = H5Scopy (dataspace);
+ CHECK(pnt1_space, FAIL, "H5Scopy");
+
+ points[0][0]=3;
+ points[0][1]=3;
+ points[1][0]=3;
+ points[1][1]=8;
+ points[2][0]=8;
+ points[2][1]=3;
+ points[3][0]=8;
+ points[3][1]=8;
+ points[4][0]=1; /* In same chunk as point #0, but "earlier" in chunk */
+ points[4][1]=1;
+ points[5][0]=1; /* In same chunk as point #1, but "earlier" in chunk */
+ points[5][1]=6;
+ points[6][0]=6; /* In same chunk as point #2, but "earlier" in chunk */
+ points[6][1]=1;
+ points[7][0]=6; /* In same chunk as point #3, but "earlier" in chunk */
+ points[7][1]=6;
+ ret = H5Sselect_elements(pnt1_space,H5S_SELECT_SET,SPACE7_NPOINTS,(const hssize_t **)points);
+ CHECK(ret, FAIL, "H5Sselect_elements");
+
+ /* Create 1st hyperslab selection */
+ hyp1_space = H5Scopy (dataspace);
+ CHECK(hyp1_space, FAIL, "H5Scopy");
+
+ start[0]=2; start[1]=2;
+ count[0]=4; count[1]=2;
+ ret = H5Sselect_hyperslab(hyp1_space,H5S_SELECT_SET,start,NULL,count,NULL);
+ CHECK(ret, FAIL, "H5Sselect_hyperslab");
+
+ /* Write out data using 1st point selection for file & hyperslab for memory */
+ ret=H5Dwrite(dataset,H5T_NATIVE_UINT,hyp1_space,pnt1_space,H5P_DEFAULT,data);
+ CHECK(ret, FAIL, "H5Dwrite");
+
+ /* Create 2nd point selection */
+ pnt2_space = H5Scopy (dataspace);
+ CHECK(pnt2_space, FAIL, "H5Scopy");
+
+ points[0][0]=4;
+ points[0][1]=4;
+ points[1][0]=4;
+ points[1][1]=9;
+ points[2][0]=9;
+ points[2][1]=4;
+ points[3][0]=9;
+ points[3][1]=9;
+ points[4][0]=2; /* In same chunk as point #0, but "earlier" in chunk */
+ points[4][1]=2;
+ points[5][0]=2; /* In same chunk as point #1, but "earlier" in chunk */
+ points[5][1]=7;
+ points[6][0]=7; /* In same chunk as point #2, but "earlier" in chunk */
+ points[6][1]=2;
+ points[7][0]=7; /* In same chunk as point #3, but "earlier" in chunk */
+ points[7][1]=7;
+ ret = H5Sselect_elements(pnt2_space,H5S_SELECT_SET,SPACE7_NPOINTS,(const hssize_t **)points);
+ CHECK(ret, FAIL, "H5Sselect_elements");
+
+ /* Create 2nd hyperslab selection */
+ hyp2_space = H5Scopy (dataspace);
+ CHECK(hyp2_space, FAIL, "H5Scopy");
+
+ start[0]=2; start[1]=4;
+ count[0]=4; count[1]=2;
+ ret = H5Sselect_hyperslab(hyp2_space,H5S_SELECT_SET,start,NULL,count,NULL);
+ CHECK(ret, FAIL, "H5Sselect_hyperslab");
+
+ /* Write out data using 2nd hyperslab selection for file & point for memory */
+ ret=H5Dwrite(dataset,H5T_NATIVE_UINT,pnt2_space,hyp2_space,H5P_DEFAULT,data);
+ CHECK(ret, FAIL, "H5Dwrite");
+
+ /* Close everything (except selections) */
+ ret=H5Pclose (dcpl);
+ CHECK(ret, FAIL, "H5Pclose");
+ ret=H5Sclose (dataspace);
+ CHECK(ret, FAIL, "H5Sclose");
+ ret=H5Dclose (dataset);
+ CHECK(ret, FAIL, "H5Dclose");
+ ret=H5Fclose (file);
+ CHECK(ret, FAIL, "H5Fclose");
+
+ /* Re-open file & dataset */
+ file = H5Fopen (FILENAME, H5F_ACC_RDONLY, H5P_DEFAULT);
+ CHECK(file, FAIL, "H5Fopen");
+ dataset = H5Dopen (file, DATASETNAME);
+ CHECK(dataset, FAIL, "H5Dopen");
+
+ /* Read data using 1st point selection for file and hyperslab for memory */
+ ret=H5Dread(dataset,H5T_NATIVE_UINT,hyp1_space,pnt1_space,H5P_DEFAULT,data_out);
+ CHECK(ret, FAIL, "H5Dread");
+
+/* Verify data (later) */
+
+ /* Read data using 2nd hyperslab selection for file and point for memory */
+ ret=H5Dread(dataset,H5T_NATIVE_UINT,pnt2_space,hyp2_space,H5P_DEFAULT,data_out);
+ CHECK(ret, FAIL, "H5Dread");
+
+/* Verify data (later) */
+
+ /* Close everything (inclusing selections) */
+ ret=H5Sclose (pnt1_space);
+ CHECK(ret, FAIL, "H5Sclose");
+ ret=H5Sclose (pnt2_space);
+ CHECK(ret, FAIL, "H5Sclose");
+ ret=H5Sclose (hyp1_space);
+ CHECK(ret, FAIL, "H5Sclose");
+ ret=H5Sclose (hyp2_space);
+ CHECK(ret, FAIL, "H5Sclose");
+ ret=H5Dclose (dataset);
+ CHECK(ret, FAIL, "H5Dclose");
+ ret=H5Fclose (file);
+ CHECK(ret, FAIL, "H5Fclose");
+
+ free (data);
+ free (data_out);
+} /* test_select_point_chunk() */
+
+/****************************************************************
+**
** test_select_valid(): Test basic H5S (dataspace) selection code.
** Tests selection validity
**
@@ -6267,6 +6465,9 @@ test_select(void)
/* Test "same shape" routine */
test_shape_same();
+ /* Test point selections in chunked datasets */
+ test_select_point_chunk();
+
} /* test_select() */