# 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.25, ρ̃₀ = 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_25_10000_9000", dynrange = 1.0, Λ = 1.0)
iter	time	function_value	psnr	ssim
1	8.25000461190939e-6	557104.045973098	-0.14938477322295188	0.005295187802576954
2	0.027825208991998807	472032.55631312955	0.3312041982470264	0.007808015136724637
3	0.04130933398846537	439975.34049901913	0.698710679415212	0.009991751300485092
4	0.06960362600511871	394007.8572519722	1.0229076961508048	0.012165106308337944
5	0.08292183399316855	357427.2526143727	1.3247690829026135	0.014278659069774931
6	0.11079516800236888	318613.19779510604	1.6064739392716765	0.015015169788863799
7	0.12421116800396703	280413.57161490456	1.9131707416410153	0.016963578296895397
8	0.15206475200830027	243131.58848860345	2.2237779767987935	0.019341200192784295
9	0.16541175200836733	202784.64361420387	2.555543571772586	0.021898507889675387
10	0.19311108501278795	164774.40481589106	2.8783463795958024	0.023782940914571647
11	0.20648716803407297	128497.51980957978	3.2164942712579387	0.026848932049455036
12	0.23422908503562212	92199.7879964225	3.5523227601029745	0.02974436600379135
13	0.24766154401004314	57009.968974685966	3.900527671336635	0.03359027302190973
14	0.276375794026535	21857.100284763845	4.253260284438286	0.0386436848735255
15	0.2897641690215096	-6660.330053249684	4.586538554501873	0.041598722801951996
16	0.31765008601360023	-37881.97310471128	4.947726432137596	0.04832654693033259
17	0.3310372950218152	-54738.058230931754	5.229844554323778	0.04982974414391782
18	0.35865983602707274	-99993.47996343624	5.691026609144715	0.06630568930375105
19	0.3719451690267306	-130807.81113600837	6.079345083729724	0.0797018113720652
20	0.3997344610397704	-160766.92859874366	6.473688457809558	0.09200752652358714
30	0.6035979200387374	-372386.22021442483	10.60626727872252	0.35010810373406537
40	0.8067512120469473	-458793.27135393856	14.326132735901755	0.5181147271597553
50	1.009617545059882	-484459.6946676675	16.384844996171214	0.6119366101123619
60	1.21258758706972	-489074.94950663066	17.232033355312126	0.6784893935028062
70	1.4161872540717013	-489470.2273326881	17.652710205766486	0.7122212066038374
80	1.6198465460620355	-489782.9953867319	17.978985519208425	0.7368479408505025
90	1.8236814210831653	-490462.23346188315	18.196376929434763	0.7523853184192726
100	2.027906713075936	-490435.12695505214	18.341301117570726	0.7634223270688452
110	2.2167444210790563	-490461.4946484469	18.380104201218867	0.7717887512742079
120	2.4205890040902887	-490739.47945067723	18.46444934711335	0.777179179175655
130	2.624739754071925	-490999.847779897	18.17161127901947	0.7726550311686369
140	2.828652712079929	-491343.1088129584	17.802171790321026	0.7571092365873925
150	3.032371295063058	-491076.6105152501	18.45266198776274	0.7831681933544719
160	3.236122962058289	-490608.1575233154	17.825972345591165	0.7727788688400669
170	3.439890379057033	-490747.67871578137	18.70512032410371	0.7907928108047726
180	3.6434688370500226	-490413.5035164302	18.47116991989026	0.790458403939364
190	3.8476936700753868	-491251.6542599818	18.480513479629906	0.7912955341071015
200	4.05153342007543	-490319.46360084874	18.231869981426783	0.7903685570797601
300	6.081855462078238	-490822.16210825345	18.7524555580695	0.7978147860498512
400	8.122511878085788	-490656.1221309062	18.79892370679593	0.8018015087537103
500	10.156011211074656	-490360.57407659385	18.745675687498235	0.8012763614703408
600	12.23714462807402	-491264.0413028057	18.91835568138744	0.8094843438752218
700	14.276405711076222	-490791.74194827763	18.87397955269322	0.806794872099291
800	16.31397966906661	-490590.8382533909	18.824960427819963	0.8068693031266551
900	18.383862086076988	-490880.3335927989	18.926530688507125	0.8095033031981897
1000	20.446723836066667	-491401.07552724215	18.6550693091996	0.8008033522448385
1100	22.534541836066637	-491788.2055406168	20.54832370019883	0.8307233108375354
1200	24.60110846106545	-491524.9310839032	21.05085933900191	0.8217514568333076
1300	26.677505086059682	-491811.6529150594	21.262465781164085	0.8110656945929036
1400	28.783322586066788	-491705.60521569557	21.369532990503373	0.802033454649987
1500	30.877530253084842	-491652.8601364169	21.434589031495417	0.7960534487922497
1600	32.97039954509819	-491698.4506875924	21.475095593059685	0.7921918181188484
1700	35.04611450311495	-491590.9809307605	21.50312936806131	0.7894368759804411
1800	37.112186253099935	-491765.69135945477	21.52344465704055	0.7862002381290228
1900	39.195609837101074	-491727.3342129692	21.53799227912682	0.7831691148469722
2000	41.27421321210568	-491764.2345426281	21.54847334711805	0.7823426980235039
2100	43.34361275311676	-490759.6059225207	18.967809025672167	0.8087139188998447
2200	45.4125870441203	-491997.59679986915	18.6483721857109	0.7852088168054554
2300	47.476981961139245	-489750.89067404135	18.75597185961822	0.8014082528006686
2400	49.560746336122975	-490631.300019826	18.83510575586628	0.8062457686211124
2500	51.62614758609561	-491464.43260874826	18.572222030650718	0.7992855363812701
2600	53.709673795092385	-491065.7154281846	18.939235140345943	0.8097475721349302
2700	55.760243795084534	-490772.87515798624	18.839722369266646	0.8051772618068038
2800	57.80251121107722	-490881.92516855634	18.941806993306333	0.8087237103384954
2900	59.869123169075465	-491334.11271207366	18.89282813461027	0.809736148508977
3000	61.913092127098935	-490827.96297204	18.859556553059747	0.8058627825401694
3100	63.996678502095165	-492141.8710384871	18.3054759376084	0.7886008707364921
3200	66.03315066808136	-489827.40780479024	17.188528430883594	0.7574097912444878
3300	68.06209029309684	-491448.03543468774	18.341458837434924	0.7832237241577746
3400	70.10792346010567	-490533.4602178033	18.443732001746323	0.7818784505054792
3500	72.17749704408925	-491442.921524731	18.708593036185313	0.8003797123832219
3600	74.21204454411054	-491521.1928713146	20.547890074573402	0.8315100877054485
3700	76.2398600440938	-491780.9948699053	21.07354238941699	0.8213133681586493
3800	78.27360550209414	-491536.46981797874	21.293222048203052	0.8118991256845762
3900	80.30115916911745	-491692.612208493	21.4027285183851	0.8034071490537611
4000	82.34527358610649	-491444.8795078573	21.463080910831927	0.7981941031232542
