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