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/emap_Cwray.cpp $
5 $LastChangedRevision: 85 $
6 $Date: 2014-07-02 16:07:08 -0400 (Wed, 02 Jul 2014) $
8 ***********************************************************
22 BOOLEAN emap_class::findqpart(INTEGER np, INTEGER nr, REAL yi, REAL* recon,
25 //----- local variables
26 INTEGER t, nt, bigT, n[4096];
27 REAL l[4096], snorm, si;
29 //------ Body of outer loop in Algorithm FIND_q~ for the MAP algorithm * /
33 wray(np, nr, n, l, &bigT, &snorm);
37 Blob.bwray(np, nr, n, l, &bigT, &snorm); // l and n must have enough size for blob basis. Why 4096? Hstau 6/25/2003
43 for (t = 0; t < bigT; t++)
45 si += recon[n[t]] * l[t];
51 "\n****The pseudo raysum for np = %d, nr = %d is zero, EMAP is not applicable, EMAP aborted!\n",
54 } //added by wei 1/2005
57 for (t = 0; t < bigT; t++)
60 qtilde[nt] += yi / si * l[t];