# 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, phantom_ρ = 100, save_prefix = "img/shepplogan256x256_pdps_known_proximal_15_10000_9000", dynrange = 1.0, Λ = 1.0)
iter	time	function_value	psnr	ssim
1	8.041010005399585e-6	553620.2403760896	-0.14938477322295188	0.005295187802576954
2	0.03069729101844132	469886.6434678431	0.3325667942232765	0.007825150603900979
3	0.04551933202310465	438452.399246117	0.7003780239057795	0.01003259629975117
4	0.07478524901671335	392823.2576830987	1.0247865392072957	0.012237505258307507
5	0.08965062402421609	356406.9427192285	1.3267579181986042	0.014396657211476215
6	0.1186729570326861	317596.36829379085	1.6091304376846438	0.015181043796362167
7	0.1335239980253391	279532.1235887144	1.9159083124560805	0.017192393487540363
8	0.16267874801997095	242343.3008720884	2.22657557850146	0.019656056812340193
9	0.1774780810228549	202087.44595524625	2.558430766877264	0.022324033820729387
10	0.20645491502364166	164111.74692227016	2.881570255112655	0.024328275395201326
11	0.221251374023268	127891.82081332077	3.219819708158233	0.027549763850388202
12	0.25031370800570585	91626.74176962694	3.5558517862769916	0.030607847773516862
13	0.265114748995984	56469.39019681547	3.9042662222934488	0.03463410148700637
14	0.2944792069902178	21352.944719634957	4.257229213538639	0.03990966550249596
15	0.3093187069753185	-7135.055572150569	4.590617777389773	0.04307798443012988
16	0.3386613319744356	-38332.99172298403	4.952007449211173	0.05012495263898051
17	0.3536556239705533	-55145.80818194775	5.234341212762086	0.05167620433098248
18	0.38503833196591586	-100355.95009261483	5.695784985806601	0.06832894048073922
19	0.40068291494389996	-131158.08215612752	6.0844814799483835	0.08207336623436605
20	0.43066387294675224	-161122.7629007586	6.4792637856803195	0.09445023016776465
30	0.6564781639317516	-372655.4849988004	10.616696510497539	0.3583403380158988
40	0.8901153299375437	-458940.02552637935	14.345774349906122	0.528884106185903
50	1.1265840379346628	-484584.32377259724	16.405992701051296	0.6268745188812505
60	1.365575037925737	-489164.64922545594	17.243250543289147	0.6935134515562562
70	1.604953454923816	-489543.85793579096	17.656282658129506	0.7261412206984049
80	1.8444753719086293	-489846.0178380555	17.972689900921722	0.7497163747817716
90	2.0805412049230654	-490518.853288945	18.182478589832275	0.7641489002165422
100	2.307423038932029	-490486.946877934	18.31951620062909	0.7735246224067344
110	2.527256080938969	-490508.79500259116	18.35230427637383	0.7803013212918248
120	2.7647909969382454	-490786.37982319965	18.434112435764465	0.7843367887828216
130	3.0084882049413864	-491050.15233541146	18.149351396554263	0.7788924918569573
140	3.2481242469511926	-491408.5128021545	17.789899431495318	0.762896403289868
150	3.4865297889336944	-491119.67354952585	18.416650963946903	0.7866089369373241
160	3.7222978309437167	-490660.98128085927	17.807346350945576	0.7758036829421631
170	3.9591228729404975	-490781.4041493002	18.648003473716734	0.7920804708068319
180	4.196134956931928	-490443.560222762	18.40865996269063	0.7905649813812416
190	4.432954706921009	-491286.72366598673	18.427090149018927	0.7912083203371496
200	4.6693691239343025	-490355.05629778164	18.183193810069312	0.7897510038458442
300	6.947048665926559	-490849.89196449466	18.685572652813576	0.7955141346308495
400	9.274081915937131	-490684.9151824536	18.724885400404176	0.7995592940197578
500	11.543492831930052	-490391.23816837324	18.677831476709695	0.7991257207569236
600	13.818848956929287	-491292.2028674344	18.847793031601128	0.8070094082886927
700	16.046179623925127	-490821.0413906906	18.80757101604938	0.8046476130023548
800	18.315781081910245	-490622.9152277328	18.765082097416176	0.8048718595724964
900	20.56411433190806	-490909.6114896431	18.859481383307415	0.8073274491007594
1000	22.81077899891534	-491436.99147829035	18.602570631192414	0.7989310027301462
1100	24.969960665912367	-491943.4090950333	20.633867316258897	0.8644306227504361
1200	27.126065582910087	-491765.2531463949	21.42785531764048	0.8880819840862505
1300	29.289119873923482	-492105.789139024	21.912754376145234	0.9005450781767975
1400	31.445614457916236	-492035.7402961455	22.26127043641575	0.9080419569948897
1500	33.598773873905884	-492007.01014658605	22.54372482020749	0.9142097488947939
1600	35.75794012390543	-492069.2079078603	22.775017884054712	0.919246108330359
1700	37.90856304092449	-491973.13335274096	22.976039635490647	0.9227267734373354
1800	40.07270979092573	-492158.2003985968	23.151465783570117	0.9256890165813207
1900	42.22690699892701	-492128.89406789024	23.300572290114253	0.9266638474435216
2000	44.4074779159273	-492173.1323091511	23.436351564359917	0.9295870815230814
2100	46.57178795692744	-490816.1538875897	18.986187702482553	0.8119571543758224
2200	48.76058779092273	-492034.5738376045	18.592946098985543	0.7846447190115298
2300	50.92433741592686	-489784.29328177724	18.690933580975774	0.7996557334536801
2400	53.086628873919835	-490661.0977429429	18.77090992728449	0.8042398146895393
2500	55.244362915924285	-491496.62115471764	18.512318211434653	0.79716902392803
2600	57.40704845791333	-491093.3420714276	18.868900558367187	0.8075833074186625
2700	59.577090457925806	-490800.22291179147	18.771518668606415	0.8028138829416416
2800	61.74316166591598	-490912.1701690984	18.87267628772809	0.8064162491593027
2900	63.90556595791713	-491365.11996370833	18.825387716368926	0.8075229604420743
3000	66.06402624890325	-490859.1921361896	18.793445229240056	0.8038283588004382
3100	68.22591366490815	-492181.57488757436	18.258044416875112	0.7870093002629274
3200	70.42490649889805	-489879.38343074237	17.173957648786086	0.7563465706289124
3300	72.62534766588942	-491479.85327734065	18.282804191900496	0.7812714137605374
3400	74.78620091587072	-490559.5501416354	18.373131347910387	0.7799485514094795
3500	76.9623342078703	-491467.5324558249	18.6356367323254	0.7976755197463377
3600	79.13639554087422	-491677.1777941653	20.60537075562256	0.8640868971038267
3700	81.29112604088732	-492026.75153306033	21.41003010313372	0.8865898675892729
3800	83.44647604090278	-491836.22383070586	21.897737267830287	0.8999137214630423
3900	85.61152016589767	-492025.5594304397	22.24837219997315	0.9073386405744372
4000	87.76990991589264	-491800.9788958434	22.51729267403548	0.9145295028127571
