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_cri1.cpp $
5 $LastChangedRevision: 107 $
6 $Date: 2014-07-07 12:06:15 -0400 (Mon, 07 Jul 2014) $
8 ***********************************************************
18 REAL sec_cri1_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 = 1; i < pixelcount; i++)
27 // skip right-most column
28 if (i % (dimensions) == 0) continue;
30 tau += sqrt(pow(image[i - 1] - image[i - 1 + dimensions], 2) + pow(image[i - 1] - image[i], 2));