# 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_primalonly_20_10000_9000", dynrange = 1.0, Λ = 1.0)
iter	time	function_value	psnr	ssim
1	1.1416996130719781e-5	563529.9047468598	-0.02636674332749045	0.0041734697881560985
2	0.027969208982540295	478723.86635741394	0.472614316916481	0.007390950701247952
3	0.04138429198064841	446609.11112007604	0.8542534085296956	0.010070229902135195
4	0.06907050000154413	400212.5521294223	1.182157012069232	0.012072309200388077
5	0.08258220800780691	362898.2297963329	1.4928609753526136	0.01405688655153932
6	0.11079316600807942	323518.3653354468	1.777019464404555	0.015142925401845865
7	0.12437520798994228	284303.01439888566	2.085161258934029	0.016051413524515147
8	0.15283320800517686	246288.3373402894	2.3995274193779057	0.017804162796283263
9	0.1666248330147937	205528.49188343948	2.7264220371171715	0.019915993586429696
10	0.19550987402908504	166772.6632505268	3.0490575401786084	0.021566592844316487
11	0.2092884160229005	129659.92783665891	3.386152939236541	0.023425353347655054
12	0.2380187910166569	92655.77068662681	3.719746054441953	0.026273795955601614
13	0.25187379101407714	56441.589052720075	4.063352887172254	0.02908460055966324
14	0.28092845701030456	20864.868133207037	4.409469395742472	0.03530258368410493
15	0.2946922070113942	-10416.74148512728	4.749697831345213	0.03615008036416433
16	0.3234319160110317	-42477.84204140522	5.103266241805367	0.04407482175628345
17	0.33723516602185555	-63070.337297614016	5.420894590513371	0.038671257529147775
18	0.36596408201148733	-105523.77513876763	5.818396526111193	0.0702255818438765
19	0.37984533203416504	-135958.70458176133	6.185791004982363	0.08934450534422858
20	0.4077442900452297	-165124.66547075234	6.554334251939668	0.11627737931121589
30	0.6136708730482496	-350392.6622369389	10.259267643256749	0.45047883885783624
40	0.82341995605384	-421303.9127842722	13.759333910933595	0.5367777287636595
50	1.0327229140675627	-448448.3939904807	16.561592037409522	0.5617676915900466
60	1.243159705074504	-457322.1361262212	18.304396335723876	0.5923110638304869
70	1.45113916308037	-459386.6370526964	19.20170138056132	0.611123217089754
80	1.6614027040777728	-460545.6349042846	19.671534693366752	0.6233511840425825
90	1.8731394960777834	-460958.55090008094	19.960187325226478	0.6275585110392166
100	2.082457830081694	-460961.76024366653	20.160038878467624	0.6304646436345803
110	2.276240330073051	-460810.39763043064	20.27167216496915	0.6314250198417443
120	2.4893671220634133	-461385.3373383841	20.357658865070526	0.6297257996218593
130	2.7025286630669143	-461633.91212114727	20.28950824123649	0.618895086117552
140	2.913505455071572	-459661.6746927151	20.019524260854467	0.5989227745242516
150	3.119934205082245	-461362.1713702133	20.45767379462715	0.6258399495810821
160	3.3275364960718434	-459570.41321647767	20.0982260455574	0.605945132600692
170	3.537113955069799	-461240.16975068813	20.586131821318347	0.6326292588154661
180	3.7488110800914	-460758.0965295466	20.458896240915045	0.6197426489008615
190	3.9550359970889986	-460360.77198140696	20.45449192403316	0.620639841758876
200	4.163410580076743	-460308.34330217086	20.359917758147734	0.6141085658920069
300	6.209775413066382	-461451.4501927924	20.675740725342934	0.6311798230336816
400	8.259066413069377	-461044.5050682531	20.714637128288196	0.6315464934397176
500	10.303505705058342	-460866.784172736	20.68967098114775	0.6324426543482435
600	12.360726747050649	-461683.55796069675	20.796821035956036	0.641653406091522
700	14.412162372056628	-461143.56994406116	20.80006041460581	0.6399669860249377
800	16.443685872043716	-461226.4025474807	20.760438301274764	0.6397296194652646
900	18.480401831067866	-461328.149611983	20.844621441422824	0.6425070169465426
1000	20.52367724807118	-461348.2966200302	20.650457284365004	0.6302645736696586
1100	22.566389373067068	-461479.49426298344	21.572603653963146	0.6795144446094246
1200	24.595038664061576	-461194.79871351196	21.883096302149248	0.6948066189441205
1300	26.633042331057368	-461581.1901964803	22.051951042198898	0.7019791373281925
1400	28.705013456055894	-461415.5473180235	22.14815358569895	0.7073841379190982
1500	30.781105331057915	-461419.5991404839	22.20441470379911	0.7099318819104695
1600	32.83703845605487	-461451.3362044265	22.242254646543497	0.7128995219324443
1700	34.88081508106552	-461358.6398237299	22.269058836143945	0.7145110966039772
1800	36.909250622062245	-461467.9828940076	22.2900779294736	0.7152609055732249
1900	38.993327330070315	-461471.4868080469	22.302286832662773	0.7142766485912573
2000	41.05312824709108	-461551.35889428295	22.31381037958778	0.7160158001749969
2100	43.15472374710953	-460761.2825640375	20.901206492710344	0.6347501863446159
2200	45.2223707891244	-461591.93724485405	20.635107168704426	0.6204701826148219
2300	47.27907178911846	-460825.99434226967	20.717282867247096	0.6334835948135353
2400	49.36809083013213	-461316.82439889555	20.776420184922657	0.6393569735611241
2500	51.47236987212091	-460683.36712313414	20.588042145887595	0.6252041760278695
2600	53.55706341311452	-461517.59649860335	20.793500417700944	0.6419901996858943
2700	55.624818580108695	-460915.0611538726	20.750054167682794	0.6381432273879163
2800	57.72617374712718	-461577.6334020887	20.803494762529766	0.6402816943270919
2900	59.83240124711301	-462039.3704012938	20.776206232747988	0.6393101421847449
3000	61.93398412212264	-461907.996813662	20.771676785349328	0.6366623044999391
3100	64.044800455129	-459641.5031279542	20.3920002449596	0.6089472802926528
3200	66.15119278815109	-459477.34857467946	19.87822805619336	0.581087140357154
3300	68.25394212212996	-460750.1986670097	20.433043155520632	0.6097867528647885
3400	70.35511670613778	-460777.54665878665	20.514960820506843	0.6146702756173856
3500	72.4273487891478	-460839.07120107766	20.636751422629473	0.6292350175577288
3600	74.51519503915915	-461928.5823309728	21.538165015867573	0.6779423640740626
3700	76.5682702481572	-461972.32910919335	21.869484120631427	0.6923363014738877
3800	78.6262610811682	-461915.0873571728	22.0490145818345	0.7020058243961733
3900	80.69734162316308	-462037.6231640439	22.151664048021956	0.706868165205581
4000	82.72620991518488	-461832.3300872139	22.21569034632101	0.7104498956983064
