/*
 * Aaron Weeden, Shodor, May 2014
 */
#include <hdf5.h>

const int N_ROWS = 3;
const int N_COLS = 4;

int flatten(int row, int col) {
    return row * N_COLS + col;
}

int main() {

    hid_t file_id, dataset_id, dataspace_id;
    hsize_t dims[2] = {N_ROWS, N_COLS};
    int row, col, data[N_ROWS][N_COLS];

    for (row=0; row<N_ROWS; row++) {
        for (col=0; col<N_COLS; col++) {
            data[row][col] = flatten(row,col);
        }
    }

    file_id = H5Fcreate("sample-dataset.h5", H5F_ACC_TRUNC, H5P_DEFAULT,
            H5P_DEFAULT);

    dataspace_id = H5Screate_simple(2, dims, NULL);

    dataset_id = H5Dcreate2(file_id, "/sample-dataset", H5T_STD_I32BE,
            dataspace_id, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);

    H5Dwrite(dataset_id, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, data);
    H5Dread(dataset_id, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, data);

    H5Dclose(dataset_id);
    H5Sclose(dataspace_id);
    H5Fclose(file_id);

    return 0;
}