#ifndef DATATYPES_H
#define DATATYPES_H

#include <math.h>

#define PRECISION_SINGLE

#ifdef PRECISION_SINGLE
typedef float ptype;
#elif defined PRECISION_DOUBLE
typedef double ptype;
#endif

#define BLOCK	256
#define SQRT_BLOCK 16

#define PI	3.14159

//a very small number
#define EPSILON		0.00001

//CUDA hybrid code - complex class should run on both the CPU and GPU


typedef ptype fieldPoint;

//hybrid GPU/CPU complex data typ
#include "rts/rtsComplex.h"
#include "rts/rtsVector.h"
#include "rts/rtsPoint.h"
#include "rts/rtsQuad.h"

typedef rts::rtsComplex<ptype> bsComplex;
typedef rts::rtsVector<ptype, 3> bsVector;
typedef rts::rtsPoint<ptype, 3> bsPoint;
typedef rts::rtsQuad<ptype, 3> bsRect;


#endif