#ifndef TIPSY_INCLUDED
#define TIPSY_INCLUDED

#include <rpc/types.h>
#include <rpc/xdr.h>

#define TIPSY_TYPE_GAS	1
#define TIPSY_TYPE_DARK	2
#define TIPSY_TYPE_STAR	3

#define TIPSY_BLOCK 10000

/*
** Set DOUBLE_POS or DOUBLE_POS_VEL for double precision positions and
** velocities.
*/
//#define DOUBLE_POS
//#define DOUBLE_POS_VEL
struct gas_particle {
    float mass;
#if defined(DOUBLE_POS) || defined(DOUBLE_POS_VEL)
    double pos[3];
#else
    float pos[3];
#endif
#ifdef DOUBLE_POS_VEL
	double vel[3];
#else
    float vel[3];
#endif
    float rho;
    float temp;
    float hsmooth;
    float metals;
    float phi;
    };

struct dark_particle {
    float mass;
#if defined(DOUBLE_POS) || defined(DOUBLE_POS_VEL)
    double pos[3];
#else
    float pos[3];
#endif
#ifdef DOUBLE_POS_VEL
	double vel[3];
#else
    float vel[3];
#endif
    float eps;
    float phi;
    };

struct star_particle {
    float mass;
#if defined(DOUBLE_POS) || defined(DOUBLE_POS_VEL)
    double pos[3];
#else
    float pos[3];
#endif
#ifdef DOUBLE_POS_VEL
    double vel[3];
#else
    float vel[3];
#endif
    float metals;
    float tform;
    float eps;
    float phi;
    };

struct base_particle {
    float mass;
#if defined(DOUBLE_POS) || defined(DOUBLE_POS_VEL)
    double pos[3];
#else
    float pos[3];
#endif
#ifdef DOUBLE_POS_VEL
	double vel[3];
#else
    float vel[3];
#endif
    };

struct dump {
    double time;
    unsigned nbodies;
    unsigned ndim;
    unsigned nsph;
    unsigned ndark;
    unsigned nstar;
    };


typedef struct TipsyContext {
    int bNative;
    unsigned iCurr;
    int iCurrType;
    unsigned nRemain;
    unsigned nBlock;
    unsigned iBlock;
    struct gas_particle gpCurr;
    struct dark_particle dpCurr;
    struct star_particle spCurr;
    unsigned nMaxGas,nGas;
    unsigned nMaxDark,nDark;
    unsigned nMaxStar,nStar;
    struct gas_particle *gp;
    struct dark_particle *dp;
    struct star_particle *sp;
    FILE *fp;
    XDR xdr;
    struct dump head;
    char *pchInFile;
    } * TCTX;


int xdr_header(XDR *xdrs, struct dump *header);
int xdr_gas(XDR *xdrs,struct gas_particle *p);
int xdr_dark(XDR *xdrs,struct dark_particle *p);
int xdr_star(XDR *xdrs,struct star_particle *p);

void TipsyInitialize(TCTX *pctx,int bNative,char *pchInFile);
void TipsyFinish(TCTX ctx);
struct base_particle *pTipsyReadNative(TCTX ctx,int *piType,double *pdSoft);
struct base_particle *pTipsyRead(TCTX ctx,int *piType,double *pdSoft);
struct base_particle *pTipsyParticle(TCTX ctx,unsigned iIndex,int *piType,double *pdSoft);
void TipsyReadAll(TCTX ctx);
void TipsyAddGas(TCTX ctx,struct gas_particle *pGas);
void TipsyAddDark(TCTX ctx,struct dark_particle *pDark);
void TipsyAddStar(TCTX ctx,struct star_particle *pStar);
void TipsyWriteAll(TCTX ctx,double dTime,char *pchFileName);
unsigned iTipsyNumParticles(TCTX ctx);
double dTipsyTime(TCTX ctx);

#endif