/*  
 *  h5float.c  - Read in floating point values from text
 *               file, and write them to an HDF5 file.
 */
 
#include<stdio.h>
#include <hdf5.h>

#define DATASETNAME "FloatData" 
#define NX     7                      /* dataset dimensions */
#define NY     2 
#define RANK   2 

int
main (void)
{
    hid_t       file, dataset;         /* file and dataset handles */
    hid_t       datatype, dataspace;   /* handles */
    hsize_t     dimsf[2];              /* dataset dimensions */
    herr_t      status;                             
    float       data[NX][NY];          /* data to write */
    float       data1[NX][NY];          /* data to write */
    int         i, j, num;
    float     val;
    FILE *fp;

    fp = fopen ("float.txt", "r");

    for (i=0; i<NX; i++)
        for (j=0; j<NY; j++)
        {
           num = fscanf (fp, "%f", &val);
           data[i][j] = val;
        }

    /*
     * Create a new file using H5F_ACC_TRUNC access,
     * default file creation properties, and default file
     * access properties.
     */

    file = H5Fcreate("float.h5", H5F_ACC_TRUNC, H5P_DEFAULT, H5P_DEFAULT);


    /*
     * Describe the size of the array and create the data space for fixed
     * size dataset. 
     */
    dimsf[0] = NX;
    dimsf[1] = NY;
    dataspace = H5Screate_simple(RANK, dimsf, NULL); 

    dataset = H5Dcreate(file, DATASETNAME, H5T_IEEE_F32BE, dataspace,
			H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);

    /*
     * Write the data to the dataset using default transfer properties.
     */
    status = H5Dwrite(dataset, H5T_NATIVE_FLOAT, H5S_ALL, H5S_ALL,
		      H5P_DEFAULT, data);

    status = H5Dread(dataset, H5T_NATIVE_FLOAT, H5S_ALL, H5S_ALL, 
             H5P_DEFAULT, data1);

    printf ("\n");
    for (i=0; i<NX; i++)
    {
        for (j=0; j<NY; j++)
        { 
         printf (" %f", data1[i][j]);
        }
        printf ("\n");
   }

    /*
     * Close/release resources.
     */
    H5Sclose(dataspace);
    H5Dclose(dataset);
    H5Fclose(file);
 
    return 0;
}