src/hdf5/H5ecloud/H5PartF.c

Go to the documentation of this file.

#include "H5Part.h"
#include "Underscore.h"
#include <hdf5.h>
/*

********* Perhaps call this pfile **********

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

#if 0
#define f_h5pt_openr F77NAME(h5prt_openr_,h5prt_openr,H5PRT_OPENR)
#define f_h5pt_openw F77NAME(h5prt_openw_,h5prt_openw,H5PRT_OPENW)
#define f_h5pt_close F77NAME(h5prt_close_,h5prt_close,H5PRT_CLOSE)
#define f_h5pt_npoints F77NAME(h5prt_npoints_,h5prt_npoints,H5PRT_NPOINTS)
#define f_h5pt_ndata F77NAME(h5prt_ndata_,h5prt_ndata,H5PRT_NDATA)
#define f_h5pt_nstep F77NAME(h5prt_nstep_,h5prt_nstep,H5PRT_NSTEP)
#define f_h5pt_readstep F77NAME(h5prt_readstep_,h5prt_readstep,H5PRT_READSTEP)
#define f_h5pt_step F77NAME(h5prt_step_,h5prt_step,H5PRT_STEP)
#define f_h5pt_wreal8 F77NAME(h5prt_wreal8_,h5prt_wreal8,H5PRT_WREAL8)
#define f_h5pt_wint8 F77NAME(h5prt_wint8_,h5prt_wint8,H5PRT_WINT8)
#endif

/* open/close interface */
#define f_h5pt_openr F77NAME(h5pt_openr_,h5pt_openr,H5PT_OPENR) /* func returns INT8 */
#define f_h5pt_openw F77NAME(h5pt_openw_,h5pt_openw,H5PT_OPENW) /* func returns INT8 */

#define f_h5pt_openr_par F77NAME(h5pt_openr_par_,h5pt_openr_par,H5PT_OPENR_PAR) /* func returns INT8 */
#define f_h5pt_openw_par F77NAME(h5pt_openw_par_,h5pt_openw_par,H5PT_OPENW_PAR) /* func returns INT8 */
#define f_h5pt_close F77NAME(h5pt_close_,h5pt_close,H5PT_CLOSE) 
#define f_h5pt_isvalid F77NAME(h5pt_isvalid_,h5pt_isvalid,H5PT_ISVALID) /* func returns INT8 */

/* writing interface */
#define f_h5pt_setnpoints F77NAME(h5pt_setnpoints_,h5pt_setnpoints,H5PT_SETNPOINTS)
#define f_h5pt_setstep F77NAME(h5pt_setstep_,h5pt_setstep,H5PT_SETSTEP)
#define f_h5pt_writedata_r8 F77NAME(h5pt_writedata_r8_,h5pt_writedata_r8,H5PT_WRITEDATA_R8)
#define f_h5pt_writedata_i8 F77NAME(h5pt_writedata_i8_,h5pt_writedata_i8,H5PT_WRITEDATA_I8)

/* Reading interface  (define dataset, step, particles, attributes) */
#define f_h5pt_getnsteps F77NAME(h5pt_getnsteps_,h5pt_getnsteps,H5PT_GETNSTEPS)
#define f_h5pt_getndatasets F77NAME(h5pt_getndatasets_,h5pt_getndatasets,H5PT_GETNDATASETS)
#define f_h5pt_getnpoints F77NAME(h5pt_getnpoints_,h5pt_getnpoints,H5PT_GETNPOINTS)
#define f_h5pt_getdatasetname F77NAME(h5pt_getdatasetname_,h5pt_getdatasetname,H5PT_GETDATASETNAME)

/* Views and parallelism */
#define f_h5pt_setview F77NAME(h5pt_setview_,h5pt_setview,H5PT_SETVIEW)
#define f_h5pt_resetview F77NAME(h5pt_resetview_,h5pt_resetview,H5PT_RESETVIEW)
#define f_h5pt_hasview F77NAME(h5pt_hasview_,h5pt_hasview,H5PT_HASVIEW)

/* Reading data */
#define f_h5pt_readdata_r8 F77NAME(h5pt_readdata_r8_,h5pt_readdata_r8,H5PT_READDATA_R8)
#define f_h5pt_readdata_i8 F77NAME(h5pt_readdata_i8_,h5pt_readdata_i8,H5PT_READDATA_I8)
#define f_h5pt_readdata F77NAME(h5pt_readdata_,h5pt_readdata,H5PT_READDATA)

/* Attributes */
/* writing */
#define f_h5pt_writefileattrib_r8 F77NAME(h5pt_writefileattrib_r8_,h5pt_writefileattrib_r8,H5PT_WRITEFILEATTRIB_R8)
#define f_h5pt_writefileattrib_i8 F77NAME(h5pt_writefileattrib_i8_,h5pt_writefileattrib_i8,H5PT_WRITEFILEATTRIB_I8)
#define f_h5pt_writefileattrib_string F77NAME(h5pt_writefileattrib_string_,h5pt_writefileattrib_string,H5PT_writefileattrib_string)
#define f_h5pt_writestepattrib_r8 F77NAME(h5pt_writestepattrib_r8_,h5pt_writestepattrib_r8,H5PT_WRITESTEPATTRIB_R8)
#define f_h5pt_writestepattrib_i8 F77NAME(h5pt_writestepattrib_i8_,h5pt_writestepattrib_i8,H5PT_WRITESTEPATTRIB_I8)
#define f_h5pt_writestepattrib_string F77NAME(h5pt_writestepattrib_string_,h5pt_writestepattrib_string,H5PT_WRITESTEPATTRIB_STRING)
/* reading */
#define f_h5pt_getnstepattribs F77NAME(h5pt_getnstepattribs_,h5pt_getnstepattribs,H5PT_GETNSTEPATTRIBS)
#define f_h5pt_getnfileattribs F77NAME(h5pt_getnfileattribs_,h5pt_getnfileattribs,H5PT_GETNFILEATTRIBS)
#define f_h5pt_getstepattribinfo F77NAME(h5pt_getstepattribinfo_,h5pt_getstepattribinfo,H5PT_GETSTEPATTRIBINFO)
#define f_h5pt_getfileattribinfo F77NAME(h5pt_getfileattribinfo_,h5pt_getfileattribinfo,H5PT_GETFILEATTRIBINFO)
#define f_h5pt_readstepattrib F77NAME(h5pt_readstepattrib_,h5pt_readstepattrib,H5PT_READSTEPATTRIB)
#define f_h5pt_readfileattrib F77NAME(h5pt_readfileattrib_,h5pt_readfileattrib,H5PT_READFILEATTRIB)


/* open/close interface */
haddr_t f_h5pt_openr(char *file,int flen){ /* func returns INT8 */
  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_h5pt_openw(char *file,int flen){ /* func returns INT8 */
  haddr_t fh;
  H5PartFile* f;
  file[flen]='\0'; /* null terminate */
  f = (H5PartOpenFile(file,H5PART_WRITE));
  /*  printf("openr file=[%s] flen=%u haddr=%u\n",file,flen,f); */
  fh = (haddr_t)f;
  /* printf("FileHandle=%llu\n",fh); */
  return fh;
}

#ifdef PARALLEL_IO
haddr_t f_h5pt_openr_par(char *file,void *c,int flen){ /* func returns INT8 */
  haddr_t fh;
  H5PartFile* f;
  file[flen]='\0'; /* null terminate */
#ifdef PARALLEL_IO
  f = H5PartOpenFileParallel(file,H5PART_READ,*((MPI_Comm*)c));
  /*  printf("openr file=[%s] flen=%u haddr=%u\n",file,flen,f); */
#else
  f=H5PartOpenFile(file,H5PART_READ);
#endif
  fh = (haddr_t)f;
  /* printf("FileHandle=%llu\n",fh); */
  return fh;
}
haddr_t f_h5pt_openw_par(char *file,MPI_Comm *c,int flen){ /* func returns INT8 */
  haddr_t fh;
  H5PartFile* f;
  file[flen]='\0'; /* null terminate */
#ifdef PARALLEL_IO
  f = H5PartOpenFileParallel(file,H5PART_WRITE,*((MPI_Comm*)c));
  /*  printf("openr file=[%s] flen=%u haddr=%u\n",file,flen,f); */
#else
  f=H5PartOpenFile(file,H5PART_WRITE);
#endif
  fh = (haddr_t)f;
  /* printf("FileHandle=%llu\n",fh); */
  return fh;
}
#endif

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

int f_h5pt_isvalid(haddr_t *file){
  return H5PartFileIsValid((H5PartFile*)*file);
}

/*==============Writing and Setting Dataset info========*/
#if 0
void f_h5pt_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);
}
#endif

void f_h5pt_setnpoints(haddr_t *file,uint64_t *np){
  H5PartSetNumParticles((H5PartFile*)*file,*np);
}

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

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

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

/*==============Reading Data Characteristics============*/

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

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

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

/* probably should get index from name */
int f_h5pt_getdatasetname(haddr_t *file,int *index,char *name,int namelen){
  return H5PartGetDatasetName((H5PartFile *)*file,*index,name,namelen);
}

uint64_t f_h5pt_getnumpoints(haddr_t *file){
  return H5PartGetNumParticles((H5PartFile*)*file);
}

/*=============Setting and getting views================*/

void f_h5pt_setview(haddr_t *file,int64_t *start,int64_t *end){
  H5PartSetView((H5PartFile*)*file,*start,*end);
}

void f_h5pt_resetview(haddr_t *file){
  H5PartResetView((H5PartFile*)*file);
}

int f_h5pt_hasview(haddr_t *file){
 return H5PartHasView(((H5PartFile*)*file));
}

void f_h5pt_getview(haddr_t *file,int64_t *start,int64_t *end){
  H5PartGetView((H5PartFile*)*file,start,end);
}
/*==================Reading data ============*/
int f_h5pt_readdata_r8(haddr_t *file,char *name,void *array,int namelen){
  int r;
  char *newname = (char*)malloc(namelen+1); /* be safe? */
  strncpy(newname,name,namelen);
  newname[namelen]='\0';
  r=H5PartReadDataFloat64((H5PartFile*)*file,newname,array);
  free(newname);
  return r;
}
int f_h5pt_readdata_i8(haddr_t *file,char *name,void *array,int namelen){
  int r;
  char *newname = (char*)malloc(namelen+1); /* be safe? */
  strncpy(newname,name,namelen);
  newname[namelen]='\0';
  r=H5PartReadDataInt64((H5PartFile*)*file,newname,array);
  free(newname);
  return r;
}

/*=================== Attributes ================*/

/* Attributes */
/* writing */
int f_h5pt_writefileattrib_r8(haddr_t *f,char *name,
                                void *attrib,int *nelem,int namelen){
  int r;
  char *newname = (char*)malloc(namelen+1); /* be safe? */
  strncpy(newname,name,namelen);
  newname[namelen]='\0';
  r=H5PartWriteFileAttrib((H5PartFile *)*f,newname,
                          H5T_NATIVE_DOUBLE,attrib,*nelem);
  free(newname);
  return r;
}

int f_h5pt_writefileattrib_i8(haddr_t *f,char *name,
                                void *attrib,int *nelem,int namelen){
  int r;
  char *newname = (char*)malloc(namelen+1); /* be safe? */
  strncpy(newname,name,namelen);
  newname[namelen]='\0';
  r=H5PartWriteFileAttrib((H5PartFile *)*f,newname,
                          H5T_NATIVE_INT64,attrib,*nelem);
  free(newname);
  return r;
}

int f_h5pt_writefileattrib_string(haddr_t *f,char *name,
                                char *attrib,int namelen,int attriblen){
  int r;
  char *newattrib;
  char *newname;
  newname = (char*)malloc(namelen+1); /* be safe? */
  strncpy(newname,name,namelen);
  newname[namelen]='\0';
  newattrib = (char*)malloc(attriblen+1); /* be safe? */
  strncpy(newattrib,attrib,attriblen);
  newattrib[attriblen]='\0';
  r=H5PartWriteFileAttribString((H5PartFile *)*f,newname,newattrib);
  free(newname);
  free(newattrib);
  return r;
}

int f_h5pt_writestepattrib_r8(haddr_t *f,char *name,
                                void *attrib,int *nelem,int namelen){
  int r;
  char *newname = (char*)malloc(namelen+1); /* be safe? */
  strncpy(newname,name,namelen);
  newname[namelen]='\0';
  r=H5PartWriteStepAttrib((H5PartFile *)*f,newname,
                          H5T_NATIVE_DOUBLE,attrib,*nelem);
  free(newname);
  return r;
}

int f_h5pt_writestepattrib_i8(haddr_t *f,char *name,
                                void *attrib,int *nelem,int namelen){
  int r;
  char *newname = (char*)malloc(namelen+1); /* be safe? */
  strncpy(newname,name,namelen);
  newname[namelen]='\0';
  r=H5PartWriteStepAttrib((H5PartFile *)*f,newname,
                          H5T_NATIVE_INT64,attrib,*nelem);
  free(newname);
  return r;
}

int f_h5pt_writestepattrib_string(haddr_t *f,char *name,
                                char *attrib,int namelen,int attriblen){
  int r;
  char *newattrib;
  char *newname;
  newname = (char*)malloc(namelen+1); /* be safe? */
  strncpy(newname,name,namelen);
  newname[namelen]='\0';
  newattrib = (char*)malloc(attriblen+1); /* be safe? */
  strncpy(newattrib,attrib,attriblen);
  newattrib[attriblen]='\0';
  r=H5PartWriteStepAttribString((H5PartFile *)*f,newname,newattrib);
  free(newname);
  free(newattrib);
  return r;
}

/* reading attributes ************************* */
int f_h5pt_getnstepattribs(haddr_t *f){
  return H5PartGetNumStepAttribs((H5PartFile*)*f);
}

int f_h5pt_getnfileattribs(haddr_t *f){
  return H5PartGetNumFileAttribs((H5PartFile*)*f);
}

int f_h5pt_getstepattribinfo(haddr_t *f,int *idx,char *name,
                                int *nelem,int maxnamelen){
  hid_t type;
  H5PartGetStepAttribInfo((H5PartFile*)*f,*idx,name,maxnamelen,&type,nelem);
  return 1;
}

int f_h5pt_getfileattribinfo(haddr_t *f,int *idx,char *name,
                                int *nelem,int maxnamelen){
  hid_t type;
  H5PartGetFileAttribInfo((H5PartFile*)*f,*idx,name,maxnamelen,&type,nelem);
  return 1;
}

int f_h5pt_readstepattrib(haddr_t *f,char *name,void *data,int namelen){
  int r;
  char *newname = (char*)malloc(namelen+1); /* be safe? */
  strncpy(newname,name,namelen);
  newname[namelen]='\0';
  H5PartReadStepAttrib((H5PartFile*)*f,newname,data);
}

int f_h5pt_readfileattrib(haddr_t *f,char *name,void *data,int namelen){
  int r;
  char *newname = (char*)malloc(namelen+1); /* be safe? */
  strncpy(newname,name,namelen);
  newname[namelen]='\0';
  H5PartReadFileAttrib((H5PartFile*)*f,newname,data);
}

---