vchoi_fork

/*

 * Example of using the parallel HDF5 library to access datasets.

 * Last revised: April 24, 2001.

 *

 * This program contains two parts.  In the first part, the mpi processes

 * collectively create a new parallel HDF5 file and create two fixed

 * dimension datasets in it.  Then each process writes a hyperslab into

 * each dataset in an independent mode.  All processes collectively

 * close the datasets and the file.

 * In the second part, the processes collectively open the created file

 * and the two datasets in it.  Then each process reads a hyperslab from

 * each dataset in an independent mode and prints them out.

 * All processes collectively close the datasets and the file.

 */

​

#include <assert.h>

#include <hdf5.h>

​

#ifdef H5_HAVE_PARALLEL

/* Temporary source code */

#define FAIL -1

/* temporary code end */

​

/* Define some handy debugging shorthands, routines, ... */

/* debugging tools */

#define MESG(x)\

    if (verbose) printf("%s\n", x);\

​

#define MPI_BANNER(mesg)\

    {printf("--------------------------------\n");\

    printf("Proc %d: ", mpi_rank); \

    printf("*** %s\n", mesg);\

    printf("--------------------------------\n");}

​

#define SYNC(comm)\

    {MPI_BANNER("doing a SYNC"); MPI_Barrier(comm); MPI_BANNER("SYNC DONE");}

/* End of Define some handy debugging shorthands, routines, ... */

​

/* Constants definitions */

/* 24 is a multiple of 2, 3, 4, 6, 8, 12.  Neat for parallel tests. */

#define SPACE1_DIM1 24

#define SPACE1_DIM2 24

#define SPACE1_RANK 2

#define DATASETNAME1    "Data1"

#define DATASETNAME2    "Data2"

#define DATASETNAME3    "Data3"

/* hyperslab layout styles */

#define BYROW       1   /* divide into slabs of rows */

#define BYCOL       2   /* divide into blocks of columns */

​

​