# 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, save_prefix = "img/shepplogan256x256_pdps_known_rotation_20_10000_9000", dynrange = 1.0, Λ = 1.0)
iter	time	function_value	psnr	ssim
1	7.916009053587914e-6	555362.1431745938	-0.14938477322295188	0.005295187802576954
2	0.028921707009430975	470965.17176531116	0.3314460538266129	0.0078112081340624105
3	0.04338974901475012	439247.722726294	0.6990065571636511	0.01000280042917199
4	0.07238212402444333	393414.46786418976	1.0234748563592682	0.012188894853834335
5	0.08678862403030507	356906.48086817714	1.3253380988403776	0.014323816235725742
6	0.11576983201666735	318080.9296750219	1.6079186818975129	0.015113447028961991
7	0.13010712401592173	279921.6744559892	1.914997819459777	0.01712460792522995
8	0.1596317490329966	242682.6308663172	2.225452962423629	0.0195285316279422
9	0.17388183204457164	202376.45878938324	2.557121331359826	0.022139025685291726
10	0.2027538740367163	164385.4376645373	2.880424218627551	0.02413730817109596
11	0.2170544990513008	128118.68814002613	3.2184288538537342	0.02727091978986875
12	0.24574058305006474	91851.33001951037	3.554548498981464	0.03030492713417986
13	0.2600235410500318	56681.16955108162	3.9027757056978	0.034236322892530785
14	0.2888022070401348	21552.881138652316	4.255620029932976	0.03942275663241934
15	0.3030879150319379	-6933.842547822498	4.58867260343696	0.04236172394214785
16	0.3320818740176037	-38132.101143003	4.94972639684636	0.049097947019988195
17	0.3463603740383405	-54950.97610416152	5.2315032090791656	0.05055641217965871
18	0.37509304101695307	-100191.64330465607	5.692929102158221	0.0669595625808021
19	0.38937816602992825	-130986.15671658724	6.081354183297224	0.0804122998964257
20	0.41815712401876226	-160929.6735047109	6.476139558165176	0.09280670823502037
30	0.6326543740287889	-372453.4770328692	10.610832643010042	0.3541549342305579
40	0.8474082490429282	-458830.1150923402	14.33649197626317	0.5247926833124967
50	1.0618021650589071	-484512.2216579852	16.395580273378712	0.62145547658029
60	1.2760071230586618	-489110.8098082604	17.233446981962892	0.6860347846006175
70	1.4902914150734432	-489493.5456482296	17.648150486461812	0.7170528027234633
80	1.706353457091609	-489798.1483663213	17.966588151906507	0.7403406101920298
90	1.9217430820863228	-490472.61962010484	18.17812666849408	0.7553146457796567
100	2.1367785400652792	-490445.90931598045	18.321600948869254	0.7653702148591773
110	2.337928248074604	-490472.45087029884	18.356923153723706	0.7729860962995853
120	2.557809457066469	-490750.19311752776	18.442219497370218	0.7776725528732843
130	2.777898291067686	-491015.2652741599	18.15125769794784	0.7726343378075466
140	2.9987229580583517	-491360.9308131314	17.790461952142287	0.7567620589522677
150	3.2194779580750037	-491085.10261178087	18.423602716673837	0.7822435558244306
160	3.440022999071516	-490621.22843595606	17.804317761989388	0.7714531218223233
170	3.660403333051363	-490751.650734164	18.664923909834084	0.7891179230393812
180	3.881455208058469	-490415.39736765524	18.42854659897678	0.7884016700461268
190	4.102477624051971	-491255.52630533517	18.443951163711972	0.7891489170348913
200	4.3231108330364805	-490327.3844236469	18.199620536527735	0.7879765141075625
300	6.468131458037533	-490828.402269911	18.710463277395135	0.7949112069511671
400	8.601025000039954	-490657.30489844363	18.745774553782837	0.7985664743676746
500	10.755992167047225	-490365.9108786195	18.698391786549934	0.7984690923849465
600	12.953326209040824	-491266.5642689452	18.87257902914499	0.8065955316836935
700	15.119053418049589	-490798.2045174672	18.835028711036845	0.8041953521037447
800	17.29391883505741	-490597.3739339626	18.784284094216712	0.8042651162926697
900	19.499163210071856	-490885.490345645	18.88587254249189	0.8065952223968069
1000	21.71555000107037	-491411.34590388986	18.623623628095043	0.7980608430664407
1100	23.93544612606638	-491886.0974568922	20.649698025621746	0.8633466679064813
1200	26.143367084063357	-491705.7308923672	21.457087219661727	0.8870656136220734
1300	28.350397917063674	-492045.91125816474	21.950986809617042	0.8993776600995834
1400	30.56020525007625	-491975.9869600194	22.300140212272126	0.9066646474611746
1500	32.7633470830624	-491947.8893534837	22.579178901155288	0.912625234738903
1600	34.97749224904692	-492010.3194659754	22.800510793511574	0.917519336658416
1700	37.17977579103899	-491914.4379371507	22.990062958829952	0.9208068443475634
1800	39.384890083019855	-492098.5892206226	23.15146561439853	0.9235874718237778
1900	41.5874458750186	-492069.016144178	23.282656876182482	0.9243799038243842
2000	43.789207041001646	-492113.59831658355	23.40295639230332	0.9270847666335853
2100	46.02147124899784	-490787.5160109572	19.002162087136607	0.8105708554674146
2200	48.22677070699865	-492005.6208189675	18.618956551555733	0.7834133048603815
2300	50.40908970698365	-489757.3591864787	18.712343842199218	0.7988575973764007
2400	52.609418498992454	-490638.31456929236	18.795922590063544	0.8034116873510242
2500	54.79632858198602	-491467.1563346714	18.53591716641413	0.7963006043835451
2600	56.99749420696753	-491067.42941324477	18.890696220695837	0.8064537915675937
2700	59.18418478994863	-490776.55936448934	18.797008824983024	0.802561456361496
2800	61.35068958194461	-490885.96624323964	18.89782972661592	0.8055497272989023
2900	63.49295416593668	-491338.28397288435	18.84987507380553	0.8066656855976988
3000	65.61763249893556	-490832.6215952387	18.814497440948664	0.8029209311587783
3100	67.7472604579234	-492148.6608641416	18.26435420859942	0.7858244386848604
3200	69.8768968329241	-489847.1052282641	17.177693664862137	0.7552690595667536
3300	72.00020158290863	-491452.44934257993	18.297178784746684	0.7804132116008675
3400	74.13553483289434	-490537.3265450032	18.39632600749497	0.7794094026116936
3500	76.27203679090599	-491441.11791133235	18.663784111139805	0.7972560279756472
3600	78.41300745788612	-491618.7912805367	20.614136258939453	0.863039318462748
3700	80.5416302078811	-491965.7037865886	21.432945088394728	0.8856858203180422
3800	82.66252904088469	-491775.38126708806	21.930123786507828	0.8989322462459487
3900	84.80186437387601	-491965.0177946527	22.28129559383111	0.9060268886615894
4000	86.92848195787519	-491741.07056272787	22.54282081767927	0.913016873824646
