summaryrefslogtreecommitdiffstats
path: root/HDF5Examples/C/H5D/h5ex_d_transform.c
diff options
context:
space:
mode:
Diffstat (limited to 'HDF5Examples/C/H5D/h5ex_d_transform.c')
-rw-r--r--HDF5Examples/C/H5D/h5ex_d_transform.c154
1 files changed, 154 insertions, 0 deletions
diff --git a/HDF5Examples/C/H5D/h5ex_d_transform.c b/HDF5Examples/C/H5D/h5ex_d_transform.c
new file mode 100644
index 0000000..2d48051
--- /dev/null
+++ b/HDF5Examples/C/H5D/h5ex_d_transform.c
@@ -0,0 +1,154 @@
+/************************************************************
+
+ This example shows how to read and write data to a dataset
+ using a data transform expression. The program first
+ writes integers to a dataset using the transform
+ expression TRANSFORM, then closes the file. Next, it
+ reopens the file, reads back the data without a transform,
+ and outputs the data to the screen. Finally it reads the
+ data using the transform expression RTRANSFORM and outputs
+ the results to the screen.
+
+ ************************************************************/
+
+#include "hdf5.h"
+#include <stdio.h>
+#include <stdlib.h>
+
+#define FILE "h5ex_d_transform.h5"
+#define DATASET "DS1"
+#define DIM0 4
+#define DIM1 7
+#define TRANSFORM "x+1"
+#define RTRANSFORM "x-1"
+
+int
+main(void)
+{
+ hid_t file = H5I_INVALID_HID;
+ hid_t space = H5I_INVALID_HID;
+ hid_t dset = H5I_INVALID_HID;
+ hid_t dxpl = H5I_INVALID_HID;
+ /* Handles */
+ herr_t status;
+ hsize_t dims[2] = {DIM0, DIM1};
+ int wdata[DIM0][DIM1]; /* Write buffer */
+ int rdata[DIM0][DIM1]; /* Read buffer */
+ hsize_t i, j;
+
+ /*
+ * Initialize data.
+ */
+ for (i = 0; i < DIM0; i++)
+ for (j = 0; j < DIM1; j++)
+ wdata[i][j] = i * j - j;
+
+ /*
+ * Output the data to the screen.
+ */
+ printf("Original Data:\n");
+ for (i = 0; i < DIM0; i++) {
+ printf(" [");
+ for (j = 0; j < DIM1; j++)
+ printf(" %3d", wdata[i][j]);
+ printf("]\n");
+ }
+
+ /*
+ * Create a new file using the default properties.
+ */
+ file = H5Fcreate(FILE, H5F_ACC_TRUNC, H5P_DEFAULT, H5P_DEFAULT);
+
+ /*
+ * Create dataspace. Setting maximum size to NULL sets the maximum
+ * size to be the current size.
+ */
+ space = H5Screate_simple(2, dims, NULL);
+
+ /*
+ * Create the dataset transfer property list and define the
+ * transform expression.
+ */
+ dxpl = H5Pcreate(H5P_DATASET_XFER);
+ status = H5Pset_data_transform(dxpl, TRANSFORM);
+
+ /*
+ * Create the dataset using the default properties. Unfortunately
+ * we must save as a native type or the transform operation will
+ * fail.
+ */
+ dset = H5Dcreate(file, DATASET, H5T_NATIVE_INT, space, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
+
+ /*
+ * Write the data to the dataset using the dataset transfer
+ * property list.
+ */
+ status = H5Dwrite(dset, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, dxpl, wdata[0]);
+
+ /*
+ * Close and release resources.
+ */
+ status = H5Pclose(dxpl);
+ status = H5Dclose(dset);
+ status = H5Sclose(space);
+ status = H5Fclose(file);
+
+ /*
+ * Now we begin the read section of this example.
+ */
+
+ /*
+ * Open file and dataset using the default properties.
+ */
+ file = H5Fopen(FILE, H5F_ACC_RDONLY, H5P_DEFAULT);
+ dset = H5Dopen(file, DATASET, H5P_DEFAULT);
+
+ /*
+ * Read the data using the default properties.
+ */
+ status = H5Dread(dset, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata[0]);
+
+ /*
+ * Output the data to the screen.
+ */
+ printf("\nData as written with transform \"%s\":\n", TRANSFORM);
+ for (i = 0; i < DIM0; i++) {
+ printf(" [");
+ for (j = 0; j < DIM1; j++)
+ printf(" %3d", rdata[i][j]);
+ printf("]\n");
+ }
+
+ /*
+ * Create the dataset transfer property list and define the
+ * transform expression.
+ */
+ dxpl = H5Pcreate(H5P_DATASET_XFER);
+ status = H5Pset_data_transform(dxpl, RTRANSFORM);
+
+ /*
+ * Read the data using the dataset transfer property list.
+ */
+ status = H5Dread(dset, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, dxpl, rdata[0]);
+
+ /*
+ * Output the data to the screen.
+ */
+ printf("\nData as written with transform \"%s\" and read with transform \"%s\":\n", TRANSFORM,
+ RTRANSFORM);
+ for (i = 0; i < DIM0; i++) {
+ printf(" [");
+ for (j = 0; j < DIM1; j++)
+ printf(" %3d", rdata[i][j]);
+ printf("]\n");
+ }
+
+ /*
+ * Close and release resources.
+ */
+ status = H5Pclose(dxpl);
+ status = H5Dclose(dset);
+ status = H5Fclose(file);
+
+ return 0;
+}