# 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/brainphantom256x256_pdps_known_rotation_15_10000_9000", dynrange = 1.0, Λ = 1.0)
iter	time	function_value	psnr	ssim
1	7.708003977313638e-6	561764.4291806474	-0.02636674332749045	0.0041734697881560985
2	0.029065873997751623	477783.63857486524	0.4729199726544222	0.007393204824260974
3	0.04343904001871124	445914.98053813865	0.8540988097744042	0.010070524861940346
4	0.07216749800136313	399703.4619365499	1.1818564235542763	0.012070432378445475
5	0.08648608098155819	362481.21718638134	1.492624401470657	0.014044332917057607
6	0.11520670598838478	323114.64572887955	1.7765480664531115	0.01512062115529108
7	0.12950920598814264	283941.3190742547	2.0843937745929946	0.016003783187425232
8	0.15845974799594842	245974.10581154193	2.3988362281864912	0.017755379561047058
9	0.17286991400760598	205235.7778453111	2.725802814004962	0.019880962453011814
10	0.20257887200568803	166475.59466548232	3.0482458139683204	0.02151455748385437
11	0.21704245600267313	129375.70167079523	3.3854079596284974	0.023383357174318956
12	0.24637187301414087	92382.18068978301	3.7188332581323715	0.02621359170145969
13	0.26094516401644796	56195.716466207334	4.062216160007936	0.029050188095151058
14	0.2899840390018653	20619.578748356867	4.408076607080682	0.03511821654186564
15	0.3045761640241835	-10652.608019288964	4.748345128163269	0.03587564789601703
16	0.3339359140372835	-42704.427068743425	5.1018828134553065	0.043841042512859404
17	0.34827058104565367	-63314.16952852689	5.419469715669077	0.03851713357065999
18	0.3770917060319334	-105729.86007062733	5.816955445957763	0.06984060310111218
19	0.3913970390567556	-136155.53084407336	6.184119677870833	0.08891154772924184
20	0.41999249704531394	-165308.35368757602	6.552438632074759	0.11563014980116891
30	0.6351549140526913	-350493.8980523829	10.256451908917413	0.4491381205116422
40	0.8495419140672311	-421390.53652183624	13.75407350682556	0.5365634706141134
50	1.0641407470684499	-448516.83831450564	16.552494127924174	0.5613678120318857
60	1.2790414970659185	-457373.2113617165	18.295142435802013	0.5915914930439161
70	1.4936183720710687	-459431.06127289287	19.192911704593293	0.6103702444914032
80	1.708056247094646	-460586.60831924144	19.665002042864327	0.6229839277103291
90	1.9218420381075703	-460998.52510959347	19.95722957104912	0.627725467941526
100	2.135208621097263	-461001.4847622098	20.159025981744236	0.6305442552993215
110	2.334023830102524	-460846.63972027064	20.274919563371306	0.6317170726951721
120	2.547699955117423	-461422.19252325303	20.363637799534146	0.630169791552601
130	2.7616720381192863	-461670.980807086	20.29627611540824	0.6189122235339211
140	2.9764231211156584	-459693.47982493864	20.025667213829884	0.5992023549045603
150	3.194136287114816	-461400.7603589639	20.472366957532984	0.6265291420142822
160	3.408284537115833	-459606.170090121	20.106604477370503	0.605851947530117
170	3.6218501201074105	-461279.33634610986	20.609913881824255	0.6336855070494156
180	3.8352556200989056	-460795.7837269131	20.483926607318303	0.6208375834573114
190	4.050199995108414	-460399.2758792732	20.480865810809252	0.6216443750755236
200	4.264386786118848	-460343.7159777495	20.383175489732274	0.6147825291268055
300	6.395652245118981	-461494.7951912815	20.71131949308186	0.6320103187738491
400	8.526247745117871	-461088.5839015691	20.757681718041006	0.6329576997087213
500	10.667058328108396	-460910.3257245551	20.735439924184377	0.6339836089823129
600	12.794829120131908	-461728.48459210614	20.840615854606984	0.6429158391955876
700	14.927600787108531	-461189.57051261695	20.843436601375508	0.6412483701447755
800	17.078017662104685	-461269.24692860455	20.800323139009752	0.6409788457271608
900	19.244590662099654	-461373.6717090735	20.88739133463733	0.6438276324653864
1000	21.372524245118257	-461390.6427800393	20.685447628383994	0.6311305290457396
1100	23.50148220313713	-461532.1470327497	21.664540533534968	0.6827871503581056
1200	25.62585007812595	-461249.7700741688	22.015195898971772	0.6998969685766705
1300	27.764985328132752	-461636.5348705611	22.214623409209324	0.7082991724546306
1400	29.893507828150177	-461471.14034994907	22.337992042686974	0.715162746152892
1500	32.02450824514381	-461476.16387205996	22.41506105799241	0.7186487980739497
1600	34.15395953715779	-461507.8774804134	22.46936002263268	0.7224757001135167
1700	36.27969191217562	-461415.5557213107	22.509223478607893	0.7247166776169618
1800	38.40836320316885	-461525.8222676012	22.538912727495003	0.725926844059773
1900	40.56684641216998	-461529.3404227838	22.558536750026455	0.7252983695602662
2000	42.71517262017005	-461608.54779140104	22.57649505964899	0.7274419930344935
2100	44.860992828180315	-460804.44635363447	20.958239742268475	0.6377353132730768
2200	46.98667645320529	-461636.4251494054	20.679737077682194	0.6221960451995789
2300	49.11004324519308	-460869.8979837401	20.758811303688613	0.6350125670557403
2400	51.24471899520722	-461362.23067702865	20.820333837778797	0.6404706189575793
2500	53.36901903722901	-460724.68512689235	20.623790695962633	0.6260346711724998
2600	55.494336829229724	-461562.8800748145	20.840919527330062	0.6434417313275931
2700	57.6836062882212	-460958.53812500084	20.793347770597602	0.639247448987441
2800	59.86122245521983	-461622.2314742859	20.848371664542974	0.6416139169628532
2900	62.01120199621073	-462083.7446326095	20.82045056101171	0.6404805087686212
3000	64.1428712462075	-461953.37896732794	20.81588684437946	0.6381163191111507
3100	66.2932096622244	-459681.44512010436	20.431892616959356	0.6099637172812262
3200	68.428635287215	-459514.4532973278	19.899426968526686	0.5810634035987677
3300	70.5758698702266	-460792.0444522406	20.472351883736316	0.6110193082522237
3400	72.7093916622398	-460819.31242350716	20.556212479685463	0.616206239893945
3500	74.83704391223728	-460880.89434794354	20.67825685645735	0.6302559650044816
3600	76.96735387123772	-461981.98968116625	21.630545544292445	0.6811376842976082
3700	79.11625083023682	-462027.4541997789	21.99931693741647	0.6971489301081195
3800	81.2435840392427	-461970.1063094816	22.205709919201052	0.7079239998889968
3900	83.38119324823492	-462093.75577076443	22.331285777492965	0.7138889587497761
4000	85.50739920625347	-461888.2925844729	22.41365700160263	0.7183127441821946
