src/hdf5/H5ecloud/H5PartTest.cc

Go to the documentation of this file.

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <hdf5.h>
#include "H5Part.hh"

/*
  A simple regression test that shows how you use this API
  to write and read multi-timestep files of particle data.
*/
#ifdef PARALLEL_IO

int main(int argc,char *argv[]){
  int sz=5;
  double *x,*y,*z;
  long long *id;
  char name[64];
  H5PartFile *file;
  int i,t,nt,nds;
  int nprocs,myproc;
  hid_t gid;
  MPI_Comm comm=MPI_COMM_WORLD;

  MPI_Init(&argc,&argv);
  MPI_Comm_size(comm,&nprocs);
  MPI_Comm_rank(comm,&myproc);

  x=(double*)malloc(sz*nprocs*sizeof(double));
  y=(double*)malloc(sz*nprocs*sizeof(double));
  z=(double*)malloc(sz*nprocs*sizeof(double));
  id=(long long*)malloc(sz*nprocs*sizeof(long long));
  /* parallel file creation */
  file=H5PartOpenFileParallel("parttest.h5",H5PART_WRITE,comm);
  if(!file) {
    perror("File open failed:  exiting!");
    exit(0);
  }

  for(t=0;t<5;t++){
    MPI_Barrier(comm);
    for(i=0;i<sz;i++) {
      x[i]=(double)(i+t)+10.0*(double)myproc;
      y[i]=0.1 + (double)(i+t);
      z[i]=0.2 + (double)(i+t*10);
      id[i]=i+sz*myproc;
    }
    printf("Proc[%u] Writing timestep %u file=%u\n",myproc,t,file->file);
    H5PartSetStep(file,t); /* must set the current timestep in file */
    H5PartSetNumParticles(file,sz); /* then set number of particles to store */
    /* now write different tuples of data into this timestep of the file */
    H5PartWriteDataFloat64(file,"x",x); 
    H5PartWriteDataFloat64(file,"y",y);
    H5PartWriteDataFloat64(file,"z",z);

    H5PartWriteDataFloat64(file,"px",x); 
    H5PartWriteDataFloat64(file,"py",y);
    H5PartWriteDataFloat64(file,"pz",z);

    H5PartWriteDataInt64(file,"id",id);
  }

  unsigned int idStart = 0+sz*myproc;
  unsigned int idEnd = (sz-1)+sz*myproc;

  printf("AllDone p[%u]\n",myproc);
  H5PartCloseFile(file);
  fprintf(stderr,"Closed files p[%u]\n",myproc);
  MPI_Barrier(comm);
  
  fprintf(stderr,"p[%u:%u] : OK, close file and reopen for reading idStart %u  idEnd %u \n",myproc,nprocs,idStart,idEnd);

  file=H5PartOpenFileParallel("parttest.h5",H5PART_READ,comm);
  H5PartSetStep(file,0);
 // unsigned int np = 0;
  unsigned int np = (int)H5PartGetNumParticles(file);
  nt=H5PartGetNumSteps(file); /* get number of steps in file */
  nds=H5PartGetNumDatasets(file); /* get number of datasets in timestep 0 */

  MPI_Barrier(comm);
  
  H5PartSetView(file,idStart,idEnd);

  np = (int)H5PartGetNumParticles(file);
  printf("After SetView(%d,%d): steps= %u  datasets= %u particles= %u\n",
        (int)idStart,(int)idEnd,
        nt,nds,np);

  if(x) 
    free(x); 
  if(y) 
    free(y);
  if(z) 
    free(z);
  if(id) 
    free(id);

  H5PartCloseFile(file);
  MPI_Barrier(comm);
  fprintf(stderr,"proc[%u]:  done\n",myproc);
  return MPI_Finalize();
}

#else
int main(int argc,char *argv[]){
  int sz=10;
  double *x,*y,*z;
  long long *id;
  char name[64];
  H5PartFile *file;
  int i,t,nt,nds,np;
  hid_t gid;
  long long idStart = 0;
  long long idEnd = 0;


  x=(double*)malloc(sz*sizeof(double));
  y=(double*)malloc(sz*sizeof(double));
  z=(double*)malloc(sz*sizeof(double));
  id=(long long*)malloc(sz*sizeof(long long));
  /* parallel file creation */
  file=H5PartOpenFile("parttest.h5",H5PART_WRITE);
  if(!file) {
    perror("File open failed:  exiting!");
    exit(0);
  }

  for(t=0;t<5;t++){
    fprintf(stdout,"Writing timestep %u\n",t);
    for(i=0;i<sz;i++) {
      x[i]=(double)(i+t);
      y[i]=0.1 + (double)(i+t);
      z[i]=0.2 + (double)(i+t*10);
      id[i]=i;
      fprintf(stdout,"\tp[%u] x=%f y=%f z=%f id=%d\n",
              i,x[i],y[i],z[i],(int)id[i]);
    }
    H5PartSetStep(file,t); /* must set the current timestep in file */
    H5PartSetNumParticles(file,sz); /* then set number of particles to store */
    /* now write different tuples of data into this timestep of the file */
    H5PartWriteDataFloat64(file,"x",x); 
    H5PartWriteDataFloat64(file,"y",y);
    H5PartWriteDataFloat64(file,"z",z);

    H5PartWriteDataFloat64(file,"px",x); 
    H5PartWriteDataFloat64(file,"py",y);
    H5PartWriteDataFloat64(file,"pz",z);

    H5PartWriteDataInt64(file,"id",id);
  }

  printf("AllDone writing\n");
  H5PartCloseFile(file);
  
 
  /*+++++++++++++ Reopen File for Reading +++++++++++*/
  file=H5PartOpenFile("parttest.h5",H5PART_READ);

  
  /********************************************/
  H5PartSetStep(file,0);
  nt=H5PartGetNumSteps(file); /* get number of steps in file */
  nds=H5PartGetNumDatasets(file); /* get number of datasets in timestep 0 */
  np=H5PartGetNumParticles(file);
  

  fprintf(stdout,"OK, close file and reopen for reading\n");
  fprintf(stdout,"steps= %u\tdatasets=%u\tparticles= %u\n",
         nt,nds,np);
  
  // clear the particles
  for(i=0;i<np;i++){
    x[i]=y[i]=z[i]=0.0;
    id[i]=0;
  }

  H5PartReadDataFloat64(file,"x",x);
  H5PartReadDataFloat64(file,"y",y);
  H5PartReadDataFloat64(file,"z",z);
  H5PartReadDataInt64(file,"id",id);

  for(i=0;i<np;i++){
    fprintf(stdout,"\tp[%3u] x=%lf y=%lf z=%lf id=%lu\n",
            i,x[i],y[i],z[i],(unsigned)(id[i]));
  }
  /********************************************/
  printf("Set to last step and reload data\n");
  H5PartSetStep(file,nt-1);
  H5PartReadDataFloat64(file,"x",x);
  H5PartReadDataFloat64(file,"y",y);
  H5PartReadDataFloat64(file,"z",z);
  H5PartReadDataInt64(file,"id",id);
  for(i=0;i<np;i++){
    fprintf(stdout,"\tp[%3u] x=%lf y=%lf z=%lf id=%lu\n",
            i,x[i],y[i],z[i],(unsigned)(id[i]));
  }
  
  /********************************************/
  idEnd=np;
  printf("Old View is %d:%d\n",(int)idStart,(int)idEnd);
  H5PartSetView(file,idStart,idEnd>>1);
  printf("Set new view = %d:%d\n",(int)idStart,(int)(idEnd>>1));
  H5PartGetView(file,&idStart,&idEnd);
  np=H5PartGetNumParticles(file);
  printf("steps= %u  datasets= %u particles= %d with view %d:%d\n",
         nt,nds,(int)np,(int)idStart,(int)idEnd);
  H5PartSetStep(file,nt-1); // set to last step
  printf("Setting to last step = %u\n",nt-1);
  for(i=0;i<10;i++){ x[i]=y[i]=z[i]=0.0; id[i]=0; } /* clear the arrays */
  H5PartReadDataFloat64(file,"x",x);
  H5PartReadDataFloat64(file,"y",y);
  H5PartReadDataFloat64(file,"z",z);
  H5PartReadDataInt64(file,"id",id);

  for(i=0;i<np;i++){
    fprintf(stdout,"\tp[%3u] x=%lf y=%lf z=%lf id=%lu\n",i,x[i],y[i],z[i],(unsigned)id[i]);
  }

  /********************************************/
  printf("Now set the view to the latter half of the data in step #%u\n",nt-1);
  H5PartResetView(file);
  H5PartGetView(file,&idStart,&idEnd);
  printf("Reset view = %d:%d\nSetting to %u:%u\n",
         (int)idStart,(int)idEnd,
         (int)idEnd>>1,(int)idEnd);
  H5PartSetView(file,(idEnd>>1),idEnd);
  np=H5PartGetNumParticles(file);
  printf("Now particles in selection are %d\n",np);
  printf("doubleCheck=%d\n",H5PartGetView(file,0,0));
  
  for(i=0;i<10;i++){ x[i]=y[i]=z[i]=0.0; id[i]=0; } /* clear the arrays */
 
  H5PartReadDataFloat64(file,"x",x);
  H5PartReadDataFloat64(file,"y",y);
  H5PartReadDataFloat64(file,"z",z);
  H5PartReadDataInt64(file,"id",id);
  for(i=0;i<np;i++){
    fprintf(stdout,"\tp[%3u] x=%lf y=%lf z=%lf id=%lu\n",i,x[i],y[i],z[i],(unsigned)id[i]);
  }
  if(x) 
    free(x); 
  if(y) 
    free(y);
  if(z) 
    free(z);
  if(id) 
    free(id);

  H5PartCloseFile(file);
  fprintf(stderr,"done\n");
}

#endif

---