/*
##############################################################################
# 
# Copyright by The HDF Group.
# All rights reserved.
#
# This file is part of the hl_region High-Level HDF5 APIs.  The full copyright 
# notice, including terms governing use, modification, and redistribution, 
# is contained in the file COPYING, which can be found at the root of the
# source code distribution tree and in the documentation directory (doc/html/). 
# If you do not have access to this file, you may request a copy of 
# "the hl_region High-Level HDF5 APIs copyright and license statement" from 
# help@hdfgroup.org.
# 
##############################################################################
*/

/*
  This program shows how to create a dataset and region references associated
  with hyperslabs. It then creates a new dataset composed of data from the 
  data associated with a recursive loop over all the region references. 
  Main illustrative function: H5LRcreate_ref_to_all
*/

#include "hdf5.h"
#include "hdf5_hl.h"
#include "h5hl_region.h"

#define filename "ex_ref_to_all.h5"
#define dsetnamev "MY_DATA" /* dataset name */
#define dsetnamer "MY_REF" /* region references name */

#define DIM0  9   /* dataset dimensions */
#define DIM1  8
#define RANK  2   /* dataset rank */


int main(void)
{
  hid_t file_id;        /* file identifier */
  hid_t space_id_ref;   /* region reference dataspace identifier */
  hid_t dset_id_ref;    /* region reference dataset identifiers*/
  hsize_t dims[RANK] =  {DIM0,DIM1}; /* dataset dimensions */
  hsize_t dims_ref[1] = {2}; /* region reference dimensions */
  herr_t status;
  hdset_reg_ref_t ref[2]; /* region references */
  int data[DIM0][DIM1];   /* data */
  int i, j;
  const char *path[2];    /* paths to the data for the region references */
  hsize_t block_coord[8] ={ 0, 0, 2, 2, 3, 3, 5, 5}; /* hyperslab coordinates defining region references */
  
  /*********************************************************  
   This section writes data to the HDF5 file.  
  *********************************************************/ 
  
  path[0] ="/MY_DATA";
  path[1] ="/MY_DATA";
  /* 
   * Data initialization. 
   */
  printf("FULL 2D DATA:");
  for (i=0; i<DIM0; i++) {
    printf("\n[ ");
    for (j=0; j<DIM1; j++) {
      data[i][j] =  10*i+j;
      printf("%02d ", data[i][j]);
    }
    printf("]");
  }
  printf("\n");

  /* Open the file. */
  
  file_id = H5Fcreate(filename, H5F_ACC_TRUNC, H5P_DEFAULT, H5P_DEFAULT);

  /*
   * Create and write the dataset.
   */
  status = H5LTmake_dataset ( file_id, dsetnamev, RANK, dims, H5T_NATIVE_INT, data);
  
  /*
   * Create two region references with hyperslab coordinates: (0,0)-(2,2) and (3,3)-(5,5).
   */
  status = H5LRcreate_region_references(file_id, 2, path, block_coord, ref);
  
  /*
   * Write the region references to the file
   */
  space_id_ref = H5Screate_simple(1, dims_ref, NULL);
  dset_id_ref = H5Dcreate2(file_id, dsetnamer, H5T_STD_REF_DSETREG, space_id_ref, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
  status = H5Dwrite(dset_id_ref, H5T_STD_REF_DSETREG, H5S_ALL, H5S_ALL, H5P_DEFAULT,ref);
  
  /*
   * Close all objects.
   */
  status = H5Sclose(space_id_ref);
  status = H5Dclose(dset_id_ref);
  status = H5Fclose(file_id);
  
/*************************************************************  

  This section creates a data set from region references found
  recursively in the file starting at "/" for this example.

************************************************************/  
  
  /*
   * Reopen the file.
   */
  file_id = H5Fopen(filename, H5F_ACC_RDWR,  H5P_DEFAULT);
  
  /*
   * Create the "/MY_DATAc" dataset from regions of data.
   */
  
  status = H5LRcreate_ref_to_all(file_id, "/",
				    "/NEW_DATA", H5_INDEX_NAME, H5_ITER_INC, H5R_DATASET_REGION);
  status = H5Fclose(file_id);

  return 0;
}