1 snark14UserDefined.s170426 - A PICTURE RECONSTRUCTION PROGRAM
6 <*> Creation of a small test pattern and projection data with execution of
8 <*> a user-defined reconstruction algorithm and analysis.
14 Example 2, A simple test pattern and projection data
17 <#> SPECTRUM MONOCHROMATIC with energy at 75 KeV
18 energy spectrum is monochromatic at energy level 75
21 <#> OBJECTS cx cy u v ang den
22 description of objects
24 numb type x-coord y-coord x-length y-length angle av dens 75
26 1 elip 0.0000 0.0000 0.9000 0.9000 0.0000 2.5000 2.5000
28 2 elip 0.0000 0.0000 0.8000 0.6000 0.0000 -1.2500 -1.2500
30 3 elip 0.0000 0.0000 0.2100 0.2500 0.0000 0.1000 0.1000
32 scale factor multiplying object densities 1.0000
35 inhomogeneity set to 0.0000
37 <#> PHANTOM AVERAGE 1 point for estimating pixel density
39 this run will generate a phantom
40 density in each pixel is obtained as the average of 1 x 1 points
43 <#> number of elements is 31 Size of a pixel is 0.07
44 picture size 31 x 31, pixel size 0.0700
47 <#> RAYSUM AVERAGE 1 line for estimation of ray sum with weight
49 this run will generate projection data
50 projection data are calculated by dividing each ray interval into 1 substrips
52 with aperture (substrip) weights 1
58 <#> PARALLEL UNIFORM LINE
59 rays are parallel with uniform spacing between rays
60 data collected along lines
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
69 <#> ANGLES 30 EQUAL SPACING between the first and last angles which are
70 total number of projections 30
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
77 <#> MEASUREMENT PERFECT
78 projection data are noiseless
81 <#> BACKGROUND absorption is 0.0
84 background absorption 0.0000
88 0.004 seconds phantom creation
89 0.000 seconds projection data creation
90 0.004 seconds used for processing command crea
97 Example 2, A simple test pattern and projection data
101 energy spectrum is monochromatic at energy level 75
105 description of objects
107 numb type x-coord y-coord x-length y-length angle av dens 75
109 1 elip 0.0000 0.0000 0.9000 0.9000 0.0000 2.5000 2.5000
111 2 elip 0.0000 0.0000 0.8000 0.6000 0.0000 -1.2500 -1.2500
113 3 elip 0.0000 0.0000 0.2100 0.2500 0.0000 0.1000 0.1000
115 scale factor multiplying object densities 1.0000
118 inhomogeneity set to 0.0000
122 density in each pixel is obtained as the average of 1 x 1 points
125 <#> pixe 31 size 0.0700
126 picture size 31 x 31, pixel size 0.0700
129 Example 2, A simple test pattern and projection data
130 0.000 seconds used for processing command pict
137 Example 2, A simple test pattern and projection data
141 energy spectrum is monochromatic at energy level 75
145 description of objects
147 numb type x-coord y-coord x-length y-length angle av dens 75
149 1 elip 0.0000 0.0000 0.9000 0.9000 0.0000 2.5000 2.5000
151 2 elip 0.0000 0.0000 0.8000 0.6000 0.0000 -1.2500 -1.2500
153 3 elip 0.0000 0.0000 0.2100 0.2500 0.0000 0.1000 0.1000
155 scale factor multiplying object densities 1.0000
158 inhomogeneity set to 0.0000
162 projection data are calculated by dividing each ray interval into 1 substrips
164 with aperture (substrip) weights 1
171 rays are parallel with uniform spacing between rays
172 data collected along lines
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
183 total number of projections 30
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
192 projection data are noiseless
198 background absorption 0.0000
200 estimate of totlen = 2823.757863
201 estimate of totden = 2677.896607
202 estimate of average density = 0.9483
204 Example 2, A simple test pattern and projection data
205 0.001 seconds used for processing command proj
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
218 <#> SELECT USER RAYSEQ
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
227 <#> STOP TERMINATION VARIANCE at 0.01
228 termination test vari
230 0.000 seconds used for processing command stop
235 <#> EXECUTE AVERAGE ART SMOOTH
237 Example 2a underrelaxed art4 with variable tolerance
239 <#> smoothing threshold is 0.2 Weights are 2.25 1.5 1.0
240 reconstruction is smoothed after
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
249 threshold = 0.20000 w1 = 2.25000 w2 = 1.50000 w3 = 1.00000
251 <#> ART4 RELAXATION CONSTANT = 0.025 TOLERANCE VARIABLE decreasing
253 relaxation parameter is 0.0250
254 tolerance is variable
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
355 <#> STOP ITERATIONS = 30
357 0.000 seconds used for processing command stop
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
519 <*> Reset mode to unconstrained
522 0.000 seconds used for processing command mode
527 <#> EVALUATE RESOLUTION
529 Example 2a underrelaxed art4
532 Region cx cy u v ang t1 t2
534 wholepic -1e+20 1e+20
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
542 residual calculation using line integration
543 0.047 seconds used for processing command eval
548 <#> LINES COLUMNS 15 and 18 are to be printed
553 lines display of reconstruction using ART iter 10
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
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
595 lines display of reconstruction using ART iter 20
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
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
637 lines display of reconstruction using ART iter 21
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
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
679 lines display of reconstruction using ALP1 iter 10
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
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
721 lines display of reconstruction using ALP1 iter 20
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
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
763 lines display of reconstruction using ALP1 iter 30
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
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
806 0.002 seconds used for processing command line