Exclude examples with random number generators from regression. Add new examples...

diff --git a/examples/b12/README b/examples/b12/README
index 9127728d06eb1c592318352fe2e704e6285e243a..156db76dadb0c74e1612270c7fc717b6ab24736b 100644 (file)
--- a/examples/b12/README
+++ b/examples/b12/README
@@ -5,3 +5,9 @@ libraries.
 This version uses PERFECT data so it can be used to test regression
 across different platforms.
 
 This version uses PERFECT data so it can be used to test regression
 across different platforms.
 

+WARNING
+
+This version does not pass regression testing. There is a very small difference in value of phi between platforms:
+
+gcc 5.4 on Windows 10 WSL and Centos original version match
+gcc 7.2 on Ubuntu 2017.10 does not match

diff --git a/examples/b12/b12r-centos-orig.out b/examples/b12/b12r-centos-orig.out
new file mode 100644 (file)
index 0000000..ce900f5
--- /dev/null
+++ b/examples/b12/b12r-centos-orig.out
@@ -0,0 +1,505 @@
+   snark14.s170710 - A PICTURE RECONSTRUCTION PROGRAM
+
+
+     <*> *******EXAMPLE 11
+
+     <*>  SUPERIORIZED MAP EM ALGORITHM FOR EMISSION TOMOGRAPHY. RECONSTRUCTION
+
+     <*>  OF BRAIN PHANTOM. SIMULATING PET GEOMETRY WITH A RING OF 300 DETECTORS
+
+     <*>  WITH EACH DETECTOR IN COINCIDENCE WITH 101 DETECTORS OPPOSITE IT.
+
+     <*>  COEFFICIENT OF PENALTLY TERM IS SET TO 0
+
+     <*> 
+
+     <#> CREATE                                                                                                                                                                                                                                                           
+
+         EXAMPLE 11 Brain Phantom 
+
+
+     <#> SPECTRUM MONOCHROMATIC 511
+         energy spectrum is monochromatic at energy level   511
+
+
+     <#> OBJECTS
+         description of objects
+                                                                density at levels
+ numb type  x-coord  y-coord x-length y-length    angle  av dens       511
+
+    1 elip  -7.0000  46.0000   3.0000   6.0000  17.0000   0.9500    0.9500
+
+    2 elip   7.0000  46.0000   3.0000   6.0000 -17.0000   1.0000    1.0000
+
+    3 rect -12.0000  64.0000   7.5000   4.5000   5.0000   1.0000    1.0000
+
+    4 rect  12.0000  64.0000   7.5000   4.5000  -5.0000   0.9500    0.9500
+
+    5 rect -38.0000  51.0000   3.5000  13.0000 -39.0000   0.9500    0.9500
+
+    6 rect  38.0000  51.0000   3.5000  13.0000  39.0000   1.0000    1.0000
+
+    7 rect -46.0000  24.0000   6.5000   6.0000 -18.0000   0.9500    0.9500
+
+    8 rect  46.0000  24.0000   6.5000   6.0000  18.0000   1.0000    1.0000
+
+    9 rect -49.0000   6.0000   2.5000  10.0000  63.0000   1.0000    1.0000
+
+   10 rect  49.0000   6.0000   2.5000  10.0000 -63.0000   0.9500    0.9500
+
+   11 rect -52.0000 -14.0000   9.0000   7.0000 -14.0000   0.9500    0.9500
+
+   12 rect  52.0000 -14.0000   9.0000   7.0000  14.0000   1.0000    1.0000
+
+   13 rect -10.0000 -56.0000   5.5000  10.0000  -1.0000   0.9500    0.9500
+
+   14 rect  10.0000 -56.0000   5.5000  10.0000   1.0000   1.0000    1.0000
+
+   15 elip -40.0000 -47.0000   9.0000  22.5000  48.0000   1.0000    1.0000
+
+   16 elip  40.0000 -47.0000   9.0000  22.5000 -48.0000   0.9500    0.9500
+
+   17 elip  -8.0000 -22.0000   3.5000  15.5000  -9.0000   1.0000    1.0000
+
+   18 elip   8.0000 -22.0000   3.5000  15.5000   9.0000   0.9500    0.9500
+
+   19 elip -27.0000  -6.0000   5.5000  23.5000  -5.0000   0.9500    0.9500
+
+   20 elip  27.0000  -6.0000   5.5000  23.5000   5.0000   1.0000    1.0000
+
+   21 elip -25.0000  38.0000   6.5000  10.5000 -14.0000   1.0000    1.0000
+
+   22 elip  25.0000  38.0000   6.5000  10.5000  14.0000   0.9500    0.9500
+
+   23 rect  -8.0000  32.0000   1.5000   6.5000  38.0000   1.0000    1.0000
+
+   24 rect   8.0000  32.0000   1.5000   6.5000 -38.0000   0.9500    0.9500
+
+   25 rect  -8.0000   3.0000   1.0000   9.0000 -33.0000   0.9500    0.9500
+
+   26 rect   8.0000   3.0000   1.0000   9.0000  33.0000   1.0000    1.0000
+
+   27 elip   0.0000   0.0000  66.5000  74.0000   0.0000   1.0000    1.0000
+
+         scale factor multiplying object densities     0.5100
+
+         seed set to 1
+         inhomogeneity set to     0.0500
+
+     <#> PHANTOM AVERAGE 7
+
+         this run will generate a phantom
+         density in each pixel is obtained as the average of 7 x 7 points
+
+
+     <#> 381 PIXELS OF SIZE 0.4
+         picture size 381 x 381,  pixel size     0.4000
+
+
+     <#> RAYSUM AVERAGE 1
+
+         this run will generate projection data
+         projection data are calculated by dividing each ray interval into 1 substrips
+
+         with aperture (substrip) weights     1
+
+
+     <#> GEOMETRY
+
+
+     <#> divergent arc 153 306
+         rays are divergent from point sources
+         source to origin distance    153.0000
+         the detectors lie on an arc with source to detector distance =   306.0000
+
+
+     <#> RAYS USER 301 DETECTOR SPACING 1.1
+         number of rays per projection   301
+         at detector spacing     1.1000
+
+
+     <#> ANGLES 500 EQUAL SPACING
+         total number of projections   500
+
+         projection angles    0.0000    0.7190    1.4381    2.1571    2.8762    3.5952    4.3142    5.0333    5.7523    6.4713
+                              7.1904    7.9094    8.6285    9.3475   10.0665   10.7856   11.5046   12.2236   12.9427   13.6617
+                             14.3808   15.0998   15.8188   16.5379   17.2569   17.9760   18.6950   19.4140   20.1331   20.8521
+                             21.5711   22.2902   23.0092   23.7283   24.4473   25.1663   25.8854   26.6044   27.3234   28.0425
+                             28.7615   29.4806   30.1996   30.9186   31.6377   32.3567   33.0758   33.7948   34.5138   35.2329
+                             35.9519   36.6709   37.3900   38.1090   38.8281   39.5471   40.2661   40.9852   41.7042   42.4232
+                             43.1423   43.8613   44.5804   45.2994   46.0184   46.7375   47.4565   48.1756   48.8946   49.6136
+                             50.3327   51.0517   51.7707   52.4898   53.2088   53.9279   54.6469   55.3659   56.0850   56.8040
+                             57.5230   58.2421   58.9611   59.6802   60.3992   61.1182   61.8373   62.5563   63.2754   63.9944
+                             64.7134   65.4325   66.1515   66.8705   67.5896   68.3086   69.0277   69.7467   70.4657   71.1848
+                             71.9038   72.6228   73.3419   74.0609   74.7800   75.4990   76.2180   76.9371   77.6561   78.3752
+                             79.0942   79.8132   80.5323   81.2513   81.9703   82.6894   83.4084   84.1275   84.8465   85.5655
+                             86.2846   87.0036   87.7226   88.4417   89.1607   89.8798   90.5988   91.3178   92.0369   92.7559
+                             93.4749   94.1940   94.9130   95.6321   96.3511   97.0701   97.7892   98.5082   99.2273   99.9463
+                            100.6653  101.3844  102.1034  102.8224  103.5415  104.2605  104.9796  105.6986  106.4176  107.1367
+                            107.8557  108.5747  109.2938  110.0128  110.7319  111.4509  112.1699  112.8890  113.6080  114.3271
+                            115.0461  115.7651  116.4842  117.2032  117.9222  118.6413  119.3603  120.0794  120.7984  121.5174
+                            122.2365  122.9555  123.6745  124.3936  125.1126  125.8317  126.5507  127.2697  127.9888  128.7078
+                            129.4269  130.1459  130.8649  131.5840  132.3030  133.0220  133.7411  134.4601  135.1792  135.8982
+                            136.6172  137.3363  138.0553  138.7743  139.4934  140.2124  140.9315  141.6505  142.3695  143.0886
+                            143.8076  144.5267  145.2457  145.9647  146.6838  147.4028  148.1218  148.8409  149.5599  150.2790
+                            150.9980  151.7170  152.4361  153.1551  153.8741  154.5932  155.3122  156.0313  156.7503  157.4693
+                            158.1884  158.9074  159.6265  160.3455  161.0645  161.7836  162.5026  163.2216  163.9407  164.6597
+                            165.3788  166.0978  166.8168  167.5359  168.2549  168.9739  169.6930  170.4120  171.1311  171.8501
+                            172.5691  173.2882  174.0072  174.7263  175.4453  176.1643  176.8834  177.6024  178.3214  179.0405
+                            179.7595  180.4786  181.1976  181.9166  182.6357  183.3547  184.0737  184.7928  185.5118  186.2309
+                            186.9499  187.6689  188.3880  189.1070  189.8261  190.5451  191.2641  191.9832  192.7022  193.4212
+                            194.1403  194.8593  195.5784  196.2974  197.0164  197.7355  198.4545  199.1735  199.8926  200.6116
+                            201.3307  202.0497  202.7687  203.4878  204.2068  204.9259  205.6449  206.3639  207.0830  207.8020
+                            208.5210  209.2401  209.9591  210.6782  211.3972  212.1162  212.8353  213.5543  214.2733  214.9924
+                            215.7114  216.4305  217.1495  217.8685  218.5876  219.3066  220.0257  220.7447  221.4637  222.1828
+                            222.9018  223.6208  224.3399  225.0589  225.7780  226.4970  227.2160  227.9351  228.6541  229.3731
+                            230.0922  230.8112  231.5303  232.2493  232.9683  233.6874  234.4064  235.1255  235.8445  236.5635
+                            237.2826  238.0016  238.7206  239.4397  240.1587  240.8778  241.5968  242.3158  243.0349  243.7539
+                            244.4729  245.1920  245.9110  246.6301  247.3491  248.0681  248.7872  249.5062  250.2253  250.9443
+                            251.6633  252.3824  253.1014  253.8204  254.5395  255.2585  255.9776  256.6966  257.4156  258.1347
+                            258.8537  259.5727  260.2918  261.0108  261.7299  262.4489  263.1679  263.8870  264.6060  265.3251
+                            266.0441  266.7631  267.4822  268.2012  268.9202  269.6393  270.3583  271.0774  271.7964  272.5154
+                            273.2345  273.9535  274.6725  275.3916  276.1106  276.8297  277.5487  278.2677  278.9868  279.7058
+                            280.4248  281.1439  281.8629  282.5820  283.3010  284.0200  284.7391  285.4581  286.1772  286.8962
+                            287.6152  288.3343  289.0533  289.7723  290.4914  291.2104  291.9295  292.6485  293.3675  294.0866
+                            294.8056  295.5246  296.2437  296.9627  297.6818  298.4008  299.1198  299.8389  300.5579  301.2770
+                            301.9960  302.7150  303.4341  304.1531  304.8721  305.5912  306.3102  307.0293  307.7483  308.4673
+                            309.1864  309.9054  310.6244  311.3435  312.0625  312.7816  313.5006  314.2196  314.9387  315.6577
+                            316.3768  317.0958  317.8148  318.5339  319.2529  319.9719  320.6910  321.4100  322.1291  322.8481
+                            323.5671  324.2862  325.0052  325.7242  326.4433  327.1623  327.8814  328.6004  329.3194  330.0385
+                            330.7575  331.4766  332.1956  332.9146  333.6337  334.3527  335.0717  335.7908  336.5098  337.2289
+                            337.9479  338.6669  339.3860  340.1050  340.8240  341.5431  342.2621  342.9812  343.7002  344.4192
+                            345.1383  345.8573  346.5764  347.2954  348.0144  348.7335  349.4525  350.1715  350.8906  351.6096
+                            352.3287  353.0477  353.7667  354.4858  355.2048  355.9238  356.6429  357.3619  358.0810  358.8000
+
+
+     <#> MEASUREMENT PERFECT
+         projection data are noiseless
+
+
+     <#> BACKGROUND 0.0
+                               at levels
+                                     511 
+          background absorption   0.0000 
+ 
+
+     <#> RUN
+              0.438 seconds phantom creation
+              1.219 seconds projection data creation
+              1.656 seconds used for processing command crea
+
+
+     <*> 
+
+     <#> PICTURE TEST                                                                                                                                                                                                                                                     
+
+         EXAMPLE 11 Brain Phantom 
+
+
+     <#> spec    mono  511
+         energy spectrum is monochromatic at energy level   511
+
+
+     <#> obje
+         description of objects
+                                                                density at levels
+ numb type  x-coord  y-coord x-length y-length    angle  av dens       511
+
+    1 elip  -7.0000  46.0000   3.0000   6.0000  17.0000   0.4845    0.4845
+
+    2 elip   7.0000  46.0000   3.0000   6.0000 -17.0000   0.5100    0.5100
+
+    3 rect -12.0000  64.0000   7.5000   4.5000   5.0000   0.5100    0.5100
+
+    4 rect  12.0000  64.0000   7.5000   4.5000  -5.0000   0.4845    0.4845
+
+    5 rect -38.0000  51.0000   3.5000  13.0000 -39.0000   0.4845    0.4845
+
+    6 rect  38.0000  51.0000   3.5000  13.0000  39.0000   0.5100    0.5100
+
+    7 rect -46.0000  24.0000   6.5000   6.0000 -18.0000   0.4845    0.4845
+
+    8 rect  46.0000  24.0000   6.5000   6.0000  18.0000   0.5100    0.5100
+
+    9 rect -49.0000   6.0000   2.5000  10.0000  63.0000   0.5100    0.5100
+
+   10 rect  49.0000   6.0000   2.5000  10.0000 -63.0000   0.4845    0.4845
+
+   11 rect -52.0000 -14.0000   9.0000   7.0000 -14.0000   0.4845    0.4845
+
+   12 rect  52.0000 -14.0000   9.0000   7.0000  14.0000   0.5100    0.5100
+
+   13 rect -10.0000 -56.0000   5.5000  10.0000  -1.0000   0.4845    0.4845
+
+   14 rect  10.0000 -56.0000   5.5000  10.0000   1.0000   0.5100    0.5100
+
+   15 elip -40.0000 -47.0000   9.0000  22.5000  48.0000   0.5100    0.5100
+
+   16 elip  40.0000 -47.0000   9.0000  22.5000 -48.0000   0.4845    0.4845
+
+   17 elip  -8.0000 -22.0000   3.5000  15.5000  -9.0000   0.5100    0.5100
+
+   18 elip   8.0000 -22.0000   3.5000  15.5000   9.0000   0.4845    0.4845
+
+   19 elip -27.0000  -6.0000   5.5000  23.5000  -5.0000   0.4845    0.4845
+
+   20 elip  27.0000  -6.0000   5.5000  23.5000   5.0000   0.5100    0.5100
+
+   21 elip -25.0000  38.0000   6.5000  10.5000 -14.0000   0.5100    0.5100
+
+   22 elip  25.0000  38.0000   6.5000  10.5000  14.0000   0.4845    0.4845
+
+   23 rect  -8.0000  32.0000   1.5000   6.5000  38.0000   0.5100    0.5100
+
+   24 rect   8.0000  32.0000   1.5000   6.5000 -38.0000   0.4845    0.4845
+
+   25 rect  -8.0000   3.0000   1.0000   9.0000 -33.0000   0.4845    0.4845
+
+   26 rect   8.0000   3.0000   1.0000   9.0000  33.0000   0.5100    0.5100
+
+   27 elip   0.0000   0.0000  66.5000  74.0000   0.0000   0.5100    0.5100
+
+         scale factor multiplying object densities     0.5100
+
+         seed set to 1
+         inhomogeneity set to     0.0500
+
+     <#> phan    aver    7
+
+         density in each pixel is obtained as the average of 7 x 7 points
+
+
+     <#> pixe      381    size        0.4000
+         picture size 381 x 381,  pixel size     0.4000
+
+         test picture read
+         EXAMPLE 11 Brain Phantom 
+              0.062 seconds used for processing command pict
+
+
+     <*> 
+
+     <#> PROJECTION REAL                                                                                                                                                                                                                                                  
+
+         EXAMPLE 11 Brain Phantom 
+
+
+     <#> spec    mono  511
+         energy spectrum is monochromatic at energy level   511
+
+
+     <#> obje
+         description of objects
+                                                                density at levels
+ numb type  x-coord  y-coord x-length y-length    angle  av dens       511
+
+    1 elip  -7.0000  46.0000   3.0000   6.0000  17.0000   0.4845    0.4845
+
+    2 elip   7.0000  46.0000   3.0000   6.0000 -17.0000   0.5100    0.5100
+
+    3 rect -12.0000  64.0000   7.5000   4.5000   5.0000   0.5100    0.5100
+
+    4 rect  12.0000  64.0000   7.5000   4.5000  -5.0000   0.4845    0.4845
+
+    5 rect -38.0000  51.0000   3.5000  13.0000 -39.0000   0.4845    0.4845
+
+    6 rect  38.0000  51.0000   3.5000  13.0000  39.0000   0.5100    0.5100
+
+    7 rect -46.0000  24.0000   6.5000   6.0000 -18.0000   0.4845    0.4845
+
+    8 rect  46.0000  24.0000   6.5000   6.0000  18.0000   0.5100    0.5100
+
+    9 rect -49.0000   6.0000   2.5000  10.0000  63.0000   0.5100    0.5100
+
+   10 rect  49.0000   6.0000   2.5000  10.0000 -63.0000   0.4845    0.4845
+
+   11 rect -52.0000 -14.0000   9.0000   7.0000 -14.0000   0.4845    0.4845
+
+   12 rect  52.0000 -14.0000   9.0000   7.0000  14.0000   0.5100    0.5100
+
+   13 rect -10.0000 -56.0000   5.5000  10.0000  -1.0000   0.4845    0.4845
+
+   14 rect  10.0000 -56.0000   5.5000  10.0000   1.0000   0.5100    0.5100
+
+   15 elip -40.0000 -47.0000   9.0000  22.5000  48.0000   0.5100    0.5100
+
+   16 elip  40.0000 -47.0000   9.0000  22.5000 -48.0000   0.4845    0.4845
+
+   17 elip  -8.0000 -22.0000   3.5000  15.5000  -9.0000   0.5100    0.5100
+
+   18 elip   8.0000 -22.0000   3.5000  15.5000   9.0000   0.4845    0.4845
+
+   19 elip -27.0000  -6.0000   5.5000  23.5000  -5.0000   0.4845    0.4845
+
+   20 elip  27.0000  -6.0000   5.5000  23.5000   5.0000   0.5100    0.5100
+
+   21 elip -25.0000  38.0000   6.5000  10.5000 -14.0000   0.5100    0.5100
+
+   22 elip  25.0000  38.0000   6.5000  10.5000  14.0000   0.4845    0.4845
+
+   23 rect  -8.0000  32.0000   1.5000   6.5000  38.0000   0.5100    0.5100
+
+   24 rect   8.0000  32.0000   1.5000   6.5000 -38.0000   0.4845    0.4845
+
+   25 rect  -8.0000   3.0000   1.0000   9.0000 -33.0000   0.4845    0.4845
+
+   26 rect   8.0000   3.0000   1.0000   9.0000  33.0000   0.5100    0.5100
+
+   27 elip   0.0000   0.0000  66.5000  74.0000   0.0000   0.5100    0.5100
+
+         scale factor multiplying object densities     0.5100
+
+         seed set to 1
+         inhomogeneity set to     0.0500
+
+     <#> rays    aver    1
+
+         projection data are calculated by dividing each ray interval into 1 substrips
+
+         with aperture (substrip) weights     1
+
+
+     <#> geom
+
+
+     <#> dive    arc     source at  153.0000     det dist  306.0000
+         rays are divergent from point sources
+         source to origin distance    153.0000
+         the detectors lie on an arc with source to detector distance =   306.0000
+
+
+     <#> rays    user      301    spacing        1.1000
+         number of rays per projection   301
+         snark computed number of rays   437
+         at detector spacing     1.1000
+
+
+     <#> angl      500
+         total number of projections   500
+
+
+         projection angles    0.0000    0.7190    1.4381    2.1571    2.8762    3.5952    4.3142    5.0333    5.7523    6.4713
+                              7.1904    7.9094    8.6285    9.3475   10.0665   10.7856   11.5046   12.2236   12.9427   13.6617
+                             14.3808   15.0998   15.8188   16.5379   17.2569   17.9760   18.6950   19.4140   20.1331   20.8521
+                             21.5711   22.2902   23.0092   23.7283   24.4473   25.1663   25.8854   26.6044   27.3234   28.0425
+                             28.7615   29.4806   30.1996   30.9186   31.6377   32.3567   33.0758   33.7948   34.5138   35.2329
+                             35.9519   36.6709   37.3900   38.1090   38.8281   39.5471   40.2661   40.9852   41.7042   42.4232
+                             43.1423   43.8613   44.5804   45.2994   46.0184   46.7375   47.4565   48.1756   48.8946   49.6136
+                             50.3327   51.0517   51.7707   52.4898   53.2088   53.9279   54.6469   55.3659   56.0850   56.8040
+                             57.5230   58.2421   58.9611   59.6802   60.3992   61.1182   61.8373   62.5563   63.2754   63.9944
+                             64.7134   65.4325   66.1515   66.8705   67.5896   68.3086   69.0277   69.7467   70.4657   71.1848
+                             71.9038   72.6228   73.3419   74.0609   74.7800   75.4990   76.2180   76.9371   77.6561   78.3752
+                             79.0942   79.8132   80.5323   81.2513   81.9703   82.6894   83.4084   84.1275   84.8465   85.5655
+                             86.2846   87.0036   87.7226   88.4417   89.1607   89.8798   90.5988   91.3178   92.0369   92.7559
+                             93.4749   94.1940   94.9130   95.6321   96.3511   97.0701   97.7892   98.5082   99.2273   99.9463
+                            100.6653  101.3844  102.1034  102.8224  103.5415  104.2605  104.9796  105.6986  106.4176  107.1367
+                            107.8557  108.5747  109.2938  110.0128  110.7319  111.4509  112.1699  112.8890  113.6080  114.3271
+                            115.0461  115.7651  116.4842  117.2032  117.9222  118.6413  119.3603  120.0794  120.7984  121.5174
+                            122.2365  122.9555  123.6745  124.3936  125.1126  125.8317  126.5507  127.2697  127.9888  128.7078
+                            129.4269  130.1459  130.8649  131.5840  132.3030  133.0220  133.7411  134.4601  135.1792  135.8982
+                            136.6172  137.3363  138.0553  138.7743  139.4934  140.2124  140.9315  141.6505  142.3695  143.0886
+                            143.8076  144.5267  145.2457  145.9647  146.6838  147.4028  148.1218  148.8409  149.5599  150.2790
+                            150.9980  151.7170  152.4361  153.1551  153.8741  154.5932  155.3122  156.0313  156.7503  157.4693
+                            158.1884  158.9074  159.6265  160.3455  161.0645  161.7836  162.5026  163.2216  163.9407  164.6597
+                            165.3788  166.0978  166.8168  167.5359  168.2549  168.9739  169.6930  170.4120  171.1311  171.8501
+                            172.5691  173.2882  174.0072  174.7263  175.4453  176.1643  176.8834  177.6024  178.3214  179.0405
+                            179.7595  180.4786  181.1976  181.9166  182.6357  183.3547  184.0737  184.7928  185.5118  186.2309
+                            186.9499  187.6689  188.3880  189.1070  189.8261  190.5451  191.2641  191.9832  192.7022  193.4212
+                            194.1403  194.8593  195.5784  196.2974  197.0164  197.7355  198.4545  199.1735  199.8926  200.6116
+                            201.3307  202.0497  202.7687  203.4878  204.2068  204.9259  205.6449  206.3639  207.0830  207.8020
+                            208.5210  209.2401  209.9591  210.6782  211.3972  212.1162  212.8353  213.5543  214.2733  214.9924
+                            215.7114  216.4305  217.1495  217.8685  218.5876  219.3066  220.0257  220.7447  221.4637  222.1828
+                            222.9018  223.6208  224.3399  225.0589  225.7780  226.4970  227.2160  227.9351  228.6541  229.3731
+                            230.0922  230.8112  231.5303  232.2493  232.9683  233.6874  234.4064  235.1255  235.8445  236.5635
+                            237.2826  238.0016  238.7206  239.4397  240.1587  240.8778  241.5968  242.3158  243.0349  243.7539
+                            244.4729  245.1920  245.9110  246.6301  247.3491  248.0681  248.7872  249.5062  250.2253  250.9443
+                            251.6633  252.3824  253.1014  253.8204  254.5395  255.2585  255.9776  256.6966  257.4156  258.1347
+                            258.8537  259.5727  260.2918  261.0108  261.7299  262.4489  263.1679  263.8870  264.6060  265.3251
+                            266.0441  266.7631  267.4822  268.2012  268.9202  269.6393  270.3583  271.0774  271.7964  272.5154
+                            273.2345  273.9535  274.6725  275.3916  276.1106  276.8297  277.5487  278.2677  278.9868  279.7058
+                            280.4248  281.1439  281.8629  282.5820  283.3010  284.0200  284.7391  285.4581  286.1772  286.8962
+                            287.6152  288.3343  289.0533  289.7723  290.4914  291.2104  291.9295  292.6485  293.3675  294.0866
+                            294.8056  295.5246  296.2437  296.9627  297.6818  298.4008  299.1198  299.8389  300.5579  301.2770
+                            301.9960  302.7150  303.4341  304.1531  304.8721  305.5912  306.3102  307.0293  307.7483  308.4673
+                            309.1864  309.9054  310.6244  311.3435  312.0625  312.7816  313.5006  314.2196  314.9387  315.6577
+                            316.3768  317.0958  317.8148  318.5339  319.2529  319.9719  320.6910  321.4100  322.1291  322.8481
+                            323.5671  324.2862  325.0052  325.7242  326.4433  327.1623  327.8814  328.6004  329.3194  330.0385
+                            330.7575  331.4766  332.1956  332.9146  333.6337  334.3527  335.0717  335.7908  336.5098  337.2289
+                            337.9479  338.6669  339.3860  340.1050  340.8240  341.5431  342.2621  342.9812  343.7002  344.4192
+                            345.1383  345.8573  346.5764  347.2954  348.0144  348.7335  349.4525  350.1715  350.8906  351.6096
+                            352.3287  353.0477  353.7667  354.4858  355.2048  355.9238  356.6429  357.3619  358.0810  358.8000
+
+
+     <#> meas    perf
+         projection data are noiseless
+
+
+     <#> back      0.0000
+                               at levels
+                                     511 
+          background absorption   0.0000 
+
+         estimate of totlen =  20962469.574466
+         estimate of totden =   9803802.737080
+         estimate of average density =     0.4677
+         projection data read
+         EXAMPLE 11 Brain Phantom 
+              0.094 seconds used for processing command proj
+
+
+     <*> 
+
+     <#> STOP TERMINATION KLDS 100000 RPRT                                                                                                                                                                                                                                
+         termination test klds
+         reporting is enabled
+         reporting file: RPRTklds
+         reporting on every iteration
+         epsilon = 100000
+              0.000 seconds used for processing command stop
+
+
+     <*> 
+
+     <#> SUPERIORIZE 16 0.999 1 SMOO RPRT                                                                                                                                                                                                                                 
+         Superiorization is enabled
+         N = 16
+         a = 0.999
+         b = 1
+         secondary criterion: smoo
+         reporting is enabled
+         reporting file: RPRTsuperiorization
+         reporting on every iteration
+              0.000 seconds used for processing command supe
+
+
+     <#> EXECUTE AVERAGE EMAP                                                                                                                                                                                                                                             
+
+         Example 11 Illustrating the Superiorized MAP EM algorithm for PET
+
+     <#> gamma is 0 EVAL
+
+  -----------------------------------------------------------
+
+   maximum a-posteriori probability expectation maximization
+
+            gamma:   0.000
+            evaluation flag is set 
+
+  -----------------------------------------------------------
+
+     *** WARNING - There exist negative values in the projection data ***
+     *** Negative values set to 0.0 ***
+              value of l: 15
+              value of phi before algorithm operator: 4.42622e-28
+              value of phi after algorithm operator:  2.39806
+         algorithm executed in iteration    1
+              1.734 seconds for the execution of the algorithm
+****Error: Negative raysum in prjfil detected, KLDS is not applicable! Terminating algorithm execution.
+         reconstruction completed after iteration    1
+              1.734 seconds for this iteration
+              1.922 seconds used for processing command exec
+
+
+     <*> 
+
+     <#> END                                                                                                                                                                                                                                                              

diff --git a/examples/b12/b12r-ubuntu201710.out b/examples/b12/b12r-ubuntu201710.out
new file mode 100644 (file)
index 0000000..e8f9220
--- /dev/null
+++ b/examples/b12/b12r-ubuntu201710.out
@@ -0,0 +1,505 @@
+   snark14.s170710 - A PICTURE RECONSTRUCTION PROGRAM
+
+
+     <*> *******EXAMPLE 11
+
+     <*>  SUPERIORIZED MAP EM ALGORITHM FOR EMISSION TOMOGRAPHY. RECONSTRUCTION
+
+     <*>  OF BRAIN PHANTOM. SIMULATING PET GEOMETRY WITH A RING OF 300 DETECTORS
+
+     <*>  WITH EACH DETECTOR IN COINCIDENCE WITH 101 DETECTORS OPPOSITE IT.
+
+     <*>  COEFFICIENT OF PENALTLY TERM IS SET TO 0
+
+     <*> 
+
+     <#> CREATE                                                                                                                                                                                                                                                           
+
+         EXAMPLE 11 Brain Phantom 
+
+
+     <#> SPECTRUM MONOCHROMATIC 511
+         energy spectrum is monochromatic at energy level   511
+
+
+     <#> OBJECTS
+         description of objects
+                                                                density at levels
+ numb type  x-coord  y-coord x-length y-length    angle  av dens       511
+
+    1 elip  -7.0000  46.0000   3.0000   6.0000  17.0000   0.9500    0.9500
+
+    2 elip   7.0000  46.0000   3.0000   6.0000 -17.0000   1.0000    1.0000
+
+    3 rect -12.0000  64.0000   7.5000   4.5000   5.0000   1.0000    1.0000
+
+    4 rect  12.0000  64.0000   7.5000   4.5000  -5.0000   0.9500    0.9500
+
+    5 rect -38.0000  51.0000   3.5000  13.0000 -39.0000   0.9500    0.9500
+
+    6 rect  38.0000  51.0000   3.5000  13.0000  39.0000   1.0000    1.0000
+
+    7 rect -46.0000  24.0000   6.5000   6.0000 -18.0000   0.9500    0.9500
+
+    8 rect  46.0000  24.0000   6.5000   6.0000  18.0000   1.0000    1.0000
+
+    9 rect -49.0000   6.0000   2.5000  10.0000  63.0000   1.0000    1.0000
+
+   10 rect  49.0000   6.0000   2.5000  10.0000 -63.0000   0.9500    0.9500
+
+   11 rect -52.0000 -14.0000   9.0000   7.0000 -14.0000   0.9500    0.9500
+
+   12 rect  52.0000 -14.0000   9.0000   7.0000  14.0000   1.0000    1.0000
+
+   13 rect -10.0000 -56.0000   5.5000  10.0000  -1.0000   0.9500    0.9500
+
+   14 rect  10.0000 -56.0000   5.5000  10.0000   1.0000   1.0000    1.0000
+
+   15 elip -40.0000 -47.0000   9.0000  22.5000  48.0000   1.0000    1.0000
+
+   16 elip  40.0000 -47.0000   9.0000  22.5000 -48.0000   0.9500    0.9500
+
+   17 elip  -8.0000 -22.0000   3.5000  15.5000  -9.0000   1.0000    1.0000
+
+   18 elip   8.0000 -22.0000   3.5000  15.5000   9.0000   0.9500    0.9500
+
+   19 elip -27.0000  -6.0000   5.5000  23.5000  -5.0000   0.9500    0.9500
+
+   20 elip  27.0000  -6.0000   5.5000  23.5000   5.0000   1.0000    1.0000
+
+   21 elip -25.0000  38.0000   6.5000  10.5000 -14.0000   1.0000    1.0000
+
+   22 elip  25.0000  38.0000   6.5000  10.5000  14.0000   0.9500    0.9500
+
+   23 rect  -8.0000  32.0000   1.5000   6.5000  38.0000   1.0000    1.0000
+
+   24 rect   8.0000  32.0000   1.5000   6.5000 -38.0000   0.9500    0.9500
+
+   25 rect  -8.0000   3.0000   1.0000   9.0000 -33.0000   0.9500    0.9500
+
+   26 rect   8.0000   3.0000   1.0000   9.0000  33.0000   1.0000    1.0000
+
+   27 elip   0.0000   0.0000  66.5000  74.0000   0.0000   1.0000    1.0000
+
+         scale factor multiplying object densities     0.5100
+
+         seed set to 1
+         inhomogeneity set to     0.0500
+
+     <#> PHANTOM AVERAGE 7
+
+         this run will generate a phantom
+         density in each pixel is obtained as the average of 7 x 7 points
+
+
+     <#> 381 PIXELS OF SIZE 0.4
+         picture size 381 x 381,  pixel size     0.4000
+
+
+     <#> RAYSUM AVERAGE 1
+
+         this run will generate projection data
+         projection data are calculated by dividing each ray interval into 1 substrips
+
+         with aperture (substrip) weights     1
+
+
+     <#> GEOMETRY
+
+
+     <#> divergent arc 153 306
+         rays are divergent from point sources
+         source to origin distance    153.0000
+         the detectors lie on an arc with source to detector distance =   306.0000
+
+
+     <#> RAYS USER 301 DETECTOR SPACING 1.1
+         number of rays per projection   301
+         at detector spacing     1.1000
+
+
+     <#> ANGLES 500 EQUAL SPACING
+         total number of projections   500
+
+         projection angles    0.0000    0.7190    1.4381    2.1571    2.8762    3.5952    4.3142    5.0333    5.7523    6.4713
+                              7.1904    7.9094    8.6285    9.3475   10.0665   10.7856   11.5046   12.2236   12.9427   13.6617
+                             14.3808   15.0998   15.8188   16.5379   17.2569   17.9760   18.6950   19.4140   20.1331   20.8521
+                             21.5711   22.2902   23.0092   23.7283   24.4473   25.1663   25.8854   26.6044   27.3234   28.0425
+                             28.7615   29.4806   30.1996   30.9186   31.6377   32.3567   33.0758   33.7948   34.5138   35.2329
+                             35.9519   36.6709   37.3900   38.1090   38.8281   39.5471   40.2661   40.9852   41.7042   42.4232
+                             43.1423   43.8613   44.5804   45.2994   46.0184   46.7375   47.4565   48.1756   48.8946   49.6136
+                             50.3327   51.0517   51.7707   52.4898   53.2088   53.9279   54.6469   55.3659   56.0850   56.8040
+                             57.5230   58.2421   58.9611   59.6802   60.3992   61.1182   61.8373   62.5563   63.2754   63.9944
+                             64.7134   65.4325   66.1515   66.8705   67.5896   68.3086   69.0277   69.7467   70.4657   71.1848
+                             71.9038   72.6228   73.3419   74.0609   74.7800   75.4990   76.2180   76.9371   77.6561   78.3752
+                             79.0942   79.8132   80.5323   81.2513   81.9703   82.6894   83.4084   84.1275   84.8465   85.5655
+                             86.2846   87.0036   87.7226   88.4417   89.1607   89.8798   90.5988   91.3178   92.0369   92.7559
+                             93.4749   94.1940   94.9130   95.6321   96.3511   97.0701   97.7892   98.5082   99.2273   99.9463
+                            100.6653  101.3844  102.1034  102.8224  103.5415  104.2605  104.9796  105.6986  106.4176  107.1367
+                            107.8557  108.5747  109.2938  110.0128  110.7319  111.4509  112.1699  112.8890  113.6080  114.3271
+                            115.0461  115.7651  116.4842  117.2032  117.9222  118.6413  119.3603  120.0794  120.7984  121.5174
+                            122.2365  122.9555  123.6745  124.3936  125.1126  125.8317  126.5507  127.2697  127.9888  128.7078
+                            129.4269  130.1459  130.8649  131.5840  132.3030  133.0220  133.7411  134.4601  135.1792  135.8982
+                            136.6172  137.3363  138.0553  138.7743  139.4934  140.2124  140.9315  141.6505  142.3695  143.0886
+                            143.8076  144.5267  145.2457  145.9647  146.6838  147.4028  148.1218  148.8409  149.5599  150.2790
+                            150.9980  151.7170  152.4361  153.1551  153.8741  154.5932  155.3122  156.0313  156.7503  157.4693
+                            158.1884  158.9074  159.6265  160.3455  161.0645  161.7836  162.5026  163.2216  163.9407  164.6597
+                            165.3788  166.0978  166.8168  167.5359  168.2549  168.9739  169.6930  170.4120  171.1311  171.8501
+                            172.5691  173.2882  174.0072  174.7263  175.4453  176.1643  176.8834  177.6024  178.3214  179.0405
+                            179.7595  180.4786  181.1976  181.9166  182.6357  183.3547  184.0737  184.7928  185.5118  186.2309
+                            186.9499  187.6689  188.3880  189.1070  189.8261  190.5451  191.2641  191.9832  192.7022  193.4212
+                            194.1403  194.8593  195.5784  196.2974  197.0164  197.7355  198.4545  199.1735  199.8926  200.6116
+                            201.3307  202.0497  202.7687  203.4878  204.2068  204.9259  205.6449  206.3639  207.0830  207.8020
+                            208.5210  209.2401  209.9591  210.6782  211.3972  212.1162  212.8353  213.5543  214.2733  214.9924
+                            215.7114  216.4305  217.1495  217.8685  218.5876  219.3066  220.0257  220.7447  221.4637  222.1828
+                            222.9018  223.6208  224.3399  225.0589  225.7780  226.4970  227.2160  227.9351  228.6541  229.3731
+                            230.0922  230.8112  231.5303  232.2493  232.9683  233.6874  234.4064  235.1255  235.8445  236.5635
+                            237.2826  238.0016  238.7206  239.4397  240.1587  240.8778  241.5968  242.3158  243.0349  243.7539
+                            244.4729  245.1920  245.9110  246.6301  247.3491  248.0681  248.7872  249.5062  250.2253  250.9443
+                            251.6633  252.3824  253.1014  253.8204  254.5395  255.2585  255.9776  256.6966  257.4156  258.1347
+                            258.8537  259.5727  260.2918  261.0108  261.7299  262.4489  263.1679  263.8870  264.6060  265.3251
+                            266.0441  266.7631  267.4822  268.2012  268.9202  269.6393  270.3583  271.0774  271.7964  272.5154
+                            273.2345  273.9535  274.6725  275.3916  276.1106  276.8297  277.5487  278.2677  278.9868  279.7058
+                            280.4248  281.1439  281.8629  282.5820  283.3010  284.0200  284.7391  285.4581  286.1772  286.8962
+                            287.6152  288.3343  289.0533  289.7723  290.4914  291.2104  291.9295  292.6485  293.3675  294.0866
+                            294.8056  295.5246  296.2437  296.9627  297.6818  298.4008  299.1198  299.8389  300.5579  301.2770
+                            301.9960  302.7150  303.4341  304.1531  304.8721  305.5912  306.3102  307.0293  307.7483  308.4673
+                            309.1864  309.9054  310.6244  311.3435  312.0625  312.7816  313.5006  314.2196  314.9387  315.6577
+                            316.3768  317.0958  317.8148  318.5339  319.2529  319.9719  320.6910  321.4100  322.1291  322.8481
+                            323.5671  324.2862  325.0052  325.7242  326.4433  327.1623  327.8814  328.6004  329.3194  330.0385
+                            330.7575  331.4766  332.1956  332.9146  333.6337  334.3527  335.0717  335.7908  336.5098  337.2289
+                            337.9479  338.6669  339.3860  340.1050  340.8240  341.5431  342.2621  342.9812  343.7002  344.4192
+                            345.1383  345.8573  346.5764  347.2954  348.0144  348.7335  349.4525  350.1715  350.8906  351.6096
+                            352.3287  353.0477  353.7667  354.4858  355.2048  355.9238  356.6429  357.3619  358.0810  358.8000
+
+
+     <#> MEASUREMENT PERFECT
+         projection data are noiseless
+
+
+     <#> BACKGROUND 0.0
+                               at levels
+                                     511 
+          background absorption   0.0000 
+ 
+
+     <#> RUN
+              0.214 seconds phantom creation
+              0.604 seconds projection data creation
+              0.820 seconds used for processing command crea
+
+
+     <*> 
+
+     <#> PICTURE TEST                                                                                                                                                                                                                                                     
+
+         EXAMPLE 11 Brain Phantom 
+
+
+     <#> spec    mono  511
+         energy spectrum is monochromatic at energy level   511
+
+
+     <#> obje
+         description of objects
+                                                                density at levels
+ numb type  x-coord  y-coord x-length y-length    angle  av dens       511
+
+    1 elip  -7.0000  46.0000   3.0000   6.0000  17.0000   0.4845    0.4845
+
+    2 elip   7.0000  46.0000   3.0000   6.0000 -17.0000   0.5100    0.5100
+
+    3 rect -12.0000  64.0000   7.5000   4.5000   5.0000   0.5100    0.5100
+
+    4 rect  12.0000  64.0000   7.5000   4.5000  -5.0000   0.4845    0.4845
+
+    5 rect -38.0000  51.0000   3.5000  13.0000 -39.0000   0.4845    0.4845
+
+    6 rect  38.0000  51.0000   3.5000  13.0000  39.0000   0.5100    0.5100
+
+    7 rect -46.0000  24.0000   6.5000   6.0000 -18.0000   0.4845    0.4845
+
+    8 rect  46.0000  24.0000   6.5000   6.0000  18.0000   0.5100    0.5100
+
+    9 rect -49.0000   6.0000   2.5000  10.0000  63.0000   0.5100    0.5100
+
+   10 rect  49.0000   6.0000   2.5000  10.0000 -63.0000   0.4845    0.4845
+
+   11 rect -52.0000 -14.0000   9.0000   7.0000 -14.0000   0.4845    0.4845
+
+   12 rect  52.0000 -14.0000   9.0000   7.0000  14.0000   0.5100    0.5100
+
+   13 rect -10.0000 -56.0000   5.5000  10.0000  -1.0000   0.4845    0.4845
+
+   14 rect  10.0000 -56.0000   5.5000  10.0000   1.0000   0.5100    0.5100
+
+   15 elip -40.0000 -47.0000   9.0000  22.5000  48.0000   0.5100    0.5100
+
+   16 elip  40.0000 -47.0000   9.0000  22.5000 -48.0000   0.4845    0.4845
+
+   17 elip  -8.0000 -22.0000   3.5000  15.5000  -9.0000   0.5100    0.5100
+
+   18 elip   8.0000 -22.0000   3.5000  15.5000   9.0000   0.4845    0.4845
+
+   19 elip -27.0000  -6.0000   5.5000  23.5000  -5.0000   0.4845    0.4845
+
+   20 elip  27.0000  -6.0000   5.5000  23.5000   5.0000   0.5100    0.5100
+
+   21 elip -25.0000  38.0000   6.5000  10.5000 -14.0000   0.5100    0.5100
+
+   22 elip  25.0000  38.0000   6.5000  10.5000  14.0000   0.4845    0.4845
+
+   23 rect  -8.0000  32.0000   1.5000   6.5000  38.0000   0.5100    0.5100
+
+   24 rect   8.0000  32.0000   1.5000   6.5000 -38.0000   0.4845    0.4845
+
+   25 rect  -8.0000   3.0000   1.0000   9.0000 -33.0000   0.4845    0.4845
+
+   26 rect   8.0000   3.0000   1.0000   9.0000  33.0000   0.5100    0.5100
+
+   27 elip   0.0000   0.0000  66.5000  74.0000   0.0000   0.5100    0.5100
+
+         scale factor multiplying object densities     0.5100
+
+         seed set to 1
+         inhomogeneity set to     0.0500
+
+     <#> phan    aver    7
+
+         density in each pixel is obtained as the average of 7 x 7 points
+
+
+     <#> pixe      381    size        0.4000
+         picture size 381 x 381,  pixel size     0.4000
+
+         test picture read
+         EXAMPLE 11 Brain Phantom 
+              0.043 seconds used for processing command pict
+
+
+     <*> 
+
+     <#> PROJECTION REAL                                                                                                                                                                                                                                                  
+
+         EXAMPLE 11 Brain Phantom 
+
+
+     <#> spec    mono  511
+         energy spectrum is monochromatic at energy level   511
+
+
+     <#> obje
+         description of objects
+                                                                density at levels
+ numb type  x-coord  y-coord x-length y-length    angle  av dens       511
+
+    1 elip  -7.0000  46.0000   3.0000   6.0000  17.0000   0.4845    0.4845
+
+    2 elip   7.0000  46.0000   3.0000   6.0000 -17.0000   0.5100    0.5100
+
+    3 rect -12.0000  64.0000   7.5000   4.5000   5.0000   0.5100    0.5100
+
+    4 rect  12.0000  64.0000   7.5000   4.5000  -5.0000   0.4845    0.4845
+
+    5 rect -38.0000  51.0000   3.5000  13.0000 -39.0000   0.4845    0.4845
+
+    6 rect  38.0000  51.0000   3.5000  13.0000  39.0000   0.5100    0.5100
+
+    7 rect -46.0000  24.0000   6.5000   6.0000 -18.0000   0.4845    0.4845
+
+    8 rect  46.0000  24.0000   6.5000   6.0000  18.0000   0.5100    0.5100
+
+    9 rect -49.0000   6.0000   2.5000  10.0000  63.0000   0.5100    0.5100
+
+   10 rect  49.0000   6.0000   2.5000  10.0000 -63.0000   0.4845    0.4845
+
+   11 rect -52.0000 -14.0000   9.0000   7.0000 -14.0000   0.4845    0.4845
+
+   12 rect  52.0000 -14.0000   9.0000   7.0000  14.0000   0.5100    0.5100
+
+   13 rect -10.0000 -56.0000   5.5000  10.0000  -1.0000   0.4845    0.4845
+
+   14 rect  10.0000 -56.0000   5.5000  10.0000   1.0000   0.5100    0.5100
+
+   15 elip -40.0000 -47.0000   9.0000  22.5000  48.0000   0.5100    0.5100
+
+   16 elip  40.0000 -47.0000   9.0000  22.5000 -48.0000   0.4845    0.4845
+
+   17 elip  -8.0000 -22.0000   3.5000  15.5000  -9.0000   0.5100    0.5100
+
+   18 elip   8.0000 -22.0000   3.5000  15.5000   9.0000   0.4845    0.4845
+
+   19 elip -27.0000  -6.0000   5.5000  23.5000  -5.0000   0.4845    0.4845
+
+   20 elip  27.0000  -6.0000   5.5000  23.5000   5.0000   0.5100    0.5100
+
+   21 elip -25.0000  38.0000   6.5000  10.5000 -14.0000   0.5100    0.5100
+
+   22 elip  25.0000  38.0000   6.5000  10.5000  14.0000   0.4845    0.4845
+
+   23 rect  -8.0000  32.0000   1.5000   6.5000  38.0000   0.5100    0.5100
+
+   24 rect   8.0000  32.0000   1.5000   6.5000 -38.0000   0.4845    0.4845
+
+   25 rect  -8.0000   3.0000   1.0000   9.0000 -33.0000   0.4845    0.4845
+
+   26 rect   8.0000   3.0000   1.0000   9.0000  33.0000   0.5100    0.5100
+
+   27 elip   0.0000   0.0000  66.5000  74.0000   0.0000   0.5100    0.5100
+
+         scale factor multiplying object densities     0.5100
+
+         seed set to 1
+         inhomogeneity set to     0.0500
+
+     <#> rays    aver    1
+
+         projection data are calculated by dividing each ray interval into 1 substrips
+
+         with aperture (substrip) weights     1
+
+
+     <#> geom
+
+
+     <#> dive    arc     source at  153.0000     det dist  306.0000
+         rays are divergent from point sources
+         source to origin distance    153.0000
+         the detectors lie on an arc with source to detector distance =   306.0000
+
+
+     <#> rays    user      301    spacing        1.1000
+         number of rays per projection   301
+         snark computed number of rays   437
+         at detector spacing     1.1000
+
+
+     <#> angl      500
+         total number of projections   500
+
+
+         projection angles    0.0000    0.7190    1.4381    2.1571    2.8762    3.5952    4.3142    5.0333    5.7523    6.4713
+                              7.1904    7.9094    8.6285    9.3475   10.0665   10.7856   11.5046   12.2236   12.9427   13.6617
+                             14.3808   15.0998   15.8188   16.5379   17.2569   17.9760   18.6950   19.4140   20.1331   20.8521
+                             21.5711   22.2902   23.0092   23.7283   24.4473   25.1663   25.8854   26.6044   27.3234   28.0425
+                             28.7615   29.4806   30.1996   30.9186   31.6377   32.3567   33.0758   33.7948   34.5138   35.2329
+                             35.9519   36.6709   37.3900   38.1090   38.8281   39.5471   40.2661   40.9852   41.7042   42.4232
+                             43.1423   43.8613   44.5804   45.2994   46.0184   46.7375   47.4565   48.1756   48.8946   49.6136
+                             50.3327   51.0517   51.7707   52.4898   53.2088   53.9279   54.6469   55.3659   56.0850   56.8040
+                             57.5230   58.2421   58.9611   59.6802   60.3992   61.1182   61.8373   62.5563   63.2754   63.9944
+                             64.7134   65.4325   66.1515   66.8705   67.5896   68.3086   69.0277   69.7467   70.4657   71.1848
+                             71.9038   72.6228   73.3419   74.0609   74.7800   75.4990   76.2180   76.9371   77.6561   78.3752
+                             79.0942   79.8132   80.5323   81.2513   81.9703   82.6894   83.4084   84.1275   84.8465   85.5655
+                             86.2846   87.0036   87.7226   88.4417   89.1607   89.8798   90.5988   91.3178   92.0369   92.7559
+                             93.4749   94.1940   94.9130   95.6321   96.3511   97.0701   97.7892   98.5082   99.2273   99.9463
+                            100.6653  101.3844  102.1034  102.8224  103.5415  104.2605  104.9796  105.6986  106.4176  107.1367
+                            107.8557  108.5747  109.2938  110.0128  110.7319  111.4509  112.1699  112.8890  113.6080  114.3271
+                            115.0461  115.7651  116.4842  117.2032  117.9222  118.6413  119.3603  120.0794  120.7984  121.5174
+                            122.2365  122.9555  123.6745  124.3936  125.1126  125.8317  126.5507  127.2697  127.9888  128.7078
+                            129.4269  130.1459  130.8649  131.5840  132.3030  133.0220  133.7411  134.4601  135.1792  135.8982
+                            136.6172  137.3363  138.0553  138.7743  139.4934  140.2124  140.9315  141.6505  142.3695  143.0886
+                            143.8076  144.5267  145.2457  145.9647  146.6838  147.4028  148.1218  148.8409  149.5599  150.2790
+                            150.9980  151.7170  152.4361  153.1551  153.8741  154.5932  155.3122  156.0313  156.7503  157.4693
+                            158.1884  158.9074  159.6265  160.3455  161.0645  161.7836  162.5026  163.2216  163.9407  164.6597
+                            165.3788  166.0978  166.8168  167.5359  168.2549  168.9739  169.6930  170.4120  171.1311  171.8501
+                            172.5691  173.2882  174.0072  174.7263  175.4453  176.1643  176.8834  177.6024  178.3214  179.0405
+                            179.7595  180.4786  181.1976  181.9166  182.6357  183.3547  184.0737  184.7928  185.5118  186.2309
+                            186.9499  187.6689  188.3880  189.1070  189.8261  190.5451  191.2641  191.9832  192.7022  193.4212
+                            194.1403  194.8593  195.5784  196.2974  197.0164  197.7355  198.4545  199.1735  199.8926  200.6116
+                            201.3307  202.0497  202.7687  203.4878  204.2068  204.9259  205.6449  206.3639  207.0830  207.8020
+                            208.5210  209.2401  209.9591  210.6782  211.3972  212.1162  212.8353  213.5543  214.2733  214.9924
+                            215.7114  216.4305  217.1495  217.8685  218.5876  219.3066  220.0257  220.7447  221.4637  222.1828
+                            222.9018  223.6208  224.3399  225.0589  225.7780  226.4970  227.2160  227.9351  228.6541  229.3731
+                            230.0922  230.8112  231.5303  232.2493  232.9683  233.6874  234.4064  235.1255  235.8445  236.5635
+                            237.2826  238.0016  238.7206  239.4397  240.1587  240.8778  241.5968  242.3158  243.0349  243.7539
+                            244.4729  245.1920  245.9110  246.6301  247.3491  248.0681  248.7872  249.5062  250.2253  250.9443
+                            251.6633  252.3824  253.1014  253.8204  254.5395  255.2585  255.9776  256.6966  257.4156  258.1347
+                            258.8537  259.5727  260.2918  261.0108  261.7299  262.4489  263.1679  263.8870  264.6060  265.3251
+                            266.0441  266.7631  267.4822  268.2012  268.9202  269.6393  270.3583  271.0774  271.7964  272.5154
+                            273.2345  273.9535  274.6725  275.3916  276.1106  276.8297  277.5487  278.2677  278.9868  279.7058
+                            280.4248  281.1439  281.8629  282.5820  283.3010  284.0200  284.7391  285.4581  286.1772  286.8962
+                            287.6152  288.3343  289.0533  289.7723  290.4914  291.2104  291.9295  292.6485  293.3675  294.0866
+                            294.8056  295.5246  296.2437  296.9627  297.6818  298.4008  299.1198  299.8389  300.5579  301.2770
+                            301.9960  302.7150  303.4341  304.1531  304.8721  305.5912  306.3102  307.0293  307.7483  308.4673
+                            309.1864  309.9054  310.6244  311.3435  312.0625  312.7816  313.5006  314.2196  314.9387  315.6577
+                            316.3768  317.0958  317.8148  318.5339  319.2529  319.9719  320.6910  321.4100  322.1291  322.8481
+                            323.5671  324.2862  325.0052  325.7242  326.4433  327.1623  327.8814  328.6004  329.3194  330.0385
+                            330.7575  331.4766  332.1956  332.9146  333.6337  334.3527  335.0717  335.7908  336.5098  337.2289
+                            337.9479  338.6669  339.3860  340.1050  340.8240  341.5431  342.2621  342.9812  343.7002  344.4192
+                            345.1383  345.8573  346.5764  347.2954  348.0144  348.7335  349.4525  350.1715  350.8906  351.6096
+                            352.3287  353.0477  353.7667  354.4858  355.2048  355.9238  356.6429  357.3619  358.0810  358.8000
+
+
+     <#> meas    perf
+         projection data are noiseless
+
+
+     <#> back      0.0000
+                               at levels
+                                     511 
+          background absorption   0.0000 
+
+         estimate of totlen =  20962469.574466
+         estimate of totden =   9804758.593222
+         estimate of average density =     0.4677
+         projection data read
+         EXAMPLE 11 Brain Phantom 
+              0.070 seconds used for processing command proj
+
+
+     <*> 
+
+     <#> STOP TERMINATION KLDS 100000 RPRT                                                                                                                                                                                                                                
+         termination test klds
+         reporting is enabled
+         reporting file: RPRTklds
+         reporting on every iteration
+         epsilon = 100000
+              0.000 seconds used for processing command stop
+
+
+     <*> 
+
+     <#> SUPERIORIZE 16 0.999 1 SMOO RPRT                                                                                                                                                                                                                                 
+         Superiorization is enabled
+         N = 16
+         a = 0.999
+         b = 1
+         secondary criterion: smoo
+         reporting is enabled
+         reporting file: RPRTsuperiorization
+         reporting on every iteration
+              0.000 seconds used for processing command supe
+
+
+     <#> EXECUTE AVERAGE EMAP                                                                                                                                                                                                                                             
+
+         Example 11 Illustrating the Superiorized MAP EM algorithm for PET
+
+     <#> gamma is 0 EVAL
+
+  -----------------------------------------------------------
+
+   maximum a-posteriori probability expectation maximization
+
+            gamma:   0.000
+            evaluation flag is set 
+
+  -----------------------------------------------------------
+
+     *** WARNING - There exist negative values in the projection data ***
+     *** Negative values set to 0.0 ***
+              value of l: 15
+              value of phi before algorithm operator: 0
+              value of phi after algorithm operator:  2.39758
+         algorithm executed in iteration    1
+              1.266 seconds for the execution of the algorithm
+****Error: Negative raysum in prjfil detected, KLDS is not applicable! Terminating algorithm execution.
+         reconstruction completed after iteration    1
+              1.268 seconds for this iteration
+              1.392 seconds used for processing command exec
+
+
+     <*> 
+
+     <#> END                                                                                                                                                                                                                                                              

diff --git a/examples/b13/b13r.out b/examples/b13/b13r.out
index e419bf7a7a6ddd308af6598749b96f02084ba38d..8eeb4b8400759b9ae78614a93221b3966e57ad51 100644 (file)
--- a/examples/b13/b13r.out
+++ b/examples/b13/b13r.out
@@ -164,9 +164,9 @@
  
 
      <#> RUN
  
 
      <#> RUN

-              0.014 seconds phantom creation
-              0.078 seconds projection data creation
-              0.093 seconds used for processing command crea
+              0.031 seconds phantom creation
+              0.109 seconds projection data creation
+              0.141 seconds used for processing command crea

 
 
      <*> 
 
 
      <*> 


@@ -254,7 +254,7 @@
 
          test picture read
          EXAMPLE 7 Brain Phantom 
 
          test picture read
          EXAMPLE 7 Brain Phantom 

-              0.003 seconds used for processing command pict
+              0.000 seconds used for processing command pict

 
 
      <*> 
 
 
      <*> 


@@ -400,11 +400,11 @@
           background absorption   0.0000 
 
          estimate of totlen =   4263169.769963
           background absorption   0.0000 
 
          estimate of totlen =   4263169.769963

-         estimate of totden =   2023083.983824
-         estimate of average density =     0.4745
+         estimate of totden =   2020712.203355
+         estimate of average density =     0.4740

          projection data read
          EXAMPLE 7 Brain Phantom 
          projection data read
          EXAMPLE 7 Brain Phantom 

-              0.017 seconds used for processing command proj
+              0.016 seconds used for processing command proj

 
 
      <*> 
 
 
      <*> 


@@ -434,87 +434,87 @@
      *** WARNING - There exist negative values in the projection data ***
      *** Negative values set to 0.0 ***
          algorithm executed in iteration    1
      *** WARNING - There exist negative values in the projection data ***
      *** Negative values set to 0.0 ***
          algorithm executed in iteration    1

-              0.063 seconds for the execution of the algorithm
+              0.172 seconds for the execution of the algorithm

          iteration    1 completed
          iteration    1 completed

-              0.064 seconds for this iteration
+              0.172 seconds for this iteration

          algorithm executed in iteration    2
          algorithm executed in iteration    2

-              0.040 seconds for the execution of the algorithm
+              0.062 seconds for the execution of the algorithm

          iteration    2 completed
          iteration    2 completed

-              0.040 seconds for this iteration
+              0.062 seconds for this iteration

          algorithm executed in iteration    3
          algorithm executed in iteration    3

-              0.040 seconds for the execution of the algorithm
+              0.031 seconds for the execution of the algorithm

          iteration    3 completed
          iteration    3 completed

-              0.040 seconds for this iteration
+              0.031 seconds for this iteration

          algorithm executed in iteration    4
          algorithm executed in iteration    4

-              0.040 seconds for the execution of the algorithm
+              0.094 seconds for the execution of the algorithm

          iteration    4 completed
          iteration    4 completed

-              0.040 seconds for this iteration
+              0.094 seconds for this iteration

          algorithm executed in iteration    5
          algorithm executed in iteration    5

-              0.041 seconds for the execution of the algorithm
+              0.094 seconds for the execution of the algorithm

          iteration    5 completed
          iteration    5 completed

-              0.041 seconds for this iteration
+              0.094 seconds for this iteration

          algorithm executed in iteration    6
          algorithm executed in iteration    6

-              0.042 seconds for the execution of the algorithm
+              0.047 seconds for the execution of the algorithm

          iteration    6 completed
          iteration    6 completed

-              0.042 seconds for this iteration
+              0.047 seconds for this iteration

          algorithm executed in iteration    7
          algorithm executed in iteration    7

-              0.044 seconds for the execution of the algorithm
+              0.109 seconds for the execution of the algorithm

          iteration    7 completed
          iteration    7 completed

-              0.044 seconds for this iteration
+              0.109 seconds for this iteration

          algorithm executed in iteration    8
          algorithm executed in iteration    8

-              0.041 seconds for the execution of the algorithm
+              0.047 seconds for the execution of the algorithm

          iteration    8 completed
          iteration    8 completed

-              0.041 seconds for this iteration
+              0.047 seconds for this iteration

          algorithm executed in iteration    9
          algorithm executed in iteration    9

-              0.040 seconds for the execution of the algorithm
+              0.031 seconds for the execution of the algorithm

          iteration    9 completed
          iteration    9 completed

-              0.040 seconds for this iteration
+              0.031 seconds for this iteration

          algorithm executed in iteration   10
          algorithm executed in iteration   10

-              0.041 seconds for the execution of the algorithm
+              0.031 seconds for the execution of the algorithm

          iteration   10 completed
          iteration   10 completed

-              0.041 seconds for this iteration
+              0.031 seconds for this iteration

          algorithm executed in iteration   11
          algorithm executed in iteration   11

-              0.040 seconds for the execution of the algorithm
+              0.062 seconds for the execution of the algorithm

          iteration   11 completed
          iteration   11 completed

-              0.040 seconds for this iteration
+              0.062 seconds for this iteration

          algorithm executed in iteration   12
          algorithm executed in iteration   12

-              0.040 seconds for the execution of the algorithm
+              0.094 seconds for the execution of the algorithm

          iteration   12 completed
          iteration   12 completed

-              0.040 seconds for this iteration
+              0.094 seconds for this iteration

          algorithm executed in iteration   13
          algorithm executed in iteration   13

-              0.040 seconds for the execution of the algorithm
+              0.062 seconds for the execution of the algorithm

          iteration   13 completed
          iteration   13 completed

-              0.040 seconds for this iteration
+              0.062 seconds for this iteration

          algorithm executed in iteration   14
          algorithm executed in iteration   14

-              0.040 seconds for the execution of the algorithm
+              0.031 seconds for the execution of the algorithm

          iteration   14 completed
          iteration   14 completed

-              0.040 seconds for this iteration
+              0.031 seconds for this iteration

          algorithm executed in iteration   15
          algorithm executed in iteration   15

-              0.040 seconds for the execution of the algorithm
+              0.062 seconds for the execution of the algorithm

          iteration   15 completed
          iteration   15 completed

-              0.040 seconds for this iteration
+              0.062 seconds for this iteration

          algorithm executed in iteration   16
          algorithm executed in iteration   16

-              0.040 seconds for the execution of the algorithm
+              0.078 seconds for the execution of the algorithm

          iteration   16 completed
          iteration   16 completed

-              0.041 seconds for this iteration
+              0.078 seconds for this iteration

          algorithm executed in iteration   17
          algorithm executed in iteration   17

-              0.040 seconds for the execution of the algorithm
+              0.078 seconds for the execution of the algorithm

          iteration   17 completed
          iteration   17 completed

-              0.040 seconds for this iteration
+              0.078 seconds for this iteration

          algorithm executed in iteration   18
          algorithm executed in iteration   18

-              0.040 seconds for the execution of the algorithm
+              0.062 seconds for the execution of the algorithm

          iteration   18 completed
          iteration   18 completed

-              0.040 seconds for this iteration
+              0.062 seconds for this iteration

          algorithm executed in iteration   19
          algorithm executed in iteration   19

-              0.041 seconds for the execution of the algorithm
+              0.078 seconds for the execution of the algorithm

          iteration   19 completed
          iteration   19 completed

-              0.041 seconds for this iteration
+              0.078 seconds for this iteration

          algorithm executed in iteration   20
          algorithm executed in iteration   20

-              0.041 seconds for the execution of the algorithm
+              0.062 seconds for the execution of the algorithm

          reconstruction completed after iteration   20
          reconstruction completed after iteration   20

-              0.041 seconds for this iteration
-              0.836 seconds for all iterations
-              0.883 seconds used for processing command exec
+              0.062 seconds for this iteration
+              1.391 seconds for all iterations
+              1.453 seconds used for processing command exec

 
 
      <*> 
 
 
      <*> 

diff --git a/examples/regression b/examples/regression
index fb4c2f1db5e9416e222b520e70b1454d9ed58e1f..524fbc64196d5c042ecdb99bddc902f725b84f7d 100755 (executable)
--- a/examples/regression
+++ b/examples/regression
@@ -1,8 +1,9 @@
 #!/bin/bash
 
 # Exclude examples 7 and 11 as they use random number generator not compatible across platforms
 #!/bin/bash
 
 # Exclude examples 7 and 11 as they use random number generator not compatible across platforms

+# Exclude example 12 as small differences in phi across gcc5.4 vs gcc7.2

 
 

-for i in 1 2 3 4 5 6 8 9 10 12 13; 
+for i in 1 2 3 4 5 6 8 9 10 13; 

 do
        echo Computing the diff for example b"$i".
        cd b"$i"
 do
        echo Computing the diff for example b"$i".
        cd b"$i"