# 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_zerodual_20_10000_9000", dynrange = 1.0, Λ = 1.0)
iter	time	function_value	psnr	ssim
1	8.749979315325618e-6	563529.9047468598	-0.02636674332749045	0.0041734697881560985
2	0.027876249951077625	478856.08122394385	0.4727406757380256	0.007392494550781689
3	0.0412792079441715	446629.7517429079	0.8539140607374829	0.010069665864595012
4	0.06899879095726646	400269.21028804214	1.1814775051209048	0.01206703098275115
5	0.08268387493444607	362962.4296867711	1.4922690080185277	0.014042451163532112
6	0.11036995792528614	323591.6316846677	1.7754944575834577	0.01508482569419906
7	0.12366241592098959	284365.7473419522	2.0830435884394536	0.015941529087566995
8	0.15149541592109017	246354.00136495073	2.397592669599506	0.01769874438581332
9	0.1650198329007253	205571.0802674675	2.7246366424854465	0.019819222915664732
10	0.19340866690617986	166785.36886898818	3.046645922327711	0.021419885750702854
11	0.20701108290813863	129674.12709935356	3.383895663760555	0.023311111836831093
12	0.23496912492555566	92653.28690360574	3.7170275054469064	0.026093821873056477
13	0.24838299994007684	56452.42106017006	4.060330686248957	0.02889385552418672
14	0.27641983394278213	20857.101703946475	4.406071274468021	0.034948265482444334
15	0.28997666694340296	-10438.770138025126	4.7464343299667675	0.03570432755419898
16	0.31794804194942117	-42505.71385518951	5.09999542278511	0.04368141902332796
17	0.33136675096466206	-63126.3486904598	5.417650255125764	0.03836780958154917
18	0.35913070995593444	-105540.98876235927	5.815058755545951	0.06961002429483774
19	0.37249391796649434	-135974.1889506175	6.182079030583943	0.08845381424367746
20	0.4002704179729335	-165136.5459505766	6.549971469614685	0.11495759048821691
30	0.6057929589878768	-350415.0314217395	10.252427224594552	0.4471126464336167
40	0.8108901259838603	-421351.0510403995	13.745084438160491	0.5359227561921107
50	1.0156440009886865	-448479.0073656686	16.53690856952951	0.5606411527744863
60	1.2197823339665774	-457336.8096845315	18.27931489930829	0.5904973995238473
70	1.4234844589664135	-459398.0840489778	19.179624185637728	0.6093611914294952
80	1.627001124958042	-460559.5591591289	19.6561073243845	0.6226485701401986
90	1.8303712499618996	-460974.94356359495	19.952200148056296	0.6280188369990043
100	2.033651915960945	-460976.9330205219	20.158006520071332	0.6310120146415426
110	2.223528998962138	-460823.35601861565	20.28045260848213	0.6323698991507318
120	2.4310180399625096	-461399.4134079207	20.371377467926813	0.6309711207896302
130	2.6399422899703495	-461647.8006028427	20.30448183839954	0.6193954448838336
140	2.8482378729677293	-459669.0638747321	20.033771146525563	0.5999480920533595
150	3.056630247971043	-461382.3516274438	20.488673631379672	0.6275619497138916
160	3.261942831974011	-459584.76990198955	20.117378070081273	0.6064410771422082
170	3.4653351229499094	-461262.8949247521	20.632698623775173	0.6349720698087877
180	3.6690817899652757	-460776.93593953113	20.508040251631897	0.6220042768785026
190	3.873239622975234	-460380.741208267	20.506010644425544	0.6229426538602285
200	4.0777439149678685	-460321.5768809585	20.405097040206886	0.6158879739860437
300	6.1171862069750205	-461481.0754801085	20.74744573932619	0.6339098676634498
400	8.15905512298923	-461076.03712357854	20.80750624836734	0.6350963001174728
500	10.19770728997537	-460897.45434653264	20.785116329988476	0.6361090967152319
600	12.221155372972135	-461715.21774791036	20.8885812245061	0.6447264460010231
700	14.252476414956618	-461174.46579402883	20.889219459912386	0.6433037109515745
800	16.28088845696766	-461252.9004520723	20.843357967092256	0.6427375585920718
900	18.369778872962343	-461358.4785052773	20.93078477800351	0.6456148736715817
1000	20.434952913958114	-461372.8941923818	20.718958180955326	0.6327190681059426
1100	22.510358997969888	-461414.9701224553	21.766102959638015	0.673821653762958
1200	24.568871122959536	-461041.9166435804	22.003929647010448	0.6807280440248624
1300	26.620753414958017	-461372.63751805713	22.090299527481303	0.6834402398892566
1400	28.689906873943983	-461170.01896055706	22.128976681629965	0.68542545689322
1500	30.758403623942286	-461160.2126787338	22.14814323034254	0.686039807025556
1600	32.84527916595107	-461181.97184875305	22.158852066450976	0.6873869724378071
1700	34.89148724995903	-461081.1592314979	22.166042041787154	0.68789123111372
1800	36.95888554197154	-461185.26824670244	22.17030504773458	0.6875923962828542
1900	39.00933399997302	-461184.5751019353	22.173225250094085	0.6859667063517108
2000	41.08802108396776	-461258.5980084426	22.17592448055875	0.6878784679003866
2100	43.132536291959696	-460771.53124172956	20.927094120595704	0.6373595114078217
2200	45.182526332937414	-461620.79945889197	20.701856244574948	0.6236782213101687
2300	47.241347124945605	-460856.5242060274	20.79806687347351	0.6369615329273195
2400	49.273692499962635	-461346.24083914765	20.86233673773571	0.6422327739324779
2500	51.30414429094526	-460707.61172871926	20.662297654560877	0.6274252244047018
2600	53.32980708294781	-461549.6827994127	20.89159447838363	0.64551323986934
2700	55.377858207939425	-460944.5930327483	20.840278730844854	0.6412001573412833
2800	57.3990448739205	-461607.95922741835	20.89508058479119	0.6435383338540954
2900	59.44118787391926	-462069.6194101821	20.868371095578603	0.6423312207735403
3000	61.509484915906796	-461939.7728975037	20.863462578975806	0.6400346980387652
3100	63.57368387392489	-459662.8110055467	20.473002640122978	0.6119351687778724
3200	65.65636933292262	-459492.1410088092	19.92206478983645	0.5818023435679106
3300	67.69076954090269	-460777.87024429097	20.512673853422765	0.6128868277123647
3400	69.73519895688514	-460807.84031042387	20.59940201450136	0.6182434157116807
3500	71.76092108190642	-460866.9791645211	20.721798843383645	0.6319543021296905
3600	73.82680079090642	-461865.55813633965	21.721787076828377	0.6718406246758974
3700	75.96071062490228	-461811.69497260876	21.96519136008503	0.6772952205049021
3800	78.05097733391449	-461697.5668382382	22.052473958616435	0.6810839610534091
3900	80.13434595891158	-461790.21643061546	22.09139141098718	0.6825147958910878
4000	82.19546791689936	-461565.1222510449	22.11075340450595	0.6839747971613462
