# 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, save_prefix = "img/brainphantom256x256_pdps_known_primalonly_10_10000_9000", dynrange = 1.0, Λ = 1.0)
iter	time	function_value	psnr	ssim
1	1.362498733215034e-5	559998.953614435	-0.02636674332749045	0.0041734697881560985
2	0.02784554098616354	476637.35912192235	0.47267920025842847	0.007391754948748003
3	0.041244540974730626	445197.7947481975	0.8540720863234396	0.010069591451451773
4	0.06904329097596928	399148.83304854366	1.1817876128730762	0.012067979032692038
5	0.08244441598071717	362057.34112770046	1.4925445085252673	0.014048663581902442
6	0.11008258195943199	322712.0276680034	1.77625688092415	0.015113954262856466
7	0.12350737297674641	283613.942531232	2.0841034085453805	0.015997034260319107
8	0.15137403999688104	245692.6887982395	2.3985638099617166	0.017752176084650434
9	0.1647218319994863	204991.2014572536	2.725536206825913	0.019868554868042693
10	0.192437498975778	166249.2978845601	3.0478886908498932	0.021496579449602808
11	0.205797206988791	129185.61702440766	3.385058820207485	0.023370880491656944
12	0.2336025409749709	92203.74685047589	3.7184541719575357	0.026191260906605833
13	0.24706212399178185	56025.27350298222	4.061921750086968	0.02899817938262596
14	0.27482458198210225	20461.21758546986	4.407828903050753	0.035128934545417295
15	0.2881886239920277	-10803.835451415036	4.748127372773057	0.03593677250031736
16	0.31598945800215006	-42847.47513523618	5.101695794560708	0.04389137899125017
17	0.32935162499779835	-63443.03086125404	5.419338240571236	0.038525443219608387
18	0.35718479100614786	-105849.34578828198	5.8167690134808305	0.06991085698562982
19	0.3704934990091715	-136270.28729403834	6.18398950960363	0.08887967661497398
20	0.39825449901400134	-165423.9498583278	6.5522365427921105	0.11559522628463147
30	0.6036216650390998	-350573.8956128067	10.256017481392474	0.4488015606315464
40	0.8078587070631329	-421450.4933523644	13.752475400285118	0.5363463569506601
50	1.0131478320690803	-448565.9078563955	16.549782624811098	0.5612579634887314
60	1.219132957077818	-457413.77144058904	18.29261865964377	0.5914725586549482
70	1.4242339570773765	-459467.16001480375	19.191430899016197	0.610261019676593
80	1.628582707082387	-460619.61631071015	19.66434193789005	0.6229644137988117
90	1.834326957090525	-461031.0386310068	19.956398383022776	0.627803601691895
100	2.0403210820804816	-461031.13346871175	20.159135779817227	0.630741857396114
110	2.2310876240953803	-460875.74331417785	20.27523472254238	0.6318715477288588
120	2.436164790095063	-461450.2475751854	20.36413428421689	0.6303173218704683
130	2.6420754150894936	-461697.6980984421	20.29720805229376	0.6191836439919975
140	2.847974540083669	-459724.2421457958	20.027421845526327	0.5994916723929916
150	3.053096415096661	-461427.2499585363	20.47259842608105	0.6267595721326389
160	3.257866915082559	-459633.46127603634	20.107912575821327	0.6063536965936422
170	3.4632560820900835	-461305.8989554447	20.60744054689911	0.6337971522556158
180	3.669635540107265	-460823.3785132904	20.48138582818285	0.6208835347608487
190	3.874398415122414	-460426.2850959753	20.47889956978143	0.6217922978560577
200	4.078688207111554	-460370.9770649099	20.380778268882427	0.6150258549859013
300	6.118503416102612	-461520.1376735285	20.70880961871906	0.6325514839724236
400	8.155699208116857	-461116.24131294305	20.75656769628761	0.6333643365930558
500	10.192763667117106	-460936.683348056	20.733395906811815	0.6342673215053536
600	12.246076500130584	-461755.4454213588	20.83835424322887	0.6431341639349641
700	14.302084750146605	-461216.4835107925	20.841205597404322	0.6416098140941502
800	16.340288000152213	-461295.9977181666	20.79793714920456	0.6411942690991013
900	18.398094333155314	-461400.0677499938	20.883949797762426	0.6440362068034247
1000	20.44767337513622	-461416.8657634884	20.68173142246717	0.6314792278695214
1100	22.49863491713768	-461588.1847619136	21.739560272029053	0.6853997285028702
1200	24.571823583130026	-461309.6783841132	22.144193399933002	0.7043563936452136
1300	26.638367000123253	-461696.888455064	22.38344112962733	0.714285377264023
1400	28.6956044161343	-461531.70410181576	22.540770399427323	0.7224340698755756
1500	30.74302745811292	-461538.05449086736	22.645864804974526	0.7269549321108091
1600	32.788170124113094	-461569.84094938164	22.72407969885808	0.7318112774335912
1700	34.84249912411906	-461477.7612231863	22.78539312208317	0.7349524412817143
1800	36.89821208210196	-461588.9884935569	22.8329513456417	0.7369587678028388
1900	38.959518665098585	-461592.7915392147	22.867055422906244	0.7369801819508492
2000	41.02157912310213	-461671.28494592634	22.89815904497654	0.7398011114615928
2100	43.08711653909995	-460836.5429331339	20.97280086513432	0.6383892738854783
2200	45.140766581112985	-461662.63731138944	20.679595343746477	0.6225348046167729
2300	47.18438574808533	-460896.4318020762	20.75675276002409	0.6353095491499551
2400	49.25304383109324	-461388.4795449203	20.817256212980976	0.6408292491449734
2500	51.30997299711453	-460752.1507607577	20.622606972880625	0.62638038478049
2600	53.35904570511775	-461590.6152649233	20.838921663380848	0.6436822455613906
2700	55.413077538105426	-460984.9326845244	20.791739415696547	0.6396629252393182
2800	57.451073121104855	-461649.6280823049	20.844840605857016	0.6418022359694735
2900	59.52025820410927	-462111.7301477763	20.81824818388245	0.6407750592454109
3000	61.57589291210752	-461980.35959904495	20.813427975481243	0.6383245494698039
3100	63.62265437009046	-459709.25462360156	20.430860715066178	0.6104528798986557
3200	65.6799219120876	-459542.04862854624	19.90164077519999	0.5816407949088012
3300	67.73195774509804	-460816.1629657935	20.469083956654607	0.6113765177168107
3400	69.78927820309764	-460843.63436738896	20.55242485248078	0.6164338205648605
3500	71.83943599509075	-460907.08966585406	20.67517585158097	0.6306037845625683
3600	73.88942099511041	-462038.7018278711	21.703228135665636	0.6836132840214901
3700	75.97925149509683	-462087.4079842227	22.120189477145065	0.7011464160108352
3800	78.0274011620786	-462030.12746266456	22.364756970963285	0.7133359233565696
3900	80.08979612108669	-462154.87912725785	22.52437596511985	0.7205670193348057
4000	82.16952887107618	-461949.3815315335	22.634326310508794	0.7260937418260256
