# 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_noprediction_15_10000_9000", dynrange = 1.0, Λ = 1.0)
iter	time	function_value	psnr	ssim
1	7.499998901039362e-6	561764.4291806474	-0.02636674332749045	0.0041734697881560985
2	0.027392791991587728	520530.870463891	0.21206741550356306	0.004011763030691591
3	0.04032858400023542	479985.8368318991	0.4663308950867685	0.004574691264828214
4	0.06776616699062288	438673.3897938882	0.7195881726740667	0.004248131432570155
5	0.0806483750056941	398096.90506932477	0.9572820306501764	0.0033185673268244804
6	0.10795720800524577	358753.0812446605	1.213803188498757	0.0032617276795959844
7	0.12083804100984707	319803.76206136314	1.4782130573365877	0.0034637155035860167
8	0.1480742910061963	281054.7065820969	1.7433168966793522	0.0023884058621807137
9	0.16094287400483154	242364.87745857172	2.023565721070737	0.0035113667305756587
10	0.18817570700775832	203882.28699976287	2.3053111621379934	0.0037153834569130044
11	0.20107454003300518	167102.53971932502	2.5907101054391504	0.0014215566657160852
12	0.22827495701494627	130708.37467803882	2.878719647103553	0.0017473641654045932
13	0.24116879000212066	94928.59062117501	3.177622121775929	0.0029341038537563048
14	0.26850012299837545	60839.54148341778	3.4630272317860786	0.0018744720016096936
15	0.2814390400017146	27264.721034679984	3.796196521798074	0.006773964975820862
16	0.3086632490158081	-4245.77969848381	4.0873750298498	0.005206999592546855
17	0.3218367070076056	-33747.833814727004	4.323337831540398	0.001778486612477689
18	0.34975608202512376	-64174.24768509745	4.706084450788128	0.00899564575588497
19	0.36271604005014524	-92737.60331175382	5.010938028871835	0.009513206426156972
20	0.3901843320345506	-120067.41039135329	5.318981358069066	0.011085212922876408
30	0.5954704570467584	-315746.1942847589	8.297337779321772	0.09271699619493276
40	0.8014624150528107	-407017.78266570024	11.057719738918998	0.27657124828767965
50	1.0065846650395542	-441171.17660338327	12.815058139593404	0.3847467355526497
60	1.2133112070441712	-452820.6992331579	14.243712397270961	0.4446997422275714
70	1.4188539160531946	-454846.17996091873	14.571914396643828	0.4569011983005828
80	1.6254644160508178	-456079.1841640098	14.926062713688141	0.478912216100132
90	1.8248776660475414	-457996.92442521313	15.719345891666343	0.49308884382510343
100	2.0252518330235034	-458302.09490080207	16.14845397475363	0.5094411970675139
110	2.2120233330351766	-456174.2969518729	15.720073020529586	0.4866280177574167
120	2.411957167030778	-457784.37145470426	16.159784473973577	0.49485216593751336
130	2.611702917027287	-457572.927780753	16.28179346685444	0.4845204733035019
140	2.8123948340071365	-456164.3567261998	16.497253044474736	0.48147982156166935
150	3.0121579999977257	-459055.40283440665	17.163089661412094	0.5159470579455723
160	3.211642375012161	-457772.57691079704	17.306091185525187	0.5039505458241422
170	3.4124848330102395	-459222.3024026213	17.533392399313154	0.5238901302409257
180	3.6120526659942698	-459187.9366404675	17.743624642992405	0.5252546515104065
190	3.8126501250080764	-458767.1021743071	17.92834096382889	0.5159041227220037
200	4.0126757500111125	-459062.42995950405	18.10017688821087	0.5250623402993604
300	5.997012458014069	-461008.24518443976	19.246291231428334	0.5694446652674443
400	7.991836375003913	-460323.1141507816	19.288888405423986	0.5755707032804991
500	9.97516450000694	-459119.922754881	18.496593174892368	0.5594615645316354
600	11.969736291997833	-460835.6919411478	18.857198610534386	0.5763808707785094
700	14.007171541976277	-460141.2237245469	18.852426092329438	0.5635060351168601
800	16.01030924997758	-458932.4253019166	18.188664222443556	0.5545149403755596
900	18.021704541984946	-460280.2986701431	18.951522931381522	0.5732961462398796
1000	20.007491416996345	-460792.6303104337	19.257243133077104	0.5776464432253751
1100	22.00448800000595	-461445.45187062793	20.94584751352056	0.652914360096078
1200	24.03432158299256	-461219.4533854995	21.585505313131698	0.6802418506733003
1300	26.02518279198557	-461621.7266319394	21.93391484717701	0.6946204964403497
1400	28.024397874978604	-461462.5889903576	22.148021441484936	0.7052851212796574
1500	30.019189082988305	-461471.1301179847	22.28247506711291	0.711652816715506
1600	32.04167058298481	-461504.40488798975	22.372154872072798	0.7172670708068869
1700	34.0250668329827	-461413.33191349363	22.436085556141244	0.7207620160858309
1800	36.023497999995016	-461524.29217828397	22.482218179933252	0.7228304670673408
1900	38.017824916983955	-461528.33272297896	22.51357582323644	0.722919525775275
2000	40.00033483398147	-461607.7135888349	22.539840804414904	0.7255484806685252
2100	41.98530279198894	-459281.75414725277	18.025835527211747	0.5464400202477273
2200	43.97091537597589	-461093.2110290225	19.047972603532827	0.5614701056493945
2300	45.95654020897928	-459975.84019264486	18.88360716262732	0.5737296007615795
2400	47.96479183397605	-460386.8684691104	18.81290258577932	0.5614262394430742
2500	49.95945604197914	-459725.04777559603	18.857897959802806	0.5605885276996079
2600	51.96820933397976	-459702.98100995633	18.35265797538008	0.5415867387243142
2700	53.9674101669807	-460297.45896365546	19.513762378699656	0.5872493201873015
2800	55.96272387597128	-460161.3174441715	18.71458728512711	0.5536405086532002
2900	57.97251075098757	-459466.8969991279	18.197186236192326	0.5454479520012276
3000	59.978143541986356	-459482.3468816678	18.016772273041106	0.5508689542807944
3100	62.09218191698892	-459200.7747675715	19.154788752992758	0.5628239646535916
3200	64.234035417001	-459242.3011468117	18.975492336324386	0.5555391872418787
3300	66.412120708992	-460096.73195519444	18.974177162604498	0.5605069155977577
3400	68.61117045898573	-457729.72478505666	17.837461748261795	0.5367819087628174
3500	70.82053733398789	-460108.88446760486	18.884348621496684	0.5710911279681508
3600	73.02428891698946	-461878.0293521962	20.806836105938388	0.6476695327307788
3700	75.23888608397101	-461992.3075567792	21.521403780006082	0.6751528404434636
3800	77.45018454198726	-461953.156448564	21.907216521285097	0.6929851687313527
3900	79.66850587597582	-462084.3403574424	22.136274337106062	0.7035520007911169
4000	81.87820100097451	-461882.6361696941	22.282260647962147	0.7112144403647047
