# 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_dualscaling_15_10000_9000", dynrange = 1.0, Λ = 1.0)
iter	time	function_value	psnr	ssim
1	8.625007467344403e-6	561764.4291806474	-0.02636674332749045	0.0041734697881560985
2	0.02840716700302437	477750.0871983116	0.4727944903244874	0.007393534638103009
3	0.04356154199922457	445911.2963937415	0.8540761861508194	0.010072131606553826
4	0.07201504200929776	399705.91422077495	1.1817166219636182	0.01207129706203814
5	0.0864079590246547	362507.87416993175	1.4924757006186684	0.014049008034708325
6	0.11467904201708734	323146.952461721	1.7760151729737184	0.015108191661339605
7	0.12865683404379524	283981.67516653036	2.0837512221407133	0.015982012387161158
8	0.1570541670662351	246010.802260521	2.398239993397958	0.01773937484023549
9	0.1709365830756724	205268.3251055606	2.7252256812324744	0.019857151238669303
10	0.19913341707433574	166501.99727205865	3.047411863843154	0.021472489654892877
11	0.21327937507885508	129413.58672650585	3.384603232498194	0.023352236216542446
12	0.241477542091161	92410.31701395214	3.7178591125103737	0.026158875672678664
13	0.25553854211466387	56221.660540963356	4.0612291987287445	0.02897193262492952
14	0.284448167134542	20641.91327636402	4.407020555203639	0.03507073833314676
15	0.298305584146874	-10633.162622630509	4.747331213012094	0.03585109036815751
16	0.32691204216098413	-42690.724347838564	5.100870202104618	0.043820181517713094
17	0.341421334160259	-63293.16919669225	5.418487743986934	0.03845894035862168
18	0.3696674171660561	-105699.52673235018	5.815890285724544	0.06974261899971085
19	0.38370391714852303	-136128.26475905208	6.1829851376140645	0.08861108429456918
20	0.4125696261471603	-165285.10254382424	6.5510663945352725	0.11528193096295947
30	0.623694292124128	-350488.9246768555	10.254053791141626	0.4478387327161209
40	0.8346277501259465	-421396.5083914074	13.748604212213252	0.5361229300594774
50	1.0617425411182921	-448517.8580109811	16.543825691839906	0.5610124726378126
60	1.299742041126592	-457372.81727921445	18.287352838534773	0.591124422685422
70	1.5165846661257092	-459432.3112424644	19.187470363554205	0.6099445765490726
80	1.7306366241245996	-460589.5988086035	19.662381238928617	0.6229227630874126
90	1.942063333117403	-461003.41738116206	19.957162646634792	0.6280429966942783
100	2.1524315410933923	-461004.8909787913	20.162022203742367	0.6310478894422439
110	2.3515221250709146	-460850.62078651955	20.281595255731574	0.6323648965502002
120	2.5671906250645407	-461426.1974177158	20.372799003272437	0.6309610451422232
130	2.783018084068317	-461673.8603926775	20.304970007188384	0.6196304308251852
140	2.9988964170624968	-459699.0461199116	20.03580610448562	0.6001695688593818
150	3.214583751076134	-461406.2403825748	20.4868577030074	0.627488316749595
160	3.431308293074835	-459609.1535149404	20.11434759094803	0.6066298058817097
170	3.646896877064137	-461286.43607217993	20.628511757093793	0.6347422544957879
180	3.857774168049218	-460802.2913088794	20.502640642406437	0.621737742522334
190	4.070336751057766	-460406.00013474294	20.50147430597046	0.6228203940282614
200	4.282350584049709	-460347.1611084251	20.397013168009096	0.615491504578176
300	6.374203251034487	-461503.52585578227	20.737065766805536	0.6335788331163876
400	8.49777462601196	-461099.6348919455	20.793392282841687	0.6345943938694065
500	10.597963042993797	-460920.7718803054	20.77067362396887	0.6355283943315537
600	12.713830293010687	-461738.75319005456	20.87461086201666	0.6443512483049748
700	14.808132668025792	-461198.6329400556	20.87438915046812	0.6427418796790549
800	16.92705858402769	-461277.2193066326	20.829539535896217	0.6423713706259029
900	19.035732417018153	-461382.7984948262	20.91751954825242	0.6453027011901213
1000	21.13279316702392	-461397.8120509304	20.70820984514509	0.6324916124832246
1100	23.244127126032254	-461533.4191340724	21.675241874244673	0.6835042308478766
1200	25.390301543026	-461250.28330985014	22.021284279442167	0.7003487854977832
1300	27.516353209037334	-461636.7660385329	22.217963191247296	0.7086081032129683
1400	29.621415876055835	-461471.2447380908	22.33991887275551	0.7153724500545252
1500	31.72809520905139	-461476.2070636908	22.41629887952014	0.7188018501233393
1600	33.83117258406128	-461507.9124822948	22.470197444652477	0.722588523964396
1700	35.935097375069745	-461415.58269320626	22.509752159711816	0.7247975242025798
1800	38.03706204108312	-461525.8351032552	22.53926033022498	0.7259852880905031
1900	40.16495241608936	-461529.35733083595	22.558797635656433	0.7253440495260162
2000	42.26639445807086	-461608.5567739438	22.576704054300833	0.72747869645522
2100	44.38388979109004	-460810.4115249227	20.97969507545388	0.6389346488319273
2200	46.53480920707807	-461646.22775586694	20.7054667621733	0.623591716490285
2300	48.67784141609445	-460879.6905067876	20.788398589358728	0.6364680896330827
2400	50.776830375078134	-461370.7807381112	20.850228208998608	0.641957348655893
2500	52.87676416707109	-460732.16502670216	20.650786969276396	0.6271794367148912
2600	54.986830584093696	-461573.63497209456	20.87684247239335	0.6449957945648095
2700	57.114487376093166	-460967.705382693	20.82614778806271	0.6406696388412856
2800	59.26857270911569	-461631.85177789757	20.880100760881596	0.6430523709454216
2900	61.42454541812185	-462093.9048533554	20.85362992059981	0.6420011994041429
3000	63.53704600111814	-461963.6302101359	20.85016969647489	0.6396479290160653
3100	65.6513670431159	-459689.07805887266	20.461220747620317	0.6116302680721639
3200	67.75288691811147	-459518.1966020551	19.9163069387429	0.5818585444708159
3300	69.86483650110313	-460800.6698643648	20.500424380145862	0.6124995684458916
3400	71.97946508409223	-460829.6756734714	20.587154846762964	0.6177623182404888
3500	74.1186725840962	-460890.05914300407	20.707927199703583	0.6315562730973766
3600	76.26610145907034	-461983.697791248	21.64389244743553	0.682005402584309
3700	78.37587508407887	-462028.0299582434	22.006297864768637	0.6977169133997414
3800	80.53468683408573	-461970.3789106687	22.209404415105134	0.7082865232360976
3900	82.65739766709157	-462093.8825830858	22.33347019014792	0.7141336197253889
4000	84.76434708409943	-461888.364981896	22.415116241503554	0.7184822682376056
