From 93394341ab3f11455fc247c8e4e5050b965d8c49 Mon Sep 17 00:00:00 2001 From: "Kevin M. Rosenberg" Date: Mon, 12 Feb 2018 15:26:08 -0700 Subject: [PATCH] Exclude examples with random number generators from regression. Add new examples for same test, but with MEASUREMENT PERFECT. --- examples/b12/README | 6 + examples/b12/b12r-centos-orig.out | 505 +++++++++++++++++++++++++++++ examples/b12/b12r-ubuntu201710.out | 505 +++++++++++++++++++++++++++++ examples/b13/b13r.out | 98 +++--- examples/regression | 3 +- 5 files changed, 1067 insertions(+), 50 deletions(-) create mode 100644 examples/b12/b12r-centos-orig.out create mode 100644 examples/b12/b12r-ubuntu201710.out diff --git a/examples/b12/README b/examples/b12/README index 9127728..156db76 100644 --- 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. +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 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 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 e419bf7..8eeb4b8 100644 --- a/examples/b13/b13r.out +++ b/examples/b13/b13r.out @@ -164,9 +164,9 @@ <#> 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 - 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 - 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 - 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 - 0.063 seconds for the execution of the algorithm + 0.172 seconds for the execution of the algorithm iteration 1 completed - 0.064 seconds for this iteration + 0.172 seconds for this iteration 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 - 0.040 seconds for this iteration + 0.062 seconds for this iteration 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 - 0.040 seconds for this iteration + 0.031 seconds for this iteration 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 - 0.040 seconds for this iteration + 0.094 seconds for this iteration 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 - 0.041 seconds for this iteration + 0.094 seconds for this iteration 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 - 0.042 seconds for this iteration + 0.047 seconds for this iteration 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 - 0.044 seconds for this iteration + 0.109 seconds for this iteration 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 - 0.041 seconds for this iteration + 0.047 seconds for this iteration 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 - 0.040 seconds for this iteration + 0.031 seconds for this iteration 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 - 0.041 seconds for this iteration + 0.031 seconds for this iteration 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 - 0.040 seconds for this iteration + 0.062 seconds for this iteration 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 - 0.040 seconds for this iteration + 0.094 seconds for this iteration 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 - 0.040 seconds for this iteration + 0.062 seconds for this iteration 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 - 0.040 seconds for this iteration + 0.031 seconds for this iteration 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 - 0.040 seconds for this iteration + 0.062 seconds for this iteration 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 - 0.041 seconds for this iteration + 0.078 seconds for this iteration 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 - 0.040 seconds for this iteration + 0.078 seconds for this iteration 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 - 0.040 seconds for this iteration + 0.062 seconds for this iteration 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 - 0.041 seconds for this iteration + 0.078 seconds for this iteration 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 - 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 fb4c2f1..524fbc6 100755 --- 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 +# 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" -- 2.34.1