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/nav3.cpp $
5 $LastChangedRevision: 49 $
6 $Date: 2014-06-16 21:49:04 -0400 (Mon, 16 Jun 2014) $
8 ***********************************************************
10 *************************************************************************************
11 Dummy class for the calculation of the non-ascending vector corresponding to sec_cri3
12 *************************************************************************************
20 void nav3_class::Run(REAL* xkn, REAL* vkn)
26 pixelcount = GeoPar.area;
30 pixelcount = Blob.area;
33 REAL gradientSquareNorm = 0;
35 // calculate the gradient of your phi function here
36 // calcGradient(xkn, vkn);
38 for (INTEGER i = 0; i < pixelcount; i++) gradientSquareNorm += pow(vkn[i], 2);
39 REAL gradientNorm = sqrt(gradientSquareNorm);
40 for (INTEGER i = 0; i < pixelcount; i++) vkn[i] = (-1.0) * vkn[i] / gradientNorm;