summaryrefslogtreecommitdiffstats
path: root/test/dtransform.c
diff options
context:
space:
mode:
Diffstat (limited to 'test/dtransform.c')
-rw-r--r--test/dtransform.c110
1 files changed, 110 insertions, 0 deletions
diff --git a/test/dtransform.c b/test/dtransform.c
new file mode 100644
index 0000000..0f2e4e8
--- /dev/null
+++ b/test/dtransform.c
@@ -0,0 +1,110 @@
+#include "h5test.h"
+
+#define ROWS 12
+#define COLS 18
+#define TOL 2
+
+int compare(int* a, int* b);
+
+int main()
+{
+
+ hid_t file_id, dxpl_id_f_to_c, dxpl_id_c_to_f, dset_id, datatype, dataspace;
+ hsize_t dim[2] = {ROWS, COLS};
+ char* f_to_c = "(5/9.0)*(x-32)";
+ char* c_to_f = "(9/5.0)*x + 32";
+ int windchillF[ROWS][COLS] =
+ { {36, 31, 25, 19, 13, 7, 1, -5, -11, -16, -22, -28, -34, -40, -46, -52, -57, -63 },
+ {34, 27, 21, 15, 9, 3, -4, -10, -16, -22, -28, -35, -41, -47, -53, -59, -66, -72 } ,
+ {32, 25, 19, 13, 6, 0, -7, -13, -19, -26, -32, -39, -45, -51, -58, -64, -71, -77 },
+ {30, 24, 17, 11, 4, -2, -9, -15, -22, -29, -35, -42, -48, -55, -61, -68, -74, -81 },
+ {29, 23, 16, 9, 3, -4, -11, -17, -24, -31, -37, -44, -51, -58, -64, -71, -78, -84 },
+ {28, 22, 15, 8, 1, -5, -12, -19, -26, -33, -39, -46, -53, -60, -67, -73, -80, -87 },
+ {28, 21, 14, 7, 0, -7, -14, -21, -27, -34, -41, -48, -55, -62, -69, -76, -82, -89 },
+ {27, 20, 13, 6, -1, -8, -15, -22, -29, -36, -43, -50, -57, -64, -71, -78, -84, -91 },
+ {26, 19, 12, 5, -2, -9, -16, -23, -30, -37, -44, -51, -58, -65, -72, -79, -86, -93 },
+ {26, 19, 12, 4, -3, -10, -17, -24, -31, -38, -45, -52, -60, -67, -74, -81, -88, -95},
+ {25, 18, 11, 4, -3, -11, -18, -25, -32, -39, -46, -54, -61, -68, -75, -82, -89, -97},
+ {25, 17, 10, 3, -4, -11, -19, -26, -33, -40, -48, -55, -62, -69, -76, -84, -91, -98}
+ };
+ int windchillCread[ROWS][COLS];
+ int windchillFread[ROWS][COLS];
+ int windchillCcalc[ROWS][COLS];
+ herr_t err;
+ int row, col;
+
+
+ TESTING("data transform");
+
+ dxpl_id_f_to_c = H5Pcreate(H5P_DATASET_XFER);
+ dxpl_id_c_to_f = H5Pcreate(H5P_DATASET_XFER);
+
+ err= H5Pset_data_transform(dxpl_id_f_to_c, f_to_c);
+ err= H5Pset_data_transform(dxpl_id_c_to_f, c_to_f);
+
+
+ file_id = H5Fcreate("dtransform.h5", H5F_ACC_TRUNC, H5P_DEFAULT, H5P_DEFAULT);
+
+ dataspace = H5Screate_simple(2, dim, NULL);
+
+ datatype = H5Tcopy(H5T_NATIVE_INT);
+
+ dset_id = H5Dcreate(file_id, "/transformtest", datatype, dataspace, H5P_DEFAULT);
+
+ err = H5Dwrite(dset_id, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, dxpl_id_f_to_c, windchillF);
+ err = H5Dread(dset_id, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, windchillCread);
+
+ for(row = 0; row<ROWS; row++)
+ {
+ for(col = 0; col<COLS; col++)
+ windchillCcalc[row][col] = (5/9.0)*(windchillF[row][col] - 32);
+ }
+
+ if( (compare(windchillCread, windchillCcalc)) == 0)
+ {
+ fprintf(stderr, "ERROR: Conversion failed to match computed data\n");
+ TEST_ERROR;
+ }
+ else
+ {
+ err = H5Dread(dset_id, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, dxpl_id_c_to_f, windchillFread);
+ if( (compare(windchillF, windchillFread)) == 0)
+ {
+ fprintf(stderr, "ERROR: Conversion failed to match computed data\n");
+ TEST_ERROR;
+ }
+ }
+
+
+ /* Close the objects we opened/created */
+ err = H5Sclose(dataspace);
+ err = H5Tclose(datatype);
+ err = H5Dclose(dset_id);
+ err = H5Fclose(file_id);
+ err = H5Pclose(dxpl_id_f_to_c);
+ err = H5Pclose(dxpl_id_c_to_f);
+
+ PASSED();
+
+ return 0;
+
+error:
+ return -1;
+}
+
+int compare(int* a, int* b)
+{
+ int row, col;
+ for(row = 0; row<ROWS; row++)
+ {
+ for(col = 0; col<COLS; col++)
+ {
+ if( !((a[row*COLS + col] <= (b[row*COLS + col] + TOL)) && (a[row*COLS + col] >= (b[row*COLS + col] - TOL))) )
+ return 0;
+ }
+ }
+ return 1;
+}
+
+
+