--- a/img/brainphantom256x256_pdps_known_dualscaling_20_10000_9000.txt Wed Apr 24 16:55:04 2024 +0300 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,78 +0,0 @@ -# params = (ρ = 0, verbose_iter = 100, maxiter = 4000, save_results = true, save_images = true, save_images_iters = Set([500, 1800, 2300, 3000, 1300, 3800, 1000, 2800, 300, 4000, 3500, 2000, 2500, 3300, 800, 100, 1500]), pixelwise_displacement = false, dual_flow = true, prox_predict = true, handle_interrupt = true, init = :zero, plot_movement = false, stable_interval = Set([1144, 3634, 1175, 1953, 1546, 1703, 1956, 3601, 1028, 1299, 1438, 1074, 3672, 3705, 3631, 1812, 1199, 1090, 1437, 1662, 1370, 1991, 3930, 3608, 1704, 1663, 1957, 1266, 1646, 1823, 1015, 1373, 3981, 3719, 3646, 3982, 3655, 1168, 3855, 1039, 1085, 3914, 3794, 1050, 1276, 1285, 3540, 3524, 1356, 3762, 1466, 3711, 1036, 3549, 1171, 1739, 1625, 1669, 1058, 1814, 3963, 1102, 1388, 1905, 3765, 3943, 1570, 3651, 1745, 1574, 1864, 1117, 1859, 1338, 1636, 3915, 1861, 3754, 3570, 1896, 3510, 1075, 1342, 3557, 1725, 1873, 1526, 1769, 1185, 1928, 3937, 1922, 1939, 1504, 1052, 1374, 1807, 1724, 3587, 1245, 3791, 3652, 3628, 1138, 1095, 1604, 1848, 1155, 1111, 1916, 1906, 1251, 3938, 1341, 1753, 1278, 1221, 1621, 3731, 1497, 1819, 1845, 3728, 3809, 3677, 3617, 1521, 3713, 1539, 1670, 3803, 1672, 1808, 1661, 3729, 1056, 1610, 1236, 3504, 3571, 3743, 3722, 3507, 1070, 1994, 1548, 1357, 1903, 1580, 1678, 1004, 1091, 1147, 1045, 1255, 1635, 3526, 1001, 3786, 1972, 1274, 1336, 3940, 3792, 1157, 1836, 1073, 1612, 1282, 3835, 1448, 3856, 3757, 3740, 3648, 3815, 1578, 3691, 1304, 2000, 1372, 1849, 3877, 3979, 1475, 1585, 1783, 3859, 1378, 3602, 1719, 3614, 1259, 1528, 3956, 1325, 1287, 1935, 3574, 1019, 1913, 3770, 1354, 1238, 1559, 1249, 1818, 1021, 1519, 3844, 1201, 1313, 3851, 1362, 1326, 3696, 1449, 1054, 3825, 1120, 1908, 3657, 1923, 1302, 1637, 3690, 1232, 1395, 1344, 3965, 1698, 3530, 1208, 1811, 1333, 1941, 3817, 1458, 1865, 3886, 1116, 3931, 3926, 1268, 1732, 3568, 3828, 1714, 1494, 1846, 3702, 3923, 1642, 1398, 3616, 1883, 1079, 1506, 1186, 1582, 3773, 3806, 1043, 1535, 3858, 1872, 1103, 1319, 3972, 1246, 1152, 1409, 1560, 1996, 3903, 1145, 1767, 1680, 1792, 1141, 1009, 1566, 1137, 1360, 1565, 1854, 1588, 1538, 3519, 1315, 1681, 1723, 3681, 3917, 1283, 3509, 1825, 3868, 3870, 3819, 1770, 3593, 1063, 1235, 1614, 1166, 1064, 3641, 1355, 1151, 1169, 3700, 3973, 3986, 4000, 3784, 1452, 1479, 3962, 3508, 1576, 1830, 3644, 1183, 3599, 1842, 1377, 3566, 1046, 1564, 1493, 1180, 3550, 1007, 3837, 1256, 1993, 3850, 1759, 1998, 3922, 3697, 1909, 3609, 1040, 1937, 1728, 3525, 3990, 3759, 1384, 1562, 3658, 3925, 3971, 3649, 1675, 3928, 1600, 3693, 1543, 1172, 1135, 1607, 3737, 3977, 3852, 1966, 1471, 1618, 3680, 1927, 1309, 1529, 1797, 3764, 3909, 1664, 1999, 3801, 1476, 1403, 1041, 1839, 1813, 1595, 1900, 3620, 1442, 1296, 3751, 3505, 1123, 3832, 1743, 1942, 1228, 1381, 3501, 1899, 1901, 3703, 1727, 1177, 1850, 3613, 1393, 1886, 1434, 1469, 3880, 1591, 1673, 1027, 1581, 1002, 1113, 3882, 3686, 1391, 1260, 1109, 1383, 3622, 3576, 1429, 1435, 3556, 1634, 3861, 1676, 1012, 1987, 1730, 1425, 1347, 1620, 1831, 3821, 3905, 1198, 3952, 3625, 1140, 1404, 3955, 1720, 1414, 3901, 1875, 3684, 3933, 1544, 3692, 1394, 1761, 1006, 1748, 3541, 3788, 1207, 3862, 1523, 1067, 3579, 3727, 1284, 1044, 1066, 1488, 1520, 1961, 1016, 1617, 1763, 1706, 1159, 3590, 1337, 1781, 1771, 1766, 3607, 1912, 1483, 3747, 3920, 1252, 1167, 1324, 1793, 1294, 1127, 1192, 1352, 3753, 1944, 3775, 1988, 3744, 1156, 1895, 1921, 1917, 1170, 1631, 3582, 1068, 1701, 1542, 3793, 3569, 1705, 1345, 1597, 1321, 1290, 1666, 1189, 3538, 3633, 1623, 1505, 1744, 1622, 1047, 1834, 3769, 3534, 3707, 1779, 1218, 3584, 3807, 1272, 1958, 1639, 1982, 1132, 1816, 1667, 1887, 1379, 3567, 3585, 1173, 1051, 1729, 1303, 1507, 1650, 1205, 1862, 1871, 3954, 1463, 1656, 1777, 1918, 3547, 3627, 1240, 1334, 1653, 3708, 1129, 3823, 1416, 1482, 1624, 1413, 3799, 1933, 1665, 1086, 1496, 1555, 1243, 3748, 3522, 1866, 3863, 3883, 3983, 1525, 1204, 1380, 1225, 1919, 1550, 1035, 1470, 1184, 1649, 1216, 3725, 1194, 1693, 3517, 3621, 1963, 1048, 1514, 1660, 1440, 1231, 1545, 1286, 3798, 1758, 1125, 1584, 1700, 1852, 1332, 3826, 1220, 1415, 1734, 1615, 3812, 1951, 3527, 3985, 1776, 1267, 1472, 1013, 3756, 1323, 3957, 1754, 3726, 1485, 1930, 1212, 1645, 1997, 1134, 1847, 1322, 3787, 1311, 1551, 1078, 1081, 1101, 1490, 1110, 1764, 3660, 1003, 1697, 1230, 3736, 1197, 1643, 1911, 1298, 1527, 3811, 1242, 1420, 3839, 1696, 1474, 3561, 1806, 3797, 1568, 1269, 1348, 3624, 3539, 1209, 1153, 3642, 3988, 3704, 1626, 3896, 1737, 1178, 1486, 3712, 1423, 1353, 1890, 3518, 3885, 3996, 3721, 1014, 1694, 1860, 3818, 1952, 1387, 1879, 1124, 3829, 3664, 3595, 3610, 1695, 3667, 1461, 1804, 1567, 3827, 1389, 1587, 3881, 1328, 1182, 1349, 1386, 1426, 1445, 3654, 1306, 1112, 1782, 1084, 1513, 1467, 1736, 1133, 1262, 1931, 1340, 3942, 1794, 1060, 3558, 1264, 1271, 1596, 3783, 1962, 1162, 3560, 3929, 1833, 1508, 3647, 1316, 3813, 1291, 1206, 1613, 3531, 1837, 1281, 1740, 3840, 1799, 3932, 1417, 1428, 1258, 1254, 1358, 1707, 1577, 3750, 3968, 1072, 1686, 3678, 3573, 1108, 3589, 1644, 1750, 1838, 3565, 3854, 1973, 1853, 3715, 3594, 3842, 1751, 1950, 3669, 1011, 3998, 1008, 1396, 1093, 1983, 1148, 1450, 1314, 1176, 1583, 3867, 1757, 3846, 3626, 1843, 1010, 1422, 3830, 3710, 1503, 3778, 1457, 1722, 3910, 1965, 3911, 1468, 1616, 1191, 1293, 1351, 1094, 3913, 1515, 3665, 1844, 1401, 1193, 1037, 1092, 1932, 1877, 3761, 3875, 1477, 3884, 1366, 1335, 1297, 3777, 3969, 1702, 1361, 1088, 1473, 1594, 3895, 3572, 3802, 1510, 1163, 1301, 1142, 1265, 1684, 1975, 1022, 1224, 1139, 1436, 3772, 1752, 3833, 3675, 1233, 1214, 1237, 1465, 1652, 1367, 3987, 1690, 1893, 3597, 1431, 1130, 1055, 1870, 3745, 1000, 1305, 1970, 3831, 1257, 1629, 1460, 1196, 1363, 1876, 1327, 3536, 1190, 1885, 1762, 3716, 1924, 1029, 1627, 1891, 3662, 1076, 1034, 1061, 1462, 1599, 3637, 1682, 1119, 1765, 1534, 1026, 1312, 1948, 1943, 1598, 1146, 1300, 1484, 1042, 1784, 1331, 1097, 1310, 1824, 1203, 1499, 1685, 1810, 1082, 3876, 3580, 3653, 1721, 3537, 1914, 1747, 1657, 1986, 1715, 3709, 1981, 3992, 1934, 1980, 1679, 1571, 3927, 1089, 3623, 3533, 3650, 1279, 1018, 3553, 3741, 3984, 1775, 1925, 1495, 3735, 1023, 3600, 1211, 3611, 1822, 1524, 1829, 3892, 1926, 3699, 1418, 1031, 1080, 1840, 1263, 3908, 1270, 1756, 3659, 1248, 3562, 1350, 3544, 1537, 3874, 3894, 3897, 3532, 1478, 1215, 1411, 1553, 1638, 1038, 1364, 3959, 1406, 3635, 1541, 1976, 1441, 1407, 1710, 3845, 1516, 3780, 3666, 3824, 1057, 1857, 3720, 1408, 3888, 1330, 1894, 1114, 3763, 1531, 1969, 3934, 3906, 1785, 1558, 1726, 1960, 3790, 1532, 3695, 3995, 1795, 1716, 1668, 3596, 3732, 3618, 1907, 3563, 1502, 1511, 1633, 3918, 1791, 3630, 3766, 1371, 3948, 1802, 1946, 1187, 1421, 1552, 3723, 1889, 1755, 1397, 3640, 1005, 1995, 1165, 1659, 3834, 1400, 3890, 1971, 1164, 1104, 1234, 1881, 3991, 1244, 1647, 3871, 1222, 3964, 3506, 3774, 3916, 1121, 1498, 1878, 1547, 1025, 1640, 1424, 3742, 1863, 3724, 3836, 3843, 1540, 3632, 1405, 3670, 3980, 1343, 1059, 1671, 3873, 1049, 1735, 3767, 1630, 3779, 1821, 1464, 3588, 3978, 1786, 1455, 1851, 1979, 3545, 1856, 3760, 3717, 1832, 3639, 3535, 1805, 3528, 3872, 1892, 3781, 3730, 3543, 1365, 1053, 3796, 3899, 1788, 3752, 1920, 3866, 1122, 3936, 1945, 1586, 1984, 1898, 3701, 3604, 3838, 3921, 3583, 1512, 1107, 1247, 1787, 1033, 1855, 1651, 3638, 1087, 1882, 1188, 3941, 1077, 3513, 3999, 1692, 3820, 1619, 1359, 1742, 1253, 1609, 1593, 1590, 3970, 3516, 1131, 1569, 1501, 1017, 1772, 1874, 3912, 1654, 3864, 1223, 1241, 3661, 3668, 1150, 1768, 3961, 1250, 1632, 1280, 1884, 1712, 1575, 1105, 3521, 1433, 1202, 1481, 3656, 3689, 1938, 1954, 3860, 1213, 1641, 3902, 1557, 1158, 3629, 1796, 3552, 1320, 3945, 3575, 1688, 1803, 3958, 1790, 3810, 1154, 1292, 3603, 3785, 3967, 1992, 1800, 3694, 1382, 1308, 3949, 3994, 1275, 1674, 3869, 1940, 1369, 3673, 1691, 1518, 3889, 1071, 1897, 1985, 1603, 1967, 1024, 3887, 1867, 3739, 1446, 1841, 1427, 3768, 3577, 1317, 3808, 3619, 1974, 1798, 3891, 1815, 3546, 1959, 1746, 1443, 3822, 1915, 3733, 1128, 1454, 3947, 1329, 1227, 3879, 1554, 3679, 1828, 1376, 3795, 3520, 3953, 3612, 1717, 3848, 1563, 1536, 1447, 3800, 1868, 1160, 1989, 3771, 3898, 1990, 1556, 1683, 1880, 3893, 1273, 1605, 1749, 1492, 1936, 1738, 1573, 1606, 1219, 1709, 1589, 1713, 1339, 3636, 1239, 1677, 1459, 1731, 1968, 1977, 1118, 3924, 1020, 1910, 1390, 1100, 3904, 3758, 1106, 1419, 1307, 3682, 1161, 1385, 3849, 3976, 1820, 3688, 3746, 1126, 1530, 3605, 3551, 1500, 1444, 1689, 3559, 1195, 1229, 3663, 1522, 1432, 1453, 1099, 3789, 3542, 1602, 3975, 3841, 1655, 1083, 3512, 1947, 1592, 3878, 1733, 1835, 3564, 3685, 3997, 1489, 3529, 1136, 1062, 3907, 1289, 1069, 3615, 1611, 1904, 1741, 3503, 1773, 1929, 3586, 1978, 1412, 1699, 3865, 3714, 1658, 3554, 1098, 1517, 1760, 1096, 1949, 1708, 3950, 1491, 3951, 1277, 1858, 1509, 1456, 1030, 1346, 1179, 1572, 1817, 3993, 1780, 1210, 3674, 1648, 1801, 1964, 3782, 3814, 1392, 3804, 1809, 3944, 3581, 3749, 3776, 3598, 3671, 3935, 1826, 3755, 3989, 1955, 1410, 3939, 3853, 1439, 1628, 1288, 3974, 3738, 1149, 3645, 1200, 1181, 3523, 3511, 3718, 1778, 1295, 3606, 3578, 1174, 1789, 3946, 3687, 1601, 1430, 1718, 1902, 3919, 1487, 1261, 3500, 1888, 3816, 1399, 3706, 3548, 3676, 1533, 1869, 3900, 1368, 3643, 1032, 1318, 1579, 3592, 1226, 3805, 3502, 1143, 3960, 3857, 1065, 1451, 1774, 1375, 1549, 1561, 3698, 1480, 1687, 3515, 3966, 3734, 1608, 1217, 3683, 3555, 3591, 3847, 1711, 3514, 1827, 1402, 1115]), noise_level = 0.5, shake_noise_level = 0.1, shake = 1.0, rotation_factor = 0.075, rotation_noise_level = 0.0075, α = 0.2, ρ̃₀ = 1.0, σ̃₀ = 1.0, δ = 0.9, σ₀ = 1.0, τ₀ = 0.9, λ = 1, radondims = [128, 64], sz = (256, 256), scale = 1, c = 1.0, sino_sparsity = 0.5, L = 300.0, L_experiment = false, save_prefix = "img/brainphantom256x256_pdps_known_dualscaling_20_10000_9000", dynrange = 1.0, Λ = 1.0) -iter time function_value psnr ssim -1 7.458991603925824e-6 563529.9047468598 -0.02636674332749045 0.0041734697881560985 -2 0.029424667009152472 478815.8877180156 0.4728121726799157 0.007393873939648631 -3 0.04358695898554288 446621.01861306693 0.8541320904970705 0.010072996745564798 -4 0.07331191797857173 400246.5976206233 1.1818011527566534 0.012073039664698679 -5 0.08814312596223317 362939.42699044186 1.4925472412903251 0.014051465408188281 -6 0.1169954169599805 323561.56825124007 1.7761864497222082 0.015115982742480386 -7 0.13146695794421248 284334.7445670818 2.083984724712945 0.015995479575748504 -8 0.1617940409341827 246317.43252711956 2.3984532997457633 0.01775300802462601 -9 0.17621249894727953 205541.18394848844 2.7254194078883742 0.01987000906395938 -10 0.2056500819453504 166760.72379936458 3.047659401876086 0.021490143444487752 -11 0.22073262394405901 129647.68974389756 3.38483165103416 0.02336657830950552 -12 0.25023349895491265 92628.68937525783 3.7181241614148997 0.026179082006162938 -13 0.2646568319760263 56425.64864111258 4.061515258347494 0.028993688389069808 -14 0.2947926649649162 20836.433660134353 4.40732646249738 0.03510883949506924 -15 0.3092432899575215 -10447.174209328285 4.747618196325444 0.03589697170627342 -16 0.338872456952231 -42515.32141731963 5.101149517948931 0.04386117251882236 -17 0.3533309979538899 -63119.92731434457 5.418751879823237 0.038486425229260726 -18 0.38261895696632564 -105541.5590402533 5.816154608941648 0.0697799959817118 -19 0.3972433319722768 -135978.0067528888 6.1832696742224735 0.08865849326642813 -20 0.42656941595487297 -165140.42845348633 6.551405598593208 0.11537748725480909 -30 0.6460802909568883 -350400.964625847 10.254553895713274 0.4480472948524359 -40 0.8662906249519438 -421329.19411009795 13.749689090507713 0.5361863559230877 -50 1.0854535839462187 -448461.76957347745 16.545919845403713 0.561111758201957 -60 1.3040628759481478 -457328.0313995502 18.289752616418966 0.5912700910568013 -70 1.5243164589628577 -459393.6253244579 19.189672613018946 0.6100606982283885 -80 1.7422732089762576 -460554.6697854429 19.664051519895654 0.6229758318346524 -90 1.9607078339613508 -460969.91655252344 19.958358234686596 0.6280277094389275 -100 2.1794170839712024 -460972.8429374658 20.16286433750038 0.631024707112307 -110 2.383851208956912 -460819.9756663424 20.28182889835419 0.6323596563905578 -120 2.6032532089739107 -461396.1571078476 20.372963420122762 0.6309302960109029 -130 2.8211496669682674 -461644.35279775906 20.30496029353259 0.6196621176027335 -140 3.039623666991247 -459669.82433511456 20.03605616783858 0.6001705381214962 -150 3.2599434169824235 -461377.35270378046 20.486447655698687 0.6274620176056287 -160 3.4781825419631787 -459580.0507116893 20.11379030880716 0.6066532546411252 -170 3.696168583963299 -461257.50699902576 20.62777850610058 0.6346811376387271 -180 3.9152012509584893 -460773.2431908192 20.501751581510604 0.6216592463738658 -190 4.1331271259405185 -460377.1688911911 20.50051644746936 0.6227511068280037 -200 4.351693667966174 -460318.23402201524 20.395684534611476 0.6153784219961699 -300 6.529022000962868 -461474.8130140242 20.7352801847634 0.6335093814404125 -400 8.694981916982215 -461069.68202049204 20.790708368997876 0.6344496899833714 -500 10.850913000962464 -460891.5380552021 20.767851301424326 0.6353979422863717 -600 13.026992833969416 -461708.6926510307 20.87160501616868 0.6442133578777001 -700 15.19843595896964 -461168.18074641726 20.871303131511645 0.6425844047058101 -800 17.365949917963007 -461247.0025275965 20.82680637872876 0.642275216497501 -900 19.540130876965122 -461352.5949219489 20.91470643832399 0.6451834225975633 -1000 21.701169793959707 -461368.207073563 20.705871796946923 0.6324081353467753 -1100 23.877596876962343 -461482.94301367877 21.59782831134988 0.6804547983325663 -1200 26.046273167972686 -461196.07190656563 21.89529875580126 0.695351975793228 -1300 28.20146191798267 -461581.75374825916 22.05797647498557 0.7022800214554449 -1400 30.343181126983836 -461415.7796921233 22.15131349974166 0.7075633892084601 -1500 32.502722001983784 -461419.7571235548 22.20611169905932 0.7100437766381722 -1600 34.653629543987336 -461451.41976794874 22.243285163550297 0.7129767882990969 -1700 36.76414687698707 -461358.69571003475 22.269733554172724 0.7145654958270448 -1800 38.89133045999915 -461468.02711054875 22.29053732340955 0.7153006248654192 -1900 41.00445283501176 -461471.51086336945 22.302629393744308 0.7143084066198495 -2000 43.117388710030355 -461551.3816673975 22.314105859243234 0.7160416940397969 -2100 45.25618321003276 -460776.03475141205 20.954106532528847 0.6378267364151339 -2200 47.389108710049186 -461616.8578864591 20.698179671294298 0.6232929064419316 -2300 49.50521750203916 -460849.94480730133 20.783907101253245 0.6362590030038515 -2400 51.59733458503615 -461340.719244005 20.84683821598235 0.6418304291650572 -2500 53.68617733503925 -460702.2744652891 20.64804326972326 0.6270398956831558 -2600 55.80447171005653 -461543.3489333875 20.87364614197037 0.6448347575445377 -2700 57.89097816805588 -460938.2183221819 20.823092840713187 0.6404984989843274 -2800 59.991070001065964 -461601.446279374 20.876892628119208 0.6429010223436173 -2900 62.09843408406596 -462063.3652945971 20.850453553284535 0.6418678044493948 -3000 64.20349112607073 -461933.3975917351 20.847361158103244 0.6395210676695535 -3100 66.31828466808656 -459659.56728269195 20.458863264811484 0.6115348679572872 -3200 68.4254074601049 -459489.04598748975 19.914914736645635 0.5817828364836369 -3300 70.54199137710384 -460772.89982956165 20.498095371031873 0.612402626074702 -3400 72.69053458608687 -460802.34275263146 20.584735637276054 0.6176447599026642 -3500 74.80276500308537 -460861.31976160425 20.705016237476364 0.6314018325413513 -3600 76.94524450309109 -461932.851046751 21.568440605773766 0.6791966812035373 -3700 79.10879846208263 -461973.8230159011 21.884764527642258 0.6931355199714218 -3800 81.28451871208381 -461915.7062584949 22.056968039369096 0.7024661733798362 -3900 83.46225454605883 -462037.9705417774 22.156093295047306 0.7071359325770814 -4000 85.64281308805221 -461832.4826507634 22.218266177302247 0.710610200112095