# 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, phantom_ρ = 100, save_prefix = "img/brainphantom256x256_pdps_known_proximal_10_10000_9000", dynrange = 1.0, Λ = 1.0)
iter	time	function_value	psnr	ssim
1	8.25000461190939e-6	559998.953614435	-0.02636674332749045	0.0041734697881560985
2	0.029035667015705258	476702.3999876519	0.4736890264940851	0.0074078534592577835
3	0.04385950000141747	445151.4481718742	0.8550663013597882	0.01010294007468283
4	0.07300341699738055	399133.27899208653	1.1827770219107538	0.012120605310117936
5	0.08786016699741594	362008.4855032445	1.4936496364619263	0.01411402225321938
6	0.11680654197698459	322654.3503474456	1.7773415133048265	0.015180179866191456
7	0.13160895896726288	283565.40985823533	2.0849445745208834	0.016056959011475136
8	0.16060195898171514	245646.0334588395	2.399551420737132	0.017831241019007488
9	0.17540287598967552	204949.45601712886	2.7266725854195832	0.01996743703556237
10	0.2043553759867791	166198.75281774256	3.048934201364329	0.021586229639314482
11	0.21914800099330023	129139.83227267055	3.38629931788309	0.023463913526352458
12	0.24898762599332258	92146.48981368991	3.7196063990592005	0.02628638442116191
13	0.2637879589747172	55969.75909607799	4.063119787839123	0.0291595348765065
14	0.2927376259758603	20402.351838801358	4.409062290691955	0.03525977460462214
15	0.3076033339602873	-10873.969340725318	4.749566106401687	0.036084415880783025
16	0.3369116669637151	-42928.91544535525	5.103328697163116	0.044076540386587454
17	0.35172399994917214	-63525.556235502176	5.421253993373715	0.038745691661673715
18	0.3809049579722341	-105948.6933970444	5.818765235944866	0.0703554550624216
19	0.39570637498400174	-136377.61731976288	6.186038498219217	0.08956216080503844
20	0.42473954198067077	-165541.97554229252	6.554205172526113	0.11615938447927442
30	0.6460121249838267	-350648.2549096943	10.258824175687966	0.4501141373552548
40	0.8665228329773527	-421468.089218312	13.756716788849403	0.5364913055768086
50	1.0872767489636317	-448573.47643583163	16.555396138498402	0.5611676253524074
60	1.3082574989530258	-457417.8536520089	18.297517821711008	0.5911608168747564
70	1.5291949149686843	-459470.92348895833	19.193238748369296	0.6096673157502177
80	1.7499865819700062	-460621.6345669247	19.661466388442868	0.6223294402338269
90	1.971698123961687	-461032.3905852317	19.949734982905074	0.6271998581882117
100	2.1924066659703385	-461029.7169623344	20.146783952248583	0.629716192384679
110	2.4004639579798095	-460872.4185433884	20.260791383825442	0.6308848234729219
120	2.621559499966679	-461446.6328197424	20.346392625203297	0.6294244730535569
130	2.8433849579596426	-461693.2754711352	20.28424894327953	0.618749345596139
140	3.066550749965245	-459726.78951810824	20.019602933954314	0.5995944228894794
150	3.2910687499679625	-461423.2056702522	20.451826635780144	0.6259477122072603
160	3.515915415948257	-459631.40511938627	20.098040243485084	0.6062835235383089
170	3.7397873329464346	-461300.0178988426	20.578317709443557	0.6326277546942092
180	3.9613202079490293	-460819.9929916063	20.45515235649507	0.6200189619946171
190	4.181550499954028	-460421.6844373105	20.45285834142375	0.620967116884689
200	4.399071540974546	-460366.3009629885	20.354009723480267	0.6141567349698136
300	6.580437123979209	-461510.6340248639	20.670162065272955	0.6311798885101705
400	8.775280373956775	-461109.31488016684	20.719260579863196	0.6320365045963248
500	10.952228248934262	-460929.44942193036	20.69657023395072	0.6330079251828638
600	13.169707081950037	-461746.79410263704	20.796054883620165	0.641493099766703
700	15.361864498932846	-461206.74202251434	20.796769048212926	0.6397848748850262
800	17.56121866594185	-461286.10636824416	20.757068661184558	0.639635109316252
900	19.772845624916954	-461390.8764714184	20.83882793997118	0.6422888450320935
1000	21.9738249159127	-461407.52806499985	20.643298173687565	0.6302184334096915
1100	24.19308987489785	-461586.3464608615	21.714727552295344	0.6847989458556734
1200	26.371729499893263	-461309.2493362168	22.128736288275782	0.7040388033245408
1300	28.56248712487286	-461697.02369169856	22.373184078092287	0.7140875046042392
1400	30.739504958881298	-461532.060447717	22.53327339722152	0.7223438466197112
1500	32.94195437588496	-461538.6328266802	22.640227457251804	0.7269184999698559
1600	35.147573417867534	-461570.39979150286	22.719460443849783	0.7318081471035895
1700	37.32171300187474	-461478.2880634861	22.781911836912393	0.7350081597576864
1800	39.529310834856005	-461589.6973758689	22.83036748998034	0.7369972607890444
1900	41.72529770983965	-461593.3157013211	22.864685101755235	0.7370869252314568
2000	43.929937501845416	-461672.062019953	22.895891268729272	0.7399227241661679
2100	46.11643445986556	-460831.29395805806	20.936662399543998	0.6369015618657763
2200	48.3090147098701	-461655.8055382073	20.64523631901625	0.6214557969374217
2300	50.505537125864066	-460888.7439540354	20.719644310526093	0.6340752176739488
2400	52.70069054188207	-461378.7613495	20.773491251355875	0.6390824101918089
2500	54.89738570788177	-460745.52557162923	20.59091357192212	0.6253200337119199
2600	57.099454790877644	-461582.4252614616	20.799306139559324	0.6422328118940932
2700	59.29250366589986	-460975.75993916084	20.751926230670588	0.6381800903861549
2800	61.52068008290371	-461640.6205145294	20.801231183914133	0.6400815578204478
2900	63.72112591590849	-462102.25989439106	20.77663659142264	0.6391663323752093
3000	65.95480183290783	-461972.59659505886	20.774923493991242	0.6367666652511604
3100	68.1544514579	-459706.3833287126	20.406713505413517	0.6097266340231027
3200	70.3588437079161	-459540.3579830785	19.891873668526927	0.5817929157252311
3300	72.58158487491892	-460809.57105874363	20.438339788572712	0.6107432619285195
3400	74.79197470791405	-460835.99837281246	20.519542876967876	0.6156733969525865
3500	76.97564924889593	-460900.0054587607	20.63829058443242	0.6292179164119205
3600	79.18227833291166	-462037.07721628685	21.681415099018988	0.6830625002904731
3700	81.40557187490049	-462087.17023677495	22.106245075618002	0.7008438743463233
3800	83.59339433291461	-462030.6164721507	22.354595356880775	0.7131423049886483
3900	85.83075062491116	-462155.1903594393	22.51666247477873	0.7205039944188413
4000	88.00754766690079	-461949.80336541944	22.628432061240154	0.7260688222455439
