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/quad_matrix.cpp $
5 $LastChangedRevision: 80 $
6 $Date: 2014-07-01 21:01:54 -0400 (Tue, 01 Jul 2014) $
8 ***********************************************************
11 THIS ROUTINE CALCULATES
13 Y = D*M#*A*M*X IF COPT .LE. 3
17 Y = C**P*(C*(SA+M#*A*M + SB*B**C) + SC)*X OTHERWISE
31 void quad_class::matrix(REAL* y, REAL* x, REAL* d, INTEGER* list, REAL* weight,
47 // COMPUTE SA*M#*A*M*X STORE IN Y
50 mtamx(x, y, list, weight);
60 for (i = 0; i < area; i++)
69 adsmos(j, l, bw1, bw2, bw3, bbcon);
73 badsmos(j, l, bw1, bw2, bbcon);
76 for (i = 0; i < area; i++)
80 delete[] j; // bug 92 - Lajos - 03/02/2005
86 adsmos(y, 1, cw1, cw2, cw3, cbcon);
90 badsmos(y, 1, cw1, cw2, cbcon);
97 for (i = 0; i < area; i++)
103 if ((copt != 5) && ((t - 1) > 0))
107 adsmos(y, t - 1, cw1, cw2, cw3, cbcon);
111 badsmos(y, t - 1, cw1, cw2, cbcon);
117 for (i = 0; i < area; i++)
123 if ((copt != 5) && (t > 0))
127 adsmos(y, t, cw1, cw2, cw3, cbcon);
131 badsmos(y, t, cw1, cw2, cbcon);
138 for (i = 0; i < area; i++)