# 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, phantom_ρ = 100, save_prefix = "img/brainphantom256x256_pdps_known_proximal_15_10000_9000", dynrange = 1.0, Λ = 1.0)
iter	time	function_value	psnr	ssim
1	7.208000170066953e-6	561764.4291806474	-0.02636674332749045	0.0041734697881560985
2	0.029386874986812472	477772.60352856905	0.47416210950828774	0.007415531622632549
3	0.04420633299741894	445831.54201752564	0.8556395922189581	0.010119589144049099
4	0.07331924999016337	399652.6986019082	1.1834220640538975	0.0121473521051133
5	0.08836183298262767	362400.4389922261	1.494333619204808	0.014149753113820093
6	0.11764654098078609	323024.03425067704	1.778254249533942	0.01522774635184404
7	0.13247116599814035	283882.8640395682	2.0858838868164566	0.016114624787739226
8	0.1616962069820147	245919.64670498838	2.4005183696836827	0.017897547454711907
9	0.17658429098082706	205195.3352253661	2.7276762982458527	0.020041397138472895
10	0.2055951239890419	166431.96416691126	3.0500610179333396	0.021669226677280964
11	0.22041212397743948	129350.44266024753	3.387480496295	0.023541047843757906
12	0.24945612397277728	92339.67907049721	3.7208727047749894	0.02638172254501852
13	0.26427566597703844	56147.77029549168	4.064487757284027	0.02929538506547372
14	0.29331962496507913	20572.118110416806	4.410528453062292	0.03540857231991382
15	0.3081306659732945	-10717.21148112784	4.7510965050743	0.03626512680470327
16	0.3384145409800112	-42786.43826475472	5.10495633423716	0.04426236117452481
17	0.3533202069811523	-63384.87202863484	5.4230117332797025	0.03893011155319675
18	0.38246474898187444	-105837.90759620132	5.8205728962821714	0.07071721320744478
19	0.39733216498279944	-136277.43848340202	6.18796899455673	0.09007742956462768
20	0.4263431649596896	-165454.1526176333	6.556268797608462	0.1167238847931111
30	0.645438123960048	-350596.48582448537	10.261929469722748	0.4515319317673516
40	0.8636804159614258	-421404.5116904354	13.762296794871032	0.5367584291902869
50	1.0820691239496227	-448518.3978258499	16.564121685319297	0.5614030823582099
60	1.3005964159383439	-457373.9085048995	18.30588057221346	0.5914100338476127
70	1.5189499579137191	-459432.86409579293	19.19921858592214	0.609679583569413
80	1.7374591249099467	-460585.59002685547	19.66333186496991	0.6220467099606226
90	1.9560254999087192	-460996.87652548985	19.947649623108553	0.6266827743602128
100	2.173608874902129	-460994.3826579701	20.140385260812828	0.6290103992314371
110	2.378357624897035	-460837.70738044876	20.25010008894009	0.6300870357362052
120	2.5969945419055875	-461412.0924345887	20.332815310297896	0.6285638546521161
130	2.8144480418995954	-461658.89207825245	20.272447037818814	0.6182707705242505
140	3.03174116689479	-459697.5887098162	20.010457281062187	0.5991754925400398
150	3.248936666874215	-461388.0745708661	20.432047592212626	0.6249815949447837
160	3.4659844998677727	-459598.46030530456	20.086587765554206	0.6060105063583754
170	3.682812124839984	-461263.03743456514	20.55083790223796	0.6314303911125991
180	3.8998824578593485	-460785.1983080981	20.428615008101904	0.6190426413310206
190	4.116731790883932	-460386.24951323256	20.425940217084047	0.6199531030423062
200	4.333990249899216	-460332.37921665737	20.328014325985503	0.6132154651987625
300	6.510508208884858	-461470.7027605368	20.63198556687535	0.6297499279720288
400	8.689072999870405	-461068.78618832177	20.676103651816092	0.6304578785624024
500	10.872397541883402	-460889.52990054723	20.65336294609731	0.6313994099699308
600	13.04934995790245	-461705.69892908004	20.7511004631057	0.6399725002028397
700	15.216933290910674	-461165.2126459676	20.7517820270491	0.6380474790019816
800	17.381680123886326	-461245.2000950096	20.714654637638752	0.638138517691965
900	19.5476342489128	-461349.7385930514	20.79399066338013	0.6406764080286608
1000	21.74007312388858	-461368.10407306807	20.605447117399546	0.6289342386097674
1100	23.924216457875445	-461527.83026946115	21.616291773315723	0.6817806732254383
1200	26.09317041587201	-461248.647813059	21.98721232568007	0.6993059160637243
1300	28.260992207884556	-461636.5369824387	22.197349355010175	0.7079662332492418
1400	30.429287915903842	-461471.7551362793	22.326098831792347	0.714994318112899
1500	32.5962968749227	-461477.0226392162	22.405944486933382	0.7185020441894621
1600	34.76356783290976	-461508.7001345785	22.46188308603088	0.7223598761565722
1700	36.94523666691384	-461416.29635482014	22.503393805210333	0.7246931297489942
1800	39.1128853749251	-461526.79680164054	22.535319903902668	0.7259730883269055
1900	41.295438707922585	-461530.28649795457	22.555507666566527	0.7254438064507713
2000	43.46189366694307	-461609.84374154016	22.57423108919398	0.7276674269361936
2100	45.65475829193019	-460790.5980782592	20.882845551605985	0.6342910307319234
2200	47.82176862593042	-461616.32627589116	20.604385686471993	0.6197987399008086
2300	49.99434533395106	-460848.1997743013	20.67723909016766	0.6324957287172485
2400	52.1647558749537	-461337.6764989341	20.728983824634383	0.6374979506212082
2500	54.331151999969734	-460707.1907771798	20.55463800214148	0.6241803058581371
2600	56.500439666968305	-461541.28551366477	20.754409053805468	0.6406513388218529
2700	58.66917033394566	-460935.4430990912	20.708950402739337	0.6366894086388818
2800	60.88020183396293	-461599.3130413398	20.755683434617033	0.638431925728777
2900	63.04865087597864	-462060.66440661007	20.732091954982465	0.6376257718833535
3000	65.21420025097905	-461931.5317004886	20.73139839092415	0.6351303235745192
3100	67.39173433399992	-459670.6183546943	20.373206920737836	0.6084590027724093
3200	69.56313725101063	-459507.5774207035	19.87401427579497	0.5815126734639297
3300	71.73393629200291	-460772.1549017727	20.401709350018148	0.6095554652807589
3400	73.92400345901842	-460797.7804659181	20.480097357608685	0.614274134915554
3500	76.1067618750385	-460861.64256642124	20.597212522349935	0.6278432623094587
3600	78.2815276660258	-461978.03524543456	21.585960742251714	0.6802457280868743
3700	80.45292812405387	-462026.5804320709	21.973028053232863	0.6966235300338117
3800	82.62317233305657	-461970.6727958898	22.188329138762313	0.7076371472070837
3900	84.7941199990455	-462094.1608076924	22.318997047023824	0.7137776150322142
4000	86.97895983204944	-461888.7879098893	22.40499723244197	0.7182709936334599
