From 5d9f7ebd4c1fb7a008aaa5190fa212dbca5c711c Mon Sep 17 00:00:00 2001
From: Leon Arber <larber@ncsa.uiuc.edu>
Date: Tue, 5 Oct 2004 17:44:17 -0500
Subject: [svn-r9369] Purpose:

Added tests for most of the types that were added to the data transform code.

Further tests will also be added shortly.


Platforms tested:
sol + eirene

Misc. update:
---
 test/dtransform.c | 516 ++++++++++++++++++++++++++++++++++++++++++------------
 1 file changed, 402 insertions(+), 114 deletions(-)

diff --git a/test/dtransform.c b/test/dtransform.c
index 1d50df9..37b2a77 100644
--- a/test/dtransform.c
+++ b/test/dtransform.c
@@ -3,18 +3,44 @@
 #define ROWS    12
 #define COLS    18
 
-int compare_int(int* a, int* b, int tol);
-int compare_float(float* a, float* b, double tol);
+int init_test(hid_t* dxpl_id_c_to_f, hid_t* dxpl_id_simple, hid_t* dxpl_id_polynomial, hid_t* dxpl_id_polynomial_copy, hid_t* dxpl_id_c_to_f_copy);
+int test_int(const hid_t dxpl_id_c_to_f, const hid_t dxpl_id_polynomial);
+int test_float(const hid_t dxpl_id_c_to_f, const hid_t dxpl_id_polynomial);
+int test_char(const hid_t dxpl_id_c_to_f);
+int test_double(const hid_t dxpl_id_c_to_f);
+int test_short(const hid_t dxpl_id_c_to_f);
+int test_copy(const hid_t dxpl_id_c_to_f_copy, const hid_t dxpl_id_polynomial_copy);
+int test_trivial(const hid_t dxpl_id_simple);
+
+/* These are used everywhere and are init'ed in init_test() */
+hid_t dset_id_int;
+hid_t dset_id_float;
+hid_t dset_id_char;
+hid_t dset_id_uchar;
+hid_t dset_id_schar;
+hid_t dset_id_short;
+hid_t dset_id_long;
+hid_t dset_id_ulong;
+hid_t dset_id_uint;
+hid_t dset_id_double;
+hid_t dset_id_ushort;
+
+const float windchillFfloat[ROWS][COLS] = 
+    {   {36.0, 31.0, 25.0, 19.0, 13.0, 7.0, 1.0, -5.0, -11.0, -16.0, -22.0, -28.0, -34.0, -40.0, -46.0, -52.0, -57.0, -63.0 },
+	{34.0, 27.0, 21.0, 15.0, 9.0, 3.0, -4.0, -10.0, -16.0, -22.0, -28.0, -35.0, -41.0, -47.0, -53.0, -59.0, -66.0, -72.0 } ,  
+	{32.0, 25.0, 19.0, 13.0, 6.0, 0.0, -7.0, -13.0, -19.0, -26.0, -32.0, -39.0, -45.0, -51.0, -58.0, -64.0, -71.0, -77.0 },  
+	{30.0, 24.0, 17.0, 11.0, 4.0, -2.0, -9.0, -15.0, -22.0, -29.0, -35.0, -42.0, -48.0, -55.0, -61.0, -68.0, -74.0, -81.0 },  
+	{29.0, 23.0, 16.0, 9.0, 3.0, -4.0, -11.0, -17.0, -24.0, -31.0, -37.0, -44.0, -51.0, -58.0, -64.0, -71.0, -78.0, -84.0 },  
+	{28.0, 22.0, 15.0, 8.0, 1.0, -5.0, -12.0, -19.0, -26.0, -33.0, -39.0, -46.0, -53.0, -60.0, -67.0, -73.0, -80.0, -87.0 },  
+	{28.0, 21.0, 14.0, 7.0, 0.0, -7.0, -14.0, -21.0, -27.0, -34.0, -41.0, -48.0, -55.0, -62.0, -69.0, -76.0, -82.0, -89.0 },  
+	{27.0, 20.0, 13.0, 6.0, -1.0, -8.0, -15.0, -22.0, -29.0, -36.0, -43.0, -50.0, -57.0, -64.0, -71.0, -78.0, -84.0, -91.0 },  
+	{26.0, 19.0, 12.0, 5.0, -2.0, -9.0, -16.0, -23.0, -30.0, -37.0, -44.0, -51.0, -58.0, -65.0, -72.0, -79.0, -86.0, -93.0 },
+	{26.0, 19.0, 12.0, 4.0, -3.0, -10.0, -17.0, -24.0, -31.0, -38.0, -45.0, -52.0, -60.0, -67.0, -74.0, -81.0, -88.0, -95.0},
+	{25.0, 18.0, 11.0, 4.0, -3.0, -11.0, -18.0, -25.0, -32.0, -39.0, -46.0, -54.0, -61.0, -68.0, -75.0, -82.0, -89.0, -97.0},  
+	{25.0, 17.0, 10.0, 3.0, -4.0, -11.0, -19.0, -26.0, -33.0, -40.0, -48.0, -55.0, -62.0, -69.0, -76.0, -84.0, -91.0, -98.0} 
+    };
 
-int main(void)
-{
-    hid_t file_id, dxpl_id_f_to_c, dxpl_id_c_to_f, dxpl_id_c_to_f_copy, dxpl_id_simple, dxpl_id_polynomial, dxpl_id_polynomial_copy, dset_id_int, dset_id_float, datatype, dataspace;
-    hsize_t dim[2] = {ROWS, COLS};
-    const char* f_to_c = "(5/9.0)*(x-32)";
-    const char* c_to_f = "(9/5.0)*x + 32";
-    const char* polynomial = "(2+x)* ((x-8)/2)";
-    const char* simple = "(4/2) * ( (2 + 4)/(5 - 2.5))"; /* this equals 4.8 */ 
-    int windchillFint[ROWS][COLS] =
+const int windchillFint[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 }, 
@@ -28,69 +54,82 @@ int main(void)
         {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 windchillFintread[ROWS][COLS];
-    float windchillFfloat[ROWS][COLS] = 
-    {   {36.0, 31.0, 25.0, 19.0, 13.0, 7.0, 1.0, -5.0, -11.0, -16.0, -22.0, -28.0, -34.0, -40.0, -46.0, -52.0, -57.0, -63.0 },
-	{34.0, 27.0, 21.0, 15.0, 9.0, 3.0, -4.0, -10.0, -16.0, -22.0, -28.0, -35.0, -41.0, -47.0, -53.0, -59.0, -66.0, -72.0 } ,  
-	{32.0, 25.0, 19.0, 13.0, 6.0, 0.0, -7.0, -13.0, -19.0, -26.0, -32.0, -39.0, -45.0, -51.0, -58.0, -64.0, -71.0, -77.0 },  
-	{30.0, 24.0, 17.0, 11.0, 4.0, -2.0, -9.0, -15.0, -22.0, -29.0, -35.0, -42.0, -48.0, -55.0, -61.0, -68.0, -74.0, -81.0 },  
-	{29.0, 23.0, 16.0, 9.0, 3.0, -4.0, -11.0, -17.0, -24.0, -31.0, -37.0, -44.0, -51.0, -58.0, -64.0, -71.0, -78.0, -84.0 },  
-	{28.0, 22.0, 15.0, 8.0, 1.0, -5.0, -12.0, -19.0, -26.0, -33.0, -39.0, -46.0, -53.0, -60.0, -67.0, -73.0, -80.0, -87.0 },  
-	{28.0, 21.0, 14.0, 7.0, 0.0, -7.0, -14.0, -21.0, -27.0, -34.0, -41.0, -48.0, -55.0, -62.0, -69.0, -76.0, -82.0, -89.0 },  
-	{27.0, 20.0, 13.0, 6.0, -1.0, -8.0, -15.0, -22.0, -29.0, -36.0, -43.0, -50.0, -57.0, -64.0, -71.0, -78.0, -84.0, -91.0 },  
-	{26.0, 19.0, 12.0, 5.0, -2.0, -9.0, -16.0, -23.0, -30.0, -37.0, -44.0, -51.0, -58.0, -65.0, -72.0, -79.0, -86.0, -93.0 },
-	{26.0, 19.0, 12.0, 4.0, -3.0, -10.0, -17.0, -24.0, -31.0, -38.0, -45.0, -52.0, -60.0, -67.0, -74.0, -81.0, -88.0, -95.0},
-	{25.0, 18.0, 11.0, 4.0, -3.0, -11.0, -18.0, -25.0, -32.0, -39.0, -46.0, -54.0, -61.0, -68.0, -75.0, -82.0, -89.0, -97.0},  
-	{25.0, 17.0, 10.0, 3.0, -4.0, -11.0, -19.0, -26.0, -33.0, -40.0, -48.0, -55.0, -62.0, -69.0, -76.0, -84.0, -91.0, -98.0} 
-    };
-    float windchillFfloatread[ROWS][COLS];
-    float polyflres[ROWS][COLS];
-    int   polyintres[ROWS][COLS];
-    float polyflread[ROWS][COLS];
-    int   polyintread[ROWS][COLS];
-    int windchillCint[ROWS][COLS];
+
+int compare_int(const int* a, const int* b, int tol);
+int compare_float(const float* a,const  float* b, double tol);
+int compare_char(const char* a,const  float* b, int tol);
+int compare_double(const double* a, const  float* b, int tol);
+int compare_short(const short* a, const  float* b, int tol);
+
+
+int main(void)
+{
+    hid_t dxpl_id_c_to_f, dxpl_id_c_to_f_copy, dxpl_id_simple, dxpl_id_polynomial, dxpl_id_polynomial_copy;
+
     herr_t err;
-    int row, col;
 
-    TESTING("data transform, setting up")
+    /* Run all the tests */
+    if((err = init_test(&dxpl_id_c_to_f, &dxpl_id_simple, &dxpl_id_polynomial, &dxpl_id_polynomial_copy, &dxpl_id_c_to_f_copy)) < 0) TEST_ERROR;
+ 
+    if((err = test_int(dxpl_id_c_to_f, dxpl_id_polynomial)) < 0) TEST_ERROR;
+ 
+    if((err = test_float(dxpl_id_c_to_f, dxpl_id_polynomial)) < 0) TEST_ERROR;
+  
+    if((err = test_char(dxpl_id_c_to_f)) < 0) TEST_ERROR;
+
+    if((err = test_short(dxpl_id_c_to_f)) < 0) TEST_ERROR;
+    if((err = test_double(dxpl_id_c_to_f)) < 0) TEST_ERROR;
+ 
+    if((err = test_copy(dxpl_id_c_to_f_copy, dxpl_id_polynomial_copy)) < 0) TEST_ERROR;
+ 
+    if((err = test_trivial(dxpl_id_simple)) < 0) TEST_ERROR;
 
-    for(row = 0; row<ROWS; row++)
-    {
-	for(col = 0; col<COLS; col++)
-	    windchillCint[row][col] = (5/9.0)*(windchillFint[row][col] - 32);
-    }
+    /* Close the objects we opened/created */
+   if((err = H5Dclose(dset_id_int))<0) TEST_ERROR;
+   if((err = H5Dclose(dset_id_float))<0) TEST_ERROR;
+   if((err = H5Dclose(dset_id_char))<0) TEST_ERROR;
+   if((err = H5Dclose(dset_id_double))<0) TEST_ERROR;
+   if((err = H5Pclose(dxpl_id_c_to_f))<0) TEST_ERROR;
+   if((err = H5Pclose(dxpl_id_c_to_f_copy))<0) TEST_ERROR;
+   if((err = H5Pclose(dxpl_id_polynomial))<0) TEST_ERROR;
+   if((err = H5Pclose(dxpl_id_polynomial_copy))<0) TEST_ERROR;
+   if((err = H5Pclose(dxpl_id_simple))<0) TEST_ERROR;
 
-    for(row = 0; row<ROWS; row++)
-    {
-	for(col = 0; col<COLS; col++)
-	    polyflres[row][col] = (2.0+windchillCint[row][col])*((windchillCint[row][col]-8.0)/2.0);
-    }
 
-    for(row = 0; row<ROWS; row++)
-    {
-	for(col = 0; col<COLS; col++)
-	    polyintres[row][col] = (2+windchillCint[row][col])*((windchillCint[row][col]-8)/2);
-    }
+   return 0;
+
+error:
+   return -1;
+}
+
+int init_test(hid_t* dxpl_id_c_to_f, hid_t* dxpl_id_simple, hid_t* dxpl_id_polynomial, hid_t* dxpl_id_polynomial_copy, hid_t* dxpl_id_c_to_f_copy)
+{
+    const char* f_to_c = "(5/9.0)*(x-32)";
+    const char* c_to_f = "(9/5.0)*x + 32";
+    const char* polynomial = "(2+x)* ((x-8)/2)";
+    const char* simple = "(4/2) * ( (2 + 4)/(5 - 2.5))"; /* this equals 4.8 */ 
+
+    herr_t err;
+    hid_t file_id, datatype, dataspace, dxpl_id_f_to_c;
+    hsize_t dim[2] = {ROWS, COLS};
 
     if((dxpl_id_f_to_c = H5Pcreate(H5P_DATASET_XFER))<0) TEST_ERROR;
-    if((dxpl_id_c_to_f = H5Pcreate(H5P_DATASET_XFER))<0) TEST_ERROR;
-    if((dxpl_id_simple = H5Pcreate(H5P_DATASET_XFER))<0) TEST_ERROR;
-    if((dxpl_id_polynomial =  H5Pcreate(H5P_DATASET_XFER))<0) TEST_ERROR;
+    if((*dxpl_id_c_to_f = H5Pcreate(H5P_DATASET_XFER))<0) TEST_ERROR;
+    if((*dxpl_id_simple = H5Pcreate(H5P_DATASET_XFER))<0) TEST_ERROR;
+    if((*dxpl_id_polynomial =  H5Pcreate(H5P_DATASET_XFER))<0) TEST_ERROR;
     
     if((err= H5Pset_data_transform(dxpl_id_f_to_c, f_to_c))<0) TEST_ERROR;
-    if((err= H5Pset_data_transform(dxpl_id_c_to_f, c_to_f))<0) TEST_ERROR;
-    if((err = H5Pset_data_transform(dxpl_id_simple, simple))<0) TEST_ERROR;
-    if((err = H5Pset_data_transform(dxpl_id_polynomial, polynomial))<0) TEST_ERROR;
-    if((dxpl_id_polynomial_copy = H5Pcopy(dxpl_id_polynomial)) < 0) TEST_ERROR;
-    if((dxpl_id_c_to_f_copy = H5Pcopy(dxpl_id_c_to_f)) < 0) TEST_ERROR;
+    if((err= H5Pset_data_transform(*dxpl_id_c_to_f, c_to_f))<0) TEST_ERROR;
+    if((err = H5Pset_data_transform(*dxpl_id_simple, simple))<0) TEST_ERROR;
+    if((err = H5Pset_data_transform(*dxpl_id_polynomial, polynomial))<0) TEST_ERROR;
+    if((*dxpl_id_polynomial_copy = H5Pcopy(*dxpl_id_polynomial)) < 0) TEST_ERROR;
+    if((*dxpl_id_c_to_f_copy = H5Pcopy(*dxpl_id_c_to_f)) < 0) TEST_ERROR;
 
     if((file_id = H5Fcreate("dtransform.h5", H5F_ACC_TRUNC, H5P_DEFAULT, H5P_DEFAULT))<0) TEST_ERROR;
 
     if((dataspace = H5Screate_simple(2, dim, NULL))<0) TEST_ERROR;
 
-    PASSED();
-
-    TESTING("data transform, writing integer data")
+    TESTING("data transform, writing out data file")
 
     /* Write out the integer dataset to the file, converting it to c */
     if((datatype = H5Tcopy(H5T_NATIVE_INT))<0) TEST_ERROR;
@@ -98,9 +137,13 @@ int main(void)
     if((err = H5Dwrite(dset_id_int, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, dxpl_id_f_to_c, windchillFint))<0) TEST_ERROR;
     if((err = H5Tclose(datatype))<0) TEST_ERROR;
 
-    PASSED();
 
-    TESTING("data transform, writing float data")
+    /* Write out the character dataset to the file, converting it to c */
+    if((datatype = H5Tcopy(H5T_NATIVE_CHAR))<0) TEST_ERROR;
+    if((dset_id_char = H5Dcreate(file_id, "/transformtest_char", datatype, dataspace, H5P_DEFAULT))<0) TEST_ERROR;
+    if((err = H5Dwrite(dset_id_char, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, dxpl_id_f_to_c, windchillFint))<0) TEST_ERROR;
+    if((err = H5Tclose(datatype))<0) TEST_ERROR;
+
     
     /* Write out a floating point version to the file, converting it to c */
     if((datatype = H5Tcopy(H5T_NATIVE_FLOAT))<0) TEST_ERROR;
@@ -108,26 +151,155 @@ int main(void)
     if((err = H5Dwrite(dset_id_float, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, dxpl_id_f_to_c, windchillFint))<0) TEST_ERROR;
     if((err = H5Tclose(datatype))<0) TEST_ERROR;
 
+    if((datatype = H5Tcopy(H5T_NATIVE_UCHAR))<0) TEST_ERROR;
+    if((dset_id_uchar = H5Dcreate(file_id, "/transformtest_uchar", datatype, dataspace, H5P_DEFAULT))<0) TEST_ERROR;
+    if((err = H5Dwrite(dset_id_uchar, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, dxpl_id_f_to_c, windchillFint))<0) TEST_ERROR;
+    if((err = H5Tclose(datatype))<0) TEST_ERROR;
+
+
+    if((datatype = H5Tcopy(H5T_NATIVE_SCHAR))<0) TEST_ERROR;
+    if((dset_id_schar = H5Dcreate(file_id, "/transformtest_schar", datatype, dataspace, H5P_DEFAULT))<0) TEST_ERROR;
+    if((err = H5Dwrite(dset_id_schar, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, dxpl_id_f_to_c, windchillFint))<0) TEST_ERROR;
+    if((err = H5Tclose(datatype))<0) TEST_ERROR;
+
+
+    if((datatype = H5Tcopy(H5T_NATIVE_UINT))<0) TEST_ERROR;
+    if((dset_id_uint = H5Dcreate(file_id, "/transformtest_uint", datatype, dataspace, H5P_DEFAULT))<0) TEST_ERROR;
+    if((err = H5Dwrite(dset_id_uint, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, dxpl_id_f_to_c, windchillFint))<0) TEST_ERROR;
+    if((err = H5Tclose(datatype))<0) TEST_ERROR;
+
+    if((datatype = H5Tcopy(H5T_NATIVE_USHORT))<0) TEST_ERROR;
+    if((dset_id_ushort = H5Dcreate(file_id, "/transformtest_ushort", datatype, dataspace, H5P_DEFAULT))<0) TEST_ERROR;
+    if((err = H5Dwrite(dset_id_ushort, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, dxpl_id_f_to_c, windchillFint))<0) TEST_ERROR;
+    if((err = H5Tclose(datatype))<0) TEST_ERROR;
+
+    if((datatype = H5Tcopy(H5T_NATIVE_ULONG))<0) TEST_ERROR;
+    if((dset_id_ulong = H5Dcreate(file_id, "/transformtest_ulong", datatype, dataspace, H5P_DEFAULT))<0) TEST_ERROR;
+    if((err = H5Dwrite(dset_id_ulong, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, dxpl_id_f_to_c, windchillFint))<0) TEST_ERROR;
+    if((err = H5Tclose(datatype))<0) TEST_ERROR;
+
+    if((datatype = H5Tcopy(H5T_NATIVE_LONG))<0) TEST_ERROR;
+    if((dset_id_long = H5Dcreate(file_id, "/transformtest_long", datatype, dataspace, H5P_DEFAULT))<0) TEST_ERROR;
+    if((err = H5Dwrite(dset_id_long, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, dxpl_id_f_to_c, windchillFint))<0) TEST_ERROR;
+    if((err = H5Tclose(datatype))<0) TEST_ERROR;
+
+    if((datatype = H5Tcopy(H5T_NATIVE_SHORT))<0) TEST_ERROR;
+    if((dset_id_short = H5Dcreate(file_id, "/transformtest_short", datatype, dataspace, H5P_DEFAULT))<0) TEST_ERROR;
+    if((err = H5Dwrite(dset_id_short, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, dxpl_id_f_to_c, windchillFint))<0) TEST_ERROR;
+    if((err = H5Tclose(datatype))<0) TEST_ERROR;
+
+    if((datatype = H5Tcopy(H5T_NATIVE_DOUBLE))<0) TEST_ERROR;
+    if((dset_id_double = H5Dcreate(file_id, "/transformtest_double", datatype, dataspace, H5P_DEFAULT))<0) TEST_ERROR;
+    if((err = H5Dwrite(dset_id_double, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, dxpl_id_f_to_c, windchillFint))<0) TEST_ERROR;
+    if((err = H5Tclose(datatype))<0) TEST_ERROR;
+
     PASSED();
 
+
+
+   if((err = H5Fclose(file_id))<0) TEST_ERROR;
+   if((err = H5Pclose(dxpl_id_f_to_c))<0) TEST_ERROR;
+   if((err = H5Sclose(dataspace))<0) TEST_ERROR;
+
+   return 0;
+error:
+   return -1;
+}
+ 
+int test_int(const hid_t dxpl_id_c_to_f, const hid_t dxpl_id_polynomial)
+{
+    int windchillFintread[ROWS][COLS];
+    int   polyintres[ROWS][COLS];
+    int   polyintread[ROWS][COLS];
+    int windchillCint[ROWS][COLS];
+
+    herr_t err;
+
+    int row, col;
+
+    for(row = 0; row<ROWS; row++)
+    {
+	for(col = 0; col<COLS; col++)
+	    windchillCint[row][col] = (5/9.0)*(windchillFint[row][col] - 32);
+    }
+
+    for(row = 0; row<ROWS; row++)
+    {
+	for(col = 0; col<COLS; col++)
+	    polyintres[row][col] = (2+windchillCint[row][col])*((windchillCint[row][col]-8)/2);
+    }
+
+  
+
     TESTING("data transform, no data type conversion (int->int)")
     /* Read in the integer data with a data transform. */
     if((err = H5Dread(dset_id_int, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, dxpl_id_c_to_f, windchillFintread))<0) TEST_ERROR;
     
     if( (compare_int(*windchillFintread, *windchillFint, 2)) == 0)
     {
-        fprintf(stderr, "ERROR: Conversion failed to match computed data\n");
-        TEST_ERROR;
+	H5_FAILED();
+        fprintf(stderr, "    ERROR: Conversion failed to match computed data\n");
     }
     else
 	PASSED();
 
-    TESTING("data transform, no data type conversion (float->float)")
+   TESTING("data transform, with data type conversion (float->int)")
+    if((err = H5Dread(dset_id_float, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, dxpl_id_c_to_f, windchillFintread))<0) TEST_ERROR;
+    if( (compare_int(*windchillFint, *windchillFintread, 2)) == 0)
+         {   
+	     H5_FAILED();
+             fprintf(stderr, "    ERROR: Conversion failed to match computed data\n");
+    }
+    else
+       PASSED();
+
+    TESTING("data transform, polynomial transform (float->int)")
+    if((err = H5Dread(dset_id_float, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, dxpl_id_polynomial, polyintread))<0) TEST_ERROR;
+    if( (compare_int(*polyintres, *polyintread, 2)) == 0)
+    {
+	H5_FAILED();	
+	fprintf(stderr, "    ERROR: Conversion failed to match computed data\n");
+    }
+    else
+       PASSED();
+
+    return 0;
+
+error:
+    return -1;
+
+}
+
+int test_float(const hid_t dxpl_id_c_to_f, const hid_t dxpl_id_polynomial)
+{
+    float polyflres[ROWS][COLS];
+    float windchillFfloatread[ROWS][COLS];
+    float polyflread[ROWS][COLS];
+    int windchillCint[ROWS][COLS];
+    
+    herr_t err;
+    int row, col;
+
+    for(row = 0; row<ROWS; row++)
+    {
+	for(col = 0; col<COLS; col++)
+	    windchillCint[row][col] = (5/9.0)*(windchillFint[row][col] - 32);
+    }
+   
+    for(row = 0; row<ROWS; row++)
+    {
+	for(col = 0; col<COLS; col++)
+	    polyflres[row][col] = (2.0+windchillCint[row][col])*((windchillCint[row][col]-8.0)/2.0);
+    }
+
+
+
+   TESTING("data transform, no data type conversion (float->float)")
     if((err = H5Dread(dset_id_float, H5T_NATIVE_FLOAT, H5S_ALL, H5S_ALL, dxpl_id_c_to_f, windchillFfloatread))<0) TEST_ERROR;
     if( (compare_float(*windchillFfloat, *windchillFfloatread, 2.0)) == 0)
     {
-	fprintf(stderr, "ERROR: Conversion failed to match computed data\n");
-	TEST_ERROR;
+	H5_FAILED();	
+	fprintf(stderr, "    ERROR: Conversion failed to match computed data\n");
     }
     else
       PASSED();
@@ -136,65 +308,156 @@ int main(void)
     if((err = H5Dread(dset_id_int, H5T_NATIVE_FLOAT, H5S_ALL, H5S_ALL, dxpl_id_c_to_f, windchillFfloatread))<0) TEST_ERROR;
     if( (compare_float(*windchillFfloat, *windchillFfloatread, 2.0)) == 0)
          {   
-             fprintf(stderr, "ERROR: Conversion failed to match computed data\n");
-        TEST_ERROR;
+	H5_FAILED();	
+             fprintf(stderr, "    ERROR: Conversion failed to match computed data\n");
     }
     else
        PASSED();
 
-   TESTING("data transform, with data type conversion (float->int)")
-    if((err = H5Dread(dset_id_float, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, dxpl_id_c_to_f, windchillFintread))<0) TEST_ERROR;
-    if( (compare_int(*windchillFint, *windchillFintread, 2)) == 0)
-         {   
-             fprintf(stderr, "ERROR: Conversion failed to match computed data\n");
-        TEST_ERROR;
+    TESTING("data transform, polynomial transform (int-float)")
+    if((err = H5Dread(dset_id_int, H5T_NATIVE_FLOAT, H5S_ALL, H5S_ALL, dxpl_id_polynomial, polyflread))<0) TEST_ERROR;
+    if( (compare_float(*polyflres, *polyflread, 2.0)) == 0)
+    {   
+	H5_FAILED();	
+	fprintf(stderr, "    ERROR: Conversion failed to match computed data\n");
     }
     else
        PASSED();
 
-    TESTING("data transform, linear transform w/ copied property")
-    if((err = H5Dread(dset_id_float, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, dxpl_id_c_to_f_copy, windchillFintread))<0) TEST_ERROR;
-    if( (compare_int(*windchillFint, *windchillFintread, 2)) == 0)
+    return 0;
+
+error:
+    return -1;
+
+}
+
+
+int test_char(const hid_t dxpl_id_c_to_f)
+{
+    char  windchillFchar[ROWS][COLS];
+
+    herr_t err;
+    
+    TESTING("data transform, no data type conversion (char->char)")
+    if((err = H5Dread(dset_id_char, H5T_NATIVE_CHAR, H5S_ALL, H5S_ALL, dxpl_id_c_to_f, windchillFchar))<0) TEST_ERROR;
+    if( (compare_char(*windchillFchar, *windchillFfloat, 2)) == 0)
     {   
-	fprintf(stderr, "ERROR: Conversion failed to match computed data\n");
-        TEST_ERROR;
+	H5_FAILED();
+	fprintf(stderr, "    ERROR: Conversion failed to match computed data\n");
     }
     else
        PASSED();
 
-    TESTING("data transform, polynomial transform (float->int)")
-    if((err = H5Dread(dset_id_float, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, dxpl_id_polynomial, polyintread))<0) TEST_ERROR;
-    if( (compare_int(*polyintres, *polyintread, 2)) == 0)
+    return 0;
+
+error:
+    return -1;
+}
+
+int test_double(const hid_t dxpl_id_c_to_f)
+{
+    double windchillFdouble[ROWS][COLS];
+
+    herr_t err;
+    
+    TESTING("data transform, no data type conversion (double->double)")
+    if((err = H5Dread(dset_id_double, H5T_NATIVE_DOUBLE, H5S_ALL, H5S_ALL, dxpl_id_c_to_f, windchillFdouble))<0) TEST_ERROR;
+    if( (compare_double(*windchillFdouble, *windchillFfloat, 2)) == 0)
     {   
-	fprintf(stderr, "ERROR: Conversion failed to match computed data\n");
-        TEST_ERROR;
+	H5_FAILED();
+	fprintf(stderr, "    ERROR: Conversion failed to match computed data\n");
     }
     else
        PASSED();
 
-    TESTING("data transform, polynomial transform w/ copied property")
-    if((err = H5Dread(dset_id_float, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, dxpl_id_polynomial_copy, polyintread))<0) TEST_ERROR;
-    if( (compare_int(*polyintres, *polyintread, 2)) == 0)
+    return 0;
+
+error:
+    return -1;
+}
+
+
+int test_short(const hid_t dxpl_id_c_to_f)
+{
+    short windchillFshort[ROWS][COLS];
+
+    herr_t err;
+    
+    TESTING("data transform, no data type conversion (short->short)")
+    if((err = H5Dread(dset_id_short, H5T_NATIVE_SHORT, H5S_ALL, H5S_ALL, dxpl_id_c_to_f, windchillFshort))<0) TEST_ERROR;
+    if( (compare_short(*windchillFshort, *windchillFfloat, 2)) == 0)
     {   
-	fprintf(stderr, "ERROR: Conversion failed to match computed data\n");
-        TEST_ERROR;
+	H5_FAILED();
+	fprintf(stderr, "    ERROR: Conversion failed to match computed data\n");
     }
     else
        PASSED();
 
+    return 0;
 
-    TESTING("data transform, polynomial transform (int-float)")
-    if((err = H5Dread(dset_id_int, H5T_NATIVE_FLOAT, H5S_ALL, H5S_ALL, dxpl_id_polynomial, polyflread))<0) TEST_ERROR;
-    if( (compare_float(*polyflres, *polyflread, 2.0)) == 0)
+error:
+    return -1;
+}
+
+
+int test_copy(const hid_t dxpl_id_c_to_f_copy, const hid_t dxpl_id_polynomial_copy)
+{
+    int windchillCint[ROWS][COLS];
+    int polyintres[ROWS][COLS];
+    int polyintread[ROWS][COLS];
+    int windchillFintread[ROWS][COLS];
+
+    herr_t err;
+
+    int row, col;
+
+    for(row = 0; row<ROWS; row++)
+    {
+	for(col = 0; col<COLS; col++)
+	    windchillCint[row][col] = (5/9.0)*(windchillFint[row][col] - 32);
+    }
+
+    for(row = 0; row<ROWS; row++)
+    {
+	for(col = 0; col<COLS; col++)
+	    polyintres[row][col] = (2+windchillCint[row][col])*((windchillCint[row][col]-8)/2);
+    }
+
+  
+   TESTING("data transform, linear transform w/ copied property")
+    if((err = H5Dread(dset_id_float, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, dxpl_id_c_to_f_copy, windchillFintread))<0) TEST_ERROR;
+    if( (compare_int(*windchillFint, *windchillFintread, 2)) == 0)
     {   
-	fprintf(stderr, "ERROR: Conversion failed to match computed data\n");
-        TEST_ERROR;
+	H5_FAILED();	
+	fprintf(stderr, "    ERROR: Conversion failed to match computed data\n");
     }
     else
        PASSED();
 
+    TESTING("data transform, polynomial transform w/ copied property")
+    if((err = H5Dread(dset_id_float, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, dxpl_id_polynomial_copy, polyintread))<0) TEST_ERROR;
+    if( (compare_int(*polyintres, *polyintread, 2)) == 0)
+    {   
+	H5_FAILED();	
+	fprintf(stderr, "    ERROR: Conversion failed to match computed data\n");
+    }
+    else
+       PASSED();
 
-    
+    return 0;
+
+error:
+    return -1;
+}
+   
+int test_trivial(const hid_t dxpl_id_simple)
+{
+    float windchillFfloatread[ROWS][COLS];
+    int windchillFintread[ROWS][COLS];
+
+    herr_t err;
+    int row, col;
+ 
     TESTING("data transform, trivial transform, without type conversion")
     if((err = H5Dread(dset_id_float, H5T_NATIVE_FLOAT, H5S_ALL, H5S_ALL, dxpl_id_simple, windchillFfloatread))<0) TEST_ERROR;
     for(row = 0; row<ROWS; row++)
@@ -203,8 +466,8 @@ int main(void)
 	{
 	    if((windchillFfloatread[row][col] - 4.8) > 0.0001)
 	    {
-		fprintf(stderr, "ERROR: Conversion failed to match computed data\n");
-		TEST_ERROR;
+		H5_FAILED();
+		fprintf(stderr, "    ERROR: Conversion failed to match computed data\n");
 	    }
 	}
     }
@@ -218,29 +481,19 @@ int main(void)
 	{
 	    if(windchillFintread[row][col] != 4)
 	    {
-		fprintf(stderr, "ERROR: Conversion failed to match computed data\n");
-		TEST_ERROR;
+		H5_FAILED();
+		fprintf(stderr, "    ERROR: Conversion failed to match computed data\n");
 	    }
 	}
     }
     PASSED();
 
-
-    /* Close the objects we opened/created */
-   if((err = H5Sclose(dataspace))<0) TEST_ERROR;
-   if((err = H5Dclose(dset_id_int))<0) TEST_ERROR;
-   if((err = H5Dclose(dset_id_float))<0) TEST_ERROR;
-   if((err = H5Fclose(file_id))<0) TEST_ERROR;
-   if((err = H5Pclose(dxpl_id_f_to_c))<0) TEST_ERROR;
-   if((err = H5Pclose(dxpl_id_c_to_f))<0) TEST_ERROR;
-
-   return 0;
-
+    return 0;
 error:
-   return -1;
+    return -1;
 }
 
-int compare_int(int* a, int* b, int tol)
+int compare_int(const int* a, const int* b, int tol)
 {
     int i;
     for(i=0; i<ROWS*COLS; i++)
@@ -252,7 +505,7 @@ int compare_int(int* a, int* b, int tol)
    
 }
     
-int compare_float(float* a, float* b, double tol)
+int compare_float(const float* a,const  float* b, double tol)
 {
     int i;
     for(i=0; i<ROWS*COLS; i++)
@@ -263,4 +516,39 @@ int compare_float(float* a, float* b, double tol)
     return 1;
 }
 
+ 
+int compare_char(const char* a,const float* b, int tol)
+{
+    int i;
+    for(i=0; i<ROWS*COLS; i++)
+    {
+	if( !((a[i] <= ((char)b[i] + tol)) && (a[i] >= (char)(b[i] - tol))) )
+	    return 0;
+    }
+    return 1;
+}
+
+int compare_double(const double* a, const float* b, int tol)
+{
+    int i;
+    for(i=0; i<ROWS*COLS; i++)
+    {
+	if( !((a[i] <= ((double)b[i] + tol)) && (a[i] >= (double)(b[i] - tol))) )
+	    return 0;
+    }
+    return 1;
+}
+
+
+int compare_short(const short* a, const float* b, int tol)
+{
+    int i;
+    for(i=0; i<ROWS*COLS; i++)
+    {
+	if( !((a[i] <= ((short)b[i] + tol)) && (a[i] >= (short)(b[i] - tol))) )
+	    return 0;
+    }
+    return 1;
+}
+
 
-- 
cgit v0.12