#include <stdio.h>
#include <math.h>
#include <time.h>
#include <stdlib.h>
#define pi 3.1415926
#define Resolution 256
#define R Resolution
#define TubeNum 256
#define TN TubeNum
#define ProjNum 200
#define PN ProjNum
#define PixelSize 2.0/256
#define PS PixelSize
#define RelativeTubeSpace 1 /* Compared to PixelSize */
#define RTS RelativeTubeSpace
#define AbsoluteTubeSpace PS*RTS
#define ATS AbsoluteTubeSpace
#define IterationNum 999
#define IN IterationNum
int main ( )
{
double begin=time(NULL);
typedef struct {
double l;
int index;
}Coef;
Coef CoefArray[2*R];
FILE * sino; //Read sinogram in
sino=fopen("Sinogram","rb");
if (sino==NULL)
{
printf("Open file failure!\n");
}
double * Sinogram=new double[PN*TN];
fread(Sinogram,sizeof(double),PN*TN,sino);
if (feof(sino)==0)
{
printf("Error!Reading sinogram has not been finished.\n");
}
fclose(sino);
FILE * ei;
ei=fopen("EI","rb");
if (ei==NULL)
{
printf("Open file failure!\n");
}
double * EI=new double[PN*TN];
fread(EI,sizeof(double),PN*TN,ei);
if (feof(ei)==0)
{
printf("Error!Reading EI has not been finished.\n");
}