--- /dev/null
+ snark14.s170426 - A PICTURE RECONSTRUCTION PROGRAM
+
+
+ <*> *******EXAMPLE 10
+
+ <*> ROTATED RECTANGLE WITH PSEUDO PROJECTION DATA
+
+ <*> RECONSTRUCTION BY ART WITH BLOB BASIS
+
+ <*>
+
+ <*> TRACE 1
+
+ <#> CREATE
+
+ EXAMPLE 10 ROTATED RECTANGLE
+
+
+ <#> SPECTRUM MONOCHROMATIC 10
+ energy spectrum is monochromatic at energy level 10
+
+
+ <#> OBJECTS
+ description of objects
+ density at levels
+ numb type x-coord y-coord x-length y-length angle av dens 10
+
+ 1 rect 0.0000 0.0000 2.0000 10.0000 -30.0000 1.0000 1.0000
+
+ scale factor multiplying object densities 1.0000
+
+ seed set to 0
+ inhomogeneity set to 0.0000
+
+ <#> PHANTOM AVERAGE 1
+
+ this run will generate a phantom
+ density in each pixel is obtained as the average of 1 x 1 points
+
+
+ <#> 25 PIXELS OF SIZE 1.0
+ picture size 25 x 25, pixel size 1.0000
+
+
+ <#> RAYSUM
+
+ 0.000 seconds phantom creation
+ 0.000 seconds used for processing command crea
+
+
+ <*> RUN
+
+ <*>
+
+ <#> PICTURE TEST
+
+ EXAMPLE 10 ROTATED RECTANGLE
+
+
+ <#> spec mono 10
+ energy spectrum is monochromatic at energy level 10
+
+
+ <#> obje
+ description of objects
+ density at levels
+ numb type x-coord y-coord x-length y-length angle av dens 10
+
+ 1 rect 0.0000 0.0000 2.0000 10.0000 -30.0000 1.0000 1.0000
+
+ 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 25 size 1.0000
+ picture size 25 x 25, pixel size 1.0000
+
+ test picture read
+ EXAMPLE 10 ROTATED RECTANGLE
+ 0.000 seconds used for processing command pict
+
+
+ <*>
+
+ <#> PROJECTION PSEUDO
+
+ EXAMPLE 10 PSEUDO PROJECTION DATA FOR TWELVE VIEWS
+
+
+
+ <#> GEOMETRY
+
+
+ <#> PARALLEL UNIFORM LINE
+ rays are parallel with uniform spacing between rays
+ data collected along lines
+
+
+
+ <#> RAYS PROGRAM 25 2.0 DETECTOR SPACING AT 2.0
+ number of rays per projection 37
+ snark computed number of rays 19
+ at detector spacing 2.0000
+
+
+ <#> ANGLES 12 EQUALLY SPACED
+ total number of projections 12
+
+ projection angles 0.0000 15.0000 30.0000 45.0000 60.0000 75.0000 90.0000 105.0000 120.0000 135.0000
+ 150.0000 165.0000
+
+
+ <#> MEASUREMENT PERFECT
+ projection data are noiseless
+
+
+ <#> BACKGROUND 50.0
+ at levels
+ 10
+ background absorption 50.0000
+
+ estimate of totlen = 3776.330760
+ estimate of totden = 498.725717
+ estimate of average density = 0.1321
+ projection data read
+ EXAMPLE 10 PSEUDO PROJECTION DATA FOR TWELVE VIEWS
+ 0.000 seconds used for processing command proj
+
+
+ <*>
+
+ <#> BASIS BLOB
+ hexagonal grid: there are 25 rows, the middle row has 21 grid points and the row next to the middle one has 22 grid points
+ blob basis with parameters: support = 2.0629 shape = 11.2829 delta = 1.1547
+ 0.083 seconds used for processing command basi
+
+
+ <*>
+
+ <#> SELECT SNARK RAYSEQ
+
+ <#> STEP 3 1
+ sequential ray selection with ray-step 1 and projection-step 3
+ 0.000 seconds used for processing command sele
+
+
+ <*>
+
+ <#> STOP ITERATION 10
+ 10 iterations
+ 0.000 seconds used for processing command stop
+
+
+ <*>
+
+ <#> EXECUTE AVERAGE ART
+
+ ART WITH BLOB BASIS
+
+ <#> ART3
+ art3 method
+ relaxation parameter is 1.0
+ norm is 2
+ tolerance is 0.0
+
+ <#> CONSTRAINT ART2
+ constraint art2 way
+ relaxation constraints with cr = 1.0
+ picture is not normalized
+ 228 rays are used for each iteration
+ 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.000 seconds for the execution of the algorithm
+ iteration 2 completed
+ 0.001 seconds for this iteration
+ algorithm executed in iteration 3
+ 0.000 seconds for the execution of the algorithm
+ iteration 3 completed
+ 0.001 seconds for this iteration
+ algorithm executed in iteration 4
+ 0.000 seconds for the execution of the algorithm
+ iteration 4 completed
+ 0.001 seconds for this iteration
+ algorithm executed in iteration 5
+ 0.000 seconds for the execution of the algorithm
+ iteration 5 completed
+ 0.001 seconds for this iteration
+ algorithm executed in iteration 6
+ 0.000 seconds for the execution of the algorithm
+ iteration 6 completed
+ 0.001 seconds for this iteration
+ algorithm executed in iteration 7
+ 0.000 seconds for the execution of the algorithm
+ iteration 7 completed
+ 0.001 seconds for this iteration
+ algorithm executed in iteration 8
+ 0.000 seconds for the execution of the algorithm
+ iteration 8 completed
+ 0.001 seconds for this iteration
+ algorithm executed in iteration 9
+ 0.000 seconds for the execution of the algorithm
+ iteration 9 completed
+ 0.001 seconds for this iteration
+ algorithm executed in iteration 10
+ 0.000 seconds for the execution of the algorithm
+ reconstruction completed after iteration 10
+ 0.001 seconds for this iteration
+ 0.007 seconds for all iterations
+ 0.008 seconds used for processing command exec
+
+
+ <*>
+
+ <#> BASIS PIXEL
+ pixel basis
+ 0.000 seconds used for processing command basi
+
+
+ <#> EXECUTE AVERAGE ART
+
+ ART WITH PIXEL BASIS
+
+ <#> ART3
+ art3 method
+ relaxation parameter is 1.0
+ norm is 2
+ tolerance is 0.0
+
+ <#> CONSTRAINT ART2
+ constraint art2 way
+ relaxation constraints with cr = 1.0
+ picture is not normalized
+ 228 rays are used for each iteration
+ algorithm executed in iteration 1
+ 0.000 seconds for the execution of the algorithm
+ iteration 1 completed
+ 0.000 seconds for this iteration
+ algorithm executed in iteration 2
+ 0.000 seconds for the execution of the algorithm
+ iteration 2 completed
+ 0.000 seconds for this iteration
+ algorithm executed in iteration 3
+ 0.000 seconds for the execution of the algorithm
+ iteration 3 completed
+ 0.000 seconds for this iteration
+ algorithm executed in iteration 4
+ 0.000 seconds for the execution of the algorithm
+ iteration 4 completed
+ 0.000 seconds for this iteration
+ algorithm executed in iteration 5
+ 0.000 seconds for the execution of the algorithm
+ iteration 5 completed
+ 0.000 seconds for this iteration
+ algorithm executed in iteration 6
+ 0.000 seconds for the execution of the algorithm
+ iteration 6 completed
+ 0.000 seconds for this iteration
+ algorithm executed in iteration 7
+ 0.000 seconds for the execution of the algorithm
+ iteration 7 completed
+ 0.000 seconds for this iteration
+ algorithm executed in iteration 8
+ 0.000 seconds for the execution of the algorithm
+ iteration 8 completed
+ 0.000 seconds for this iteration
+ algorithm executed in iteration 9
+ 0.000 seconds for the execution of the algorithm
+ iteration 9 completed
+ 0.000 seconds for this iteration
+ algorithm executed in iteration 10
+ 0.000 seconds for the execution of the algorithm
+ reconstruction completed after iteration 10
+ 0.000 seconds for this iteration
+ 0.001 seconds for all iterations
+ 0.002 seconds used for processing command exec
+
+
+ <*>
+
+ <#> EVALUATE BOTH
+
+ COMPARING BLOB AND PIXEL BASIS
+
+ <#> WHOLEPIC
+ Region cx cy u v ang t1 t2
+
+ wholepic -1e+20 1e+20
+
+ iterations 1 2 3 4 5 6 7 8 9
+ last iteration
+
+ residual calculation using line integration
+ 0.004 seconds used for processing command eval
+
+
+ <#> END