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_CReadInp.cpp $
5 $LastChangedRevision: 85 $
6 $Date: 2014-07-02 16:07:08 -0400 (Wed, 02 Jul 2014) $
8 ***********************************************************
10 ===========================================================================
11 Read information of the input file and allocates the work space
12 ===========================================================================
24 void emap_class::ReadInp(INTEGER area, INTEGER nelem, INTEGER usrays,
25 INTEGER prjnum, BOOLEAN* flag) // nelem is not used. hstau 6/9/2003, flag added 1/04
29 INTEGER doit = CHAR2INT('e', 'v', 'a', 'l');
31 //------ read parameters from input file
32 ///MAPEM.gamma1 = getnum(input, &eol);
35 MAPEM.gamma1 = InFile.getnum(TRUE, &eol);
38 fprintf(output, "\n **** error - must specify gamma ***");
48 eval = InFile.getwrd(FALSE, &eol, exp, 1);
55 MAPEM.sclfactor = 1.0;
58 //------ write the parameters into the output file
60 "\n\n -----------------------------------------------------------\n");
63 "\n maximum a-posteriori probability expectation maximization");
65 fprintf(output, "\n\n gamma:%8.3f", MAPEM.gamma1);
69 fprintf(output, "\n evaluation flag is set ");
74 fprintf(output, "\n evaluation flag is not set");
78 "\n\n -----------------------------------------------------------\n");
80 //------ allocate space for the intermediate results
82 MAPEM.qj = new REAL[area];
86 MAPEM.Kx = new REAL[area];
87 MAPEM.Sx = new REAL[area];
88 MAPEM.Sjjinv = new REAL[area];
91 MAPEM.Wjj = new REAL[area];
92 MAPEM.prj = new REAL[prjnum * usrays];