#include <stdio.h>
#include <string.h>
#include <math.h>

#define TINY 0.000001
#define GIANT 100000

#define TRUE 1
#define FALSE 0

#define SEPARABLE  1
#define NON_SEPARABLE  2
#define NON_SYMMETRIC 3

#define KERNEL_1D_NX 500
#define KERNEL_3D_NX 100

#define Swap(x,y,z) z=x; x=y; y=z;
#define Ceil(x) (int)(x+1.0-TINY)
#define Floor(x) (int)x
#define Max(a,b) (a>b ? a : b)
#define Min(a,b) (a<b ? a : b)
#define Rad(a) (a*M_PI/180.0)
#define Sqr(x) (x*x)

typedef struct {
  double x,y,z;
} Vector;

typedef double Matrix[2][2];
typedef double Matrix3[3][3];

typedef struct {
  double r,g,b;
} Color;

typedef struct {
  double kd,ks,ka,phongExp;
  Color color; 
} ObjectProperties;

typedef struct {
  Vector pos;
  Color color;
} LightSource;

typedef struct {
  Vector VRP,VPN,VUP;
  Vector EYE,VPU,VPV;
  Matrix3 volRotationMat;
  Vector rotAngle;
  double imgWidth,imgHeight;
  int viewAngle;
  int nxImg,nyImg;
  int stripIntegral;
  char volumeSetName[300];
  char splatFileName[300];
  char splatMethod[20];
  char sliceFileName[300];
  int firstSlice,lastSlice;
  int numberOfRuns;
  int floatingSplat;
  double rayStepSize;
  char intKernel[300];
  double intKernelCoeffB,intKernelCoeffC;
  char derivIntKernel[300];
  double derivIntKernelCoeffB,derivIntKernelCoeffC;
  char derivKernel[300];
  double derivKernelCoeffB,derivKernelCoeffC;
  char derivMethod[100];
  char normalEstimationMethod[20];
  double isovalue;
  LightSource lightSource;
  Color ambientColor;
  ObjectProperties objProps;
  int colImgFlag;
  char intersectionImageFlag[40];
  char gradientVolumeFlag[40];
  char outputFileName[300];
  char inter_class[3]; 
  char deriv_class[3];
  int smoothIntFilterCode, smoothDerivFilterCode;
  float smoothIntFilterCoeff, smoothDerivFilterCoeff;
} ViewParameters;

typedef struct {
  unsigned char *data;
  signed char *sdata;
  float *fdata;
  signed short *stdata;
  unsigned short *ustdata; 
  int nx,ny,nz,n;
  Vector bbLow,bbHi;
} Volume;

typedef struct {
  double *data;
  double extent,extent2,scale;
  int nx,xc;
  double extentX,extentY,extentX2,extentY2;
  int ny;
} SplatTable;

typedef struct{
  double *data;
  double extent,scale;
  int nx,ny,nz;
  int type;
} InterpolationKernel;

typedef struct {
  int nx,ny;
  double *data;
  Color *colorData;
} Image;

typedef struct {
  Vector *pt;
  int nx,ny;
  int vMin,vMax;
  int *uMin,*uMax;
} VolumeIntersectionPoints;

typedef struct {
  int vertIndx[4];
} Polygon;  

typedef struct {
  double u,v;
} ImageSpaceVector;


void SplatParallelLine(ViewParameters *vpars,Volume **volArr,SplatTable *splat,double *imageData);    
void SplatParallelStrip(ViewParameters *vpars,Volume *vol,SplatTable *splat,double *imageData);
void SplatPerspectiveLine(ViewParameters *vpars,Volume **volArr,SplatTable *splat,double *imageData);
void SplatPerspectiveStrip(ViewParameters *vpars,Volume *vol,SplatTable *splat,double *imageData);
void  Raycast(ViewParameters *vpars,Volume **volArr,Color *imageData,InterpolationKernel **kernelArr);

int ReadViewParameterFile(ViewParameters *vpars,char *fileName);
void ComputeViewPlaneVectors(ViewParameters *vpars);
void NormalizeV(Vector *v);
double Dot(Vector *v1,Vector *v2);
Vector Cross(Vector *v1,Vector *v2);
void ComputeEyePoint(ViewParameters *vpars);
void SetRotationMatrix(Matrix3 mat,double angle,char axis);
void InitMat3(Matrix3 mat);
void Mat3Mult(Matrix3 mat1,Matrix3 mat2,Matrix3 matout);
Vector VecMat3Mult(Vector *v,Matrix3 m);
int ComputeVolumeIntersectionPoints(ViewParameters *vpars,Volume *vol,VolumeIntersectionPoints *volIntersectPts);
int CreateGradientVolumes(Volume *vol,Volume *dxVol,Volume *dyVol,Volume *dzVol,Volume *magVol,ViewParameters *vpars);
int getFiltWeights(int filtCode,float a,float t,float *coeff);