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/lines.cpp $
5 $LastChangedRevision: 80 $
6 $Date: 2014-07-01 21:01:54 -0400 (Tue, 01 Jul 2014) $
8 ***********************************************************
17 void lines(REAL* a, REAL* b, INTEGER nrow, INTEGER mcol, REAL scale, INTEGER lineno, INTEGER* linumb, FILE* file)
29 fprintf(file, "\n scaling factor = %10.5f\n", scale);
30 fprintf(file, "\n column = %4i", linumb[0]);
31 for (i = 1; i < lineno; i++)
33 fprintf(file, " %4i", linumb[i]);
36 fprintf(file, "\n row ");
37 for (i = 0; i < lineno; i++)
39 fprintf(file, " original estimate difference ");
43 for (irow = 0; irow < nrow; irow++)
46 for (jcolum = 0; jcolum < lineno; jcolum++)
48 index = (irow) * mcol + linumb[jcolum];
51 linout[one] = scale * a[index];
52 linout[one + 1] = scale * b[index];
53 linout[one + 2] = linout[one + 1] - linout[one];
57 fprintf(file, "\n%4i ", irow);
58 for (i = 0; i < linno3; i++)
60 fprintf(file, " %9.5f", linout[i]);
63 fprintf(file, "\n\n");