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

#define CHUNK_SIZE 1024
#define SIZE 1024*1024*CHUNK_SIZE

int main() {
   hid_t file, dataset, dcpl, memspace, dataspace;
   hsize_t size = SIZE;
   hsize_t chunk_size = CHUNK_SIZE;
   herr_t status;

   float t[CHUNK_SIZE];
   for(int j=0; j<chunk_size; j++)
     t[j] = (float)rand()/(float)(RAND_MAX/200000.) - 100000.;

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

   dataspace = H5Screate_simple(1, &size, NULL);
   dcpl = H5Pcreate(H5P_DATASET_CREATE);
   status = H5Pset_chunk(dcpl, 1, &chunk_size);
   status = H5Pset_scaleoffset(dcpl, H5Z_SO_FLOAT_DSCALE, 5);

   dataset = H5Dcreate2(file, "data", H5T_INTEL_F32, dataspace,
			H5P_DEFAULT, dcpl, H5P_DEFAULT);

   status = H5Sclose(dataspace);
   status = H5Pclose(dcpl);

   memspace = H5Screate_simple(1, &chunk_size, NULL);
   dataspace = H5Dget_space(dataset);
   for(int i=0; i<size/chunk_size; i++)
     {
       hsize_t offset = i*chunk_size;
       status = H5Sselect_hyperslab(dataspace, H5S_SELECT_SET, &offset, NULL, &chunk_size, NULL);
       status = H5Dwrite(dataset, H5T_NATIVE_FLOAT, memspace, dataspace, H5P_DEFAULT, t);
       status = H5Fflush(file, H5F_SCOPE_GLOBAL);
     }
   status = H5Sclose(dataspace);
   status = H5Sclose(memspace);

   status = H5Dclose(dataset);
   status = H5Fclose(file);
}