# 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.1, ρ̃₀ = 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_10_10000_9000", dynrange = 1.0, Λ = 1.0)
iter	time	function_value	psnr	ssim
1	7.291993824765086e-6	559998.953614435	-0.02636674332749045	0.0041734697881560985
2	0.028488707990618423	476684.2614162629	0.472776685077411	0.007393190719904383
3	0.04260420799255371	445201.0668526283	0.8540211734194507	0.010071288856104837
4	0.07170529098948464	399164.7156773093	1.1816343381859216	0.012069717808264213
5	0.08577174998936243	362075.92599906976	1.4924054250201109	0.014046724314610785
6	0.114293666993035	322731.28062713816	1.775842037312815	0.015100394823069997
7	0.12821533397072926	283627.39264009707	2.0835159093802877	0.015968523150916077
8	0.15672383399214596	245703.1173269409	2.398025040044546	0.01772582219707191
9	0.170821249979781	204994.36217082603	2.7250301301548956	0.01984448665596076
10	0.19918158298241906	166241.65344840795	3.047159370975063	0.021455002809732168
11	0.2130916239693761	129178.45129480983	3.3843701954098955	0.02333859391619731
12	0.24252858295221813	92190.42031359406	3.7175865032323996	0.02613786139749342
13	0.2564456659310963	56016.17486070832	4.060935003692993	0.02894865292149698
14	0.28496229092706926	20445.83049671856	4.406706896393983	0.03503060016232723
15	0.29965408294810914	-10820.539848524004	4.747035004401772	0.03580267150248734
16	0.32841912494041026	-42867.34069739799	5.100581166784598	0.04377450625718033
17	0.3424750829290133	-63467.280042452985	5.418211303063193	0.03842882592373254
18	0.3715502499253489	-105858.3205435936	5.815615648726185	0.06969672909321349
19	0.3856304999208078	-136279.4728838505	6.182687539054198	0.08855388825449675
20	0.4140370829263702	-165430.88657077836	6.550709492174072	0.11517468834708523
30	0.626232332928339	-350577.72715674265	10.253524197073357	0.44760542531021974
40	0.8376356239314191	-421464.7013663732	13.74745725372146	0.5360560430763901
50	1.0500765819451772	-448574.87274658383	16.541589427607462	0.5608994395049307
60	1.2621018319623545	-457418.322494328	18.28477094083371	0.5909416800508989
70	1.4744321229518391	-459471.6155837879	19.18505037053402	0.6097926747664869
80	1.685690789949149	-460625.18758905394	19.66050792728736	0.6228569197811515
90	1.897452747944044	-461037.5443990291	19.95577229624207	0.6280569007798111
100	2.108945622952888	-461037.5207565714	20.161057934645868	0.6310742150259324
110	2.3078746229584794	-460881.7701156385	20.28141537288178	0.6323770758874923
120	2.5200255399686284	-461456.70834743505	20.372638156867207	0.6309965005761459
130	2.7317270819912665	-461703.8242034038	20.30505330878801	0.619584016297692
140	2.943301249004435	-459728.37701272825	20.035420502280772	0.6001467123566366
150	3.1541914580157027	-461435.68642196327	20.48756438673483	0.6275228018008654
160	3.3668655000219587	-459638.80849275563	20.11527384764813	0.6066016564567359
170	3.5787867500039283	-461315.9426392966	20.62983860385259	0.6348387993768903
180	3.790702083002543	-460831.8463109983	20.50425346089	0.6218545770260777
190	4.003484333021333	-460435.2778542152	20.50291582184228	0.6228931522555905
200	4.215550333028659	-460376.55963348807	20.399248560738794	0.6156412709409311
300	6.324471042025834	-461532.8741688879	20.73996545384783	0.6337011675055899
400	8.428760209033499	-461130.2716086842	20.797400794223776	0.6348023915588052
500	10.535301209019963	-460950.615897312	20.774770297452587	0.6357214470329543
600	12.643214251031168	-461769.44049307716	20.87882560938038	0.6445277826958751
700	14.752102043043124	-461229.68424212997	20.87875988326106	0.6429535998025804
800	16.86368837603368	-461308.0360071368	20.83354574650479	0.6425115734229906
900	18.974047042021994	-461413.6071498565	20.92143595322821	0.6454438791167731
1000	21.075754917023005	-461427.9340122402	20.71148889818161	0.6326010398143926
1100	23.191856125049526	-461589.7919045007	21.754500860634813	0.6860226720897828
1200	25.301543500070693	-461310.2692298159	22.153301932102906	0.7048069711117542
1300	27.415782625059364	-461697.1773838522	22.389356779062716	0.7146134932110892
1400	29.53252366706147	-461531.86624387663	22.5446862423866	0.7226796158576442
1500	31.640135458059376	-461538.1497565673	22.648577381417777	0.7271443265226708
1600	33.7538428750413	-461569.9148373997	22.726057057842123	0.7319564490890298
1700	35.86513366701547	-461477.8169952978	22.78689098711219	0.7350630854441553
1800	37.97419250101666	-461589.0225341856	22.83404717689828	0.7370434833758722
1900	40.08935075101908	-461592.82003534556	22.867928299851947	0.7370530680661427
2000	42.20293112602667	-461671.304628238	22.898872440248077	0.739863872208026
2100	44.31266687603784	-460844.99766735715	21.002591493876587	0.6400089021705762
2200	46.42025608403492	-461676.1197119175	20.7127377895339	0.6239113282943877
2300	48.52856441703625	-460910.1361872444	20.793885647439275	0.6367373830296205
2400	50.63286970803165	-461401.42875568825	20.8547186753131	0.6421138474156618
2500	52.7411459160503	-460762.59986547916	20.654467323538505	0.6273314320795518
2600	54.84565312502673	-461604.54509966687	20.881244443314944	0.6452001073451017
2700	56.94960941703175	-460997.8186055902	20.83037839998733	0.6408870577834642
2800	59.05034845904447	-461662.8778691268	20.884527555373396	0.6432390171542983
2900	61.15820329205599	-462125.03635027097	20.857985036241185	0.6421555679117441
3000	63.26027033405262	-461994.49033973477	20.85406140414086	0.6398115785096935
3100	65.36187537506339	-459718.9778441281	20.46450719719097	0.6117481907128877
3200	67.47232058405643	-459547.7820764949	19.91813442670685	0.5819166675007728
3300	69.5906347920536	-460828.9564960108	20.503892102808233	0.6126426879805269
3400	71.70533354204963	-460857.5454388147	20.59068686949729	0.6179350327246754
3500	73.81569958405453	-460919.3641461849	20.712002186426005	0.6317390455407006
3600	75.93325733405072	-462040.63060978893	21.720829930773018	0.6843092471409415
3700	78.04634229204385	-462088.0623416647	22.13097096604516	0.7016704430643532
3800	80.15978345903568	-462030.4444034288	22.371578781944304	0.7137042283131156
3900	82.26893700004439	-462155.06194524525	22.52892320151577	0.7208378602463097
4000	84.37998716704897	-461949.50193975907	22.637395647085448	0.726296856645753
