# 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, save_prefix = "img/brainphantom256x256_pdps_known_zerodual_25_10000_9000", dynrange = 1.0, Λ = 1.0)
iter	time	function_value	psnr	ssim
1	9.542010957375169e-6	565295.3803130721	-0.02636674332749045	0.0041734697881560985
2	0.028729458019370213	479931.97125734057	0.4727406757380256	0.007392494550781689
3	0.04260549999889918	447342.41558703786	0.8539140607374829	0.010069665864595012
4	0.07128829200519249	400816.24284094496	1.1814775051209048	0.01206703098275115
5	0.08510683398344554	363400.2826415747	1.4922690080185277	0.014042451163532112
6	0.11325941697577946	324015.2879919201	1.7754944575834577	0.01508482569419906
7	0.12700449998374097	284728.41281393566	2.0830435884394536	0.015941529087566995
8	0.1556558329903055	246671.37771481354	2.397592669599506	0.01769874438581332
9	0.16950949898455292	205852.8536592255	2.7246366424854465	0.019819222915664732
10	0.19899603998055682	167052.39398645406	3.046645922327711	0.021419885750702854
11	0.2128069989848882	129916.2943854115	3.383895663760555	0.023311111836831093
12	0.24161258197273128	92879.81278921828	3.7170275054469064	0.026093821873056477
13	0.255499915976543	56665.16300268605	4.060330686248957	0.02889385552418672
14	0.2841892069845926	21058.857743178218	4.406071274468021	0.034948265482444334
15	0.29799845698289573	-10248.710840869799	4.7464343299667675	0.03570432755419898
16	0.32684108198736794	-42326.13060090378	5.09999542278511	0.04368141902332796
17	0.34065179096069187	-62953.14958522582	5.417650255125764	0.03836780958154917
18	0.3692843329627067	-105381.42446500185	5.815058755545951	0.06961002429483774
19	0.38316295796539634	-135821.4159877795	6.182079030583943	0.08845381424367746
20	0.4118146659457125	-164989.10123231474	6.549971469614685	0.11495759048821691
30	0.6243634569400456	-350329.2756322142	10.252427224594552	0.4471126464336167
40	0.8359275819384493	-421287.9999416141	13.745084438160491	0.5359227561921107
50	1.047335040930193	-448425.878613977	16.53690856952951	0.5606411527744863
60	1.2586926239309832	-457293.0598113892	18.27931489930829	0.5904973995238473
70	1.4692769989487715	-459359.5191783487	19.179624185637728	0.6093611914294952
80	1.6797901649551932	-460524.78842545766	19.6561073243845	0.6226485701401986
90	1.8904208319727331	-460941.6996627215	19.952200148056296	0.6280188369990043
100	2.1015561649692245	-460944.9885588711	20.158006520071332	0.6310120146415426
110	2.298138747981284	-460792.69892125804	20.28045260848213	0.6323698991507318
120	2.508908497984521	-461369.4368202822	20.371377467926813	0.6309711207896302
130	2.7218427479674574	-461618.380353424	20.30448183839954	0.6193954448838336
140	2.932817955967039	-459639.46096367826	20.033771146525563	0.5999480920533595
150	3.143063830968458	-461353.7696261568	20.488673631379672	0.6275619497138916
160	3.353338039974915	-459555.9381346951	20.117378070081273	0.6064410771422082
170	3.5637405819725245	-461234.35802084935	20.632698623775173	0.6349720698087877
180	3.775089164962992	-460748.0182927501	20.508040251631897	0.6220042768785026
190	3.9865128309756983	-460352.0384182738	20.506010644425544	0.6229426538602285
200	4.198075580963632	-460292.7151331865	20.405097040206886	0.6158879739860437
300	6.2715210389578715	-461452.82747868955	20.74744573932619	0.6339098676634498
400	8.346876788971713	-461046.54919730016	20.80750624836734	0.6350963001174728
500	10.408307038975181	-460868.65850542416	20.785116329988476	0.6361090967152319
600	12.491840038972441	-461685.74166015215	20.8885812245061	0.6447264460010231
700	14.57032862197957	-461144.6162491798	20.889219459912386	0.6433037109515745
800	16.655404454970267	-461223.13732504484	20.843357967092256	0.6427375585920718
900	18.746445038967067	-461328.76109906565	20.93078477800351	0.6456148736715817
1000	20.830292288970668	-461343.6144738312	20.718958180955326	0.6327190681059426
1100	22.93371766398195	-461329.9723444426	21.766102959638015	0.673821653762958
1200	25.048686413967516	-460931.14547009394	22.003929647010448	0.6807280440248624
1300	27.156958705978468	-461248.7530300743	22.090299527481303	0.6834402398892566
1400	29.27226199800498	-461038.61573070096	22.128976681629965	0.68542545689322
1500	31.38666170599754	-461024.4011086992	22.14814323034254	0.686039807025556
1600	33.50347795599373	-461043.73982142366	22.158852066450976	0.6873869724378071
1700	35.56376903899945	-460941.3538038143	22.166042041787154	0.68789123111372
1800	37.6392162059783	-461043.864176173	22.17030504773458	0.6875923962828542
1900	39.72243437299039	-461042.0543372998	22.173225250094085	0.6859667063517108
2000	41.776579664991004	-461115.6913654657	22.17592448055875	0.6878784679003866
2100	43.83591399798752	-460741.0142337783	20.927094120595704	0.6373595114078217
2200	45.860857205989305	-461593.00806690013	20.701856244574948	0.6236782213101687
2300	47.88142287300434	-460827.7110163703	20.79806687347351	0.6369615329273195
2400	49.91123612300726	-461316.7164819276	20.86233673773571	0.6422327739324779
2500	51.946085331001086	-460678.17355112656	20.662297654560877	0.6274252244047018
2600	53.979718831018545	-461519.95789703523	20.89159447838363	0.64551323986934
2700	56.0107506220229	-460915.6854972832	20.840278730844854	0.6412001573412833
2800	58.04343562203576	-461578.11697167065	20.89508058479119	0.6435383338540954
2900	60.09698612202192	-462039.6366551784	20.868371095578603	0.6423312207735403
3000	62.135679539031116	-461910.07843237615	20.863462578975806	0.6400346980387652
3100	64.20576449803775	-459633.40664824564	20.473002640122978	0.6119351687778724
3200	66.256937165017	-459463.02883933426	19.92206478983645	0.5818023435679106
3300	68.32724399800645	-460750.44349251135	20.512673853422765	0.6128868277123647
3400	70.35896403901279	-460780.951911309	20.59940201450136	0.6182434157116807
3500	72.3981922470266	-460838.69866366533	20.721798843383645	0.6319543021296905
3600	74.43291195502388	-461779.569859792	21.721787076828377	0.6718406246758974
3700	76.46570887204143	-461698.951156383	21.96519136008503	0.6772952205049021
3800	78.51025233004475	-461571.7922103962	22.052473958616435	0.6810839610534091
3900	80.57786433002912	-461656.8739047999	22.09139141098718	0.6825147958910878
4000	82.61989341402659	-461427.76087241067	22.11075340450595	0.6839747971613462
