# 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.25, ρ̃₀ = 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_25_10000_9000", dynrange = 1.0, Λ = 1.0)
iter	time	function_value	psnr	ssim
1	6.541988113895059e-6	565295.3803130721	-0.02636674332749045	0.0041734697881560985
2	0.028598583972780034	479767.79081317846	0.4725805895185553	0.007390520849807125
3	0.04230470897164196	447311.88610633253	0.8543522605730434	0.010070897206036938
4	0.06994341796962544	400742.03151101933	1.1823611122076518	0.012075488166981932
5	0.08372083396534435	363316.2480551129	1.4930326960382905	0.01406172802791799
6	0.11247062595793977	323918.00531252107	1.777390888186576	0.015157019493246729
7	0.12630066796555184	284643.8389803837	2.0856803877157795	0.016077804592049043
8	0.1549420009541791	246582.53273664557	2.4000004410784816	0.01782941618479579
9	0.1687113759398926	205790.80121382014	2.7268564503497235	0.01993796577127766
10	0.19735900091473013	167027.836289112	3.0496076691978447	0.021596967446818265
11	0.2111072089173831	129893.72926449643	3.3866680613952687	0.023449157768083775
12	0.23970758391078562	92877.66677149477	3.7203351620662306	0.026307411083642
13	0.2535049169091508	56645.530108664265	4.06400426887502	0.029123101768130856
14	0.2831251669267658	21062.368154960026	4.410222965270169	0.03537991282516976
15	0.296880125941243	-10227.384113214423	4.750416276959662	0.03624109127636443
16	0.3256809589220211	-42298.339176093075	5.103981258611529	0.04415354908579572
17	0.33945825093542226	-62890.40751351866	5.421603828570053	0.0387326100760103
18	0.368043042923091	-105365.8652120626	5.819140795198181	0.07036937219571023
19	0.3818204599083401	-135807.52609026985	6.1866143661415185	0.08956558299708117
20	0.41045916790608317	-164979.24441912354	6.555290473337654	0.11660888271017758
30	0.6179573339177296	-350303.053664494	10.260709569594326	0.4512442992315437
40	0.8287337509100325	-421233.12819506414	13.762461846680212	0.5369751775495705
50	1.0394433758920059	-448392.3134153097	16.56700518009705	0.5619915416853828
60	1.2465871258755215	-457278.07657181035	18.30967557015488	0.5926371094112792
70	1.4552344178664498	-459347.41557694564	19.20617865524298	0.6114364903122191
80	1.663843126851134	-460509.7602532953	19.67447557227896	0.6234767477715365
90	1.8729684608406387	-460923.4502357309	19.961570168162265	0.6274014644055955
100	2.083736876840703	-460928.09482368454	20.160077234752535	0.6302747623139757
110	2.278633460839046	-460778.70285707404	20.27019910183446	0.6311934933004957
120	2.491168752836529	-461353.8328189858	20.354750655460748	0.6294490886183985
130	2.7036972528439946	-461602.86620190006	20.285653267478445	0.618725209185867
140	2.9144578778359573	-459630.9286635343	20.01567626877283	0.5986882956222734
150	3.1221000858349726	-461330.88436953386	20.450705043380907	0.625393192209597
160	3.330150835827226	-459539.8278798046	20.09336322600251	0.6057010818704177
170	3.5397459608211648	-461208.6111581518	20.57649419191543	0.6320641957346145
180	3.750879293831531	-460726.638801538	20.44852464241106	0.6191806728278418
190	3.9583590848196764	-460329.0786615409	20.44299259217176	0.6200424586286135
200	4.16462387581123	-460277.87689042266	20.350019138528552	0.613639290922346
300	6.206045417813584	-461418.106759617	20.659653758170982	0.6304591204911283
400	8.298484042839846	-461010.0580590193	20.69521874805192	0.6306104493244393
500	10.378294084832305	-460833.0573997396	20.669307424012135	0.631553497975994
600	12.468085209839046	-461648.9030931444	20.7775013946519	0.640877561577347
700	14.513377000839682	-461108.29876310064	20.780746627739216	0.6391411348817186
800	16.563749750814168	-461192.8371295044	20.743000576121005	0.6390165859235604
900	18.623746542813024	-461293.60741329315	20.826498084162214	0.6417255881046122
1000	20.696580042800633	-461315.0626324684	20.635910262128863	0.6296310075226695
1100	22.743162876809947	-461431.55619929446	21.495522508427776	0.6762915372541085
1200	24.784542209818028	-461143.8493024875	21.766195668476445	0.6898492696700691
1300	26.82335179281654	-461529.7337628505	21.908019792840744	0.6958213463662213
1400	28.872307917801663	-461363.7378072219	21.985466043932774	0.7001549163878023
1500	30.915947208821308	-461367.0790040521	22.02879601807448	0.7020280334240361
1600	32.972215249814326	-461398.71288982074	22.05733400596702	0.704409387690822
1700	35.03863012479269	-461305.6964517738	22.077634214629377	0.7056217498616107
1800	37.077928916784	-461414.1496443342	22.0952511597853	0.7061438431323386
1900	39.13322824976058	-461417.79597357713	22.103256960819348	0.7049404932898304
2000	41.19128145877039	-461498.1667958232	22.110531904095193	0.7064510865478301
2100	43.244163374765776	-460724.8011355644	20.86518149913903	0.6330539134983187
2200	45.286666791769676	-461557.93326025194	20.613012825804844	0.6194190535262459
2300	47.3439266247733	-460792.32572446193	20.699447538582806	0.6325953499032073
2400	49.403115041786805	-461282.20131972525	20.75727781625556	0.6386173267876236
2500	51.45802879179246	-460650.05413786555	20.571840003283537	0.6245720602062248
2600	53.53432891677949	-461482.50517709996	20.77257998556233	0.6411705180258984
2700	55.58641237579286	-460881.30462796876	20.730469297979557	0.6373600832398384
2800	57.64059404178988	-461542.94555254147	20.784462867918347	0.639509080358093
2900	59.73259312479058	-462004.4712121872	20.756828049660314	0.6385450365098464
3000	61.808393749786774	-461873.182586866	20.752377344284575	0.6358161845736128
3100	63.90586199981044	-459608.70847869565	20.372987919750294	0.6081622080033958
3200	65.9708335828036	-459445.6620445786	19.865687198187874	0.5806973019702616
3300	68.0329429168196	-460718.58563566656	20.416755690070588	0.6090340837739474
3400	70.12093508383259	-460745.79147253657	20.497851917166493	0.6138606095394454
3500	72.19942287582671	-460806.30597851065	20.61931626092776	0.6285310180337031
3600	74.33442329184618	-461880.2292630959	21.461280634218667	0.6748338256146297
3700	76.49264133386896	-461921.4581927941	21.755091184629805	0.6877631201881942
3800	78.61692683387082	-461863.6752295942	21.90965323238462	0.6961454177667153
3900	80.77627791688428	-461985.2919894206	21.993361693132535	0.70003920450728
4000	82.94164525088854	-461780.120005241	22.044008210925256	0.703023104648399
