Added snark14m distribution examples
[snark14.git] / examples / b2 / b2r.out
1    snark14UserDefined.s170426 - A PICTURE RECONSTRUCTION PROGRAM
2
3
4      <*> *** EXAMPLE 2
5
6      <*>     Creation of a small test pattern and projection data with execution of
7
8      <*>     a user-defined reconstruction algorithm and analysis.
9
10      <*>     
11
12      <#> CREATE                                                                                                                                                                                                                                                           
13
14          Example 2,  A simple test pattern and projection data
15
16
17      <#> SPECTRUM MONOCHROMATIC with energy at 75 KeV
18          energy spectrum is monochromatic at energy level    75
19
20
21      <#> OBJECTS   cx    cy    u     v     ang   den
22          description of objects
23                                                                 density at levels
24  numb type  x-coord  y-coord x-length y-length    angle  av dens        75
25
26     1 elip   0.0000   0.0000   0.9000   0.9000   0.0000   2.5000    2.5000
27
28     2 elip   0.0000   0.0000   0.8000   0.6000   0.0000  -1.2500   -1.2500
29
30     3 elip   0.0000   0.0000   0.2100   0.2500   0.0000   0.1000    0.1000
31
32          scale factor multiplying object densities     1.0000
33
34          seed set to 0
35          inhomogeneity set to     0.0000
36
37      <#> PHANTOM AVERAGE 1 point for estimating pixel density
38
39          this run will generate a phantom
40          density in each pixel is obtained as the average of 1 x 1 points
41
42
43      <#> number of elements is 31  Size of a pixel is 0.07
44          picture size 31 x 31,  pixel size     0.0700
45
46
47      <#> RAYSUM AVERAGE 1 line for estimation of ray sum with weight
48
49          this run will generate projection data
50          projection data are calculated by dividing each ray interval into 1 substrips
51
52          with aperture (substrip) weights     1
53
54
55      <#> GEOMETRY
56
57
58      <#> PARALLEL UNIFORM LINE
59          rays are parallel with uniform spacing between rays
60          data collected along lines
61
62
63
64      <#> RAYS USER-defined 63 rays per projection at spacing 0.05
65          number of rays per projection    63
66          at detector spacing     0.0500
67
68
69      <#> ANGLES 30  EQUAL SPACING between the first and last angles which are
70          total number of projections    30
71
72          projection angles    0.0000    6.0000   12.0000   18.0000   24.0000   30.0000   36.0000   42.0000   48.0000   54.0000
73                              60.0000   66.0000   72.0000   78.0000   84.0000   90.0000   96.0000  102.0000  108.0000  114.0000
74                             120.0000  126.0000  132.0000  138.0000  144.0000  150.0000  156.0000  162.0000  168.0000  174.0000
75
76
77      <#> MEASUREMENT PERFECT
78          projection data are noiseless
79
80
81      <#> BACKGROUND absorption is 0.0
82                                at levels
83                                       75 
84           background absorption   0.0000 
85  
86
87      <#> RUN
88               0.004 seconds phantom creation
89               0.000 seconds projection data creation
90               0.004 seconds used for processing command crea
91
92
93      <*> 
94
95      <#> PICTURE TEST                                                                                                                                                                                                                                                     
96
97          Example 2,  A simple test pattern and projection data
98
99
100      <#> spec    mono   75
101          energy spectrum is monochromatic at energy level    75
102
103
104      <#> obje
105          description of objects
106                                                                 density at levels
107  numb type  x-coord  y-coord x-length y-length    angle  av dens        75
108
109     1 elip   0.0000   0.0000   0.9000   0.9000   0.0000   2.5000    2.5000
110
111     2 elip   0.0000   0.0000   0.8000   0.6000   0.0000  -1.2500   -1.2500
112
113     3 elip   0.0000   0.0000   0.2100   0.2500   0.0000   0.1000    0.1000
114
115          scale factor multiplying object densities     1.0000
116
117          seed set to 0
118          inhomogeneity set to     0.0000
119
120      <#> phan    aver    1
121
122          density in each pixel is obtained as the average of 1 x 1 points
123
124
125      <#> pixe       31    size        0.0700
126          picture size 31 x 31,  pixel size     0.0700
127
128          test picture read
129          Example 2,  A simple test pattern and projection data
130               0.000 seconds used for processing command pict
131
132
133      <*> 
134
135      <#> PROJECTION REAL                                                                                                                                                                                                                                                  
136
137          Example 2,  A simple test pattern and projection data
138
139
140      <#> spec    mono   75
141          energy spectrum is monochromatic at energy level    75
142
143
144      <#> obje
145          description of objects
146                                                                 density at levels
147  numb type  x-coord  y-coord x-length y-length    angle  av dens        75
148
149     1 elip   0.0000   0.0000   0.9000   0.9000   0.0000   2.5000    2.5000
150
151     2 elip   0.0000   0.0000   0.8000   0.6000   0.0000  -1.2500   -1.2500
152
153     3 elip   0.0000   0.0000   0.2100   0.2500   0.0000   0.1000    0.1000
154
155          scale factor multiplying object densities     1.0000
156
157          seed set to 0
158          inhomogeneity set to     0.0000
159
160      <#> rays    aver    1
161
162          projection data are calculated by dividing each ray interval into 1 substrips
163
164          with aperture (substrip) weights     1
165
166
167      <#> geom
168
169
170      <#> para unif line
171          rays are parallel with uniform spacing between rays
172          data collected along lines
173
174
175
176      <#> rays    user       63    spacing        0.0500
177          number of rays per projection    63
178          snark computed number of rays    63
179          at detector spacing     0.0500
180
181
182      <#> angl       30
183          total number of projections    30
184
185
186          projection angles    0.0000    6.0000   12.0000   18.0000   24.0000   30.0000   36.0000   42.0000   48.0000   54.0000
187                              60.0000   66.0000   72.0000   78.0000   84.0000   90.0000   96.0000  102.0000  108.0000  114.0000
188                             120.0000  126.0000  132.0000  138.0000  144.0000  150.0000  156.0000  162.0000  168.0000  174.0000
189
190
191      <#> meas    perf
192          projection data are noiseless
193
194
195      <#> back      0.0000
196                                at levels
197                                       75 
198           background absorption   0.0000 
199
200          estimate of totlen =      2823.757863
201          estimate of totden =      2677.896607
202          estimate of average density =     0.9483
203          projection data read
204          Example 2,  A simple test pattern and projection data
205               0.001 seconds used for processing command proj
206
207
208      <*> 
209
210      <#> MODE LOWER constraint set at 0.0  UPPER constraint set at 2.5                                                                                                                                                                                                    
211       lower constraint set to    0.0000
212       upper constraint set to    2.5000
213               0.000 seconds used for processing command mode
214
215
216      <*> 
217
218      <#> SELECT USER RAYSEQ                                                                                                                                                                                                                                               
219
220      <#> STEP projections by 10 and rays by 1
221          sequential ray selection with ray-step     1 and projection-step    10
222               0.000 seconds used for processing command sele
223
224
225      <*> 
226
227      <#> STOP TERMINATION VARIANCE at 0.01                                                                                                                                                                                                                                
228          termination test vari
229          epsilon =   0.010000
230               0.000 seconds used for processing command stop
231
232
233      <*> 
234
235      <#> EXECUTE AVERAGE ART SMOOTH                                                                                                                                                                                                                                       
236
237          Example 2a  underrelaxed art4 with variable tolerance
238
239      <#> smoothing threshold is 0.2  Weights are 2.25 1.5 1.0
240          reconstruction is smoothed after
241
242          iterations    1   2   3   4   5   6   7   8   9  10
243                       11  12  13  14  15  16  17  18  19  20
244                       21  22  23  24  25  26  27  28  29  30
245                       31  32  33  34  35  36  37  38  39  40
246                       41  42  43  44  45  46  47  48  49  50
247          last iteration
248
249          threshold =   0.20000      w1 =   2.25000      w2 =   1.50000      w3 =   1.00000
250
251      <#> ART4  RELAXATION CONSTANT = 0.025  TOLERANCE VARIABLE decreasing
252           art4 method
253           relaxation parameter is    0.0250
254           tolerance is variable
255
256      <#> CONSTRAINT BOUND CONRELAX CONSTANT = 0.5 NOMLIZE at each iteration
257           constraint set to boundaries
258           relax constraints with cr=    0.5000
259           picture is normalized to    0.9483 after each iteration
260            1890 rays are used for each iteration
261           lower constraint =    0.0000
262           upper constraint =    2.5000
263          algorithm executed in iteration    1
264               0.001 seconds for the execution of the algorithm
265          iteration    1 completed
266               0.001 seconds for this iteration
267          algorithm executed in iteration    2
268               0.001 seconds for the execution of the algorithm
269          iteration    2 completed
270               0.001 seconds for this iteration
271          algorithm executed in iteration    3
272               0.002 seconds for the execution of the algorithm
273          iteration    3 completed
274               0.002 seconds for this iteration
275          algorithm executed in iteration    4
276               0.001 seconds for the execution of the algorithm
277          iteration    4 completed
278               0.001 seconds for this iteration
279          algorithm executed in iteration    5
280               0.003 seconds for the execution of the algorithm
281          iteration    5 completed
282               0.003 seconds for this iteration
283          algorithm executed in iteration    6
284               0.002 seconds for the execution of the algorithm
285          iteration    6 completed
286               0.002 seconds for this iteration
287          algorithm executed in iteration    7
288               0.001 seconds for the execution of the algorithm
289          iteration    7 completed
290               0.001 seconds for this iteration
291          algorithm executed in iteration    8
292               0.001 seconds for the execution of the algorithm
293          iteration    8 completed
294               0.001 seconds for this iteration
295          algorithm executed in iteration    9
296               0.001 seconds for the execution of the algorithm
297          iteration    9 completed
298               0.001 seconds for this iteration
299          algorithm executed in iteration   10
300               0.002 seconds for the execution of the algorithm
301          iteration   10 completed
302               0.002 seconds for this iteration
303          algorithm executed in iteration   11
304               0.001 seconds for the execution of the algorithm
305          iteration   11 completed
306               0.001 seconds for this iteration
307          algorithm executed in iteration   12
308               0.002 seconds for the execution of the algorithm
309          iteration   12 completed
310               0.002 seconds for this iteration
311          algorithm executed in iteration   13
312               0.001 seconds for the execution of the algorithm
313          iteration   13 completed
314               0.001 seconds for this iteration
315          algorithm executed in iteration   14
316               0.001 seconds for the execution of the algorithm
317          iteration   14 completed
318               0.001 seconds for this iteration
319          algorithm executed in iteration   15
320               0.001 seconds for the execution of the algorithm
321          iteration   15 completed
322               0.001 seconds for this iteration
323          algorithm executed in iteration   16
324               0.000 seconds for the execution of the algorithm
325          iteration   16 completed
326               0.000 seconds for this iteration
327          algorithm executed in iteration   17
328               0.002 seconds for the execution of the algorithm
329          iteration   17 completed
330               0.002 seconds for this iteration
331          algorithm executed in iteration   18
332               0.001 seconds for the execution of the algorithm
333          iteration   18 completed
334               0.001 seconds for this iteration
335          algorithm executed in iteration   19
336               0.001 seconds for the execution of the algorithm
337          iteration   19 completed
338               0.001 seconds for this iteration
339          algorithm executed in iteration   20
340               0.001 seconds for the execution of the algorithm
341          iteration   20 completed
342               0.001 seconds for this iteration
343          algorithm executed in iteration   21
344               0.002 seconds for the execution of the algorithm
345           iterative process stops at iteration    21
346           the change in variance is less than   0.009852 of the variance
347          reconstruction completed after iteration   21
348               0.002 seconds for this iteration
349               0.031 seconds for all iterations
350               0.031 seconds used for processing command exec
351
352
353      <*> 
354
355      <#> STOP ITERATIONS = 30                                                                                                                                                                                                                                             
356            30 iterations
357               0.000 seconds used for processing command stop
358
359
360      <*> 
361
362      <#> EXECUTE ALP1                                                                                                                                                                                                                                                     
363
364          Example 2b  underrelaxed art with extrapolation
365          algorithm executed in iteration    1
366               0.001 seconds for the execution of the algorithm
367          iteration    1 completed
368               0.001 seconds for this iteration
369          algorithm executed in iteration    2
370               0.001 seconds for the execution of the algorithm
371          iteration    2 completed
372               0.001 seconds for this iteration
373  EXTRAPOLATION PARAMETER IN EXTRAP IS   5.062718796778e-01
374          algorithm executed in iteration    3
375               0.001 seconds for the execution of the algorithm
376          iteration    3 completed
377               0.001 seconds for this iteration
378  EXTRAPOLATION PARAMETER IN EXTRAP IS   5.062718796778e-01
379          algorithm executed in iteration    4
380               0.000 seconds for the execution of the algorithm
381          iteration    4 completed
382               0.000 seconds for this iteration
383  EXTRAPOLATION PARAMETER IN EXTRAP IS   6.560964044719e-01
384          algorithm executed in iteration    5
385               0.001 seconds for the execution of the algorithm
386          iteration    5 completed
387               0.001 seconds for this iteration
388  EXTRAPOLATION PARAMETER IN EXTRAP IS   6.560964044719e-01
389          algorithm executed in iteration    6
390               0.000 seconds for the execution of the algorithm
391          iteration    6 completed
392               0.000 seconds for this iteration
393  EXTRAPOLATION PARAMETER IN EXTRAP IS   7.480594958675e-01
394          algorithm executed in iteration    7
395               0.001 seconds for the execution of the algorithm
396          iteration    7 completed
397               0.001 seconds for this iteration
398  EXTRAPOLATION PARAMETER IN EXTRAP IS   7.480594958675e-01
399          algorithm executed in iteration    8
400               0.000 seconds for the execution of the algorithm
401          iteration    8 completed
402               0.000 seconds for this iteration
403  EXTRAPOLATION PARAMETER IN EXTRAP IS   7.951514183812e-01
404          algorithm executed in iteration    9
405               0.001 seconds for the execution of the algorithm
406          iteration    9 completed
407               0.001 seconds for this iteration
408  EXTRAPOLATION PARAMETER IN EXTRAP IS   7.951514183812e-01
409          algorithm executed in iteration   10
410               0.000 seconds for the execution of the algorithm
411          iteration   10 completed
412               0.000 seconds for this iteration
413  EXTRAPOLATION PARAMETER IN EXTRAP IS   8.250655739411e-01
414          algorithm executed in iteration   11
415               0.002 seconds for the execution of the algorithm
416          iteration   11 completed
417               0.002 seconds for this iteration
418  EXTRAPOLATION PARAMETER IN EXTRAP IS   8.250655739411e-01
419          algorithm executed in iteration   12
420               0.000 seconds for the execution of the algorithm
421          iteration   12 completed
422               0.000 seconds for this iteration
423  EXTRAPOLATION PARAMETER IN EXTRAP IS   8.467169748437e-01
424          algorithm executed in iteration   13
425               0.001 seconds for the execution of the algorithm
426          iteration   13 completed
427               0.001 seconds for this iteration
428  EXTRAPOLATION PARAMETER IN EXTRAP IS   8.467169748437e-01
429          algorithm executed in iteration   14
430               0.000 seconds for the execution of the algorithm
431          iteration   14 completed
432               0.000 seconds for this iteration
433  EXTRAPOLATION PARAMETER IN EXTRAP IS   8.628058734349e-01
434          algorithm executed in iteration   15
435               0.002 seconds for the execution of the algorithm
436          iteration   15 completed
437               0.002 seconds for this iteration
438  EXTRAPOLATION PARAMETER IN EXTRAP IS   8.628058734349e-01
439          algorithm executed in iteration   16
440               0.000 seconds for the execution of the algorithm
441          iteration   16 completed
442               0.000 seconds for this iteration
443  EXTRAPOLATION PARAMETER IN EXTRAP IS   8.749346145383e-01
444          algorithm executed in iteration   17
445               0.001 seconds for the execution of the algorithm
446          iteration   17 completed
447               0.001 seconds for this iteration
448  EXTRAPOLATION PARAMETER IN EXTRAP IS   8.749346145383e-01
449          algorithm executed in iteration   18
450               0.000 seconds for the execution of the algorithm
451          iteration   18 completed
452               0.000 seconds for this iteration
453  EXTRAPOLATION PARAMETER IN EXTRAP IS   8.843623393369e-01
454          algorithm executed in iteration   19
455               0.001 seconds for the execution of the algorithm
456          iteration   19 completed
457               0.001 seconds for this iteration
458  EXTRAPOLATION PARAMETER IN EXTRAP IS   8.843623393369e-01
459          algorithm executed in iteration   20
460               0.000 seconds for the execution of the algorithm
461          iteration   20 completed
462               0.000 seconds for this iteration
463  EXTRAPOLATION PARAMETER IN EXTRAP IS   8.920147280181e-01
464          algorithm executed in iteration   21
465               0.001 seconds for the execution of the algorithm
466          iteration   21 completed
467               0.001 seconds for this iteration
468  EXTRAPOLATION PARAMETER IN EXTRAP IS   8.920147280181e-01
469          algorithm executed in iteration   22
470               0.000 seconds for the execution of the algorithm
471          iteration   22 completed
472               0.000 seconds for this iteration
473  EXTRAPOLATION PARAMETER IN EXTRAP IS   8.985122210336e-01
474          algorithm executed in iteration   23
475               0.001 seconds for the execution of the algorithm
476          iteration   23 completed
477               0.001 seconds for this iteration
478  EXTRAPOLATION PARAMETER IN EXTRAP IS   8.985122210336e-01
479          algorithm executed in iteration   24
480               0.000 seconds for the execution of the algorithm
481          iteration   24 completed
482               0.000 seconds for this iteration
483  EXTRAPOLATION PARAMETER IN EXTRAP IS   9.042443038619e-01
484          algorithm executed in iteration   25
485               0.001 seconds for the execution of the algorithm
486          iteration   25 completed
487               0.001 seconds for this iteration
488  EXTRAPOLATION PARAMETER IN EXTRAP IS   9.042443038619e-01
489          algorithm executed in iteration   26
490               0.000 seconds for the execution of the algorithm
491          iteration   26 completed
492               0.000 seconds for this iteration
493  EXTRAPOLATION PARAMETER IN EXTRAP IS   9.094448804097e-01
494          algorithm executed in iteration   27
495               0.001 seconds for the execution of the algorithm
496          iteration   27 completed
497               0.001 seconds for this iteration
498  EXTRAPOLATION PARAMETER IN EXTRAP IS   9.094448804097e-01
499          algorithm executed in iteration   28
500               0.000 seconds for the execution of the algorithm
501          iteration   28 completed
502               0.000 seconds for this iteration
503  EXTRAPOLATION PARAMETER IN EXTRAP IS   9.142500229503e-01
504          algorithm executed in iteration   29
505               0.001 seconds for the execution of the algorithm
506          iteration   29 completed
507               0.001 seconds for this iteration
508  EXTRAPOLATION PARAMETER IN EXTRAP IS   9.142500229503e-01
509          algorithm executed in iteration   30
510               0.000 seconds for the execution of the algorithm
511          reconstruction completed after iteration   30
512               0.000 seconds for this iteration
513               0.017 seconds for all iterations
514               0.018 seconds used for processing command exec
515
516
517      <*> 
518
519      <*>  Reset mode to unconstrained
520
521      <#> MODE                                                                                                                                                                                                                                                             
522               0.000 seconds used for processing command mode
523
524
525      <*> 
526
527      <#> EVALUATE RESOLUTION                                                                                                                                                                                                                                              
528
529          Example 2a  underrelaxed art4
530
531      <#> WHOLEPIC
532  Region        cx        cy         u         v        ang      t1        t2
533
534  wholepic                                                        -1e+20     1e+20
535
536          iterations    1   2   3   4   5   6   7   8   9  10
537                       11  12  13  14  15  16  17  18  19  20
538                       21  22  23  24  25  26  27  28  29  30
539                       31  32  33  34  35  36  37  38  39
540          last iteration
541
542     residual calculation using line integration
543               0.047 seconds used for processing command eval
544
545
546      <*> 
547
548      <#> LINES COLUMNS 15 and 18 are to be printed                                                                                                                                                                                                                        
549
550          iterations   10  20  30
551          last iteration
552
553      lines display of reconstruction using ART iter 10
554
555      phantom name:    Example 2,  A simple test pattern and projection data
556      projection data: Example 2,  A simple test pattern and projection data
557      execution name:  Example 2a  underrelaxed art4 with variable tolerance
558      scaling factor =    1.00000
559
560  column =              15                               18
561  row    original   estimate difference   original   estimate difference 
562    0     0.00000    0.18210    0.18210    0.00000    0.00909    0.00909
563    1     0.00000    0.20450    0.20450    0.00000    0.24753    0.24753
564    2     0.00000    1.10539    1.10539    0.00000    0.52569    0.52569
565    3     2.50000    2.13754   -0.36246    2.50000    1.94767   -0.55233
566    4     2.50000    2.16437   -0.33563    2.50000    2.18812   -0.31188
567    5     2.50000    2.16394   -0.33606    2.50000    2.20758   -0.29242
568    6     2.50000    2.13317   -0.36683    2.50000    2.17610   -0.32390
569    7     1.25000    1.49667    0.24667    1.25000    1.68287    0.43287
570    8     1.25000    1.44121    0.19121    1.25000    1.41529    0.16529
571    9     1.25000    1.37200    0.12200    1.25000    1.38374    0.13374
572   10     1.25000    1.34328    0.09328    1.25000    1.34276    0.09276
573   11     1.25000    1.33396    0.08396    1.25000    1.32222    0.07222
574   12     1.35000    1.34036   -0.00964    1.25000    1.31411    0.06411
575   13     1.35000    1.34676   -0.00324    1.25000    1.31521    0.06521
576   14     1.35000    1.35058    0.00058    1.25000    1.31568    0.06568
577   15     1.35000    1.35053    0.00053    1.35000    1.31710   -0.03290
578   16     1.35000    1.35058    0.00058    1.25000    1.31607    0.06607
579   17     1.35000    1.34676   -0.00324    1.25000    1.31587    0.06587
580   18     1.35000    1.34035   -0.00965    1.25000    1.31417    0.06417
581   19     1.25000    1.33395    0.08395    1.25000    1.32131    0.07131
582   20     1.25000    1.34326    0.09326    1.25000    1.34021    0.09021
583   21     1.25000    1.37196    0.12196    1.25000    1.38008    0.13008
584   22     1.25000    1.44114    0.19114    1.25000    1.41090    0.16090
585   23     1.25000    1.49658    0.24658    1.25000    1.68040    0.43040
586   24     2.50000    2.13333   -0.36667    2.50000    2.17405   -0.32595
587   25     2.50000    2.16414   -0.33586    2.50000    2.20535   -0.29465
588   26     2.50000    2.16455   -0.33545    2.50000    2.18655   -0.31345
589   27     2.50000    2.13774   -0.36226    2.50000    1.94456   -0.55544
590   28     0.00000    1.10534    1.10534    0.00000    0.44393    0.44393
591   29     0.00000    0.20439    0.20439    0.00000    0.29112    0.29112
592   30     0.00000    0.18198    0.18198    0.00000    0.01326    0.01326
593
594
595      lines display of reconstruction using ART iter 20
596
597      phantom name:    Example 2,  A simple test pattern and projection data
598      projection data: Example 2,  A simple test pattern and projection data
599      execution name:  Example 2a  underrelaxed art4 with variable tolerance
600      scaling factor =    1.00000
601
602  column =              15                               18
603  row    original   estimate difference   original   estimate difference 
604    0     0.00000    0.00164    0.00164    0.00000    0.02169    0.02169
605    1     0.00000    0.00132    0.00132    0.00000    0.04480    0.04480
606    2     0.00000    0.94940    0.94940    0.00000    0.13226    0.13226
607    3     2.50000    2.43151   -0.06849    2.50000    2.25420   -0.24580
608    4     2.50000    2.43097   -0.06903    2.50000    2.44794   -0.05206
609    5     2.50000    2.37936   -0.12064    2.50000    2.44383   -0.05617
610    6     2.50000    2.31917   -0.18083    2.50000    2.39051   -0.10949
611    7     1.25000    1.46751    0.21751    1.25000    1.65881    0.40881
612    8     1.25000    1.31935    0.06935    1.25000    1.33707    0.08707
613    9     1.25000    1.29770    0.04770    1.25000    1.31051    0.06051
614   10     1.25000    1.27744    0.02744    1.25000    1.27762    0.02762
615   11     1.25000    1.28043    0.03043    1.25000    1.26720    0.01720
616   12     1.35000    1.30399   -0.04601    1.25000    1.26843    0.01843
617   13     1.35000    1.32095   -0.02905    1.25000    1.27985    0.02985
618   14     1.35000    1.32945   -0.02055    1.25000    1.28645    0.03645
619   15     1.35000    1.32847   -0.02153    1.35000    1.29098   -0.05902
620   16     1.35000    1.32944   -0.02056    1.25000    1.28714    0.03714
621   17     1.35000    1.32093   -0.02907    1.25000    1.28136    0.03136
622   18     1.35000    1.30399   -0.04601    1.25000    1.27000    0.02000
623   19     1.25000    1.28044    0.03044    1.25000    1.26871    0.01871
624   20     1.25000    1.27744    0.02744    1.25000    1.27815    0.02815
625   21     1.25000    1.29767    0.04767    1.25000    1.30973    0.05973
626   22     1.25000    1.31930    0.06930    1.25000    1.33504    0.08504
627   23     1.25000    1.46737    0.21737    1.25000    1.71927    0.46927
628   24     2.50000    2.31941   -0.18059    2.50000    2.39181   -0.10819
629   25     2.50000    2.37969   -0.12031    2.50000    2.44490   -0.05510
630   26     2.50000    2.43125   -0.06875    2.50000    2.44917   -0.05083
631   27     2.50000    2.43179   -0.06821    2.50000    2.25444   -0.24556
632   28     0.00000    0.94928    0.94928    0.00000    0.12809    0.12809
633   29     0.00000    0.00158    0.00158    0.00000    0.04576    0.04576
634   30     0.00000    0.00199    0.00199    0.00000    0.02628    0.02628
635
636
637      lines display of reconstruction using ART iter 21
638
639      phantom name:    Example 2,  A simple test pattern and projection data
640      projection data: Example 2,  A simple test pattern and projection data
641      execution name:  Example 2a  underrelaxed art4 with variable tolerance
642      scaling factor =    1.00000
643
644  column =              15                               18
645  row    original   estimate difference   original   estimate difference 
646    0     0.00000    0.00100    0.00100    0.00000    0.01742    0.01742
647    1     0.00000    0.00070    0.00070    0.00000    0.03672    0.03672
648    2     0.00000    0.93506    0.93506    0.00000    0.11759    0.11759
649    3     2.50000    2.44949   -0.05051    2.50000    2.27648   -0.22352
650    4     2.50000    2.44615   -0.05385    2.50000    2.45838   -0.04162
651    5     2.50000    2.39029   -0.10971    2.50000    2.45205   -0.04795
652    6     2.50000    2.32845   -0.17155    2.50000    2.39929   -0.10071
653    7     1.25000    1.46233    0.21233    1.25000    1.65516    0.40516
654    8     1.25000    1.31367    0.06367    1.25000    1.33085    0.08085
655    9     1.25000    1.29279    0.04279    1.25000    1.30546    0.05546
656   10     1.25000    1.27400    0.02400    1.25000    1.27412    0.02412
657   11     1.25000    1.27844    0.02844    1.25000    1.26496    0.01496
658   12     1.35000    1.30374   -0.04626    1.25000    1.26714    0.01714
659   13     1.35000    1.32159   -0.02841    1.25000    1.27948    0.02948
660   14     1.35000    1.33036   -0.01964    1.25000    1.28660    0.03660
661   15     1.35000    1.32919   -0.02081    1.35000    1.29140   -0.05860
662   16     1.35000    1.33034   -0.01966    1.25000    1.28726    0.03726
663   17     1.35000    1.32157   -0.02843    1.25000    1.28096    0.03096
664   18     1.35000    1.30372   -0.04628    1.25000    1.26874    0.01874
665   19     1.25000    1.27844    0.02844    1.25000    1.26657    0.01657
666   20     1.25000    1.27400    0.02400    1.25000    1.27479    0.02479
667   21     1.25000    1.29277    0.04277    1.25000    1.30483    0.05483
668   22     1.25000    1.31362    0.06362    1.25000    1.32898    0.07898
669   23     1.25000    1.46218    0.21218    1.25000    1.65781    0.40781
670   24     2.50000    2.32869   -0.17131    2.50000    2.40065   -0.09935
671   25     2.50000    2.39063   -0.10937    2.50000    2.45318   -0.04682
672   26     2.50000    2.44644   -0.05356    2.50000    2.45967   -0.04033
673   27     2.50000    2.44977   -0.05023    2.50000    2.27692   -0.22308
674   28     0.00000    0.93489    0.93489    0.00000    0.11431    0.11431
675   29     0.00000    0.00120    0.00120    0.00000    0.03840    0.03840
676   30     0.00000    0.00161    0.00161    0.00000    0.02249    0.02249
677
678
679      lines display of reconstruction using ALP1 iter 10
680
681      phantom name:    Example 2,  A simple test pattern and projection data
682      projection data: Example 2,  A simple test pattern and projection data
683      execution name:  Example 2b  underrelaxed art with extrapolation
684      scaling factor =    1.00000
685
686  column =              15                               18
687  row    original   estimate difference   original   estimate difference 
688    0     0.00000    0.12670    0.12670    0.00000   -0.02567   -0.02567
689    1     0.00000    0.17187    0.17187    0.00000    0.16958    0.16958
690    2     0.00000    1.20151    1.20151    0.00000    0.48350    0.48350
691    3     2.50000    2.23551   -0.26449    2.50000    2.01480   -0.48520
692    4     2.50000    2.39456   -0.10544    2.50000    2.31508   -0.18492
693    5     2.50000    2.34966   -0.15034    2.50000    2.46624   -0.03376
694    6     2.50000    2.31124   -0.18876    2.50000    2.29509   -0.20491
695    7     1.25000    1.63106    0.38106    1.25000    1.81482    0.56482
696    8     1.25000    1.35183    0.10183    1.25000    1.44203    0.19203
697    9     1.25000    1.42885    0.17885    1.25000    1.45015    0.20015
698   10     1.25000    1.31859    0.06859    1.25000    1.36973    0.11973
699   11     1.25000    1.42731    0.17731    1.25000    1.34954    0.09954
700   12     1.35000    1.31939   -0.03061    1.25000    1.31390    0.06390
701   13     1.35000    1.37323    0.02323    1.25000    1.32415    0.07415
702   14     1.35000    1.46590    0.11590    1.25000    1.37584    0.12584
703   15     1.35000    1.13749   -0.21251    1.35000    1.25214   -0.09786
704   16     1.35000    1.46594    0.11594    1.25000    1.36902    0.11902
705   17     1.35000    1.37317    0.02317    1.25000    1.30699    0.05699
706   18     1.35000    1.31944   -0.03056    1.25000    1.34043    0.09043
707   19     1.25000    1.42727    0.17727    1.25000    1.33982    0.08982
708   20     1.25000    1.31861    0.06861    1.25000    1.35654    0.10654
709   21     1.25000    1.42880    0.17880    1.25000    1.44218    0.19218
710   22     1.25000    1.35172    0.10172    1.25000    1.45375    0.20375
711   23     1.25000    1.63100    0.38100    1.25000    1.80751    0.55751
712   24     2.50000    2.31161   -0.18839    2.50000    2.27967   -0.22033
713   25     2.50000    2.34937   -0.15063    2.50000    2.50118    0.00118
714   26     2.50000    2.39545   -0.10455    2.50000    2.30942   -0.19058
715   27     2.50000    2.23526   -0.26474    2.50000    1.99542   -0.50458
716   28     0.00000    1.20133    1.20133    0.00000    0.49478    0.49478
717   29     0.00000    0.17173    0.17173    0.00000    0.15734    0.15734
718   30     0.00000    0.12647    0.12647    0.00000   -0.02832   -0.02832
719
720
721      lines display of reconstruction using ALP1 iter 20
722
723      phantom name:    Example 2,  A simple test pattern and projection data
724      projection data: Example 2,  A simple test pattern and projection data
725      execution name:  Example 2b  underrelaxed art with extrapolation
726      scaling factor =    1.00000
727
728  column =              15                               18
729  row    original   estimate difference   original   estimate difference 
730    0     0.00000   -0.09663   -0.09663    0.00000   -0.18986   -0.18986
731    1     0.00000   -0.12331   -0.12331    0.00000   -0.05402   -0.05402
732    2     0.00000    1.04579    1.04579    0.00000    0.12373    0.12373
733    3     2.50000    2.51540    0.01540    2.50000    2.25540   -0.24460
734    4     2.50000    2.59753    0.09753    2.50000    2.56869    0.06869
735    5     2.50000    2.55529    0.05529    2.50000    2.67846    0.17846
736    6     2.50000    2.40841   -0.09159    2.50000    2.45046   -0.04954
737    7     1.25000    1.52857    0.27857    1.25000    1.73851    0.48851
738    8     1.25000    1.19882   -0.05118    1.25000    1.29003    0.04003
739    9     1.25000    1.28220    0.03220    1.25000    1.30029    0.05029
740   10     1.25000    1.18525   -0.06475    1.25000    1.24281   -0.00719
741   11     1.25000    1.27819    0.02819    1.25000    1.24380   -0.00620
742   12     1.35000    1.25011   -0.09989    1.25000    1.22237   -0.02763
743   13     1.35000    1.30966   -0.04034    1.25000    1.24476   -0.00524
744   14     1.35000    1.38630    0.03630    1.25000    1.28583    0.03583
745   15     1.35000    1.11801   -0.23199    1.35000    1.20027   -0.14973
746   16     1.35000    1.38636    0.03636    1.25000    1.28050    0.03050
747   17     1.35000    1.30959   -0.04041    1.25000    1.23298   -0.01702
748   18     1.35000    1.25014   -0.09986    1.25000    1.24848   -0.00152
749   19     1.25000    1.27817    0.02817    1.25000    1.23597   -0.01403
750   20     1.25000    1.18526   -0.06474    1.25000    1.23189   -0.01811
751   21     1.25000    1.28211    0.03211    1.25000    1.29998    0.04998
752   22     1.25000    1.19865   -0.05135    1.25000    1.30195    0.05195
753   23     1.25000    1.52839    0.27839    1.25000    1.73195    0.48195
754   24     2.50000    2.40899   -0.09101    2.50000    2.44156   -0.05844
755   25     2.50000    2.55477    0.05477    2.50000    2.70939    0.20939
756   26     2.50000    2.59890    0.09890    2.50000    2.56774    0.06774
757   27     2.50000    2.51498    0.01498    2.50000    2.24565   -0.25435
758   28     0.00000    1.04557    1.04557    0.00000    0.14374    0.14374
759   29     0.00000   -0.12338   -0.12338    0.00000   -0.05331   -0.05331
760   30     0.00000   -0.09670   -0.09670    0.00000   -0.17428   -0.17428
761
762
763      lines display of reconstruction using ALP1 iter 30
764
765      phantom name:    Example 2,  A simple test pattern and projection data
766      projection data: Example 2,  A simple test pattern and projection data
767      execution name:  Example 2b  underrelaxed art with extrapolation
768      scaling factor =    1.00000
769
770  column =              15                               18
771  row    original   estimate difference   original   estimate difference 
772    0     0.00000   -0.10295   -0.10295    0.00000   -0.15511   -0.15511
773    1     0.00000   -0.21068   -0.21068    0.00000   -0.09236   -0.09236
774    2     0.00000    0.92855    0.92855    0.00000   -0.06300   -0.06300
775    3     2.50000    2.64511    0.14511    2.50000    2.37419   -0.12581
776    4     2.50000    2.61843    0.11843    2.50000    2.63291    0.13291
777    5     2.50000    2.61277    0.11277    2.50000    2.68002    0.18002
778    6     2.50000    2.39593   -0.10407    2.50000    2.48350   -0.01650
779    7     1.25000    1.47574    0.22574    1.25000    1.66721    0.41721
780    8     1.25000    1.11475   -0.13525    1.25000    1.19787   -0.05213
781    9     1.25000    1.24880   -0.00120    1.25000    1.25371    0.00371
782   10     1.25000    1.15550   -0.09450    1.25000    1.22544   -0.02456
783   11     1.25000    1.25545    0.00545    1.25000    1.24429   -0.00571
784   12     1.35000    1.26430   -0.08570    1.25000    1.22826   -0.02174
785   13     1.35000    1.32865   -0.02135    1.25000    1.25347    0.00347
786   14     1.35000    1.40845    0.05845    1.25000    1.28571    0.03571
787   15     1.35000    1.14519   -0.20481    1.35000    1.22750   -0.12250
788   16     1.35000    1.40851    0.05851    1.25000    1.28178    0.03178
789   17     1.35000    1.32857   -0.02143    1.25000    1.24374   -0.00626
790   18     1.35000    1.26433   -0.08567    1.25000    1.25314    0.00314
791   19     1.25000    1.25545    0.00545    1.25000    1.23609   -0.01391
792   20     1.25000    1.15550   -0.09450    1.25000    1.21469   -0.03531
793   21     1.25000    1.24869   -0.00131    1.25000    1.25626    0.00626
794   22     1.25000    1.11456   -0.13544    1.25000    1.20898   -0.04102
795   23     1.25000    1.47547    0.22547    1.25000    1.65907    0.40907
796   24     2.50000    2.39664   -0.10336    2.50000    2.47619   -0.02381
797   25     2.50000    2.61204    0.11204    2.50000    2.70609    0.20609
798   26     2.50000    2.62009    0.12009    2.50000    2.63048    0.13048
799   27     2.50000    2.64455    0.14455    2.50000    2.36445   -0.13555
800   28     0.00000    0.92830    0.92830    0.00000   -0.04433   -0.04433
801   29     0.00000   -0.21071   -0.21071    0.00000   -0.09436   -0.09436
802   30     0.00000   -0.10286   -0.10286    0.00000   -0.14676   -0.14676
803
804
805  
806               0.002 seconds used for processing command line
807
808
809      <*> 
810
811      <#> END