# 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_greedy_10_10000_9000", dynrange = 1.0, Λ = 1.0)
iter	time	function_value	psnr	ssim
1	1.0124989785254002e-5	559998.953614435	-0.02636674332749045	0.0041734697881560985
2	0.028456333006033674	476714.9039130535	0.47419961112481246	0.0074293444055148
3	0.04253620802774094	445128.80524157255	0.8555773119792093	0.010117635374366305
4	0.07085545803420246	399119.0265387717	1.183185286743629	0.012138577095164683
5	0.08493841704330407	361976.3871954779	1.4939285503174462	0.014134659999384753
6	0.11347358403145336	322627.32993720873	1.7776067933002384	0.015204251033975652
7	0.12742704202537425	283561.2095291894	2.0851947788971428	0.016076129699492257
8	0.1559694170136936	245640.74287148128	2.399724183184713	0.017848973471370128
9	0.1700773340126034	204944.31492415656	2.7268266076360193	0.019981535056110368
10	0.19843304302776232	166204.82329292136	3.049091390483105	0.021600439722430954
11	0.2125674590060953	129131.66918678714	3.3864384263270493	0.0234670609958266
12	0.2419587509939447	92153.35125215643	3.7197457324172896	0.026296682295172413
13	0.25608212599763647	55980.350223110145	4.0632327083948985	0.029177251014784966
14	0.2846047929779161	20414.43899301579	4.40912246327517	0.03529293231154393
15	0.2992723769857548	-10865.320420659278	4.749525580097424	0.03610480637480214
16	0.32790187699720263	-42914.74753744428	5.103208252509242	0.04408064658947526
17	0.3421555859968066	-63503.642030642855	5.421039385187347	0.03873905932023678
18	0.37069950299337506	-105925.71251661291	5.818535722270607	0.07028321196124329
19	0.38496950300759636	-136340.82754410175	6.185766472388677	0.08943234056679061
20	0.41369196100276895	-165492.98913340628	6.553986795865843	0.1160945157188977
30	0.6284502939961385	-350607.9566811528	10.258381672396597	0.44975587467296746
40	0.8430542529968079	-421455.1094278642	13.755871058279377	0.5365026429189454
50	1.0582801699929405	-448566.81943580584	16.5540860792106	0.5613623978464408
60	1.2746377119910903	-457413.8366350067	18.296820994561813	0.5915810156920336
70	1.489753378991736	-459467.6153464054	19.19456492428044	0.6103024740124721
80	1.704980920010712	-460620.2627929652	19.666480418461557	0.6230555522293311
90	1.9215493790106848	-461031.50964780897	19.95820442999396	0.627933744856654
100	2.138093004003167	-461031.8424340099	20.160484720933987	0.630892361685286
110	2.339766087010503	-460875.65069425764	20.275572380483492	0.6319260627216279
120	2.556092962011462	-461450.8268232071	20.366024227517787	0.6304421622427594
130	2.771441004035296	-461697.66234173946	20.298281173257813	0.6192288326068742
140	2.986109796038363	-459724.60763528803	20.028281527576294	0.5995325019845947
150	3.2021610460360534	-461427.4745267136	20.47421246279386	0.6269108679929178
160	3.4172441710252315	-459633.35672340955	20.107750215728586	0.6063945308059313
170	3.631832671002485	-461306.29593016626	20.610385769975736	0.6340282713095655
180	3.8467695870203897	-460824.0001435241	20.483574954373175	0.6210254904913959
190	4.062380337010836	-460426.9314087984	20.481710757588623	0.6219418600523348
200	4.277491129003465	-460371.4625908314	20.382295007031885	0.6151340753375436
300	6.4155170879967045	-461521.24392232316	20.713289372359203	0.6329277105669839
400	8.560102587973233	-461117.1499236976	20.76212788577028	0.6336981411711837
500	10.70064979599556	-460937.6660544096	20.74026569771879	0.6346327592962684
600	12.836019711976405	-461757.09033250285	20.845194861969816	0.6436000929456203
700	14.979126211983385	-461218.25804680365	20.848031675385712	0.642081326126763
800	17.11680712897214	-461296.7975302419	20.80315125490138	0.6415732703409688
900	19.25951458697091	-461401.82398155244	20.890627127268033	0.644588909883586
1000	21.408868461992824	-461418.13290174503	20.68650544353488	0.6318155609865101
1100	23.555133379006293	-461588.4116930191	21.742185293210547	0.6856197256492125
1200	25.70246262900764	-461309.7702063667	22.145585890930807	0.704492507880724
1300	27.84878812902025	-461696.93431776424	22.38419843913778	0.7143716356048007
1400	29.991105170018272	-461531.7217991018	22.541121783262508	0.7224906802918574
1500	32.144328004011186	-461538.0614080549	22.645989855367677	0.7269938156145002
1600	34.29991604600218	-461569.84670783684	22.72405044198458	0.7318389694902085
1700	36.43773196300026	-461477.762231504	22.78529593879343	0.7349779065459532
1800	38.57670933799818	-461588.98926797777	22.83279043926155	0.7369783010670452
1900	40.718937713012565	-461592.79268081445	22.86683515558497	0.736993493713237
2000	42.850999380025314	-461671.28307708265	22.897934782273587	0.7398119530594036
2100	44.98484883803758	-460836.3584065706	20.972127115484568	0.6385517321211391
2200	47.12289483801578	-461664.00439356174	20.683488347961372	0.6228742493116741
2300	49.254885588015895	-460897.0670831059	20.760601872557586	0.6356200812247147
2400	51.39105317101348	-461389.90432435717	20.823722311351954	0.6412977101410646
2500	53.526180921006016	-460753.2818393501	20.6272269319427	0.6266599890458691
2600	55.68305346299894	-461592.08589081984	20.845801809382	0.6441385829280855
2700	57.81881158798933	-460986.5921469672	20.799360315116036	0.6400885055403782
2800	59.9526721299917	-461650.8325669377	20.851064564061623	0.64221955528771
2900	62.089295296987984	-462112.63277143234	20.823177067921073	0.6411281683587434
3000	64.2309704220097	-461981.88237973215	20.820188541248243	0.638711407128069
3100	66.36894879699685	-459710.28660409804	20.434908976795416	0.6106464752175316
3200	68.51102775597246	-459542.6214941145	19.90208433268498	0.5816938138897785
3300	70.64657067297958	-460817.2954930551	20.47385543180102	0.6116535543089809
3400	72.79222050696262	-460844.37911794917	20.55675801727522	0.6166955307407066
3500	74.930579672975	-460908.4546394005	20.68153579855936	0.6309499829941464
3600	77.06690683896886	-462038.97884171107	21.70659560780871	0.6838584074844449
3700	79.20272575496347	-462087.5149081729	22.122292083500987	0.7013168630469673
3800	81.3425307549769	-462030.1803452781	22.36618805839025	0.7134581551546889
3900	83.47710637998534	-462154.90952704934	22.525354656631357	0.720654286287044
4000	85.61558771398268	-461949.3980230221	22.63502280567526	0.7261597895852216
