# 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.25, ρ̃₀ = 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, phantom_ρ = 100, save_prefix = "img/brainphantom256x256_pdps_known_proximal_25_10000_9000", dynrange = 1.0, Λ = 1.0)
iter	time	function_value	psnr	ssim
1	1.1208001524209976e-5	565295.3803130721	-0.02636674332749045	0.0041734697881560985
2	0.030958916002418846	479903.53604721784	0.47510578422212846	0.007430833192381343
3	0.04700304099242203	447179.6544496697	0.856780041675176	0.01015289692431718
4	0.07816062399069779	400678.2039598508	1.1847022392051572	0.012200709401226719
5	0.09410695798578672	363171.05617481883	1.4956881684350867	0.01422103454855555
6	0.12495074997423217	323745.6521852108	1.7800574992295406	0.015322665468890606
7	0.1403349999745842	284502.0333980168	2.0877391453491843	0.016229768192292267
8	0.17041012496338226	246452.36558857936	2.402426401602466	0.018030107393668022
9	0.18575683294329792	205673.27599040815	2.7296551174665145	0.02018915538354634
10	0.21652779093710706	166882.72325325175	3.0522796697062047	0.021835106188528505
11	0.23186587495729327	129757.9154328813	3.3898026259319227	0.023697056842494933
12	0.26188408394227736	92712.1952099948	3.723359963088117	0.02657024794773831
13	0.2772247929533478	56490.39312252364	4.067170659660131	0.02956731997949069
14	0.30726920993765816	20897.623079374418	4.413401240585629	0.03569257954416139
15	0.32279116893187165	-10415.341546443657	4.754083941460403	0.03660643575869345
16	0.353325251955539	-42511.70124290191	5.108130081393974	0.044615240959897745
17	0.3686891689430922	-63114.046320068366	5.426428898743453	0.03929203389102268
18	0.39873658496071585	-105624.81236603255	5.8240923296133325	0.0714202242999051
19	0.4140553349570837	-136085.00803839392	6.19172594361288	0.0910253181673298
20	0.44407387694809586	-165287.60969345222	6.560283452188966	0.11777774574797134
30	0.6720482939563226	-350499.8250435336	10.267948050733066	0.45421645159562163
40	0.8989842529699672	-421285.1588677082	13.772972541784924	0.5372425252799405
50	1.1295653779816348	-448415.533189213	16.58055274874455	0.56183291190998
60	1.3630736279883422	-457291.7259013016	18.321235703711658	0.5917968601281205
70	1.5945903779938817	-459361.59287753206	19.209667423636546	0.6095118918619052
80	1.8267059609934222	-460517.87694625935	19.665544794125193	0.6212708656680402
90	2.060221751977224	-460929.84522401367	19.941926301721526	0.62544810592906
100	2.293308376974892	-460927.4092791235	20.126251135021032	0.6274873757648477
110	2.510934959980659	-460771.70691945165	20.22744016236203	0.6284073125584029
120	2.7436250429600477	-461346.20160266	20.304244604540944	0.6267378958725822
130	2.9766162089654244	-461593.1420969592	20.246227462992195	0.6171235590841404
140	3.210202833986841	-459641.2621304863	19.989200846766334	0.598095202999539
150	3.4430841259891167	-461320.5286556743	20.39080936566059	0.6229226829905515
160	3.67378725099843	-459535.10978562984	20.060752238040905	0.6051829648810342
170	3.907811958983075	-461191.8001953786	20.496070294511107	0.6290585069479175
180	4.142536000988912	-460718.07765946904	20.375882188818192	0.6171498819523789
190	4.374882542993873	-460318.0033058518	20.37213834468943	0.6178932228264563
200	4.606875125988154	-460266.89521381905	20.27602904533054	0.6112377255195048
300	6.893242750986246	-461394.020528838	20.556573438130986	0.6267953947894078
400	9.197771250968799	-460990.8242873533	20.59166744328965	0.6272159234000577
500	11.473635833972367	-460812.77938580804	20.569222878331956	0.6281504723265039
600	13.804388791962992	-461626.9258708312	20.663768149340605	0.6370211369614724
700	16.1070880839834	-461085.51035547035	20.664338117622982	0.6346530519446999
800	18.37919237598544	-461166.7498599184	20.63163459979315	0.6352540046103498
900	20.678265958995325	-461270.9600552488	20.706579507588945	0.6375416513832664
1000	22.957112042000517	-461292.48613861925	20.530283814551517	0.6263470260036232
1100	25.246895876014605	-461425.7384422477	21.429317690186103	0.6745805844494788
1200	27.484348876023432	-461142.4996844493	21.727087676488637	0.6887513020672977
1300	29.67595179300406	-461529.6730147425	21.883869318815062	0.6950370179526842
1400	31.908784418017603	-461364.8612647675	21.969621017679458	0.6996601090400293
1500	34.124499001016375	-461368.4331619449	22.01760943112035	0.701767975512315
1600	36.33485954301432	-461399.93457576144	22.048466933265523	0.7041135222978161
1700	38.51358100102516	-461306.8988329608	22.07174432409296	0.7054619284085104
1800	40.69175675100996	-461415.9786913696	22.09176402326003	0.7062362082071996
1900	42.897233876021346	-461419.37612963066	22.100842298306244	0.7051402060767507
2000	45.079986293043476	-461500.050922566	22.109227342013465	0.7067727022034154
2100	47.26271108404035	-460711.1290473605	20.7722439695809	0.62915600666235
2200	49.46900125104003	-461539.8731731929	20.521494265865726	0.6162712183848869
2300	51.640878667036304	-460770.20039362134	20.593524547265048	0.6293046145362986
2400	53.82582779205404	-461258.7713221692	20.642274067437217	0.6344051645492341
2500	55.99229808404925	-460633.4954779533	20.48195393511243	0.6217774835279559
2600	58.15646787604783	-461462.54010336084	20.667080646688547	0.6375393988888594
2700	60.33018970905687	-460857.98851499916	20.624998247714377	0.6337219845376894
2800	62.50038266705815	-461520.20830195426	20.667564617534513	0.6353174192408367
2900	64.72236491707736	-461981.03334545734	20.645681798285885	0.6346383655078722
3000	66.88592625007732	-461852.6056311644	20.645873298080172	0.6318464801961852
3100	69.08854420908028	-459601.3740893914	20.305859589688016	0.6057187631120082
3200	71.29720779208583	-459444.2603333349	19.834509696456937	0.5805778759011913
3300	73.4757776671031	-460700.28659691743	20.32970332165072	0.6070298631321207
3400	75.67857858311618	-460724.2710578523	20.402470566314005	0.6113320096676764
3500	77.86582979213563	-460787.8917146806	20.51684554915498	0.6250769750078768
3600	80.08224875113228	-461875.3802847129	21.40308004915039	0.6735847201978864
3700	82.29431779313018	-461920.3732888907	21.721629643523855	0.6869987711659937
3800	84.48202604314429	-461864.6104053248	21.88785495621234	0.6957126012637542
3900	86.6943894601427	-461986.07644671825	21.978947514844723	0.6998395679111479
4000	88.85404791813926	-461781.12469112524	22.033869729054594	0.7028359751319844
