# 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.2, ρ̃₀ = 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_20_10000_9000", dynrange = 1.0, Λ = 1.0)
iter	time	function_value	psnr	ssim
1	7.875001756474376e-6	563529.9047468598	-0.02636674332749045	0.0041734697881560985
2	0.029420332983136177	478856.52972927096	0.47553123008044235	0.00746334511562588
3	0.04398308298550546	446456.1698713447	0.857056633466054	0.010156750840240049
4	0.07343649998074397	400134.653898339	1.1847149335411644	0.01219941698346098
5	0.08804262499324977	362724.0506312396	1.4954225614096608	0.014213209434067682
6	0.11763291698298417	323332.6916715397	1.779512102114193	0.015304028879733264
7	0.13228591700317338	284181.4019033931	2.0871919419532015	0.016199236739971548
8	0.16205141702084802	246173.9081278084	2.40171053917672	0.017985536022351516
9	0.17692675100988708	205422.03273018854	2.7288557061029284	0.020129350372848416
10	0.20649912601220421	166664.75403478963	3.051355171452679	0.021764930355656343
11	0.2219384180207271	129539.66796503487	3.3887848073887357	0.023593823118362577
12	0.25146112602669746	92538.5331700636	3.722250067000072	0.02646400577020073
13	0.2662775420176331	56338.932816047	4.06589959445086	0.02936028326026014
14	0.2965838329982944	20759.33732465922	4.411926057404891	0.03558525927754312
15	0.3112339159997646	-10546.09594906843	4.7523383872255724	0.036456900025501726
16	0.34083058300893754	-42616.5584264603	5.106106116536711	0.04442797376539451
17	0.35552391601959243	-63196.77126752385	5.42406630191185	0.039050442114268016
18	0.38511141602066346	-105678.29158091792	5.821642200678499	0.07085768807185569
19	0.399873666028725	-136100.87231172985	6.189054472765175	0.09026602945884304
20	0.4295062910532579	-165265.27503337548	6.557571037039256	0.11705816487098963
30	0.6497541660501156	-350459.1052074419	10.263677809406587	0.45196697676297254
40	0.8701001660665497	-421315.5408109985	13.76538921776681	0.5370176120459556
50	1.0918750830751378	-448451.4607188481	16.56901140104352	0.5619567466136817
60	1.3132858750759624	-457322.9664667305	18.311816014986	0.592493907510022
70	1.5355329170706682	-459388.7482454775	19.20729796623867	0.6111183167142602
80	1.7574523760704324	-460547.68311163905	19.67556977695091	0.6234218249025315
90	1.980398251063889	-460960.15876565117	19.96335086049688	0.6277242665177947
100	2.203165418060962	-460963.9403229104	20.163027545012216	0.6307859072373466
110	2.412259001052007	-460810.23109799705	20.272011631690493	0.6314944537130358
120	2.6360802510462236	-461386.5596528724	20.360989160103475	0.629849438725833
130	2.8596857510565314	-461634.0852239002	20.29192612297582	0.6189027260464761
140	3.0841016670747194	-459662.7802044287	20.020881412018614	0.598840886718133
150	3.3061878330772743	-461362.8767448185	20.461695030523394	0.626168381075465
160	3.528311291069258	-459570.18614364887	20.09845296318888	0.6060206512056002
170	3.7517256660794374	-461241.0770906194	20.59330120900023	0.6331153084830328
180	3.974027791089611	-460759.62783713185	20.465318293775432	0.6201638144146355
190	4.195782208087621	-460362.6346943802	20.46227517603526	0.6210262938608411
200	4.419565624091774	-460309.87794869835	20.365305166177677	0.6144848719267565
300	6.615538249083329	-461454.3989659228	20.687707393973728	0.6320828003885747
400	8.806463541084668	-461046.55484121304	20.72694525754615	0.6322775534408969
500	10.983859374071471	-460868.81441417994	20.70432722185794	0.633231038645795
600	13.16800929105375	-461687.30916277634	20.811891248585237	0.6426846289138661
700	15.327008832042338	-461148.3944773968	20.81698088596638	0.6410939263414377
800	17.477941998047754	-461228.0413439345	20.772439806057847	0.640565206474623
900	19.670127832039725	-461332.18766837864	20.86046452718952	0.6437093670818264
1000	21.824930499016773	-461351.53729155974	20.661975088680848	0.631029733281532
1100	23.987918999016983	-461480.077315633	21.57732401449882	0.6798403933252852
1200	26.11867858201731	-461195.0240296152	21.884982133277294	0.6949719284391163
1300	28.245098123006755	-461581.28444969904	22.0524776356407	0.7020521459783599
1400	30.393228706001537	-461415.58696261013	22.148231226300194	0.7074331831992058
1500	32.52483849701821	-461419.6234351667	22.204303578443554	0.7099575128004725
1600	34.69412449700758	-461451.346180693	22.242067722554417	0.7129140957164517
1700	36.83553016302176	-461358.6422149206	22.268901632251747	0.7145213781908641
1800	38.968934955017176	-461467.9818824466	22.289942234161707	0.7152701435633481
1900	41.093564039008925	-461471.4883424488	22.302142165838745	0.7142827416404962
2000	43.214659871999174	-461551.3563582529	22.31369827695516	0.7160215729571676
2100	45.36305016401457	-460761.72144077887	20.90477433867687	0.6353347216936578
2200	47.51401012201677	-461595.25106679124	20.64596472327316	0.6213206910845973
2300	49.66083462201641	-460827.305332433	20.7250589895571	0.6341908250367772
2400	51.83246512201731	-461320.60668415134	20.791563352198978	0.6404086661015876
2500	53.963754622032866	-460686.30632810015	20.598544019329616	0.6258294985922319
2600	56.120275414024945	-461521.15457940917	20.80914010809256	0.6429616483329281
2700	58.26160641401657	-460918.94595747034	20.76721003401343	0.6390688385773298
2800	60.404673705023015	-461580.1563828461	20.817198302530834	0.6411232627668964
2900	62.60801483003888	-462041.8211935099	20.788238455463887	0.6400787726798186
3000	64.81355003803037	-461910.79632258654	20.78596283042895	0.637461774879393
3100	67.02240216304199	-459644.14088962536	20.402604517199304	0.6094179230566128
3200	69.22561316302745	-459479.324780953	19.881444661107906	0.5812879030181715
3300	71.42698816303164	-460752.4928340355	20.443660873729158	0.6103814870619948
3400	73.63581128802616	-460778.8332091064	20.523322505607474	0.615209785979431
3500	75.83700362101081	-460842.0959077565	20.65024688714841	0.6300087017851805
3600	78.04689133001375	-461929.2216788807	21.544038462609155	0.6783383451343845
3700	80.25459953799145	-461972.56277914497	21.872431151513197	0.6925700871448494
3800	82.46007112099323	-461915.1818968644	22.05071845033968	0.7021442550532446
3900	84.66701970400754	-462037.6661677271	22.152686783603404	0.7069575464143085
4000	86.86832674601465	-461832.3576502744	22.216302159055125	0.7105067171079469
