/* Compile with HDF5 1.8.3 and -DH5_USE_16_API */

#include "hdf5.h"
#include <stdio.h>

void print_ints(int x, int y, int* buf){

    int i, j;
    for(i=0;i<x;i++){
        for(j=0;j<y;j++){
            fprintf(stdout, "%8d ", buf[i*y+j]);
        }
        fprintf(stdout,"\n");
    }
}

void wipe_buf(int n, int* buf, int val){
    int i;
    for(i=0;i<n;i++){
        buf[i] = val;
    }
}

static int buf[16*16];

int main(){

    hid_t fid, dset, space, dcpl;
    hsize_t chunks[] = {2,2};
    hsize_t size_small[] = {4,4};
    hsize_t size_big[] = {16,16};
    hsize_t limits[] = {H5S_UNLIMITED, H5S_UNLIMITED};
    /* int fill_value = 41; */

    fid = H5Fcreate("test.hdf5", H5F_ACC_TRUNC, H5P_DEFAULT, H5P_DEFAULT);
    
    dcpl = H5Pcreate(H5P_DATASET_CREATE);
    H5Pset_chunk(dcpl, 2, chunks);
    /* H5Pset_fill_value(dcpl, H5T_NATIVE_INT, &fill_value); */

    space = H5Screate_simple(2, size_small, limits);

    /* Initialize 4x4 dataset */
    dset = H5Dcreate(fid, "ds", H5T_NATIVE_INT, space, dcpl);

    wipe_buf(16*16,buf, 1);
    H5Dwrite(dset, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, buf);
    
    wipe_buf(16*16,buf, 2);
    H5Dread(dset, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, buf);

    fprintf(stdout, "Initial contents (4 x 4):\n");
    print_ints(4,4,buf);    
   
    /* Expand to 16 x 16 */
    H5Dset_extent(dset, size_big);

    wipe_buf(16*16,buf, 3);
    H5Dread(dset, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, buf);
   
    fprintf(stdout, "Contents after expansion (16 x 16):\n");
    print_ints(16,16,buf);

    /* Write a value */
    wipe_buf(16*16,buf, 4);
    H5Dwrite(dset, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, buf);

    wipe_buf(16*16,buf, 5);
    H5Dread(dset, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, buf);

    fprintf(stdout, "Contents after write (16 x 16):\n");
    print_ints(16,16,buf);

    /* Resize it down to original size */
    H5Dset_extent(dset,size_small);

    wipe_buf(16*16,buf, 6);
    H5Dread(dset, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, buf);

    fprintf(stdout, "Contents after shrinking (4 x 4):\n");
    print_ints(4,4,buf);

    /* Re-expand to 16 x 16.  The fill value is not correctly applied. */
    H5Dset_extent(dset,size_big);

    wipe_buf(16*16,buf, 7);
    H5Dread(dset, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, buf);

    fprintf(stdout, "Contents after resize (16 x 16):\n");
    print_ints(16,16,buf);   

    return(0);
}