--- a/img/shepplogan256x256_pdps_known_rotation_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/shepplogan256x256_pdps_known_rotation_20_10000_9000", dynrange = 1.0, Λ = 1.0) -iter time function_value psnr ssim -1 7.916009053587914e-6 555362.1431745938 -0.14938477322295188 0.005295187802576954 -2 0.028921707009430975 470965.17176531116 0.3314460538266129 0.0078112081340624105 -3 0.04338974901475012 439247.722726294 0.6990065571636511 0.01000280042917199 -4 0.07238212402444333 393414.46786418976 1.0234748563592682 0.012188894853834335 -5 0.08678862403030507 356906.48086817714 1.3253380988403776 0.014323816235725742 -6 0.11576983201666735 318080.9296750219 1.6079186818975129 0.015113447028961991 -7 0.13010712401592173 279921.6744559892 1.914997819459777 0.01712460792522995 -8 0.1596317490329966 242682.6308663172 2.225452962423629 0.0195285316279422 -9 0.17388183204457164 202376.45878938324 2.557121331359826 0.022139025685291726 -10 0.2027538740367163 164385.4376645373 2.880424218627551 0.02413730817109596 -11 0.2170544990513008 128118.68814002613 3.2184288538537342 0.02727091978986875 -12 0.24574058305006474 91851.33001951037 3.554548498981464 0.03030492713417986 -13 0.2600235410500318 56681.16955108162 3.9027757056978 0.034236322892530785 -14 0.2888022070401348 21552.881138652316 4.255620029932976 0.03942275663241934 -15 0.3030879150319379 -6933.842547822498 4.58867260343696 0.04236172394214785 -16 0.3320818740176037 -38132.101143003 4.94972639684636 0.049097947019988195 -17 0.3463603740383405 -54950.97610416152 5.2315032090791656 0.05055641217965871 -18 0.37509304101695307 -100191.64330465607 5.692929102158221 0.0669595625808021 -19 0.38937816602992825 -130986.15671658724 6.081354183297224 0.0804122998964257 -20 0.41815712401876226 -160929.6735047109 6.476139558165176 0.09280670823502037 -30 0.6326543740287889 -372453.4770328692 10.610832643010042 0.3541549342305579 -40 0.8474082490429282 -458830.1150923402 14.33649197626317 0.5247926833124967 -50 1.0618021650589071 -484512.2216579852 16.395580273378712 0.62145547658029 -60 1.2760071230586618 -489110.8098082604 17.233446981962892 0.6860347846006175 -70 1.4902914150734432 -489493.5456482296 17.648150486461812 0.7170528027234633 -80 1.706353457091609 -489798.1483663213 17.966588151906507 0.7403406101920298 -90 1.9217430820863228 -490472.61962010484 18.17812666849408 0.7553146457796567 -100 2.1367785400652792 -490445.90931598045 18.321600948869254 0.7653702148591773 -110 2.337928248074604 -490472.45087029884 18.356923153723706 0.7729860962995853 -120 2.557809457066469 -490750.19311752776 18.442219497370218 0.7776725528732843 -130 2.777898291067686 -491015.2652741599 18.15125769794784 0.7726343378075466 -140 2.9987229580583517 -491360.9308131314 17.790461952142287 0.7567620589522677 -150 3.2194779580750037 -491085.10261178087 18.423602716673837 0.7822435558244306 -160 3.440022999071516 -490621.22843595606 17.804317761989388 0.7714531218223233 -170 3.660403333051363 -490751.650734164 18.664923909834084 0.7891179230393812 -180 3.881455208058469 -490415.39736765524 18.42854659897678 0.7884016700461268 -190 4.102477624051971 -491255.52630533517 18.443951163711972 0.7891489170348913 -200 4.3231108330364805 -490327.3844236469 18.199620536527735 0.7879765141075625 -300 6.468131458037533 -490828.402269911 18.710463277395135 0.7949112069511671 -400 8.601025000039954 -490657.30489844363 18.745774553782837 0.7985664743676746 -500 10.755992167047225 -490365.9108786195 18.698391786549934 0.7984690923849465 -600 12.953326209040824 -491266.5642689452 18.87257902914499 0.8065955316836935 -700 15.119053418049589 -490798.2045174672 18.835028711036845 0.8041953521037447 -800 17.29391883505741 -490597.3739339626 18.784284094216712 0.8042651162926697 -900 19.499163210071856 -490885.490345645 18.88587254249189 0.8065952223968069 -1000 21.71555000107037 -491411.34590388986 18.623623628095043 0.7980608430664407 -1100 23.93544612606638 -491886.0974568922 20.649698025621746 0.8633466679064813 -1200 26.143367084063357 -491705.7308923672 21.457087219661727 0.8870656136220734 -1300 28.350397917063674 -492045.91125816474 21.950986809617042 0.8993776600995834 -1400 30.56020525007625 -491975.9869600194 22.300140212272126 0.9066646474611746 -1500 32.7633470830624 -491947.8893534837 22.579178901155288 0.912625234738903 -1600 34.97749224904692 -492010.3194659754 22.800510793511574 0.917519336658416 -1700 37.17977579103899 -491914.4379371507 22.990062958829952 0.9208068443475634 -1800 39.384890083019855 -492098.5892206226 23.15146561439853 0.9235874718237778 -1900 41.5874458750186 -492069.016144178 23.282656876182482 0.9243799038243842 -2000 43.789207041001646 -492113.59831658355 23.40295639230332 0.9270847666335853 -2100 46.02147124899784 -490787.5160109572 19.002162087136607 0.8105708554674146 -2200 48.22677070699865 -492005.6208189675 18.618956551555733 0.7834133048603815 -2300 50.40908970698365 -489757.3591864787 18.712343842199218 0.7988575973764007 -2400 52.609418498992454 -490638.31456929236 18.795922590063544 0.8034116873510242 -2500 54.79632858198602 -491467.1563346714 18.53591716641413 0.7963006043835451 -2600 56.99749420696753 -491067.42941324477 18.890696220695837 0.8064537915675937 -2700 59.18418478994863 -490776.55936448934 18.797008824983024 0.802561456361496 -2800 61.35068958194461 -490885.96624323964 18.89782972661592 0.8055497272989023 -2900 63.49295416593668 -491338.28397288435 18.84987507380553 0.8066656855976988 -3000 65.61763249893556 -490832.6215952387 18.814497440948664 0.8029209311587783 -3100 67.7472604579234 -492148.6608641416 18.26435420859942 0.7858244386848604 -3200 69.8768968329241 -489847.1052282641 17.177693664862137 0.7552690595667536 -3300 72.00020158290863 -491452.44934257993 18.297178784746684 0.7804132116008675 -3400 74.13553483289434 -490537.3265450032 18.39632600749497 0.7794094026116936 -3500 76.27203679090599 -491441.11791133235 18.663784111139805 0.7972560279756472 -3600 78.41300745788612 -491618.7912805367 20.614136258939453 0.863039318462748 -3700 80.5416302078811 -491965.7037865886 21.432945088394728 0.8856858203180422 -3800 82.66252904088469 -491775.38126708806 21.930123786507828 0.8989322462459487 -3900 84.80186437387601 -491965.0177946527 22.28129559383111 0.9060268886615894 -4000 86.92848195787519 -491741.07056272787 22.54282081767927 0.913016873824646