# 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_primalonly_15_10000_9000", dynrange = 1.0, Λ = 1.0)
iter	time	function_value	psnr	ssim
1	8.415983757004142e-6	561764.4291806474	-0.02636674332749045	0.0041734697881560985
2	0.030749166005989537	477680.397345021	0.4726471684408935	0.007391363917083196
3	0.045588458015117794	445904.3860085457	0.8541599527624693	0.01006979171856906
4	0.07592187402769923	399681.3677583085	1.181965605751772	0.012069788261175638
5	0.09077079003327526	362478.51260571275	1.492698133779293	0.014052620637906929
6	0.12160779003170319	323116.3309646931	1.7766406376920196	0.01512863007554914
7	0.1364782480231952	283959.62693412084	2.084634767165375	0.01602461735201811
8	0.16718812301405706	245991.7949915964	2.399048048802981	0.01777848865112338
9	0.18201504001626745	205261.68116779646	2.7259819146624027	0.0198928671550234
10	0.21358712401706725	166512.76720928695	3.048484582012936	0.021532990727437577
11	0.2284047490102239	129423.52195965033	3.385616645509007	0.02339917084037043
12	0.25908779102610424	92430.72145725321	3.719120102676909	0.026235138135556363
13	0.2739251240273006	56234.45633965816	4.06265979353628	0.02904285143897236
14	0.30462654001894407	20664.24094726918	4.408670596923166	0.035219079051320976
15	0.31957158201839775	-10609.381395861095	4.748935426912668	0.03604844848787439
16	0.3504191660322249	-42661.59035355207	5.102505281821897	0.04398840820235446
17	0.36524299904704094	-63255.65798127333	5.420140839248237	0.038603570951246294
18	0.3960326650412753	-105685.50779722535	5.817603557628814	0.07007435392424803
19	0.4108435400412418	-136113.33182701285	6.184913540940782	0.08911463336812551
20	0.4416858730255626	-165273.28286469943	6.553314331399049	0.11593650820431345
30	0.6691826650057919	-350483.0019432482	10.257701767457563	0.44966262724791584
40	0.8977807480259798	-421376.3696287504	13.755998333379935	0.5365639566579641
50	1.126997373008635	-448506.1817156033	16.55585110058356	0.5615250996648822
60	1.3573274560039863	-457367.278568895	18.298713480854506	0.5919224753440784
70	1.587885664019268	-459426.4882717257	19.19678185696968	0.6107271206636833
80	1.8165665390260983	-460582.16714712046	19.668164221450816	0.6231754868361841
90	2.045430206024321	-460994.3581329876	19.95844973973615	0.6276943480195302
100	2.273259956011316	-460996.05744529155	20.159705389263966	0.6306193343645952
110	2.4837922060105484	-460842.6671235732	20.273275564218256	0.63164131061065
120	2.7119761219946668	-461417.385214917	20.36072863701677	0.630007391405238
130	2.940203121979721	-461665.4660065767	20.293345111364772	0.6190438098421256
140	3.1677167049783748	-459692.7824968839	20.02349060492928	0.599184680141652
150	3.3947454129811376	-461394.24818732595	20.465030921427505	0.6263037242655456
160	3.6232799549761694	-459601.54683266685	20.103084853065205	0.6061660433587956
170	3.851793704990996	-461272.5478199133	20.596458900839565	0.6332099258788542
180	4.079346496990183	-460790.29285982205	20.46982947996957	0.6203140590949733
190	4.305637455981923	-460393.1500874451	20.466489672479852	0.6212235905681449
200	4.530728997982806	-460339.3084401261	20.370181767563107	0.6145768096913627
300	6.754097539989743	-461485.3614342907	20.692062460221745	0.6318770427444349
400	8.979927039996255	-461079.8071521986	20.734956392948103	0.6324710297825391
500	11.190726123022614	-460901.2263000063	20.710913093423684	0.6333472593914683
600	13.401504164008657	-461718.9832284724	20.81695904698924	0.6424006746895787
700	15.594332956010476	-461179.56224915	20.820102321469275	0.6407905476791653
800	17.774851040012436	-461260.68717863254	20.77865508062578	0.6404572206467483
900	19.95233458201983	-461363.57467234455	20.86372709664447	0.6432793041810991
1000	22.133819666021736	-461382.1738991737	20.66568004660494	0.6308858166157971
1100	24.305347123998217	-461530.9182575682	21.655390865122865	0.6827350108798891
1200	26.45305345702218	-461249.3433554145	22.010792114065772	0.6998542961654229
1300	28.594300457043573	-461636.34264641884	22.212349649320466	0.7082882509972074
1400	30.740633540030103	-461471.0447201709	22.33671066164062	0.7151614332208986
1500	32.8759546240326	-461476.0966191066	22.414328932732875	0.7186539924000441
1600	35.008345916023245	-461507.84034008486	22.468877132163282	0.7224816721872414
1700	37.13606066603097	-461415.52910419164	22.50890141212637	0.7247210937097316
1800	39.255182750028325	-461525.8074746094	22.538688852399122	0.7259340964098775
1900	41.374793250026414	-461529.332214506	22.558363223291927	0.7253043709097429
2000	43.4898561250302	-461608.5364723388	22.576345183073926	0.7274453508996068
2100	45.600828500028	-460798.63694839977	20.938046288504218	0.63656683137989
2200	47.66588137502549	-461626.8012585868	20.657439279230168	0.6215289075965023
2300	49.73420629100292	-460860.60337970284	20.736301980115265	0.6343957598347814
2400	51.817351582983974	-461352.2100266137	20.796374803916144	0.6400986819368212
2500	53.87817212397931	-460717.3587848333	20.604913499875668	0.6258168449440267
2600	55.92532595797093	-461553.5809320311	20.815509553363807	0.6428260520826519
2700	57.95500404096674	-460949.52043074305	20.77040390160846	0.6389166545074605
2800	59.981950040964875	-461613.11011076963	20.8234747171629	0.6410440047273354
2900	62.01221695798449	-462075.06543931685	20.79658819941795	0.6400555106379036
3000	64.04199791597784	-461943.6376375339	20.791911500339147	0.6375032694531415
3100	66.0718586659641	-459674.9842983623	20.411271271429342	0.6097147817219313
3200	68.10194095794577	-459509.4832587055	19.89026562792709	0.581414393751555
3300	70.13330870796926	-460782.64825368853	20.450413379545985	0.610574634570523
3400	72.1690729999682	-460810.0942876778	20.533007992925235	0.615526647958307
3500	74.21155837495462	-460872.63050879457	20.65528584387166	0.6299293023453928
3600	76.24969191694981	-461980.6096635467	21.620343000726994	0.6810593228375359
3700	78.31283991693635	-462026.97517798195	21.993853277190112	0.6970740522457501
3800	80.37395679194015	-461969.90536607994	22.202584849812126	0.7078705158244831
3900	82.40997599993716	-462093.621134614	22.329357916724145	0.7138514902973426
4000	84.44232845795341	-461888.22066056996	22.412510465316995	0.7182951186526492
