# 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/shepplogan256x256_pdps_known_zerodual_15_10000_9000", dynrange = 1.0, Λ = 1.0)
iter	time	function_value	psnr	ssim
1	7.624999852851033e-6	553620.2403760896	-0.14938477322295188	0.005295187802576954
2	0.02942737500416115	469899.7479824107	0.3312041982470264	0.007808015136724637
3	0.042885333008598536	438549.7730106998	0.698710679415212	0.009991751300485092
4	0.07075416599400342	392909.9688060973	1.0229076961508048	0.012165106308337944
5	0.0843872909899801	356548.02595477295	1.3247690829026135	0.014278659069774931
6	0.11219974898267537	317763.9767181882	1.6064739392716765	0.015015169788863799
7	0.12559458299074322	279683.3860222197	1.9131707416410153	0.016963578296895397
8	0.15342070799670182	242494.13902908965	2.2237779767987935	0.019341200192784295
9	0.16673554098815657	202220.8121917766	2.555543571772586	0.021898507889675387
10	0.1944517910014838	164242.9172240287	2.8783463795958024	0.023782940914571647
11	0.20781275001354516	128018.00095327853	3.2164942712579387	0.026848932049455036
12	0.23537595802918077	91753.96118939866	3.5523227601029745	0.02974436600379135
13	0.24872858301387168	56594.06746394796	3.900527671336635	0.03359027302190973
14	0.27645679100533016	21464.406994873985	4.253260284438286	0.0386436848735255
15	0.29060920802294277	-7028.47668920877	4.586538554501873	0.041598722801951996
16	0.31833320803707466	-38229.0781658136	4.947726432137596	0.04832654693033259
17	0.3316813740239013	-55070.525776004804	5.229844554323778	0.04982974414391782
18	0.3594583740050439	-100303.40744750749	5.691026609144715	0.06630568930375105
19	0.37276504101464525	-131107.26799339868	6.079345083729724	0.0797018113720652
20	0.40055237399064936	-161060.33193163955	6.473688457809558	0.09200752652358714
30	0.6064442490169313	-372584.5382026095	10.60626727872252	0.35010810373406537
40	0.8118365410191473	-458942.78423151973	14.326132735901755	0.5181147271597553
50	1.0167014580219984	-484575.7244275047	16.384844996171214	0.6119366101123619
60	1.2215233330207411	-489167.56967221183	17.232033355312126	0.6784893935028062
70	1.426127291022567	-489552.6904533368	17.652710205766486	0.7122212066038374
80	1.6307185830373783	-489860.3033033049	17.978985519208425	0.7368479408505025
90	1.8353290830273181	-490536.7676690472	18.196376929434763	0.7523853184192726
100	2.039830375026213	-490507.6350876701	18.341301117570726	0.7634223270688452
110	2.2313825840246864	-490531.8041237684	18.380104201218867	0.7717887512742079
120	2.4358892510354053	-490808.97256387543	18.46444934711335	0.777179179175655
130	2.6414862090314273	-491068.6066394439	18.17161127901947	0.7726550311686369
140	2.846199667052133	-491413.10319973924	17.802171790321026	0.7571092365873925
150	3.050683959037997	-491142.88125271857	18.45266198776274	0.7831681933544719
160	3.2553627930465154	-490675.0499261058	17.825972345591165	0.7727788688400669
170	3.4599318760447204	-490813.400978039	18.70512032410371	0.7907928108047726
180	3.6648374590440653	-490479.63442853355	18.47116991989026	0.790458403939364
190	3.869611334055662	-491317.9945099525	18.480513479629906	0.7912955341071015
200	4.074904626060743	-490385.64177324454	18.231869981426783	0.7903685570797601
300	6.099365043075522	-490886.4053594548	18.7524555580695	0.7978147860498512
400	8.127935084077762	-490722.57435009064	18.79892370679593	0.8018015087537103
500	10.152374292054446	-490426.22326843656	18.745675687498235	0.8012763614703408
600	12.214421126060188	-491330.67152543	18.91835568138744	0.8094843438752218
700	14.2376894170593	-490859.0737796257	18.87397955269322	0.806794872099291
800	16.267143000033684	-490658.1237398772	18.824960427819963	0.8068693031266551
900	18.323979375039926	-490947.4622502024	18.926530688507125	0.8095033031981897
1000	20.34873104203143	-491467.72312072385	18.6550693091996	0.8008033522448385
1100	22.378781292034546	-491932.99937610974	20.54832370019883	0.8307233108375354
1200	24.4092865420389	-491705.0290614756	21.05085933900191	0.8217514568333076
1300	26.44045487503172	-492012.43177625106	21.262465781164085	0.8110656945929036
1400	28.479363500053296	-491919.8444682135	21.369532990503373	0.802033454649987
1500	30.51917829204467	-491875.75777684164	21.434589031495417	0.7960534487922497
1600	32.55334804204176	-491926.7925433009	21.475095593059685	0.7921918181188484
1700	34.61764129201765	-491823.4225520114	21.50312936806131	0.7894368759804411
1800	36.649005667015444	-492002.1883006002	21.52344465704055	0.7862002381290228
1900	38.70444350002799	-491967.1751713612	21.53799227912682	0.7831691148469722
2000	40.75979050001479	-492005.7709270947	21.54847334711805	0.7823426980235039
2100	42.789726292015985	-490828.1435625837	18.967809025672167	0.8087139188998447
2200	44.8222225420177	-492062.5588837857	18.6483721857109	0.7852088168054554
2300	46.85406450001756	-489817.2469090344	18.75597185961822	0.8014082528006686
2400	48.885665166017134	-490698.37238350714	18.83510575586628	0.8062457686211124
2500	50.90876212401781	-491531.7500327265	18.572222030650718	0.7992855363812701
2600	52.938994290016126	-491132.8423503824	18.939235140345943	0.8097475721349302
2700	54.967459665000206	-490839.15821827715	18.839722369266646	0.8051772618068038
2800	57.00081758201122	-490949.4581836679	18.941806993306333	0.8087237103384954
2900	59.027221707016	-491401.72262651287	18.89282813461027	0.809736148508977
3000	61.05666304103215	-490895.4008628537	18.859556553059747	0.8058627825401694
3100	63.09268562501529	-492208.35946609575	18.3054759376084	0.7886008707364921
3200	65.12069870901178	-489895.0072907784	17.188528430883594	0.7574097912444878
3300	67.14783741702558	-491512.0125533371	18.341458837434924	0.7832237241577746
3400	69.18350991702755	-490595.103092527	18.443732001746323	0.7818784505054792
3500	71.22357466703397	-491507.91736736655	18.708593036185313	0.8003797123832219
3600	73.25976937502855	-491667.740716191	20.547890074573402	0.8315100877054485
3700	75.29769454203779	-491964.79100958916	21.07354238941699	0.8213133681586493
3800	77.32794487505453	-491739.97731807263	21.293222048203052	0.8118991256845762
3900	79.35899608404725	-491909.591070261	21.4027285183851	0.8034071490537611
4000	81.40973395903711	-491669.6773881238	21.463080910831927	0.7981941031232542
