2 ***********************************************************
3 $SNARK_Header: S N A R K 1 4 - A PICTURE RECONSTRUCTION PROGRAM $
4 $HeadURL: svn://dig.cs.gc.cuny.edu/snark/trunk/src/snark/foru_frfilt.cpp $
5 $LastChangedRevision: 85 $
6 $Date: 2014-07-02 16:07:08 -0400 (Wed, 02 Jul 2014) $
8 ***********************************************************
21 void foru_class::frfilt(REAL* vector, INTEGER size1, INTEGER nfiltr, REAL cutoff, REAL phi)
28 wvqntm = (REAL) 1.0 / ((REAL) (size1 / 2));
30 wvqntm /= MAX0((REAL) fabs(cos(phi)), (REAL) fabs(sin(phi)));
35 for (k = 0; k <= size1 - 2; k += 2)
36 { // changed "<" to "<=" to match Fortran code. Lajos, Jan 13, 2005
47 for (k = 0; k <= size1 - 2; k += 2)
48 { // changed "<" to "<=" to match Fortran code. Lajos, Jan 13, 2005
56 fact = sinc(wvnum / cutoff);
58 vector[k + 1] *= fact;
65 for (k = 0; k <= size1 - 2; k += 2)
66 { // changed "<" to "<=" to match Fortran code. Lajos, Jan 13, 2005
74 fact = (REAL) cos(Consts.pid2 * wvnum / cutoff);
76 vector[k + 1] *= fact;
83 for (k = 0; k <= size1 - 2; k += 2)
84 { // changed "<" to "<=" to match Fortran code. Lajos, Jan 13, 2005
93 + ((REAL) 1.0 - cutoff) * (REAL) cos(Consts.pi * wvnum);
95 vector[k + 1] *= fact;