# 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, save_prefix = "img/brainphantom256x256_pdps_known_greedy_15_10000_9000", dynrange = 1.0, Λ = 1.0)
iter	time	function_value	psnr	ssim
1	8.875009370967746e-6	561764.4291806474	-0.02636674332749045	0.0041734697881560985
2	0.028772582998499274	477787.90502537775	0.474881350145489	0.007446706784916812
3	0.043068708007922396	445796.5139604893	0.8563392045892859	0.010140014244119741
4	0.07137016599881463	399630.28648164694	1.1839687719052918	0.012169042560036428
5	0.08544658299069852	362352.9295885187	1.4946921806668705	0.014176765028011653
6	0.11379362500156276	322983.2978668837	1.7785798513742865	0.015259610391056472
7	0.1279981670086272	283874.7378333874	2.086210696939494	0.01614214576207136
8	0.15658200002508238	245909.9330664781	2.4007337208906057	0.017924087250776433
9	0.17085262504406273	205186.36580690258	2.727860624352012	0.02006209687919481
10	0.19961779104778543	166438.22057161957	3.0502438476329297	0.021687995707120124
11	0.21401970804436132	129339.94833053659	3.387633794214822	0.02353053242183624
12	0.24268387403571978	92349.51308382273	3.7210213554643183	0.026384096271173545
13	0.25733637405210175	56163.41793451092	4.064592449017077	0.029247068451172257
14	0.28702920806244947	20589.327268742687	4.4105536018208955	0.035445692870818765
15	0.3012140000646468	-10704.308801423429	4.750962771988851	0.03629711181209189
16	0.3302844160643872	-42764.06247240811	5.104691776752636	0.044275663916027164
17	0.34499633306404576	-63348.33489266113	5.422588788893109	0.03890137032980774
18	0.37359633308369666	-105799.94272462094	5.820131808626229	0.0705633449993735
19	0.38774462408036925	-136219.93050198903	6.18745582902326	0.08984477426423787
20	0.4168417900800705	-165379.1810066226	6.555828557513723	0.11660197924943244
30	0.6295492480858229	-350533.31146477	10.261104862521691	0.4509090990291388
40	0.8435672071063891	-421384.37341525377	13.760770905572059	0.5367715299178065
50	1.0572596231068019	-448508.12264578167	16.561831324473154	0.5616768507185164
60	1.2715385401097592	-457367.6518744297	18.304591354340534	0.5920602039193352
70	1.48640320712002	-459427.58662295796	19.201122996102878	0.6107329458456907
80	1.7004540411289781	-460583.3876703253	19.67117432682113	0.6232513944879996
90	1.9142562491470017	-460995.31601619645	19.960899880824773	0.6278455379725596
100	2.129859665146796	-460997.4171955656	20.161812976259313	0.6308444099755677
110	2.3310014571470674	-460842.5104306236	20.273715923006172	0.6317122492470526
120	2.5472612491284963	-461418.32396805997	20.36350388847638	0.6301413226922864
130	2.763978999108076	-461665.5095108739	20.295144041408278	0.6190815788505846
140	2.9779436651151627	-459693.4910713445	20.024659838166382	0.5991804751950489
150	3.1917565821204334	-461394.7099176501	20.467816699169603	0.626542691905718
160	3.407289749098709	-459601.4032299323	20.103114896270085	0.6062325839342748
170	3.6221600410935935	-461273.1945508822	20.60150488796047	0.6335706081285106
180	3.836550625099335	-460791.3863588986	20.474005080338188	0.6205869136196456
190	4.052754792122869	-460394.3546733801	20.471688083746127	0.6214866584838913
200	4.274022083118325	-460340.324049501	20.37343142719761	0.6148152349785981
300	6.416236791119445	-461487.31775097514	20.699927060435993	0.632506964105056
400	8.546468457119772	-461081.2496257523	20.743720339823092	0.6329976533720505
500	10.686527332145488	-460902.66365448193	20.721534338639707	0.6339208835597449
600	12.837084915139712	-461721.6348957348	20.827796703748472	0.6431411757030019
700	14.97894949914189	-461182.8015565137	20.831756176070986	0.6415821453682623
800	17.109473874152172	-461261.86635108304	20.787046265294773	0.6410629777595157
900	19.250070540147135	-461366.4325537632	20.87470480206253	0.6441467415198197
1000	21.38179258213495	-461384.3091481187	20.6735478614814	0.6314221363358903
1100	23.507064416160574	-461531.3016470175	21.659047902927426	0.6830312744612064
1200	25.651739874185296	-461249.49397636566	22.01240931903927	0.700009353647822
1300	27.783712040167302	-461636.41076265025	22.21298748941351	0.708383811106849
1400	29.9184551651706	-461471.0704424349	22.33684776863136	0.7152175292915564
1500	32.05576808218029	-461476.10790114774	22.414239936995422	0.7186887693034332
1600	34.1870132071781	-461507.8457360435	22.468612106908346	0.7225008365174865
1700	36.34431049917475	-461415.529303938	22.508594360056517	0.7247342974912077
1800	38.48694324915414	-461525.80424103193	22.53837772544449	0.7259455326735014
1900	40.62843091515242	-461529.33064613893	22.558067945415885	0.7253127134012363
2000	42.75649941514712	-461608.53368312016	22.57607681952402	0.7274514956573607
2100	44.88783854016219	-460798.6626148674	20.93890722065451	0.6369279333915906
2200	47.032078332151286	-461629.13094870944	20.66451176871251	0.6221104726139777
2300	49.18128049917868	-460861.52685131674	20.741966689932532	0.634888432255575
2400	51.3207911241916	-461354.738371584	20.80694714035132	0.6408562930599891
2500	53.47735774918692	-460719.3454033566	20.612344535486613	0.6262552528658131
2600	55.61178058219957	-461556.0486742443	20.826643050538753	0.6435454592027079
2700	57.767710749205435	-460952.2348473655	20.782639380804945	0.6395885297032143
2800	59.90078287420329	-461614.96056364366	20.83330525314662	0.6416727439916047
2900	62.031497665215284	-462076.6502610843	20.80487093422091	0.6406065418681653
3000	64.1616682492313	-461945.8214475742	20.802314332537087	0.6380888722146421
3100	66.31465737422695	-459676.7851599706	20.41831786021174	0.6100362621558151
3200	68.44827629122301	-459510.7110227139	19.891844542947524	0.5815312947300172
3300	70.57656454123207	-460784.3799847545	20.45797855599076	0.6110074974419297
3400	72.70946012422792	-460811.06134112016	20.53920670283435	0.6159155840971617
3500	74.84022737422492	-460874.7963219757	20.66515182685627	0.6304795714452138
3600	76.97278391523287	-461981.05753031553	21.62498620950087	0.6813867016186097
3700	79.14901954022935	-462027.1430635164	21.9964082789435	0.6972840450614574
3800	81.33977549822885	-461969.9767263539	22.2041353439557	0.7080076102668783
3900	83.46733112321817	-462093.6585456898	22.330338312300704	0.713943886996943
4000	85.59347474822425	-461888.24090691196	22.413180101821773	0.7183597960602148
