# 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.15, ρ̃₀ = 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_15_10000_9000", dynrange = 1.0, Λ = 1.0)
iter	time	function_value	psnr	ssim
1	8.417002391070127e-6	561764.4291806474	-0.02636674332749045	0.0041734697881560985
2	0.027820457995403558	477780.1911905471	0.4727406757380256	0.007392494550781689
3	0.041207000002032146	445917.087898778	0.8539140607374829	0.010069665864595012
4	0.0689165000221692	399722.1777351393	1.1814775051209048	0.01206703098275115
5	0.08229975003632717	362524.57673196757	1.4922690080185277	0.014042451163532112
6	0.11006766703212634	323167.97537741537	1.7754944575834577	0.01508482569419906
7	0.1234048750193324	284003.08186996885	2.0830435884394536	0.015941529087566995
8	0.15121916600037366	246036.62501508795	2.397592669599506	0.01769874438581332
9	0.1645486659836024	205289.30687570947	2.7246366424854465	0.019819222915664732
10	0.1933784159773495	166518.3437515223	3.046645922327711	0.021419885750702854
11	0.20668962397030555	129431.9598132956	3.383895663760555	0.023311111836831093
12	0.23449379098019563	92426.76101799321	3.7170275054469064	0.026093821873056477
13	0.24785520698060282	56239.67911765407	4.060330686248957	0.02889385552418672
14	0.27553083197562955	20655.345664714732	4.406071274468021	0.034948265482444334
15	0.2888694569701329	-10628.829435180452	4.7464343299667675	0.03570432755419898
16	0.3165709989552852	-42685.29710947524	5.09999542278511	0.04368141902332796
17	0.3298920829547569	-63299.54779569378	5.417650255125764	0.03836780958154917
18	0.3576078739424702	-105700.55305971668	5.815058755545951	0.06961002429483774
19	0.37095554094412364	-136126.96191345548	6.182079030583943	0.08845381424367746
20	0.39862887395429425	-165283.99066883844	6.549971469614685	0.11495759048821691
30	0.6042374569515232	-350500.7872112648	10.252427224594552	0.4471126464336167
40	0.8090007909631822	-421414.1021391849	13.745084438160491	0.5359227561921107
50	1.0137211659748573	-448532.13611736026	16.53690856952951	0.5606411527744863
60	1.218491624982562	-457380.55955767375	18.27931489930829	0.5904973995238473
70	1.4233064169820864	-459436.64891960693	19.179624185637728	0.6093611914294952
80	1.6267317089659628	-460594.3298928002	19.6561073243845	0.6226485701401986
90	1.8300031249818858	-461008.1874644685	19.952200148056296	0.6280188369990043
100	2.0339155419787858	-461008.8774821727	20.158006520071332	0.6310120146415426
110	2.227088333980646	-460854.0131159732	20.28045260848213	0.6323698991507318
120	2.434869250981137	-461429.3899955591	20.371377467926813	0.6309711207896302
130	2.6399893339839764	-461677.22085226135	20.30448183839954	0.6193954448838336
140	2.8437034169910476	-459698.666785786	20.033771146525563	0.5999480920533595
150	3.0473470419819932	-461410.9336287309	20.488673631379672	0.6275619497138916
160	3.250990374974208	-459613.60166928393	20.117378070081273	0.6064410771422082
170	3.454701791983098	-461291.43182865495	20.632698623775173	0.6349720698087877
180	3.657859374972759	-460805.8535863122	20.508040251631897	0.6220042768785026
190	3.8631147499836516	-460409.44399826025	20.506010644425544	0.6229426538602285
200	4.068308249989059	-460350.43862873054	20.405097040206886	0.6158879739860437
300	6.095595083985245	-461509.3234815274	20.74744573932619	0.6339098676634498
400	8.12688424999942	-461105.5250498569	20.80750624836734	0.6350963001174728
500	10.154118792008376	-460926.25018764107	20.785116329988476	0.6361090967152319
600	12.179935666994425	-461744.6938356686	20.8885812245061	0.6447264460010231
700	14.2103984580026	-461204.3153388779	20.889219459912386	0.6433037109515745
800	16.257800416002283	-461282.6635790998	20.843357967092256	0.6427375585920718
900	18.28467612399254	-461388.19591148896	20.93078477800351	0.6456148736715817
1000	20.399386873992626	-461402.1739109324	20.718958180955326	0.6327190681059426
1100	22.543398081994383	-461499.96790046804	21.766102959638015	0.673821653762958
1200	24.705900123983156	-461152.6878170669	22.003929647010448	0.6807280440248624
1300	26.865277248987695	-461496.5220060399	22.090299527481303	0.6834402398892566
1400	29.021003290981753	-461301.42219041317	22.128976681629965	0.68542545689322
1500	31.17654395799036	-461296.0242487684	22.14814323034254	0.686039807025556
1600	33.329985499964096	-461320.2038760824	22.158852066450976	0.6873869724378071
1700	35.47856583294924	-461220.9646591815	22.166042041787154	0.68789123111372
1800	37.62647145794472	-461326.67231723195	22.17030504773458	0.6875923962828542
1900	39.78949629093404	-461327.0958665708	22.173225250094085	0.6859667063517108
2000	41.94430170691339	-461401.5046514195	22.17592448055875	0.6878784679003866
2100	44.092154289915925	-460802.0482496809	20.927094120595704	0.6373595114078217
2200	46.23334453991265	-461648.59085088375	20.701856244574948	0.6236782213101687
2300	48.36546499893302	-460885.3373956845	20.79806687347351	0.6369615329273195
2400	50.50231587392045	-461375.76519636763	20.86233673773571	0.6422327739324779
2500	52.62634424891439	-460737.04990631196	20.662297654560877	0.6274252244047018
2600	54.75425203991472	-461579.40770179016	20.89159447838363	0.64551323986934
2700	56.86679974893923	-460973.5005682135	20.840278730844854	0.6412001573412833
2800	58.97865687395097	-461637.801483166	20.89508058479119	0.6435383338540954
2900	61.096633373963414	-462099.6021651859	20.868371095578603	0.6423312207735403
3000	63.20626612394699	-461969.46736263117	20.863462578975806	0.6400346980387652
3100	65.32206441593007	-459692.2153628478	20.473002640122978	0.6119351687778724
3200	67.4248589159397	-459521.25317828415	19.92206478983645	0.5818023435679106
3300	69.53400762495585	-460805.29699607065	20.512673853422765	0.6128868277123647
3400	71.63416283394326	-460834.72870953876	20.59940201450136	0.6182434157116807
3500	73.71944254194386	-460895.25966537686	20.721798843383645	0.6319543021296905
3600	75.81423754192656	-461951.5464128872	21.721787076828377	0.6718406246758974
3700	77.87708612592542	-461924.4387888345	21.96519136008503	0.6772952205049021
3800	79.95760129290284	-461823.3414660803	22.052473958616435	0.6810839610534091
3900	82.0695272088924	-461923.558956431	22.09139141098718	0.6825147958910878
4000	84.14442350086756	-461702.4836296791	22.11075340450595	0.6839747971613462
