src/hdf5/H5ecloud/ORG/H5PartF.c

Go to the documentation of this file.

#include "H5Part.h"
#include <hdf5.h>
/*
Determine underscores using perl and then send through one stage of preprocessing
to get the F77NAME fixed up.  Otherwise, must include correct flag each time.

Or just don't declare F77NAME in the header (only in source as symbols).


All pointers are cast to haddr_t (a unint64_t) within the C/C++ code
In fortran, these pointers (handles) are carried as INTEGER*8

*/

#ifdef F77_SINGLE_UNDERSCORE
#define F77NAME(a,b,c) a
#elif defined(F77_NO_UNDERSCORE)
#define F77NAME(a,b,c) b
#elif defined(F77_CRAY_UNDERSCORE)
#define F77NAME(a,b,c) c
#else
#error Error, no way to determine how to construct fortran bindings
#endif

#define f_h5prt_openr F77NAME(h5prt_openr_,h5prt_openr,H5PRT_OPENR)
#define f_h5prt_openw F77NAME(h5prt_openw_,h5prt_openw,H5PRT_OPENW)
#define f_h5prt_close F77NAME(h5prt_close_,h5prt_close,H5PRT_CLOSE)
#define f_h5prt_nprt F77NAME(h5prt_nprt_,h5prt_nprt,H5PRT_NPRT)
#define f_h5prt_ndata F77NAME(h5prt_ndata_,h5prt_ndata,H5PRT_NDATA)
#define f_h5prt_nstep F77NAME(h5prt_nstep_,h5prt_nstep,H5PRT_NSTEP)
#define f_h5prt_readstep F77NAME(h5prt_readstep_,h5prt_readstep,H5PRT_READSTEP)
#define f_h5prt_step F77NAME(h5prt_step_,h5prt_step,H5PRT_STEP)
#define f_h5prt_wreal8 F77NAME(h5prt_wreal8_,h5prt_wreal8,H5PRT_WREAL8)
#define f_h5prt_wint8 F77NAME(h5prt_wint8_,h5prt_wint8,H5PRT_WINT8)

haddr_t f_h5prt_openr (char *file,int flen){
  haddr_t fh;
  H5PartFile* f;
  file[flen]='\0'; /* null terminate */
  f = (H5PartOpenFile(file,H5PART_READ));
  /*  printf("openr file=[%s] flen=%u haddr=%u\n",file,flen,f); */
  fh = (haddr_t)f;
  /* printf("FileHandle=%llu\n",fh); */
  return fh;
}

haddr_t f_h5prt_openw (char *file,int flen){
  haddr_t f;
  file[flen]='\0'; /* null terminate */
  f=(haddr_t)(H5PartOpenFile(file,H5PART_WRITE));
  return f;
}

void f_h5prt_close(haddr_t *file){
  H5PartCloseFile((H5PartFile*)*file);
}

haddr_t f_h5prt_nprt(haddr_t *file){
  /*  printf("nprt Haddr=%llu\n",file); */
  return H5PartGetNumParticles((H5PartFile*)*file);
}

int f_h5prt_ndata(haddr_t *file){
  /* printf("ndata Haddr=%llu\n",file); */
  return H5PartGetNumDatasets((H5PartFile*)*file);
}

int f_h5prt_nstep(haddr_t *file){
  /* printf("nstep Haddr=%llu\n",file); */
   return H5PartGetNumSteps((H5PartFile*)*file);
}

void f_h5prt_readstep(haddr_t *file,int *step,
                      double *x,double *y,double *z,
                      double *px,double *py,double *pz,
                      void *id){
  H5PartReadParticleStep((H5PartFile*)*file,(*step)-1,x,y,z,px,py,pz,(long long*)id);
}

void f_h5prt_step(haddr_t *file,int *step){
  H5PartSetStep((H5PartFile*)*file,(*step)-1);
}

void f_h5prt_wreal8(haddr_t *file,
                    char *name,void *data,int flen){
  name[flen]='\0';
  H5PartWriteDataFloat64((H5PartFile*)*file,name,(double*)data);
}

void f_h5prt_wint8(haddr_t *file,
                    char *name,void *data,int flen){
  name[flen]='\0';
  H5PartWriteDataInt64((H5PartFile*)*file,name,(long long*)data);
}

---