# 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_dualscaling_20_10000_9000", dynrange = 1.0, Λ = 1.0)
iter	time	function_value	psnr	ssim
1	7.458991603925824e-6	563529.9047468598	-0.02636674332749045	0.0041734697881560985
2	0.029424667009152472	478815.8877180156	0.4728121726799157	0.007393873939648631
3	0.04358695898554288	446621.01861306693	0.8541320904970705	0.010072996745564798
4	0.07331191797857173	400246.5976206233	1.1818011527566534	0.012073039664698679
5	0.08814312596223317	362939.42699044186	1.4925472412903251	0.014051465408188281
6	0.1169954169599805	323561.56825124007	1.7761864497222082	0.015115982742480386
7	0.13146695794421248	284334.7445670818	2.083984724712945	0.015995479575748504
8	0.1617940409341827	246317.43252711956	2.3984532997457633	0.01775300802462601
9	0.17621249894727953	205541.18394848844	2.7254194078883742	0.01987000906395938
10	0.2056500819453504	166760.72379936458	3.047659401876086	0.021490143444487752
11	0.22073262394405901	129647.68974389756	3.38483165103416	0.02336657830950552
12	0.25023349895491265	92628.68937525783	3.7181241614148997	0.026179082006162938
13	0.2646568319760263	56425.64864111258	4.061515258347494	0.028993688389069808
14	0.2947926649649162	20836.433660134353	4.40732646249738	0.03510883949506924
15	0.3092432899575215	-10447.174209328285	4.747618196325444	0.03589697170627342
16	0.338872456952231	-42515.32141731963	5.101149517948931	0.04386117251882236
17	0.3533309979538899	-63119.92731434457	5.418751879823237	0.038486425229260726
18	0.38261895696632564	-105541.5590402533	5.816154608941648	0.0697799959817118
19	0.3972433319722768	-135978.0067528888	6.1832696742224735	0.08865849326642813
20	0.42656941595487297	-165140.42845348633	6.551405598593208	0.11537748725480909
30	0.6460802909568883	-350400.964625847	10.254553895713274	0.4480472948524359
40	0.8662906249519438	-421329.19411009795	13.749689090507713	0.5361863559230877
50	1.0854535839462187	-448461.76957347745	16.545919845403713	0.561111758201957
60	1.3040628759481478	-457328.0313995502	18.289752616418966	0.5912700910568013
70	1.5243164589628577	-459393.6253244579	19.189672613018946	0.6100606982283885
80	1.7422732089762576	-460554.6697854429	19.664051519895654	0.6229758318346524
90	1.9607078339613508	-460969.91655252344	19.958358234686596	0.6280277094389275
100	2.1794170839712024	-460972.8429374658	20.16286433750038	0.631024707112307
110	2.383851208956912	-460819.9756663424	20.28182889835419	0.6323596563905578
120	2.6032532089739107	-461396.1571078476	20.372963420122762	0.6309302960109029
130	2.8211496669682674	-461644.35279775906	20.30496029353259	0.6196621176027335
140	3.039623666991247	-459669.82433511456	20.03605616783858	0.6001705381214962
150	3.2599434169824235	-461377.35270378046	20.486447655698687	0.6274620176056287
160	3.4781825419631787	-459580.0507116893	20.11379030880716	0.6066532546411252
170	3.696168583963299	-461257.50699902576	20.62777850610058	0.6346811376387271
180	3.9152012509584893	-460773.2431908192	20.501751581510604	0.6216592463738658
190	4.1331271259405185	-460377.1688911911	20.50051644746936	0.6227511068280037
200	4.351693667966174	-460318.23402201524	20.395684534611476	0.6153784219961699
300	6.529022000962868	-461474.8130140242	20.7352801847634	0.6335093814404125
400	8.694981916982215	-461069.68202049204	20.790708368997876	0.6344496899833714
500	10.850913000962464	-460891.5380552021	20.767851301424326	0.6353979422863717
600	13.026992833969416	-461708.6926510307	20.87160501616868	0.6442133578777001
700	15.19843595896964	-461168.18074641726	20.871303131511645	0.6425844047058101
800	17.365949917963007	-461247.0025275965	20.82680637872876	0.642275216497501
900	19.540130876965122	-461352.5949219489	20.91470643832399	0.6451834225975633
1000	21.701169793959707	-461368.207073563	20.705871796946923	0.6324081353467753
1100	23.877596876962343	-461482.94301367877	21.59782831134988	0.6804547983325663
1200	26.046273167972686	-461196.07190656563	21.89529875580126	0.695351975793228
1300	28.20146191798267	-461581.75374825916	22.05797647498557	0.7022800214554449
1400	30.343181126983836	-461415.7796921233	22.15131349974166	0.7075633892084601
1500	32.502722001983784	-461419.7571235548	22.20611169905932	0.7100437766381722
1600	34.653629543987336	-461451.41976794874	22.243285163550297	0.7129767882990969
1700	36.76414687698707	-461358.69571003475	22.269733554172724	0.7145654958270448
1800	38.89133045999915	-461468.02711054875	22.29053732340955	0.7153006248654192
1900	41.00445283501176	-461471.51086336945	22.302629393744308	0.7143084066198495
2000	43.117388710030355	-461551.3816673975	22.314105859243234	0.7160416940397969
2100	45.25618321003276	-460776.03475141205	20.954106532528847	0.6378267364151339
2200	47.389108710049186	-461616.8578864591	20.698179671294298	0.6232929064419316
2300	49.50521750203916	-460849.94480730133	20.783907101253245	0.6362590030038515
2400	51.59733458503615	-461340.719244005	20.84683821598235	0.6418304291650572
2500	53.68617733503925	-460702.2744652891	20.64804326972326	0.6270398956831558
2600	55.80447171005653	-461543.3489333875	20.87364614197037	0.6448347575445377
2700	57.89097816805588	-460938.2183221819	20.823092840713187	0.6404984989843274
2800	59.991070001065964	-461601.446279374	20.876892628119208	0.6429010223436173
2900	62.09843408406596	-462063.3652945971	20.850453553284535	0.6418678044493948
3000	64.20349112607073	-461933.3975917351	20.847361158103244	0.6395210676695535
3100	66.31828466808656	-459659.56728269195	20.458863264811484	0.6115348679572872
3200	68.4254074601049	-459489.04598748975	19.914914736645635	0.5817828364836369
3300	70.54199137710384	-460772.89982956165	20.498095371031873	0.612402626074702
3400	72.69053458608687	-460802.34275263146	20.584735637276054	0.6176447599026642
3500	74.80276500308537	-460861.31976160425	20.705016237476364	0.6314018325413513
3600	76.94524450309109	-461932.851046751	21.568440605773766	0.6791966812035373
3700	79.10879846208263	-461973.8230159011	21.884764527642258	0.6931355199714218
3800	81.28451871208381	-461915.7062584949	22.056968039369096	0.7024661733798362
3900	83.46225454605883	-462037.9705417774	22.156093295047306	0.7071359325770814
4000	85.64281308805221	-461832.4826507634	22.218266177302247	0.710610200112095
