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/errpar.cpp $
5 $LastChangedRevision: 85 $
6 $Date: 2014-07-02 16:07:08 -0400 (Wed, 02 Jul 2014) $
8 ***********************************************************
10 THIS ROUTINE CHECKS AND PREPARES THE CONSTANTS USED BY THE ERROR
23 void errpar(BOOLEAN alg)
29 if (!(NoisePar.quanin == 0))
35 for (i = 0; i < Spctrm.nergy; i++)
37 Err.bf += Spctrm.engwt[i] * (REAL) exp(-Spctrm.backgr[i]);
40 Err.bf = (REAL) 1.0 / (Err.bf * NoisePar.quanmn);
41 Err.af = Err.bf * ((REAL) 2.0 / NoisePar.quancm + (REAL) 1.0);
47 // MULTIPLICATIVE NOISE
51 Err.mfactr = (NoisePar.ultnsd / NoisePar.ultnmn)
52 * (NoisePar.ultnsd / NoisePar.ultnmn);
61 Err.afactr = NoisePar.addnsd * NoisePar.addnsd;