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_nxtpct.cpp $
5 $LastChangedRevision: 89 $
6 $Date: 2014-07-02 17:24:53 -0400 (Wed, 02 Jul 2014) $
8 ***********************************************************
23 void quad_class::nxtpct(REAL* x, REAL* cg, REAL delta)
37 //pr = (Blob.H + Blob.M)%2;
48 for (i = 0; i < area; i++)
50 x[i] -= delta * cg[i];
53 // CONSTRAIN VALUES AND AVERAGE DENSITY ADJUSTMENT
55 prjden = GeoPar.aveden * ((REAL) (GeoPar.area));
59 for (i = 0; i < area; i++)
61 if (Blob.pix_basis || (!Blob.pix_basis && i % 2 == Blob.pr))
65 tem = MAX0(tem, Modefl.lower);
67 tem = MIN0(tem, Modefl.upper);
76 if (fabs(temden) <= Consts.zero)
79 "\n\n reconstructed picture is zero, your request for picture normalizing set the picture to aveden");
81 for (i = 0; i < area; i++)
90 ratio = prjden / temden;
93 fprintf(output, "\n\n ratio of densities ,%15.4e", ratio); //bug 147, wei, 2005/10
96 for (i = 0; i < area; i++)
106 // CALCULATE INDICATOR SIZE OF CORRECTION
109 for (i = 0; i < area; i++)
111 if (Blob.pix_basis || (!Blob.pix_basis && i % 2 == Blob.pr))
113 snorm += cg[i] * cg[i];
117 snorm = (REAL) sqrt(snorm);
118 corr = (REAL) fabs(delta) * snorm;
120 fprintf(output, "\n\n norm of correction ,%15.4e", corr); //bug 147, wei, 2005/10
121 fprintf(output, "\n delta %15.4e", delta);
122 fprintf(output, "\n norm of corr vector %15.4e", snorm);