[svn-r296] Switched prefix for dataspaces from H5P to H5S

1 files changed, 469 insertions, 0 deletions

diff --git a/src/H5Ssimp.c b/src/H5Ssimp.c
new file mode 100644
index 0000000..99900c2
--- /dev/null
+++ b/src/H5Ssimp.c
@@ -0,0 +1,469 @@
+/*
+ * Copyright (C) 1998 NCSA
+ *                    All rights reserved.
+ *
+ * Programmer:  Robb Matzke <matzke@llnl.gov>
+ *              Wednesday, January 21, 1998
+ *
+ * Purpose:	Simple data space functions.
+ */
+#include <H5private.h>
+#include <H5Eprivate.h>
+#include <H5Sprivate.h>
+#include <H5Vprivate.h>
+
+/* Interface initialization */
+#define PABLO_MASK      H5S_simp_mask
+#define INTERFACE_INIT  NULL
+static intn             interface_initialize_g = FALSE;
+
+/*-------------------------------------------------------------------------
+ * Function:	H5S_simp_init
+ *
+ * Purpose:	Generates element numbering information for the data
+ *		spaces involved in a data space conversion.
+ *
+ * Return:	Success:	Number of elements that can be efficiently
+ *				transferred at a time.
+ *
+ *		Failure:	Zero
+ *
+ * Programmer:	Robb Matzke
+ *              Wednesday, January 21, 1998
+ *
+ * Modifications:
+ *
+ *-------------------------------------------------------------------------
+ */
+size_t
+H5S_simp_init (const struct H5O_layout_t *layout, const H5S_t *mem_space,
+	       const H5S_t *file_space, H5S_number_t *numbering/*out*/)
+{
+    size_t	nelmts;
+    
+    FUNC_ENTER (H5S_simp_init, 0);
+
+    /* Check args */
+    assert (layout);
+    assert (mem_space && H5S_SIMPLE==mem_space->type);
+    assert (file_space && H5S_SIMPLE==file_space->type);
+    assert (numbering);
+
+    /* Numbering is implied by the hyperslab, C order */
+    HDmemset (numbering, 0, sizeof(H5S_number_t));
+
+    /* Data can be efficiently copied at any size */
+    nelmts = H5S_get_npoints (file_space);
+    
+    FUNC_LEAVE (nelmts);
+}
+
+/*-------------------------------------------------------------------------
+ * Function:	H5S_simp_fgath
+ *
+ * Purpose:	Gathers data points from file F and accumulates them in the
+ *		type conversion buffer BUF.  The LAYOUT argument describes
+ *		how the data is stored on disk.  ELMT_SIZE is the size in
+ *		bytes of a datum which this function treats as opaque.
+ *		FILE_SPACE describes the data space of the dataset on disk
+ *		and the elements that have been selected for reading (via
+ *		hyperslab, etc) and NUMBERING describes how those elements
+ *		are numbered (initialized by the H5S_*_init() call). This
+ *		function will copy at most NELMTS elements beginning at the
+ *		element numbered START.
+ *
+ * Return:	Success:	Number of elements copied.
+ *
+ *		Failure:	0
+ *
+ * Programmer:	Robb Matzke
+ *              Wednesday, January 21, 1998
+ *
+ * Modifications:
+ *
+ *-------------------------------------------------------------------------
+ */
+size_t
+H5S_simp_fgath (H5F_t *f, const struct H5O_layout_t *layout,
+		size_t elmt_size, const H5S_t *file_space,
+		const H5S_number_t *numbering, size_t start, size_t nelmts,
+		void *buf/*out*/)
+{
+    size_t	file_offset[H5O_LAYOUT_NDIMS];	/*offset of slab in file*/
+    size_t	hsize[H5O_LAYOUT_NDIMS];	/*size of hyperslab	*/
+    size_t	zero[H5O_LAYOUT_NDIMS];		/*zero			*/
+    size_t	sample[H5O_LAYOUT_NDIMS];	/*hyperslab sampling	*/
+#ifndef LATER
+    intn	file_offset_signed[H5O_LAYOUT_NDIMS];
+#endif
+    intn	space_ndims;			/*dimensionality of space*/
+    intn	i;				/*counters		*/
+
+    FUNC_ENTER (H5S_simp_fgath, 0);
+
+    /* Check args */
+    assert (f);
+    assert (layout);
+    assert (elmt_size>0);
+    assert (file_space);
+    assert (numbering);
+    assert (nelmts>0);
+    assert (buf);
+
+    /*
+     * The prototype doesn't support strip mining.
+     */
+    assert (0==start);
+    assert (nelmts==H5S_get_npoints (file_space));
+
+    /*
+     * Get hyperslab information to determine what elements are being
+     * selected (there might eventually be other selection methods too).
+     * We only support hyperslabs with unit sample because there's no way to
+     * currently pass sample information into H5F_arr_read() much less
+     * H5F_istore_read().
+     */
+#ifdef LATER
+    if ((space_ndims=H5S_get_hyperslab (file_space, file_offset,
+					hsize, sample))<0) {
+	HRETURN_ERROR (H5E_DATASPACE, H5E_CANTINIT, 0,
+		       "unable to retrieve hyperslab parameters");
+    }
+#else
+    if ((space_ndims=H5S_get_hyperslab (file_space, file_offset_signed,
+					hsize, sample))<0) {
+	HRETURN_ERROR (H5E_DATASPACE, H5E_CANTINIT, 0,
+		       "unable to retrieve hyperslab parameters");
+    }
+    for (i=0; i<space_ndims; i++) {
+	assert (file_offset_signed[i]>=0);
+	file_offset[i] = file_offset_signed[i];
+    }
+#endif
+    for (i=0; i<space_ndims; i++) {
+	if (sample[i]!=1) {
+	    HRETURN_ERROR (H5E_ARGS, H5E_BADVALUE, 0,
+			   "hyperslab sampling is not implemented yet");
+	}
+    }
+    file_offset[space_ndims] = 0;
+    hsize[space_ndims] = elmt_size;
+    HDmemset (zero, 0, layout->ndims*sizeof(size_t));
+
+    /*
+     * Gather from file.
+     */
+    if (H5F_arr_read (f, layout, hsize, hsize, zero, file_offset,
+		      buf/*out*/)<0) {
+	HRETURN_ERROR (H5E_DATASPACE, H5E_READERROR, 0, "read error");
+    }
+
+    FUNC_LEAVE (nelmts);
+}
+    
+/*-------------------------------------------------------------------------
+ * Function:	H5S_simp_mscat
+ *
+ * Purpose:	Scatters data points from the type conversion buffer
+ *		TCONV_BUF to the application buffer BUF.  Each element is
+ *		ELMT_SIZE bytes and they are organized in application memory
+ *		according to MEM_SPACE.  The NUMBERING information together
+ *		with START and NELMTS describe how the elements stored in
+ *		TCONV_BUF are globally numbered.
+ *
+ * Return:	Success:	SUCCEED
+ *
+ *		Failure:	FAIL
+ *
+ * Programmer:	Robb Matzke
+ *              Wednesday, January 21, 1998
+ *
+ * Modifications:
+ *
+ *-------------------------------------------------------------------------
+ */
+herr_t
+H5S_simp_mscat (const void *tconv_buf, size_t elmt_size,
+		const H5S_t *mem_space, const H5S_number_t *numbering,
+		size_t start, size_t nelmts, void *buf/*out*/)
+{
+    size_t	mem_offset[H5O_LAYOUT_NDIMS];	/*slab offset in app buf*/
+    size_t	mem_size[H5O_LAYOUT_NDIMS];	/*total size of app buf	*/
+    size_t	hsize[H5O_LAYOUT_NDIMS];	/*size of hyperslab	*/
+    size_t	zero[H5O_LAYOUT_NDIMS];		/*zero			*/
+    size_t	sample[H5O_LAYOUT_NDIMS];	/*hyperslab sampling	*/
+#ifndef LATER
+    intn	mem_offset_signed[H5O_LAYOUT_NDIMS];
+#endif
+    intn	space_ndims;			/*dimensionality of space*/
+    intn	i;				/*counters		*/
+
+    FUNC_ENTER (H5S_simp_mscat, FAIL);
+
+    /* Check args */
+    assert (tconv_buf);
+    assert (elmt_size>0);
+    assert (mem_space && H5S_SIMPLE==mem_space->type);
+    assert (numbering);
+    assert (nelmts>0);
+    assert (buf);
+
+    /*
+     * The prototype doesn't support strip mining.
+     */
+    assert (0==start);
+    assert (nelmts==H5S_get_npoints (mem_space));
+
+    /*
+     * Retrieve hyperslab information to determine what elements are being
+     * selected (there might be other selection methods in the future).  We
+     * only handle hyperslabs with unit sample because there's currently no
+     * way to pass sample information to H5V_hyper_copy().
+     */
+#ifdef LATER
+    if ((space_ndims=H5S_get_hyperslab (mem_space, mem_offset, hsize,
+					sample))<0) {
+	HRETURN_ERROR (H5E_DATASPACE, H5E_CANTINIT, FAIL,
+		       "unable to retrieve hyperslab parameters");
+    }
+#else
+    if ((space_ndims=H5S_get_hyperslab (mem_space, mem_offset_signed,
+					hsize, sample))<0) {
+	HRETURN_ERROR (H5E_DATASPACE, H5E_CANTINIT, FAIL,
+		       "unable to retrieve hyperslab parameters");
+    }
+    for (i=0; i<space_ndims; i++) {
+	assert (mem_offset_signed[i]>=0);
+	mem_offset[i] = mem_offset_signed[i];
+    }
+#endif
+    for (i=0; i<space_ndims; i++) {
+	if (sample[i]!=1) {
+	    HRETURN_ERROR (H5E_ARGS, H5E_BADVALUE, FAIL,
+			   "hyperslab sampling is not implemented yet");
+	}
+    }
+    if (H5S_get_dims (mem_space, mem_size)<0) {
+	HRETURN_ERROR (H5E_DATASPACE, H5E_CANTINIT, FAIL,
+		       "unable to retrieve data space dimensions");
+    }
+    mem_offset[space_ndims] = 0;
+    mem_size[space_ndims] = elmt_size;
+    hsize[space_ndims] = elmt_size;
+    HDmemset (zero, 0, (space_ndims+1)*sizeof(size_t));
+
+    /*
+     * Scatter from conversion buffer to application memory.
+     */
+    if (H5V_hyper_copy (space_ndims+1, hsize, mem_size, mem_offset, buf,
+			hsize, zero, tconv_buf)<0) {
+	HRETURN_ERROR (H5E_DATASPACE, H5E_CANTINIT, FAIL,
+		       "unable to scatter data to memory");
+    }
+
+    FUNC_LEAVE (SUCCEED);
+}
+
+/*-------------------------------------------------------------------------
+ * Function:	H5S_simp_mgath
+ *
+ * Purpose:	Gathers dataset elements from application memory BUF and
+ *		copies them into the data type conversion buffer TCONV_BUF.
+ *		Each element is ELMT_SIZE bytes and arranged in application
+ *		memory according to MEM_SPACE.  The elements selected from
+ *		BUF by MEM_SPACE are numbered according to NUMBERING and the
+ *		caller is requesting that at most NELMTS be gathered
+ *		beginning with number START.  The elements are packed into
+ *		TCONV_BUF in order of their NUMBERING.
+ *
+ * Return:	Success:	Number of elements copied.
+ *
+ *		Failure:	0
+ *
+ * Programmer:	Robb Matzke
+ *              Wednesday, January 21, 1998
+ *
+ * Modifications:
+ *
+ *-------------------------------------------------------------------------
+ */
+size_t
+H5S_simp_mgath (const void *buf, size_t elmt_size,
+		const H5S_t *mem_space, const H5S_number_t *numbering,
+		size_t start, size_t nelmts, void *tconv_buf/*out*/)
+{
+    size_t	mem_offset[H5O_LAYOUT_NDIMS];	/*slab offset in app buf*/
+    size_t	mem_size[H5O_LAYOUT_NDIMS];	/*total size of app buf	*/
+    size_t	hsize[H5O_LAYOUT_NDIMS];	/*size of hyperslab	*/
+    size_t	zero[H5O_LAYOUT_NDIMS];		/*zero			*/
+    size_t	sample[H5O_LAYOUT_NDIMS];	/*hyperslab sampling	*/
+#ifndef LATER
+    intn	mem_offset_signed[H5O_LAYOUT_NDIMS];
+#endif
+    intn	space_ndims;			/*dimensionality of space*/
+    intn	i;				/*counters		*/
+
+    FUNC_ENTER (H5S_simp_mgath, 0);
+
+    /* Check args */
+    assert (buf);
+    assert (elmt_size>0);
+    assert (mem_space && H5S_SIMPLE==mem_space->type);
+    assert (numbering);
+    assert (nelmts>0);
+    assert (tconv_buf);
+
+    /*
+     * The prototype doesn't support strip mining.
+     */
+    assert (0==start);
+    assert (nelmts==H5S_get_npoints (mem_space));
+
+    /*
+     * Retrieve hyperslab information to determine what elements are being
+     * selected (there might be other selection methods in the future).  We
+     * only handle hyperslabs with unit sample because there's currently no
+     * way to pass sample information to H5V_hyper_copy().
+     */
+#ifdef LATER
+    if ((space_ndims=H5S_get_hyperslab (mem_space, mem_offset, hsize,
+					sample))<0) {
+	HRETURN_ERROR (H5E_DATASPACE, H5E_CANTINIT, 0,
+		       "unable to retrieve hyperslab parameters");
+    }
+#else
+    if ((space_ndims=H5S_get_hyperslab (mem_space, mem_offset_signed,
+					hsize, sample))<0) {
+	HRETURN_ERROR (H5E_DATASPACE, H5E_CANTINIT, FAIL,
+		       "unable to retrieve hyperslab parameters");
+    }
+    for (i=0; i<space_ndims; i++) {
+	assert (mem_offset_signed[i]>=0);
+	mem_offset[i] = mem_offset_signed[i];
+    }
+#endif
+    for (i=0; i<space_ndims; i++) {
+	if (sample[i]!=1) {
+	    HRETURN_ERROR (H5E_ARGS, H5E_BADVALUE, 0,
+			   "hyperslab sampling is not implemented yet");
+	}
+    }
+    if (H5S_get_dims (mem_space, mem_size)<0) {
+	HRETURN_ERROR (H5E_DATASPACE, H5E_CANTINIT, 0,
+		       "unable to retrieve data space dimensions");
+    }
+    mem_offset[space_ndims] = 0;
+    mem_size[space_ndims] = elmt_size;
+    hsize[space_ndims] = elmt_size;
+    HDmemset (zero, 0, (space_ndims+1)*sizeof(size_t));
+
+    /*
+     * Scatter from conversion buffer to application memory.
+     */
+    if (H5V_hyper_copy (space_ndims+1, hsize, hsize, zero, tconv_buf,
+			mem_size, mem_offset, buf)<0) {
+	HRETURN_ERROR (H5E_DATASPACE, H5E_CANTINIT, 0,
+		       "unable to scatter data to memory");
+    }
+
+    FUNC_LEAVE (nelmts);
+}
+
+/*-------------------------------------------------------------------------
+ * Function:	H5S_simp_fscat
+ *
+ * Purpose:	Scatters dataset elements from the type conversion buffer BUF
+ *		to the file F where the data points are arranged according to
+ *		the file data space FILE_SPACE and stored according to
+ *		LAYOUT. Each element is ELMT_SIZE bytes and has a unique
+ *		number according to NUMBERING.  The caller is requesting that
+ *		NELMTS elements are coppied beginning with element number
+ *		START.
+ *
+ * Return:	Success:	SUCCEED
+ *
+ *		Failure:	FAIL
+ *
+ * Programmer:	Robb Matzke
+ *              Wednesday, January 21, 1998
+ *
+ * Modifications:
+ *
+ *-------------------------------------------------------------------------
+ */
+herr_t
+H5S_simp_fscat (H5F_t *f, const struct H5O_layout_t *layout,
+		size_t elmt_size, const H5S_t *file_space,
+		const H5S_number_t *numbering, size_t start, size_t nelmts,
+		const void *buf)
+{
+    size_t	file_offset[H5O_LAYOUT_NDIMS];	/*offset of hyperslab	*/
+    size_t	hsize[H5O_LAYOUT_NDIMS];	/*size of hyperslab	*/
+    size_t	zero[H5O_LAYOUT_NDIMS];		/*zero vector		*/
+    size_t	sample[H5O_LAYOUT_NDIMS];	/*hyperslab sampling	*/
+#ifndef LATER
+    intn	file_offset_signed[H5O_LAYOUT_NDIMS];
+#endif
+    intn	space_ndims;			/*space dimensionality	*/
+    intn	i;				/*counters		*/
+
+    FUNC_ENTER (H5S_simp_fscat, FAIL);
+
+    /* Check args */
+    assert (f);
+    assert (layout);
+    assert (elmt_size>0);
+    assert (file_space);
+    assert (numbering);
+    assert (nelmts>0);
+    assert (buf);
+
+    /*
+     * The prototype doesn't support strip mining.
+     */
+    assert (0==start);
+    assert (nelmts==H5S_get_npoints (file_space));
+    
+    /*
+     * Get hyperslab information to determine what elements are being
+     * selected (there might eventually be other selection methods too).
+     * We only support hyperslabs with unit sample because there's no way to
+     * currently pass sample information into H5F_arr_read() much less
+     * H5F_istore_read().
+     */
+#ifdef LATER
+    if ((space_ndims=H5S_get_hyperslab (file_space, file_offset, hsize,
+					sample))<0) {
+	HRETURN_ERROR (H5E_DATASPACE, H5E_CANTINIT, FAIL,
+		       "unable to retrieve hyperslab parameters");
+    }
+#else
+    if ((space_ndims=H5S_get_hyperslab (file_space, file_offset_signed,
+					hsize, sample))<0) {
+	HRETURN_ERROR (H5E_DATASPACE, H5E_CANTINIT, FAIL,
+		       "unable to retrieve hyperslab parameters");
+    }
+    for (i=0; i<space_ndims; i++) {
+	assert (file_offset_signed[i]>=0);
+	file_offset[i] = file_offset_signed[i];
+    }
+#endif
+    for (i=0; i<space_ndims; i++) {
+	if (sample[i]!=1) {
+	    HRETURN_ERROR (H5E_ARGS, H5E_BADVALUE, FAIL,
+			   "hyperslab sampling is not implemented yet");
+	}
+    }
+    file_offset[space_ndims] = 0;
+    hsize[space_ndims] = elmt_size;
+    HDmemset (zero, 0, layout->ndims*sizeof(size_t));
+
+    /*
+     * Scatter to file.
+     */
+    if (H5F_arr_write (f, layout, hsize, hsize, zero, file_offset, buf)<0) {
+	HRETURN_ERROR (H5E_DATASPACE, H5E_WRITEERROR, FAIL, "write error");
+    }
+
+    FUNC_LEAVE (SUCCEED);
+}