# 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, phantom_ρ = 100, save_prefix = "img/brainphantom256x256_pdps_known_proximal_20_10000_9000", dynrange = 1.0, Λ = 1.0)
iter	time	function_value	psnr	ssim
1	8.208997314795852e-6	563529.9047468598	-0.02636674332749045	0.0041734697881560985
2	0.02912275100243278	478839.64546806616	0.4746343772855547	0.007423191550118712
3	0.04398237599525601	446507.59743633814	0.8562108604554385	0.010136235965406852
4	0.07303554299869575	400167.6488076236	1.1840638534927825	0.012174050469588629
5	0.0880686680029612	362787.90545734484	1.49501333850777	0.01418542375275382
6	0.11725433502579108	323387.7239062619	1.7791597796146275	0.015275236260770611
7	0.13208383505116217	284195.00887967093	2.0868156476253907	0.016172256692559037
8	0.1610971680493094	246188.358290901	2.401476921664967	0.017963864139791758
9	0.1759393350512255	205436.58494332607	2.7286705498933532	0.02011526934129273
10	0.20559529305319302	166659.88858756062	3.0511762614112206	0.021752157646948968
11	0.2204027520492673	129556.41767937598	3.3886482584700484	0.023618208278679942
12	0.24941071003559045	92528.21914575665	3.722123873649954	0.026476223626121717
13	0.26421633502468467	56321.2701802887	4.065838015311937	0.029431835707516195
14	0.2932384190207813	20737.141250331388	4.411974773226	0.03555297734391176
15	0.3080567940196488	-10564.367461150161	4.752602212503581	0.03643909160300724
16	0.33783991902600974	-42647.30284100188	5.106556560958817	0.04444155231148236
17	0.35279108604299836	-63247.55307394409	5.4247361603015225	0.039112474808827284
18	0.38256166904466227	-105730.115467401	5.822348551451243	0.07107249226381596
19	0.39735491905594245	-136180.05411071982	6.189864661850839	0.09056709103814087
20	0.42635350304772146	-165369.44393974353	6.558294642146956	0.11726548811289159
30	0.6494821700325701	-350547.1526992722	10.264970781056025	0.4528999493197329
40	0.8708067120460328	-421343.5488377603	13.767715503548425	0.5370094185195955
50	1.0888759200461209	-448465.8042954619	16.572504762319998	0.5616250844899766
60	1.308580878045177	-457331.9084540178	18.313780398907674	0.5916182698857849
70	1.5262017530330922	-459396.46697578614	19.204691013743297	0.6096215274885372
80	1.7462172110390384	-460551.04373003234	19.66468336161593	0.62168897140696
90	1.9678123770281672	-460962.7453078999	19.945031330167488	0.6260957854721657
100	2.184723960031988	-460960.31864081806	20.133533827790423	0.6282699193730792
110	2.3947401260375045	-460804.1847973632	20.238962241936605	0.6292606279758076
120	2.617019710043678	-461378.64651609975	20.318721289669757	0.6276657769982238
130	2.8387197520351037	-461625.5584382656	20.2597131283401	0.6177231064203069
140	3.059896502032643	-459669.1011843126	20.00025015158462	0.5986654300797203
150	3.2768834180315025	-461353.87835009943	20.411624682051773	0.6239649945040456
160	3.494232335011475	-459566.38551192574	20.074082132534933	0.6056372699252903
170	3.7177616680273786	-461226.97954394517	20.52336287964841	0.6302416561646178
180	3.943570751027437	-460751.2557117686	20.402176667116315	0.6180891296163313
190	4.169205459038494	-460351.7280517628	20.3990335110627	0.6189357919991882
200	4.391255667025689	-460299.275337602	20.301991662419713	0.6122425329874197
300	6.554483459010953	-461431.8487544805	20.59407042691994	0.6282824912400324
400	8.736519250989659	-461029.3174310733	20.6335547290286	0.6288498433295564
500	10.958460417983588	-460850.637916787	20.610867563773645	0.6297752339450142
600	13.125950833957177	-461665.7863182212	20.70699242632582	0.6384892909310917
700	15.310857041942654	-461124.84484495263	20.70760218421437	0.6363379716311948
800	17.479971458931686	-461205.43928826554	20.67284589796	0.6366866648054528
900	19.653682917938568	-461309.79759255535	20.74989489491594	0.6390993715935089
1000	21.834239375952166	-461329.7766585315	20.567721906087332	0.6276404595809781
1100	24.005956458946457	-461474.983236017	21.51979659251216	0.6782025065153683
1200	26.173388083960162	-461193.75736870914	21.8514122001571	0.6939794036717877
1300	28.337121666962048	-461581.28591208084	22.03176996675282	0.7014336887725795
1400	30.52392837498337	-461416.48731623293	22.134454019977788	0.7070792591429622
1500	32.6990929999738	-461420.7851887064	22.194273230800633	0.7097229184164006
1600	34.88427820795914	-461452.393956365	22.233484670366753	0.7126639007937632
1700	37.065141707949806	-461359.65825446224	22.26266901692647	0.7143938360497855
1800	39.24909241593559	-461469.3388341785	22.286386720557203	0.7152950913531786
1900	41.4882467489515	-461472.7880771026	22.299625633679476	0.7144509495856519
2000	43.696011165942764	-461552.97970494867	22.311809774604715	0.716286368234632
2100	45.93157737393631	-460750.4497438295	20.827306301276085	0.6316830613099522
2200	48.1558433319442	-461577.67067991494	20.5629518172362	0.6180559946091699
2300	50.35970953994547	-460808.7020835183	20.635135826130323	0.6309062315403804
2400	52.60181466495851	-461297.7033023953	20.68522044602245	0.6359381909954054
2500	54.818972498964285	-460669.87516387884	20.51825934752066	0.6229951400502556
2600	57.053254790982464	-461501.36223233416	20.71032906112986	0.6390855931289572
2700	59.22159833298065	-460896.20895959245	20.666712043578187	0.6352132096314398
2800	61.44473183300579	-461559.2112633702	20.711140605293096	0.6368506127745588
2900	63.63291220800602	-462020.3009242911	20.68842811740086	0.6361190618468229
3000	65.80046170798596	-461891.55403813574	20.688323264238306	0.63348778771919
3100	67.9885538329836	-459635.61765110045	20.339523078714286	0.6071126566107001
3200	70.16799924997031	-459475.5757203825	19.854773635122264	0.5810923913122348
3300	72.33179116697283	-460735.73891212116	20.365498213880674	0.6083106151157398
3400	74.51854791698861	-460760.54374976293	20.44103467344298	0.6128114667358785
3500	76.73295445900294	-460824.3060911457	20.556684365410515	0.6264640540890752
3600	78.94666833401425	-461924.8925115768	21.491753415579446	0.6768808362557507
3700	81.10614875002648	-461971.6678571752	21.84231977498073	0.6917467933664176
3800	83.30102754200925	-461915.9762844585	22.030784364676936	0.7016760988153131
3900	85.46943495800951	-462038.296150148	22.13882784518422	0.7067336228020327
4000	87.62698287400417	-461833.0794348696	22.206571364081306	0.7103394303436182
