# 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/shepplogan256x256_pdps_known_noprediction_10_10000_9000", dynrange = 1.0, Λ = 1.0)
iter	time	function_value	psnr	ssim
1	1.2707983842119575e-5	551878.3375775854	-0.14938477322295188	0.005295187802576954
2	0.027455207979073748	511762.6762748136	0.09035124024425467	0.00519645817930506
3	0.040376540971919894	471905.12479492085	0.34074803598632303	0.0053424268543263344
4	0.06795433297520503	431561.83623518044	0.5967932959428549	0.005634381031389109
5	0.08084304098156281	392250.2538752461	0.816085117101838	0.00017454703660101464
6	0.10820512496866286	353395.3381444892	1.0723675209370906	0.0005594059015100963
7	0.12119108397746459	315382.67176386266	1.3376550145347055	0.000506423991266807
8	0.14868025097530335	277011.2142694752	1.6020717501575468	7.669629614094143e-5
9	0.16157408399158157	238782.42988323266	1.8996417748100838	0.0020841813045841644
10	0.1926674999995157	200673.84805247898	2.1823715624375124	0.0023072270285181126
11	0.20556441700318828	163796.81446780355	2.492319706825214	0.0034183215694146615
12	0.23274870900786482	127473.17183506374	2.786828382757607	0.004250110152529438
13	0.24562370899366215	91862.43047016693	3.0925617696754677	0.004532128887926853
14	0.2728040419751778	57525.40266063641	3.3869186490476206	0.004352726586196695
15	0.2858130009553861	23625.243317812627	3.7156604664606956	0.007201192596195009
16	0.313043250964256	-8528.15065962959	4.015455227778796	0.005187996518973011
17	0.32597925094887614	-38052.9506865902	4.231216063612667	0.0016292425475554623
18	0.3532711259322241	-70113.46012010332	4.663000634673421	0.00889241688230793
19	0.36620629293611273	-99694.49542750046	4.9820279440978785	0.008470538995175707
20	0.39369291794719175	-128142.39268057625	5.304982170417115	0.008799890284916043
30	0.5951338339655194	-340440.31171963277	8.387290952333986	0.04994023498655416
40	0.7977614169940352	-443061.83539941855	11.211063956251863	0.18153279511723275
50	0.998638832999859	-475476.69896889996	12.623801085069356	0.28561609430531865
60	1.2009241249761544	-484216.1414172726	13.665687098331444	0.3405936278575281
70	1.402006790973246	-484023.8835502245	13.725848126654256	0.35986891926842635
80	1.6028052489855327	-483893.03476922476	13.857422447087842	0.36388362458501855
90	1.8045765819842927	-487365.7264359214	14.74110179518759	0.3970545488510735
100	2.005351414962206	-487562.9971583935	14.964819959237392	0.42223834999309273
110	2.191947414947208	-481443.8061153458	13.779288076028873	0.3497312136374916
120	2.394672123948112	-483094.9062181984	14.215634374980732	0.35005203861815537
130	2.5955835409404244	-481465.689450985	14.013131970317136	0.3375151773531573
140	2.7959135409328155	-481721.50672382663	14.04974404502541	0.3341772172884403
150	2.9977033329196274	-485613.2367119625	14.876566668339443	0.38295964800148147
160	3.198434624937363	-487307.35702596325	15.17209916024407	0.4023929961729698
170	3.4002654999494553	-486692.83544183755	15.427045990094317	0.4243516782186473
180	3.6008963749336544	-486969.8788087038	15.623435899034705	0.4353107369267369
190	3.8017001249536406	-488046.9048776141	15.59673553606181	0.45470731173279266
200	4.0031875409476925	-488214.92731904914	15.906895449131147	0.4662133816105871
300	5.9992274579708464	-490141.55533056054	16.84130898584539	0.5989427527021232
400	7.995635582978139	-488901.5181651319	16.724408701908196	0.6681364265271069
500	9.986792624986265	-485444.78181832517	15.517815995237704	0.6413251875899091
600	11.980554375011707	-489268.42541724845	16.285093994102176	0.6531164138775719
700	13.975758875021711	-488633.5457809607	16.153627084556675	0.6468383624872704
800	15.969611458014697	-483590.43240738765	15.159762019708067	0.5849980015915978
900	17.963976875005756	-488257.49426574644	16.33443053783647	0.6729532542258758
1000	19.957228500017663	-490416.03401899495	16.814052425184073	0.7056430231861903
1100	21.94559370802017	-491852.77746538405	19.025647566825043	0.8034208528941347
1200	23.93479854101315	-491773.0237273761	19.936262100049017	0.8450551132429504
1300	25.927199666009983	-492136.17451705417	20.517600149665217	0.8662901169171147
1400	27.916584457998397	-492073.6057824432	20.958346977348967	0.8783150471962636
1500	29.90590645800694	-492048.56342839095	21.3201432721848	0.8874934898861007
1600	31.892764416028513	-492114.40826147486	21.62503347760964	0.8950477396773366
1700	33.88125370803755	-492020.3215959467	21.895793740377336	0.9009700390224764
1800	35.869459249050124	-492207.323388348	22.141744160849374	0.9060577607019956
1900	37.85484162403736	-492179.91908458516	22.365046798161508	0.9094193916209306
2000	39.8423195410287	-492223.53833846207	22.57161065820302	0.9141879637323398
2100	41.835085332015296	-486028.81118067115	15.205581539690233	0.6211119514475685
2200	43.8299174150161	-490969.46455074655	16.620091522138015	0.627249327204551
2300	45.82600674801506	-487657.40786587296	16.234350081000425	0.6518188773476752
2400	47.82188999801292	-488461.62944889977	16.075949624375536	0.6197253002974803
2500	49.8166332480032	-489734.7693556024	16.211518160882015	0.6555179571788028
2600	51.80703287301003	-486545.4182200225	15.547835555641454	0.6279845523114893
2700	53.811648706003325	-489056.1544667643	17.127782772692953	0.6633951773011924
2800	55.808713915001135	-487604.2827009419	16.039243029886165	0.662436688515945
2900	57.80081191501813	-484523.6059504403	15.166428128076614	0.6050800982934144
3000	59.7950489979994	-484141.30962730356	15.052430674491873	0.606184106868368
3100	61.794447205989854	-491323.07830380934	16.5318060178763	0.6556265598732326
3200	63.790023413981544	-489427.46878631035	16.222957283029366	0.6549220520064809
3300	65.78799203896779	-489883.7858988657	16.24462644069596	0.6390928495134239
3400	67.78219908097526	-481532.5647313757	14.796787045934447	0.5669650409822241
3500	69.77956091499072	-489727.333219283	16.250409590224013	0.6491836398145916
3600	71.77058324799873	-491551.38788407156	18.5447821523188	0.7661468091941731
3700	73.76161049801158	-492019.6047552764	19.460353255658344	0.816703661953789
3800	75.74982320601703	-491856.61371740437	20.04769880100168	0.8459420273585645
3900	77.73607729002833	-492058.48975168914	20.49741633540013	0.8644412333719713
4000	79.72538299803273	-491839.4828972337	20.867411436456944	0.8795773354700132
