# 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_noprediction_25_10000_9000", dynrange = 1.0, Λ = 1.0)
iter	time	function_value	psnr	ssim
1	9.166018571704626e-6	565295.3803130721	-0.02636674332749045	0.0041734697881560985
2	0.028202833025716245	523926.0364209317	0.21351851240197253	0.004019799987541865
3	0.04162058304063976	483245.0236139763	0.4693613135979432	0.0045938070393532
4	0.07047012401744723	441797.5439671186	0.7243328725883622	0.004274835807043116
5	0.08382316600182094	401081.61382842937	0.9638591090368822	0.0033436123404211654
6	0.11244324900326319	361593.49543775123	1.2223915488404213	0.0032934209308508955
7	0.12581237399717793	322491.80302005814	1.4889920700911525	0.003505720874476014
8	0.15483829099684954	283576.19702463347	1.7564654888021325	0.0024267177682847935
9	0.16847204099758528	244712.098438437	2.0392937187551534	0.00356432413297146
10	0.19765195800573565	206031.8814501682	2.32383032603295	0.0037781216628429057
11	0.21130716701736674	169042.26287949926	2.6122876015066017	0.0014693114336109263
12	0.24027645902242512	132424.0600731025	2.903546773061918	0.0018124603386700652
13	0.2538203760050237	96396.68710551923	3.205939337156443	0.0030261223430406033
14	0.2820467510027811	62051.5528874835	3.4950279130313904	0.0019632317355541603
15	0.2958419169881381	28215.355776214026	3.8323378836446804	0.007004601463467288
16	0.32424887499655597	-3572.419277458048	4.127650405968843	0.005448251702030085
17	0.3376936659915373	-33349.121196662876	4.367412082438091	0.0019742262157258625
18	0.3663964580046013	-64035.718353906064	4.755511332301376	0.009581874062421893
19	0.37983033299678937	-92848.01883326934	5.065189741103864	0.01034140253468932
20	0.4085233329969924	-120407.89234020142	5.378300419648627	0.012353340779787524
30	0.6187396659806836	-316938.1883672513	8.416570172403476	0.10822502834108497
40	0.8283137499820441	-407544.71817725233	11.260476634742249	0.30034296354221485
50	1.038004500005627	-441191.99001890334	13.089408153982472	0.40463055931947967
60	1.2465527090243995	-452709.0833145495	14.61256221381566	0.46136296345827116
70	1.4572892500436865	-454771.0230843682	14.967631200137738	0.4718147276934251
80	1.6666180830507074	-456085.99579523207	15.344423744224136	0.49444705319949045
90	1.8806929580459837	-458051.1503559795	16.23726307652731	0.5093003620129546
100	2.091072667040862	-458391.48151983984	16.692714376928908	0.526213998396347
110	2.2873658330354374	-456334.91419409733	16.17126850964455	0.5007382957963303
120	2.495690166048007	-457945.6759083076	16.624420836012202	0.5082653787691673
130	2.7065644570684526	-457745.92052232113	16.725142008875025	0.4982646552902247
140	2.9160279990755953	-456333.94200900494	16.91553390690286	0.493895339019698
150	3.123566707072314	-459184.4885057824	17.617281925616098	0.5309343584619322
160	3.332700123079121	-457897.3596464773	17.75571262508229	0.5172035526576583
170	3.540425164072076	-459344.19532250456	17.97560142792901	0.5375239122401061
180	3.748944747057976	-459292.8930503021	18.165609435715552	0.5393484599706178
190	3.9599639970401768	-458861.17754794285	18.33578809944114	0.5300022724484631
200	4.169113789015682	-459143.685658249	18.486398083350153	0.5409384710860784
300	6.229665455000941	-460978.92657916626	19.463567743898633	0.5818521578380125
400	8.292877454980044	-460267.6106957868	19.33949828910768	0.5799665107640387
500	10.349343579960987	-459100.1655744206	18.505684230063043	0.5617135543287441
600	12.410720954969293	-460746.18336562475	18.9019461392127	0.5803357640970563
700	14.466692246962339	-460059.76022756065	18.91767213362293	0.5678520538082183
800	16.52180520494585	-458874.9276458219	18.147799893228004	0.5553256923747589
900	18.57712691294728	-460231.911540191	18.96293317304504	0.5752616798589846
1000	20.63122049596859	-460725.3587049916	19.299442648157417	0.5803536446941447
1100	22.68277145398315	-461344.75967433315	20.918980166532847	0.6531692370708505
1200	24.712006786983693	-461115.98086155293	21.474060845908937	0.6772395075522313
1300	26.706109828985063	-461518.1343977251	21.741517049584424	0.6879840085612807
1400	28.71232095398591	-461357.6950729881	21.880740436974904	0.6948968553597735
1500	30.706591411988484	-461363.31297578855	21.958353698205173	0.6984143216169009
1600	32.70427820298937	-461396.4536541244	22.006242440450453	0.7017180589583859
1700	34.69767174497247	-461304.2502619315	22.039699676945062	0.7036513255579334
1800	36.67854778698529	-461413.2342511399	22.065871259407096	0.704593949894285
1900	38.65873020296567	-461417.15101515467	22.07997010971456	0.7037149459264208
2000	40.6440414939716	-461497.62494715495	22.091703835136695	0.7054644726065264
2100	42.62745266096317	-459199.8887442621	18.17933044315606	0.5531830186550972
2200	44.61833541095257	-461044.5990391365	19.274006400082712	0.5716812596599377
2300	46.60035936895292	-459913.7002263213	18.960543732693317	0.5768254186915276
2400	48.58351099395077	-460323.7277240528	18.903784284583857	0.5654574303711546
2500	50.56591732794186	-459671.11602936324	18.935186531076155	0.5648598911697367
2600	52.55634328594897	-459638.333330639	18.372005157503622	0.544210073715195
2700	54.54362632794073	-460239.5650984061	19.520816756905198	0.5881560151147665
2800	56.5254809949256	-460124.3223310338	18.74078751851268	0.5561916311000084
2900	58.50678353692638	-459466.5061112715	18.186799588864396	0.5469585147519052
3000	60.51412170394906	-459426.09490797337	17.983027766230343	0.5511990749529055
3100	62.51693328694091	-459132.54830025684	19.250727113524235	0.5684267196612132
3200	64.49977553694043	-459185.1319074949	19.05689781115761	0.5603675854290519
3300	66.50963503695675	-460035.2219269329	19.02234998655454	0.5630073291224362
3400	68.5725835369667	-457661.4066210271	17.778589082914902	0.5372576249785955
3500	70.61311757896328	-460032.9979074477	18.93476482963897	0.574888249118352
3600	72.66525216196897	-461776.6053701969	20.822824913343663	0.6487830120521957
3700	74.67997311998624	-461891.1723377908	21.448771823645437	0.6735708207322336
3800	76.67676437000046	-461851.0741979814	21.742302238937686	0.6877821940370992
3900	78.69431945300312	-461979.09645858954	21.89419925338237	0.694857326541838
4000	80.7099565359822	-461776.63198815036	21.981263202478658	0.699613140605402
