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/sec_cri2.cpp $
5 $LastChangedRevision: 107 $
6 $Date: 2014-07-07 12:06:15 -0400 (Mon, 07 Jul 2014) $
8 ***********************************************************
18 REAL sec_cri2_class::Run(REAL* image)
20 INTEGER dimensions = GeoPar.nelem;
21 // pixels over which to iterate (skip last line and new bottom right pixel)
22 INTEGER pixelcount = GeoPar.area - dimensions - 2;
25 for (INTEGER i = dimensions + 2; i < pixelcount; i++)
27 // skip right-most column
28 if (i % (dimensions) == 0) continue;
30 if (i % (dimensions) == 1) continue;
32 psi += pow(image[i - 1] - ((image[i - 1 - dimensions - 1] + image[i - 1 - dimensions] + image[i - dimensions] + image[i - 2] + image[i] + image[i - 1 + dimensions - 1] + image[i - 1 + dimensions] + image[i + dimensions]) / 8), 2);