snark14UserDefined.s170426 - A PICTURE RECONSTRUCTION PROGRAM <*> *** EXAMPLE 2 <*> Creation of a small test pattern and projection data with execution of <*> a user-defined reconstruction algorithm and analysis. <*> <#> CREATE Example 2, A simple test pattern and projection data <#> SPECTRUM MONOCHROMATIC with energy at 75 KeV energy spectrum is monochromatic at energy level 75 <#> OBJECTS cx cy u v ang den description of objects density at levels numb type x-coord y-coord x-length y-length angle av dens 75 1 elip 0.0000 0.0000 0.9000 0.9000 0.0000 2.5000 2.5000 2 elip 0.0000 0.0000 0.8000 0.6000 0.0000 -1.2500 -1.2500 3 elip 0.0000 0.0000 0.2100 0.2500 0.0000 0.1000 0.1000 scale factor multiplying object densities 1.0000 seed set to 0 inhomogeneity set to 0.0000 <#> PHANTOM AVERAGE 1 point for estimating pixel density this run will generate a phantom density in each pixel is obtained as the average of 1 x 1 points <#> number of elements is 31 Size of a pixel is 0.07 picture size 31 x 31, pixel size 0.0700 <#> RAYSUM AVERAGE 1 line for estimation of ray sum with weight this run will generate projection data projection data are calculated by dividing each ray interval into 1 substrips with aperture (substrip) weights 1 <#> GEOMETRY <#> PARALLEL UNIFORM LINE rays are parallel with uniform spacing between rays data collected along lines <#> RAYS USER-defined 63 rays per projection at spacing 0.05 number of rays per projection 63 at detector spacing 0.0500 <#> ANGLES 30 EQUAL SPACING between the first and last angles which are total number of projections 30 projection angles 0.0000 6.0000 12.0000 18.0000 24.0000 30.0000 36.0000 42.0000 48.0000 54.0000 60.0000 66.0000 72.0000 78.0000 84.0000 90.0000 96.0000 102.0000 108.0000 114.0000 120.0000 126.0000 132.0000 138.0000 144.0000 150.0000 156.0000 162.0000 168.0000 174.0000 <#> MEASUREMENT PERFECT projection data are noiseless <#> BACKGROUND absorption is 0.0 at levels 75 background absorption 0.0000 <#> RUN 0.004 seconds phantom creation 0.000 seconds projection data creation 0.004 seconds used for processing command crea <*> <#> PICTURE TEST Example 2, A simple test pattern and projection data <#> spec mono 75 energy spectrum is monochromatic at energy level 75 <#> obje description of objects density at levels numb type x-coord y-coord x-length y-length angle av dens 75 1 elip 0.0000 0.0000 0.9000 0.9000 0.0000 2.5000 2.5000 2 elip 0.0000 0.0000 0.8000 0.6000 0.0000 -1.2500 -1.2500 3 elip 0.0000 0.0000 0.2100 0.2500 0.0000 0.1000 0.1000 scale factor multiplying object densities 1.0000 seed set to 0 inhomogeneity set to 0.0000 <#> phan aver 1 density in each pixel is obtained as the average of 1 x 1 points <#> pixe 31 size 0.0700 picture size 31 x 31, pixel size 0.0700 test picture read Example 2, A simple test pattern and projection data 0.000 seconds used for processing command pict <*> <#> PROJECTION REAL Example 2, A simple test pattern and projection data <#> spec mono 75 energy spectrum is monochromatic at energy level 75 <#> obje description of objects density at levels numb type x-coord y-coord x-length y-length angle av dens 75 1 elip 0.0000 0.0000 0.9000 0.9000 0.0000 2.5000 2.5000 2 elip 0.0000 0.0000 0.8000 0.6000 0.0000 -1.2500 -1.2500 3 elip 0.0000 0.0000 0.2100 0.2500 0.0000 0.1000 0.1000 scale factor multiplying object densities 1.0000 seed set to 0 inhomogeneity set to 0.0000 <#> rays aver 1 projection data are calculated by dividing each ray interval into 1 substrips with aperture (substrip) weights 1 <#> geom <#> para unif line rays are parallel with uniform spacing between rays data collected along lines <#> rays user 63 spacing 0.0500 number of rays per projection 63 snark computed number of rays 63 at detector spacing 0.0500 <#> angl 30 total number of projections 30 projection angles 0.0000 6.0000 12.0000 18.0000 24.0000 30.0000 36.0000 42.0000 48.0000 54.0000 60.0000 66.0000 72.0000 78.0000 84.0000 90.0000 96.0000 102.0000 108.0000 114.0000 120.0000 126.0000 132.0000 138.0000 144.0000 150.0000 156.0000 162.0000 168.0000 174.0000 <#> meas perf projection data are noiseless <#> back 0.0000 at levels 75 background absorption 0.0000 estimate of totlen = 2823.757863 estimate of totden = 2677.896607 estimate of average density = 0.9483 projection data read Example 2, A simple test pattern and projection data 0.001 seconds used for processing command proj <*> <#> MODE LOWER constraint set at 0.0 UPPER constraint set at 2.5 lower constraint set to 0.0000 upper constraint set to 2.5000 0.000 seconds used for processing command mode <*> <#> SELECT USER RAYSEQ <#> STEP projections by 10 and rays by 1 sequential ray selection with ray-step 1 and projection-step 10 0.000 seconds used for processing command sele <*> <#> STOP TERMINATION VARIANCE at 0.01 termination test vari epsilon = 0.010000 0.000 seconds used for processing command stop <*> <#> EXECUTE AVERAGE ART SMOOTH Example 2a underrelaxed art4 with variable tolerance <#> smoothing threshold is 0.2 Weights are 2.25 1.5 1.0 reconstruction is smoothed after iterations 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 last iteration threshold = 0.20000 w1 = 2.25000 w2 = 1.50000 w3 = 1.00000 <#> ART4 RELAXATION CONSTANT = 0.025 TOLERANCE VARIABLE decreasing art4 method relaxation parameter is 0.0250 tolerance is variable <#> CONSTRAINT BOUND CONRELAX CONSTANT = 0.5 NOMLIZE at each iteration constraint set to boundaries relax constraints with cr= 0.5000 picture is normalized to 0.9483 after each iteration 1890 rays are used for each iteration lower constraint = 0.0000 upper constraint = 2.5000 algorithm executed in iteration 1 0.001 seconds for the execution of the algorithm iteration 1 completed 0.001 seconds for this iteration algorithm executed in iteration 2 0.001 seconds for the execution of the algorithm iteration 2 completed 0.001 seconds for this iteration algorithm executed in iteration 3 0.002 seconds for the execution of the algorithm iteration 3 completed 0.002 seconds for this iteration algorithm executed in iteration 4 0.001 seconds for the execution of the algorithm iteration 4 completed 0.001 seconds for this iteration algorithm executed in iteration 5 0.003 seconds for the execution of the algorithm iteration 5 completed 0.003 seconds for this iteration algorithm executed in iteration 6 0.002 seconds for the execution of the algorithm iteration 6 completed 0.002 seconds for this iteration algorithm executed in iteration 7 0.001 seconds for the execution of the algorithm iteration 7 completed 0.001 seconds for this iteration algorithm executed in iteration 8 0.001 seconds for the execution of the algorithm iteration 8 completed 0.001 seconds for this iteration algorithm executed in iteration 9 0.001 seconds for the execution of the algorithm iteration 9 completed 0.001 seconds for this iteration algorithm executed in iteration 10 0.002 seconds for the execution of the algorithm iteration 10 completed 0.002 seconds for this iteration algorithm executed in iteration 11 0.001 seconds for the execution of the algorithm iteration 11 completed 0.001 seconds for this iteration algorithm executed in iteration 12 0.002 seconds for the execution of the algorithm iteration 12 completed 0.002 seconds for this iteration algorithm executed in iteration 13 0.001 seconds for the execution of the algorithm iteration 13 completed 0.001 seconds for this iteration algorithm executed in iteration 14 0.001 seconds for the execution of the algorithm iteration 14 completed 0.001 seconds for this iteration algorithm executed in iteration 15 0.001 seconds for the execution of the algorithm iteration 15 completed 0.001 seconds for this iteration algorithm executed in iteration 16 0.000 seconds for the execution of the algorithm iteration 16 completed 0.000 seconds for this iteration algorithm executed in iteration 17 0.002 seconds for the execution of the algorithm iteration 17 completed 0.002 seconds for this iteration algorithm executed in iteration 18 0.001 seconds for the execution of the algorithm iteration 18 completed 0.001 seconds for this iteration algorithm executed in iteration 19 0.001 seconds for the execution of the algorithm iteration 19 completed 0.001 seconds for this iteration algorithm executed in iteration 20 0.001 seconds for the execution of the algorithm iteration 20 completed 0.001 seconds for this iteration algorithm executed in iteration 21 0.002 seconds for the execution of the algorithm iterative process stops at iteration 21 the change in variance is less than 0.009852 of the variance reconstruction completed after iteration 21 0.002 seconds for this iteration 0.031 seconds for all iterations 0.031 seconds used for processing command exec <*> <#> STOP ITERATIONS = 30 30 iterations 0.000 seconds used for processing command stop <*> <#> EXECUTE ALP1 Example 2b underrelaxed art with extrapolation algorithm executed in iteration 1 0.001 seconds for the execution of the algorithm iteration 1 completed 0.001 seconds for this iteration algorithm executed in iteration 2 0.001 seconds for the execution of the algorithm iteration 2 completed 0.001 seconds for this iteration EXTRAPOLATION PARAMETER IN EXTRAP IS 5.062718796778e-01 algorithm executed in iteration 3 0.001 seconds for the execution of the algorithm iteration 3 completed 0.001 seconds for this iteration EXTRAPOLATION PARAMETER IN EXTRAP IS 5.062718796778e-01 algorithm executed in iteration 4 0.000 seconds for the execution of the algorithm iteration 4 completed 0.000 seconds for this iteration EXTRAPOLATION PARAMETER IN EXTRAP IS 6.560964044719e-01 algorithm executed in iteration 5 0.001 seconds for the execution of the algorithm iteration 5 completed 0.001 seconds for this iteration EXTRAPOLATION PARAMETER IN EXTRAP IS 6.560964044719e-01 algorithm executed in iteration 6 0.000 seconds for the execution of the algorithm iteration 6 completed 0.000 seconds for this iteration EXTRAPOLATION PARAMETER IN EXTRAP IS 7.480594958675e-01 algorithm executed in iteration 7 0.001 seconds for the execution of the algorithm iteration 7 completed 0.001 seconds for this iteration EXTRAPOLATION PARAMETER IN EXTRAP IS 7.480594958675e-01 algorithm executed in iteration 8 0.000 seconds for the execution of the algorithm iteration 8 completed 0.000 seconds for this iteration EXTRAPOLATION PARAMETER IN EXTRAP IS 7.951514183812e-01 algorithm executed in iteration 9 0.001 seconds for the execution of the algorithm iteration 9 completed 0.001 seconds for this iteration EXTRAPOLATION PARAMETER IN EXTRAP IS 7.951514183812e-01 algorithm executed in iteration 10 0.000 seconds for the execution of the algorithm iteration 10 completed 0.000 seconds for this iteration EXTRAPOLATION PARAMETER IN EXTRAP IS 8.250655739411e-01 algorithm executed in iteration 11 0.002 seconds for the execution of the algorithm iteration 11 completed 0.002 seconds for this iteration EXTRAPOLATION PARAMETER IN EXTRAP IS 8.250655739411e-01 algorithm executed in iteration 12 0.000 seconds for the execution of the algorithm iteration 12 completed 0.000 seconds for this iteration EXTRAPOLATION PARAMETER IN EXTRAP IS 8.467169748437e-01 algorithm executed in iteration 13 0.001 seconds for the execution of the algorithm iteration 13 completed 0.001 seconds for this iteration EXTRAPOLATION PARAMETER IN EXTRAP IS 8.467169748437e-01 algorithm executed in iteration 14 0.000 seconds for the execution of the algorithm iteration 14 completed 0.000 seconds for this iteration EXTRAPOLATION PARAMETER IN EXTRAP IS 8.628058734349e-01 algorithm executed in iteration 15 0.002 seconds for the execution of the algorithm iteration 15 completed 0.002 seconds for this iteration EXTRAPOLATION PARAMETER IN EXTRAP IS 8.628058734349e-01 algorithm executed in iteration 16 0.000 seconds for the execution of the algorithm iteration 16 completed 0.000 seconds for this iteration EXTRAPOLATION PARAMETER IN EXTRAP IS 8.749346145383e-01 algorithm executed in iteration 17 0.001 seconds for the execution of the algorithm iteration 17 completed 0.001 seconds for this iteration EXTRAPOLATION PARAMETER IN EXTRAP IS 8.749346145383e-01 algorithm executed in iteration 18 0.000 seconds for the execution of the algorithm iteration 18 completed 0.000 seconds for this iteration EXTRAPOLATION PARAMETER IN EXTRAP IS 8.843623393369e-01 algorithm executed in iteration 19 0.001 seconds for the execution of the algorithm iteration 19 completed 0.001 seconds for this iteration EXTRAPOLATION PARAMETER IN EXTRAP IS 8.843623393369e-01 algorithm executed in iteration 20 0.000 seconds for the execution of the algorithm iteration 20 completed 0.000 seconds for this iteration EXTRAPOLATION PARAMETER IN EXTRAP IS 8.920147280181e-01 algorithm executed in iteration 21 0.001 seconds for the execution of the algorithm iteration 21 completed 0.001 seconds for this iteration EXTRAPOLATION PARAMETER IN EXTRAP IS 8.920147280181e-01 algorithm executed in iteration 22 0.000 seconds for the execution of the algorithm iteration 22 completed 0.000 seconds for this iteration EXTRAPOLATION PARAMETER IN EXTRAP IS 8.985122210336e-01 algorithm executed in iteration 23 0.001 seconds for the execution of the algorithm iteration 23 completed 0.001 seconds for this iteration EXTRAPOLATION PARAMETER IN EXTRAP IS 8.985122210336e-01 algorithm executed in iteration 24 0.000 seconds for the execution of the algorithm iteration 24 completed 0.000 seconds for this iteration EXTRAPOLATION PARAMETER IN EXTRAP IS 9.042443038619e-01 algorithm executed in iteration 25 0.001 seconds for the execution of the algorithm iteration 25 completed 0.001 seconds for this iteration EXTRAPOLATION PARAMETER IN EXTRAP IS 9.042443038619e-01 algorithm executed in iteration 26 0.000 seconds for the execution of the algorithm iteration 26 completed 0.000 seconds for this iteration EXTRAPOLATION PARAMETER IN EXTRAP IS 9.094448804097e-01 algorithm executed in iteration 27 0.001 seconds for the execution of the algorithm iteration 27 completed 0.001 seconds for this iteration EXTRAPOLATION PARAMETER IN EXTRAP IS 9.094448804097e-01 algorithm executed in iteration 28 0.000 seconds for the execution of the algorithm iteration 28 completed 0.000 seconds for this iteration EXTRAPOLATION PARAMETER IN EXTRAP IS 9.142500229503e-01 algorithm executed in iteration 29 0.001 seconds for the execution of the algorithm iteration 29 completed 0.001 seconds for this iteration EXTRAPOLATION PARAMETER IN EXTRAP IS 9.142500229503e-01 algorithm executed in iteration 30 0.000 seconds for the execution of the algorithm reconstruction completed after iteration 30 0.000 seconds for this iteration 0.017 seconds for all iterations 0.018 seconds used for processing command exec <*> <*> Reset mode to unconstrained <#> MODE 0.000 seconds used for processing command mode <*> <#> EVALUATE RESOLUTION Example 2a underrelaxed art4 <#> WHOLEPIC Region cx cy u v ang t1 t2 wholepic -1e+20 1e+20 iterations 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 last iteration residual calculation using line integration 0.047 seconds used for processing command eval <*> <#> LINES COLUMNS 15 and 18 are to be printed iterations 10 20 30 last iteration lines display of reconstruction using ART iter 10 phantom name: Example 2, A simple test pattern and projection data projection data: Example 2, A simple test pattern and projection data execution name: Example 2a underrelaxed art4 with variable tolerance scaling factor = 1.00000 column = 15 18 row original estimate difference original estimate difference 0 0.00000 0.18210 0.18210 0.00000 0.00909 0.00909 1 0.00000 0.20450 0.20450 0.00000 0.24753 0.24753 2 0.00000 1.10539 1.10539 0.00000 0.52569 0.52569 3 2.50000 2.13754 -0.36246 2.50000 1.94767 -0.55233 4 2.50000 2.16437 -0.33563 2.50000 2.18812 -0.31188 5 2.50000 2.16394 -0.33606 2.50000 2.20758 -0.29242 6 2.50000 2.13317 -0.36683 2.50000 2.17610 -0.32390 7 1.25000 1.49667 0.24667 1.25000 1.68287 0.43287 8 1.25000 1.44121 0.19121 1.25000 1.41529 0.16529 9 1.25000 1.37200 0.12200 1.25000 1.38374 0.13374 10 1.25000 1.34328 0.09328 1.25000 1.34276 0.09276 11 1.25000 1.33396 0.08396 1.25000 1.32222 0.07222 12 1.35000 1.34036 -0.00964 1.25000 1.31411 0.06411 13 1.35000 1.34676 -0.00324 1.25000 1.31521 0.06521 14 1.35000 1.35058 0.00058 1.25000 1.31568 0.06568 15 1.35000 1.35053 0.00053 1.35000 1.31710 -0.03290 16 1.35000 1.35058 0.00058 1.25000 1.31607 0.06607 17 1.35000 1.34676 -0.00324 1.25000 1.31587 0.06587 18 1.35000 1.34035 -0.00965 1.25000 1.31417 0.06417 19 1.25000 1.33395 0.08395 1.25000 1.32131 0.07131 20 1.25000 1.34326 0.09326 1.25000 1.34021 0.09021 21 1.25000 1.37196 0.12196 1.25000 1.38008 0.13008 22 1.25000 1.44114 0.19114 1.25000 1.41090 0.16090 23 1.25000 1.49658 0.24658 1.25000 1.68040 0.43040 24 2.50000 2.13333 -0.36667 2.50000 2.17405 -0.32595 25 2.50000 2.16414 -0.33586 2.50000 2.20535 -0.29465 26 2.50000 2.16455 -0.33545 2.50000 2.18655 -0.31345 27 2.50000 2.13774 -0.36226 2.50000 1.94456 -0.55544 28 0.00000 1.10534 1.10534 0.00000 0.44393 0.44393 29 0.00000 0.20439 0.20439 0.00000 0.29112 0.29112 30 0.00000 0.18198 0.18198 0.00000 0.01326 0.01326 lines display of reconstruction using ART iter 20 phantom name: Example 2, A simple test pattern and projection data projection data: Example 2, A simple test pattern and projection data execution name: Example 2a underrelaxed art4 with variable tolerance scaling factor = 1.00000 column = 15 18 row original estimate difference original estimate difference 0 0.00000 0.00164 0.00164 0.00000 0.02169 0.02169 1 0.00000 0.00132 0.00132 0.00000 0.04480 0.04480 2 0.00000 0.94940 0.94940 0.00000 0.13226 0.13226 3 2.50000 2.43151 -0.06849 2.50000 2.25420 -0.24580 4 2.50000 2.43097 -0.06903 2.50000 2.44794 -0.05206 5 2.50000 2.37936 -0.12064 2.50000 2.44383 -0.05617 6 2.50000 2.31917 -0.18083 2.50000 2.39051 -0.10949 7 1.25000 1.46751 0.21751 1.25000 1.65881 0.40881 8 1.25000 1.31935 0.06935 1.25000 1.33707 0.08707 9 1.25000 1.29770 0.04770 1.25000 1.31051 0.06051 10 1.25000 1.27744 0.02744 1.25000 1.27762 0.02762 11 1.25000 1.28043 0.03043 1.25000 1.26720 0.01720 12 1.35000 1.30399 -0.04601 1.25000 1.26843 0.01843 13 1.35000 1.32095 -0.02905 1.25000 1.27985 0.02985 14 1.35000 1.32945 -0.02055 1.25000 1.28645 0.03645 15 1.35000 1.32847 -0.02153 1.35000 1.29098 -0.05902 16 1.35000 1.32944 -0.02056 1.25000 1.28714 0.03714 17 1.35000 1.32093 -0.02907 1.25000 1.28136 0.03136 18 1.35000 1.30399 -0.04601 1.25000 1.27000 0.02000 19 1.25000 1.28044 0.03044 1.25000 1.26871 0.01871 20 1.25000 1.27744 0.02744 1.25000 1.27815 0.02815 21 1.25000 1.29767 0.04767 1.25000 1.30973 0.05973 22 1.25000 1.31930 0.06930 1.25000 1.33504 0.08504 23 1.25000 1.46737 0.21737 1.25000 1.71927 0.46927 24 2.50000 2.31941 -0.18059 2.50000 2.39181 -0.10819 25 2.50000 2.37969 -0.12031 2.50000 2.44490 -0.05510 26 2.50000 2.43125 -0.06875 2.50000 2.44917 -0.05083 27 2.50000 2.43179 -0.06821 2.50000 2.25444 -0.24556 28 0.00000 0.94928 0.94928 0.00000 0.12809 0.12809 29 0.00000 0.00158 0.00158 0.00000 0.04576 0.04576 30 0.00000 0.00199 0.00199 0.00000 0.02628 0.02628 lines display of reconstruction using ART iter 21 phantom name: Example 2, A simple test pattern and projection data projection data: Example 2, A simple test pattern and projection data execution name: Example 2a underrelaxed art4 with variable tolerance scaling factor = 1.00000 column = 15 18 row original estimate difference original estimate difference 0 0.00000 0.00100 0.00100 0.00000 0.01742 0.01742 1 0.00000 0.00070 0.00070 0.00000 0.03672 0.03672 2 0.00000 0.93506 0.93506 0.00000 0.11759 0.11759 3 2.50000 2.44949 -0.05051 2.50000 2.27648 -0.22352 4 2.50000 2.44615 -0.05385 2.50000 2.45838 -0.04162 5 2.50000 2.39029 -0.10971 2.50000 2.45205 -0.04795 6 2.50000 2.32845 -0.17155 2.50000 2.39929 -0.10071 7 1.25000 1.46233 0.21233 1.25000 1.65516 0.40516 8 1.25000 1.31367 0.06367 1.25000 1.33085 0.08085 9 1.25000 1.29279 0.04279 1.25000 1.30546 0.05546 10 1.25000 1.27400 0.02400 1.25000 1.27412 0.02412 11 1.25000 1.27844 0.02844 1.25000 1.26496 0.01496 12 1.35000 1.30374 -0.04626 1.25000 1.26714 0.01714 13 1.35000 1.32159 -0.02841 1.25000 1.27948 0.02948 14 1.35000 1.33036 -0.01964 1.25000 1.28660 0.03660 15 1.35000 1.32919 -0.02081 1.35000 1.29140 -0.05860 16 1.35000 1.33034 -0.01966 1.25000 1.28726 0.03726 17 1.35000 1.32157 -0.02843 1.25000 1.28096 0.03096 18 1.35000 1.30372 -0.04628 1.25000 1.26874 0.01874 19 1.25000 1.27844 0.02844 1.25000 1.26657 0.01657 20 1.25000 1.27400 0.02400 1.25000 1.27479 0.02479 21 1.25000 1.29277 0.04277 1.25000 1.30483 0.05483 22 1.25000 1.31362 0.06362 1.25000 1.32898 0.07898 23 1.25000 1.46218 0.21218 1.25000 1.65781 0.40781 24 2.50000 2.32869 -0.17131 2.50000 2.40065 -0.09935 25 2.50000 2.39063 -0.10937 2.50000 2.45318 -0.04682 26 2.50000 2.44644 -0.05356 2.50000 2.45967 -0.04033 27 2.50000 2.44977 -0.05023 2.50000 2.27692 -0.22308 28 0.00000 0.93489 0.93489 0.00000 0.11431 0.11431 29 0.00000 0.00120 0.00120 0.00000 0.03840 0.03840 30 0.00000 0.00161 0.00161 0.00000 0.02249 0.02249 lines display of reconstruction using ALP1 iter 10 phantom name: Example 2, A simple test pattern and projection data projection data: Example 2, A simple test pattern and projection data execution name: Example 2b underrelaxed art with extrapolation scaling factor = 1.00000 column = 15 18 row original estimate difference original estimate difference 0 0.00000 0.12670 0.12670 0.00000 -0.02567 -0.02567 1 0.00000 0.17187 0.17187 0.00000 0.16958 0.16958 2 0.00000 1.20151 1.20151 0.00000 0.48350 0.48350 3 2.50000 2.23551 -0.26449 2.50000 2.01480 -0.48520 4 2.50000 2.39456 -0.10544 2.50000 2.31508 -0.18492 5 2.50000 2.34966 -0.15034 2.50000 2.46624 -0.03376 6 2.50000 2.31124 -0.18876 2.50000 2.29509 -0.20491 7 1.25000 1.63106 0.38106 1.25000 1.81482 0.56482 8 1.25000 1.35183 0.10183 1.25000 1.44203 0.19203 9 1.25000 1.42885 0.17885 1.25000 1.45015 0.20015 10 1.25000 1.31859 0.06859 1.25000 1.36973 0.11973 11 1.25000 1.42731 0.17731 1.25000 1.34954 0.09954 12 1.35000 1.31939 -0.03061 1.25000 1.31390 0.06390 13 1.35000 1.37323 0.02323 1.25000 1.32415 0.07415 14 1.35000 1.46590 0.11590 1.25000 1.37584 0.12584 15 1.35000 1.13749 -0.21251 1.35000 1.25214 -0.09786 16 1.35000 1.46594 0.11594 1.25000 1.36902 0.11902 17 1.35000 1.37317 0.02317 1.25000 1.30699 0.05699 18 1.35000 1.31944 -0.03056 1.25000 1.34043 0.09043 19 1.25000 1.42727 0.17727 1.25000 1.33982 0.08982 20 1.25000 1.31861 0.06861 1.25000 1.35654 0.10654 21 1.25000 1.42880 0.17880 1.25000 1.44218 0.19218 22 1.25000 1.35172 0.10172 1.25000 1.45375 0.20375 23 1.25000 1.63100 0.38100 1.25000 1.80751 0.55751 24 2.50000 2.31161 -0.18839 2.50000 2.27967 -0.22033 25 2.50000 2.34937 -0.15063 2.50000 2.50118 0.00118 26 2.50000 2.39545 -0.10455 2.50000 2.30942 -0.19058 27 2.50000 2.23526 -0.26474 2.50000 1.99542 -0.50458 28 0.00000 1.20133 1.20133 0.00000 0.49478 0.49478 29 0.00000 0.17173 0.17173 0.00000 0.15734 0.15734 30 0.00000 0.12647 0.12647 0.00000 -0.02832 -0.02832 lines display of reconstruction using ALP1 iter 20 phantom name: Example 2, A simple test pattern and projection data projection data: Example 2, A simple test pattern and projection data execution name: Example 2b underrelaxed art with extrapolation scaling factor = 1.00000 column = 15 18 row original estimate difference original estimate difference 0 0.00000 -0.09663 -0.09663 0.00000 -0.18986 -0.18986 1 0.00000 -0.12331 -0.12331 0.00000 -0.05402 -0.05402 2 0.00000 1.04579 1.04579 0.00000 0.12373 0.12373 3 2.50000 2.51540 0.01540 2.50000 2.25540 -0.24460 4 2.50000 2.59753 0.09753 2.50000 2.56869 0.06869 5 2.50000 2.55529 0.05529 2.50000 2.67846 0.17846 6 2.50000 2.40841 -0.09159 2.50000 2.45046 -0.04954 7 1.25000 1.52857 0.27857 1.25000 1.73851 0.48851 8 1.25000 1.19882 -0.05118 1.25000 1.29003 0.04003 9 1.25000 1.28220 0.03220 1.25000 1.30029 0.05029 10 1.25000 1.18525 -0.06475 1.25000 1.24281 -0.00719 11 1.25000 1.27819 0.02819 1.25000 1.24380 -0.00620 12 1.35000 1.25011 -0.09989 1.25000 1.22237 -0.02763 13 1.35000 1.30966 -0.04034 1.25000 1.24476 -0.00524 14 1.35000 1.38630 0.03630 1.25000 1.28583 0.03583 15 1.35000 1.11801 -0.23199 1.35000 1.20027 -0.14973 16 1.35000 1.38636 0.03636 1.25000 1.28050 0.03050 17 1.35000 1.30959 -0.04041 1.25000 1.23298 -0.01702 18 1.35000 1.25014 -0.09986 1.25000 1.24848 -0.00152 19 1.25000 1.27817 0.02817 1.25000 1.23597 -0.01403 20 1.25000 1.18526 -0.06474 1.25000 1.23189 -0.01811 21 1.25000 1.28211 0.03211 1.25000 1.29998 0.04998 22 1.25000 1.19865 -0.05135 1.25000 1.30195 0.05195 23 1.25000 1.52839 0.27839 1.25000 1.73195 0.48195 24 2.50000 2.40899 -0.09101 2.50000 2.44156 -0.05844 25 2.50000 2.55477 0.05477 2.50000 2.70939 0.20939 26 2.50000 2.59890 0.09890 2.50000 2.56774 0.06774 27 2.50000 2.51498 0.01498 2.50000 2.24565 -0.25435 28 0.00000 1.04557 1.04557 0.00000 0.14374 0.14374 29 0.00000 -0.12338 -0.12338 0.00000 -0.05331 -0.05331 30 0.00000 -0.09670 -0.09670 0.00000 -0.17428 -0.17428 lines display of reconstruction using ALP1 iter 30 phantom name: Example 2, A simple test pattern and projection data projection data: Example 2, A simple test pattern and projection data execution name: Example 2b underrelaxed art with extrapolation scaling factor = 1.00000 column = 15 18 row original estimate difference original estimate difference 0 0.00000 -0.10295 -0.10295 0.00000 -0.15511 -0.15511 1 0.00000 -0.21068 -0.21068 0.00000 -0.09236 -0.09236 2 0.00000 0.92855 0.92855 0.00000 -0.06300 -0.06300 3 2.50000 2.64511 0.14511 2.50000 2.37419 -0.12581 4 2.50000 2.61843 0.11843 2.50000 2.63291 0.13291 5 2.50000 2.61277 0.11277 2.50000 2.68002 0.18002 6 2.50000 2.39593 -0.10407 2.50000 2.48350 -0.01650 7 1.25000 1.47574 0.22574 1.25000 1.66721 0.41721 8 1.25000 1.11475 -0.13525 1.25000 1.19787 -0.05213 9 1.25000 1.24880 -0.00120 1.25000 1.25371 0.00371 10 1.25000 1.15550 -0.09450 1.25000 1.22544 -0.02456 11 1.25000 1.25545 0.00545 1.25000 1.24429 -0.00571 12 1.35000 1.26430 -0.08570 1.25000 1.22826 -0.02174 13 1.35000 1.32865 -0.02135 1.25000 1.25347 0.00347 14 1.35000 1.40845 0.05845 1.25000 1.28571 0.03571 15 1.35000 1.14519 -0.20481 1.35000 1.22750 -0.12250 16 1.35000 1.40851 0.05851 1.25000 1.28178 0.03178 17 1.35000 1.32857 -0.02143 1.25000 1.24374 -0.00626 18 1.35000 1.26433 -0.08567 1.25000 1.25314 0.00314 19 1.25000 1.25545 0.00545 1.25000 1.23609 -0.01391 20 1.25000 1.15550 -0.09450 1.25000 1.21469 -0.03531 21 1.25000 1.24869 -0.00131 1.25000 1.25626 0.00626 22 1.25000 1.11456 -0.13544 1.25000 1.20898 -0.04102 23 1.25000 1.47547 0.22547 1.25000 1.65907 0.40907 24 2.50000 2.39664 -0.10336 2.50000 2.47619 -0.02381 25 2.50000 2.61204 0.11204 2.50000 2.70609 0.20609 26 2.50000 2.62009 0.12009 2.50000 2.63048 0.13048 27 2.50000 2.64455 0.14455 2.50000 2.36445 -0.13555 28 0.00000 0.92830 0.92830 0.00000 -0.04433 -0.04433 29 0.00000 -0.21071 -0.21071 0.00000 -0.09436 -0.09436 30 0.00000 -0.10286 -0.10286 0.00000 -0.14676 -0.14676 0.002 seconds used for processing command line <*> <#> END