/home/epics/devel/tomoRecon/tomoReconApp/src/grid.h

Go to the documentation of this file.

/*** File grid.h  12:27 PM 11/7/97 **/

#define ANSI

/**** System includes */
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>
#include <stddef.h>
#include <time.h>
#include <sys/stat.h>

#include <fftw3.h>

/**** Macros and typedefs ****/
#ifndef max
#define max(A,B) ((A)>(B)?(A):(B))
#endif
#ifndef min
#define min(A,B) ((A)<(B)?(A):(B))
#endif
#define free_matrix(A) (free(*(A)),free(A))
#define abs(A) ((A)>0 ?(A):-(A))
#define pi  3.14159265359
#define TOLERANCE 0.1   /* For comparing centers of two sinograms */
#define LTBL_DEF 512    /* Default lookup table length */

typedef struct {
    float r;  
    float i;  
} complex;

typedef struct {
   int n_ang;       
   int n_det;       
   int geom;              
   float *angles;       
   float center;        
} sg_struct;

typedef struct {
   float C;             
   int nt;          
   float lmbda;           
   float coefs[15];     
} pswf_struct;

typedef struct {
   pswf_struct *pswf;        
   float sampl;                
   float MaxPixSiz;          
   float R;                          
   float X0;                       
   float Y0;                       
   char fname[16];                   
   float (*filter)(float); 
   long ltbl;                      
   int verbose;            
} grid_struct;

#ifdef __cplusplus

class grid {  
public:
  grid(grid_struct *GP,sg_struct *SGP, long *imgsiz);
  ~grid();
  void recon(float center, float** G1,float** G2,float*** S1,float*** S2);
  void filphase_su(long pd,float fac, float(*pf)(float),complex *A);
  void pswf_su(pswf_struct *pswf,long ltbl, 
               long linv, float* wtbl,float* dwtbl,float* winv);
  
private:
  int flag;       
  int n_det;
  int n_ang;
  long pdim;
  long M;
  long M0;
  long M02;
  long ltbl;
  float sampl;
  float scale;
  float L;
  float X0;
  float Y0;
  float *SINE;
  float *COSE;
  float *wtbl; 
  float *dwtbl;
  float *work;
  float *winv;
  float previousCenter;
  float (*filter)(float);
  complex *cproj;
  complex *filphase;
  complex **H;
  fftwf_complex *HData;
  
  fftwf_plan forward_1d_plan;
  fftwf_plan backward_2d_plan;
  int verbose;   /* Debug printing flag */
};
#endif

#ifdef __cplusplus
extern "C" {
#endif

/**** Function Prototypes ****/

void get_pswf(float C, pswf_struct **P);

float (*get_filter(char *name))(float);

#ifdef __cplusplus
}
#endif