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_intrpl.cpp $
5 $LastChangedRevision: 85 $
6 $Date: 2014-07-02 16:07:08 -0400 (Wed, 02 Jul 2014) $
8 ***********************************************************
12 PURPOSE: GIVEN TWO TRANSFORMED PROJECTIONS IN ARRAY STARTING AT
13 NOLD AND NEW WITH ANGLES ANG1 AND ANG2(DANG1 AND DANG2
14 ARE ANGLES IN DEGREE ) FIND ALL POINTS THAT LIE BETWEEN
15 THE TWO RAYS AND PERFORM INTERPOLATIONS
20 ICASE = 1 IF THE TWO ANGLES LIE BETWEEN -90 AND -45
21 2 IF THE TWO ANGLES STRADDLE THE -45 DEGREE LINE
22 3 IF THE TWO ANLES LIE BETWEEN -45 AND 45 LINES
23 4 IF THE TWO ANGLES STRADDLE THE 45 DEGREE LINE
24 5 IF THE TWO ANGLES LIE BETWEEN THE 45 AND 90 LINES
25 6 IF ANG1 IS IN 1ST QUADRANT AND ANG2 IN 4TH QUADRANT
36 void foru_class::intrpl(REAL ang1, REAL ang2, REAL dang1, REAL dang2, REAL* nold, REAL* new1)
42 if ((dang1 < -45.) && (dang2 > -45.))
44 if ((dang1 >= -45.) && (dang2 <= 45.))
46 if ((dang1 < 45.) && (dang2 > 45.))
53 // IF UNI THEN CALL THE INTERPOLATION ROUTINE UNIF
54 // IF VRI THEN CALL THE ROUTINE VCASE1,...,OR VC56 DEPENDING ON ICASE
57 unif(icase, ang1, ang2, nold, new1);
67 vcase1(ang1, ang2, nold, new1);
71 vcase2(ang1, ang2, nold, new1);
75 vcase3(ang1, ang2, nold, new1);
79 vcase4(ang1, ang2, nold, new1);
84 vc56(icase, ang1, ang2, nold, new1);