summaryrefslogtreecommitdiffstats
diff options
context:
space:
mode:
authorRobb Matzke <matzke@llnl.gov>1998-01-21 20:05:56 (GMT)
committerRobb Matzke <matzke@llnl.gov>1998-01-21 20:05:56 (GMT)
commit851b17c87ad6f841efb5a5ba1b90b8c6636a2ffd (patch)
tree354fb32fa49d82ee42e0c19eef56ca2cbfc55a61
parentf5cc0548b8dacfefee4878075f568b5f124e9dc3 (diff)
downloadhdf5-851b17c87ad6f841efb5a5ba1b90b8c6636a2ffd.zip
hdf5-851b17c87ad6f841efb5a5ba1b90b8c6636a2ffd.tar.gz
hdf5-851b17c87ad6f841efb5a5ba1b90b8c6636a2ffd.tar.bz2
[svn-r162] *** empty log message ***
-rw-r--r--src/H5Psimp.c241
1 files changed, 241 insertions, 0 deletions
diff --git a/src/H5Psimp.c b/src/H5Psimp.c
new file mode 100644
index 0000000..a5dd7ec
--- /dev/null
+++ b/src/H5Psimp.c
@@ -0,0 +1,241 @@
+/*
+ * 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 <H5Pprivate.h>
+
+/* Interface initialization */
+#define PABLO_MASK H5P_simp_mask
+#define INTERFACE_INIT NULL
+static intn interface_initialize_g = FALSE;
+
+/*-------------------------------------------------------------------------
+ * Function: H5P_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
+H5P_simp_init (const struct H5O_layout_t *layout, const H5P_t *mem_space,
+ const H5P_t *file_space, H5P_number_t *numbering/*out*/)
+{
+ size_t nelmts;
+
+ FUNC_ENTER (H5P_simp_init, 0);
+
+ /* Check args */
+ assert (layout);
+ assert (mem_space && H5P_SIMPLE==mem_space->type);
+ assert (file_space && H5P_SIMPLE==file_space->type);
+ assert (numbering);
+
+ /* Numbering is implied by the hyperslab, C order */
+ HDmemset (numbering, 0, sizeof(H5P_number_t));
+
+ /* Data can be efficiently copied at any size */
+ nelmts = H5P_get_npoints (file_space);
+
+ FUNC_LEAVE (nelmts);
+}
+
+/*-------------------------------------------------------------------------
+ * Function: H5P_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 H5P_*_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
+H5P_simp_fgath (H5F_t *f, const struct H5O_layout_t *layout,
+ size_t elmt_size, const H5P_t *file_space,
+ const H5P_number_t *numbering, intn start, intn nelmts,
+ void *buf/*out*/)
+{
+ size_t offset[H5O_LAYOUT_NDIMS]; /*offset of hyperslab */
+ size_t size[H5O_LAYOUT_NDIMS]; /*size of hyperslab */
+ intn i; /*counters */
+
+ FUNC_ENTER (H5P_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==H5P_get_npoints (file_space));
+
+ /*
+ * Quincey, this is where we look at FILE_SPACE to decide what the
+ * hyperslab is to read from disk. For now, since the H5P interface
+ * doesn't support hyperslabs, we'll assume the caller is asking for the
+ * entire array. --RPM
+ */
+ assert (nelmts == H5P_get_npoints (file_space));
+ for (i=0; i<layout->ndims; i++) offset[i] = 0;
+ i = H5P_get_dims (file_space, size);
+ assert (i+1 == layout->ndims);
+ size[i] = elmt_size;
+
+ /*
+ * Gather from file.
+ */
+ if (H5F_arr_read (f, layout, size, offset, offset, size, buf/*out*/)<0) {
+ HRETURN_ERROR (H5E_DATASPACE, H5E_READERROR, 0, "read error");
+ }
+
+ FUNC_LEAVE (nelmts);
+}
+
+/*-------------------------------------------------------------------------
+ * Function: H5P_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
+H5P_simp_mscat (const void *tconv_buf, size_t elmt_size,
+ const H5P_t *mem_space, const H5P_number_t *numbering,
+ intn start, intn nelmts, void *buf/*out*/)
+{
+ FUNC_ENTER (H5P_simp_mscat, FAIL);
+
+ /* Check args */
+ assert (tconv_buf);
+ assert (elmt_size>0);
+ assert (mem_space && H5P_SIMPLE==mem_space->type);
+ assert (numbering);
+ assert (nelmts>0);
+ assert (buf);
+
+ /*
+ * The prototype doesn't support strip mining.
+ */
+ assert (0==start);
+ assert (nelmts==H5P_get_npoints (mem_space));
+
+ /*
+ * Quincey, this is where we look at the hyperslab spec of MEM_SPACE to
+ * figure out how to scatter. For now we just assume that data points
+ * are copied directly from TCONV_BUF to BUF.
+ */
+ HDmemcpy (buf, tconv_buf, nelmts*elmt_size);
+
+ FUNC_LEAVE (SUCCEED);
+}
+
+/*-------------------------------------------------------------------------
+ * Function: H5P_simp_mgath
+ *
+ * Purpose: Gathers dataset elements from application memory and copies
+ * them into the data type conversion buffer. NOT IMPLEMENTED
+ * YET.
+ *
+ * Return: Success: Number of elements copied.
+ *
+ * Failure: 0
+ *
+ * Programmer: Robb Matzke
+ * Wednesday, January 21, 1998
+ *
+ * Modifications:
+ *
+ *-------------------------------------------------------------------------
+ */
+size_t
+H5P_simp_mgath (void)
+{
+ FUNC_ENTER (H5P_simp_mgath, 0);
+
+ HRETURN_ERROR (H5E_DATASPACE, H5E_UNSUPPORTED, 0,
+ "not implemented yet");
+
+ FUNC_LEAVE (0);
+}
+
+/*-------------------------------------------------------------------------
+ * Function: H5P_simp_fscat
+ *
+ * Purpose: Scatters dataset elements from the type conversion buffer to
+ * the file. NOT IMPLEMENTED YET.
+ *
+ * Return: Success: SUCCEED
+ *
+ * Failure: FAIL
+ *
+ * Programmer: Robb Matzke
+ * Wednesday, January 21, 1998
+ *
+ * Modifications:
+ *
+ *-------------------------------------------------------------------------
+ */
+herr_t
+H5P_simp_fscat (void)
+{
+ FUNC_ENTER (H5P_simp_fscat, FAIL);
+
+ HRETURN_ERROR (H5E_DATASPACE, H5E_UNSUPPORTED, FAIL,
+ "not implemented yet");
+
+ FUNC_LEAVE (FAIL);
+}