# 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_noprediction_20_10000_9000", dynrange = 1.0, Λ = 1.0)
iter	time	function_value	psnr	ssim
1	9.749986929818988e-6	563529.9047468598	-0.02636674332749045	0.0041734697881560985
2	0.028203083988046274	522230.4121426869	0.21279333631592212	0.004015782672267686
3	0.041571167006623	481619.3631436144	0.46784761512442186	0.00458424595326561
4	0.06978720898041502	440241.3411652271	0.7219640607429758	0.004261470852131475
5	0.0831427089869976	399596.84641435475	0.9605771160065125	0.0033310742177269153
6	0.1113377089786809	360182.7230184093	1.2181080333204446	0.003277542257450899
7	0.12470233396743424	321158.8309147503	1.4836185817708007	0.0034846615050913984
8	0.15294237597845495	282327.72745036945	1.7499138574564652	0.002407495253489737
9	0.1663254179584328	243551.5715412255	2.0314605280591898	0.0035377242658122705
10	0.1945667509571649	204971.3052367428	2.3146112405493287	0.0037465857650616345
11	0.20793895894894376	168087.6077063319	2.6015506187766366	0.001445223920697655
12	0.23616591695463285	131581.95842499935	2.8911981929368222	0.0017796048358759112
13	0.24975687597179785	95678.00703163183	3.191860397111857	0.0029796778285004354
14	0.2780378349707462	61461.746567909475	3.47912397356731	0.001918121675735854
15	0.29176345994346775	27755.830565708828	3.8143838193720923	0.006887807405694318
16	0.3201446269522421	-3893.401121197153	4.107650269227293	0.0053264253518978365
17	0.33399250195361674	-33532.87057718799	4.345534238015767	0.001871913170431168
18	0.3622470019327011	-64090.021422350655	4.730985732776082	0.00927972036213337
19	0.3756835439126007	-92778.98444563332	5.038279724886831	0.009910620610891228
20	0.4044333779020235	-120223.47153768972	5.34888531669436	0.01169538908342028
30	0.6132007119012997	-316341.0958823532	8.357635065591127	0.1003391769802894
40	0.8219775448960718	-407280.41832160804	11.160242715041793	0.2887203887738666
50	1.0305977948883083	-441170.4307744398	12.953531639294782	0.39497702744138197
60	1.237591211887775	-452746.2199993733	14.429765327349909	0.45337045236288603
70	1.4473810038762167	-454786.5545682434	14.772252616086908	0.4645442965768186
80	1.6552954618819058	-456060.9388641407	15.140290546108087	0.48706626532067243
90	1.8681718778971117	-458005.4595382606	15.985904447503854	0.5016162924270622
100	2.07758825289784	-458332.28433385736	16.432883586094643	0.518465675603389
110	2.2723452939244453	-456244.8001112888	15.959814198927171	0.49431958167987833
120	2.4808666279132012	-457859.89652541955	16.409669673105153	0.5021904512040554
130	2.6883220859162975	-457658.4330316825	16.52268047581972	0.4920403244913355
140	2.8991401279054116	-456250.66166186373	16.72711190633072	0.48831333625252016
150	3.1062863779079635	-459122.6466961021	17.41561388739212	0.5241584455564836
160	3.315143377898494	-457838.8221641652	17.55702711176346	0.5113511766154799
170	3.522398460889235	-459288.09449293814	17.782500296587504	0.5316024643325038
180	3.7317859609029256	-459245.02460634016	17.98284849116627	0.5329505723636522
190	3.9397985858959146	-458819.1750288819	18.162845071003414	0.5237019940793066
200	4.146887252893066	-459107.29637756816	18.323311168994614	0.5334750342902224
300	6.203741794888629	-460997.77597588877	19.384257982696994	0.5773520879001242
400	8.266204669896979	-460297.5832697261	19.326303502041068	0.5785684566341357
500	10.320481878879946	-459112.306252804	18.50342630019859	0.5608502003769228
600	12.366652670898475	-460787.600547811	18.87865662223903	0.5786598988574547
700	14.39197358788806	-460099.6290365436	18.890632711295087	0.5664768083369967
800	16.410606670891866	-458909.6969265966	18.169656857871946	0.5551597909916686
900	18.43487117087352	-460256.09744240245	18.95909469466928	0.5744824222288933
1000	20.438566128868842	-460757.3889281838	19.28255265306516	0.5794103421405931
1100	22.447750253864797	-461392.27842845	20.941380926139786	0.6534063313853199
1200	24.44307308687712	-461165.7136774857	21.538875156253404	0.6792303756504747
1300	26.473579045879887	-461568.3221977633	21.847858494458343	0.6920818672122613
1400	28.48964442088618	-461408.6543519856	22.018547137817738	0.7006724190244622
1500	30.494022254890297	-461415.52933596604	22.1158221043029	0.7052011844930486
1600	32.49083650490502	-461448.7091727102	22.177464191585088	0.7094243711634876
1700	34.490654462890234	-461356.95833429316	22.21995122891297	0.7118903937760854
1800	36.48231975486851	-461466.86940096895	22.251813923293916	0.7132176152167666
1900	38.466594421857735	-461470.67210924195	22.27174980381048	0.7127223315498099
2000	40.453065713867545	-461550.6675872017	22.28816779282282	0.7147312803380866
2100	42.46963921387214	-459237.22029328096	18.107917370358763	0.5501779366438704
2200	44.46098704688484	-461066.44632157026	19.18100184623154	0.5675513457538064
2300	46.47050121286884	-459943.5562971577	18.931355413443796	0.5758428830839777
2400	48.52939867085661	-460354.0166790928	18.865081951625026	0.5637310439759535
2500	50.57449767086655	-459698.0172977293	18.90617159242285	0.5635546334273343
2600	52.61958746286109	-459669.880748081	18.361453882574658	0.5432155412088934
2700	54.63434887884068	-460269.61794375296	19.526725273852495	0.5881579787201829
2800	56.64962629484944	-460143.1212110824	18.728094516206617	0.5551444911398384
2900	58.68925583685632	-459469.73063015373	18.195745514717288	0.546451052601687
3000	60.71743204485392	-459455.4044809856	17.995290539777876	0.551119461453353
3100	62.750620711856754	-459164.5105180494	19.20894005690558	0.5659268156589001
3200	64.74787037886563	-459212.9589219852	19.022346322695977	0.558533388309909
3300	66.73720171186142	-460065.1604373697	19.002485865640864	0.5620524774619229
3400	68.74477333686082	-457698.43636588176	17.804722959645595	0.5371617189640575
3500	70.7534115028684	-460068.47931994154	18.90952236354945	0.5733104157508768
3600	72.80703483585967	-461823.8510147284	20.82700509179545	0.6487762017422679
3700	74.81133021085407	-461939.7948459078	21.50196277430135	0.6750930586595953
3800	76.82884762785397	-461900.771589917	21.838908940723908	0.6912482612030622
3900	78.83785746185458	-462030.3345486128	22.02280869345352	0.6999280944832776
4000	80.86260125384433	-461828.1447743712	22.13314917578639	0.7059597851425212
