# 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.1, ρ̃₀ = 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_10_10000_9000", dynrange = 1.0, Λ = 1.0)
iter	time	function_value	psnr	ssim
1	9.41698090173304e-6	559998.953614435	-0.02636674332749045	0.0041734697881560985
2	0.028081666998332366	476704.3011571503	0.4727406757380256	0.007392494550781689
3	0.04154404197470285	445204.4240546481	0.8539140607374829	0.010069665864595012
4	0.06950124999275431	399175.14518223645	1.1814775051209048	0.01206703098275115
5	0.08292037501814775	362086.723777164	1.4922690080185277	0.014042451163532112
6	0.11162025001249276	322744.319070163	1.7754944575834577	0.01508482569419906
7	0.12506225003744476	283640.4163979855	2.0830435884394536	0.015941529087566995
8	0.1532264590205159	245719.24866522517	2.397592669599506	0.01769874438581332
9	0.16680537501815706	205007.53348395147	2.7246366424854465	0.019819222915664732
10	0.19484854201436974	166251.31863405646	3.046645922327711	0.021419885750702854
11	0.20836995900026523	129189.79252723766	3.383895663760555	0.023311111836831093
12	0.236799834005069	92200.23513238067	3.7170275054469064	0.026093821873056477
13	0.2502646260254551	56026.93717513808	4.060330686248957	0.02889385552418672
14	0.27831625103135593	20453.58962548299	4.406071274468021	0.034948265482444334
15	0.2919452090281993	-10818.888732335778	4.7464343299667675	0.03570432755419898
16	0.32001745904562995	-42864.880363760974	5.09999542278511	0.04368141902332796
17	0.33341083404957317	-63472.74690092776	5.417650255125764	0.03836780958154917
18	0.36137504305224866	-105860.1173570741	5.815058755545951	0.06961002429483774
19	0.374777460063342	-136279.73487629349	6.182079030583943	0.08845381424367746
20	0.4026892520778347	-165431.4353871003	6.549971469614685	0.11495759048821691
30	0.6111022520926781	-350586.54300079006	10.252427224594552	0.4471126464336167
40	0.8264586270961445	-421477.1532379703	13.745084438160491	0.5359227561921107
50	1.0329677101108246	-448585.26486905187	16.53690856952951	0.5606411527744863
60	1.2465620441071223	-457424.309430816	18.27931489930829	0.5904973995238473
70	1.4610399191151373	-459475.21379023604	19.179624185637728	0.6093611914294952
80	1.6731328361202031	-460629.1006264714	19.6561073243845	0.6226485701401986
90	1.887541253119707	-461041.43136534194	19.952200148056296	0.6280188369990043
100	2.1189232531178277	-461040.8219438235	20.158006520071332	0.6310120146415426
110	2.3257682111288887	-460884.6702133308	20.28045260848213	0.6323698991507318
120	2.535259836120531	-461459.3665831976	20.371377467926813	0.6309711207896302
130	2.7472332111210562	-461706.64110168006	20.30448183839954	0.6193954448838336
140	2.959045545110712	-459728.26969683985	20.033771146525563	0.5999480920533595
150	3.1706759201188106	-461439.5156300179	20.488673631379672	0.6275619497138916
160	3.3780575871351175	-459642.4334365784	20.117378070081273	0.6064410771422082
170	3.5848772951285355	-461319.9687325577	20.632698623775173	0.6349720698087877
180	3.7907916291151196	-460834.7712330933	20.508040251631897	0.6220042768785026
190	3.99694358810666	-460438.14678825345	20.506010644425544	0.6229426538602285
200	4.2025808381149545	-460379.30037650256	20.405097040206886	0.6158879739860437
300	6.289980671106605	-461537.57148294634	20.74744573932619	0.6339098676634498
400	8.356046379107283	-461135.01297613536	20.80750624836734	0.6350963001174728
500	10.441382837103447	-460955.04602874955	20.785116329988476	0.6361090967152319
600	12.497597587120254	-461774.1699234268	20.8885812245061	0.6447264460010231
700	14.547355045098811	-461234.1648837269	20.889219459912386	0.6433037109515745
800	16.58932558709057	-461312.4267061273	20.843357967092256	0.6427375585920718
900	18.671138754085405	-461417.9133177006	20.93078477800351	0.6456148736715817
1000	20.72034208808327	-461431.45362948294	20.718958180955326	0.6327190681059426
1100	22.764710713090608	-461584.9656784808	21.766102959638015	0.673821653762958
1200	24.822117671079468	-461263.45899055345	22.003929647010448	0.6807280440248624
1300	26.87927008708357	-461620.4064940227	22.090299527481303	0.6834402398892566
1400	28.933095087093534	-461432.8254202693	22.128976681629965	0.68542545689322
1500	30.98696667011245	-461431.835818803	22.14814323034254	0.686039807025556
1600	33.065137462108396	-461458.4359034117	22.158852066450976	0.6873869724378071
1700	35.123556087113684	-461360.7700868651	22.166042041787154	0.68789123111372
1800	37.17862358709681	-461468.07638776145	22.17030504773458	0.6875923962828542
1900	39.23336417111568	-461469.6166312063	22.173225250094085	0.6859667063517108
2000	41.2981130880944	-461544.4112943964	22.17592448055875	0.6878784679003866
2100	43.33980583806988	-460832.56525763223	20.927094120595704	0.6373595114078217
2200	45.38649563008221	-461676.3822428756	20.701856244574948	0.6236782213101687
2300	47.44308554608142	-460914.1505853416	20.79806687347351	0.6369615329273195
2400	49.50086192108574	-461405.2895535876	20.86233673773571	0.6422327739324779
2500	51.571196713106474	-460766.4880839047	20.662297654560877	0.6274252244047018
2600	53.64862783809076	-461609.1326041676	20.89159447838363	0.64551323986934
2700	55.74452192109311	-461002.4081036787	20.840278730844854	0.6412001573412833
2800	57.83671450507245	-461667.6437389137	20.89508058479119	0.6435383338540954
2900	59.8927692550933	-462129.58492018963	20.868371095578603	0.6423312207735403
3000	61.96365771407727	-461999.1618277587	20.863462578975806	0.6400346980387652
3100	64.01866029808298	-459721.61972014885	20.473002640122978	0.6119351687778724
3200	66.07369946409017	-459550.36534775904	19.92206478983645	0.5818023435679106
3300	68.13306438107975	-460832.7237478503	20.512673853422765	0.6128868277123647
3400	70.18161571410019	-460861.61710865365	20.59940201450136	0.6182434157116807
3500	72.26864417208708	-460923.5401662326	20.721798843383645	0.6319543021296905
3600	74.36199792209663	-462037.53468943486	21.721787076828377	0.6718406246758974
3700	76.42633492211462	-462037.1826050603	21.96519136008503	0.6772952205049021
3800	78.48424167212215	-461949.1160939223	22.052473958616435	0.6810839610534091
3900	80.57924642210128	-462056.9014822466	22.09139141098718	0.6825147958910878
4000	82.63573221309343	-461839.8450083133	22.11075340450595	0.6839747971613462
