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_forwrd.cpp $
5 $LastChangedRevision: 85 $
6 $Date: 2014-07-02 16:07:08 -0400 (Wed, 02 Jul 2014) $
8 ***********************************************************
11 PURPOSE: READ IN PROJECTIONS , TRANSFORM THEN INTERPOLATE
23 void foru_class::forwrd()
44 // READ IN 1ST PROJECTION AND TRANSFORM
48 prjtrn(0, nold, &phi1, &dphi1);
52 for (np = 1; np < GeoPar.prjnum; np++)
55 // READ IN NEXT PROJECTION , TRANSFORM THEN INTERPOLATE
58 prjtrn(np, new1, &phi2, &dphi2);
63 intrpl(phi1, phi2, dphi1, dphi2, nold, new1);
75 // READ IN 1ST PROJECTION AGAIN, INTERPOLATE USING LAST AND 1ST PROJS
78 prjtrn(0, new1, &phi2, &dphi2);
83 intrpl(phi1, phi2, dphi1, dphi2, nold, new1);
89 kh = Fourie.nsize2 * (Fourie.nsize2 + 2);
91 for (k = 2; k < kh; k++)
93 Fourie.nfrpln[k] /= 2.0;
97 fprintf(output, "\ntime for projection transforms : %10.3f seconds ", tforw); // changed precision to three digits - swr 1/21/06
98 fprintf(output, "\ntime for interpolations : %10.3f seconds ", tintr); //ditto