From a95e41ac40cd2f3a4401d921618604cf33f2a904 Mon Sep 17 00:00:00 2001
From: "Kevin M. Rosenberg" <kevin@rosenberg.net>
Date: Tue, 20 Jun 2000 18:46:29 +0000
Subject: [PATCH] r116: *** empty log message ***

---
 html/simulate.html.in |  8 ++++----
 include/Makefile.am   |  3 ++-
 include/filter.h      | 38 ++++++++++++++++++++++++++++++++++++++
 3 files changed, 44 insertions(+), 5 deletions(-)
 create mode 100644 include/filter.h

diff --git a/html/simulate.html.in b/html/simulate.html.in
index 8b855a5..018a838 100644
--- a/html/simulate.html.in
+++ b/html/simulate.html.in
@@ -17,7 +17,7 @@ Phantom Name:<br>
 Nx: <input type="text" name="Phantom_Nx" size="4" value="256"><br>
 Ny: <input type="text" name="Phantom_Ny" size="4" value="256"><br>
 <p>
-Number of Samples<br>per pixel: <input type="text" name="Phantom_NSample" size="2" value="1"><br>
+Number of Samples<br>per pixel: <input type="text" name="Phantom_NSample" size="2" value="2"><br>
 <p>
 Display Window:<br>
 Minimum: <input type="text" name="Disp_Min" size=5 value="auto"><br>
@@ -30,7 +30,7 @@ MPI Supercomputing:<br>
 <h3>Simulate X-Ray acquistion</h3>
 Number of detectors: <input type="text" name="PJ_NDet" size="4" value="367"><p>
 Number of Rotations: <input type="text" name="PJ_NRot" size="4" value="320"><p>
-Number of Rays<br>(samples) per detector: <input type="text" name="PJ_NRay" size="2" value="1"><p>
+Number of Rays<br>(samples) per detector: <input type="text" name="PJ_NRay" size="2" value="3"><p>
 Rotation Angle<br>as a multiple of PI: <input type="text" name="PJ_RotAngle" size="3" value="1.0"><br>
 </td><td valign="top">
 <H3>Image Reconstruction</H3>
@@ -53,8 +53,8 @@ Backprojection Method:<br>
 <input type="radio" name="IR_Backproj" value="trig">Trig<br>
 <input type="radio" name="IR_Backproj" value="table">Trig Table<br>
 <input type="radio" name="IR_Backproj" value="diff">Difference<br>
-<input type="radio" name="IR_Backproj" value="diff2" checked>Difference Optimized<br>
-<input type="radio" name="IR_Backproj" value="idiff2">Difference Optimized (integer math)<br>
+<input type="radio" name="IR_Backproj" value="diff2">Difference Optimized<br>
+<input type="radio" name="IR_Backproj" value="idiff2" checked>Difference Optimized (integer math)<br>
 </td></tr></table>
 <p>
 <input type="submit" value="Simulate"><input type="reset" value="clear fields">
diff --git a/include/Makefile.am b/include/Makefile.am
index 087655f..fa5be04 100644
--- a/include/Makefile.am
+++ b/include/Makefile.am
@@ -1,4 +1,5 @@
-noinst_HEADERS=ct.h ezplot.h pol.h sgp.h array2d.h imagefile.h backprojectors.h mpiworld.h fnetorderstream.h phantom.h timer.h sstream scanner.h projections.h ctsupport.h
+noinst_HEADERS=ct.h ezplot.h pol.h sgp.h array2d.h imagefile.h backprojectors.h mpiworld.h fnetorderstream.h phantom.h timer.h sstream scanner.h projections.h ctsupport.h filter.h
+
 
 
 
diff --git a/include/filter.h b/include/filter.h
new file mode 100644
index 0000000..274d792
--- /dev/null
+++ b/include/filter.h
@@ -0,0 +1,38 @@
+class SignalFilter {
+ public:
+
+    SignalFilter (const FilterType filt_type, double bw, double xmin, double xmax, int n, double param, const DomainType domain, const int numInt);
+
+    ~SignalFilter (void);
+
+    double* getFilter (void) const
+      { return m_vecFilter; }
+
+    double convolve (const double f[], const double dx, const int n, const int np, const FunctionSymmetry func_type) const;
+
+    double convolve (const float f[], const double dx, const int n, const int np, const FunctionSymmetry func_type) const;
+
+    static double spatialResponseCalc (FilterType fType, double bw, double x, double param, int n);
+
+    static double spatialResponseAnalytic (FilterType fType, double bw, double x, double param);
+
+    static double frequencyResponse (FilterType fType, double bw, double u, double param);
+
+ private:
+    double *m_vecFilter;
+    double m_bw;
+    FilterType m_filterType;
+
+    double spatialResponseCalc (double x, double param, int n) const;
+
+    double spatialResponseAnalytic (double x, double param) const;
+
+    double frequencyResponse (double u, double param) const;
+
+    static double sinc (double x, double mult)
+      { return (fabs(x) > F_EPSILON ? (sin (x * mult) / x) : 1.0); }
+
+    static double integral_abscos (double u, double w);
+
+};
+
-- 
2.34.1