# 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/brainphantom256x256_pdps_known_rotation_10_10000_9000", dynrange = 1.0, Λ = 1.0)
iter	time	function_value	psnr	ssim
1	8.625007467344403e-6	559998.953614435	-0.02636674332749045	0.0041734697881560985
2	0.029205250000813976	476706.8064960277	0.47286048791687957	0.00739297165771168
3	0.043706333002774045	445203.3405791022	0.8540364686831491	0.010070180367036986
4	0.07282112500979565	399163.0669032806	1.1817254281051257	0.012068999525278803
5	0.0873190000129398	362058.3368102178	1.49250034342434	0.014043278627918097
6	0.11735020799096674	322710.5628766621	1.7761971209045668	0.015108973656727314
7	0.13202779099810869	283601.8085081453	2.0839431311672763	0.015983264844487487
8	0.1614519580034539	245679.99003893384	2.3984204752976277	0.017736447856043905
9	0.17620416701538488	204973.50929269247	2.725412001924101	0.019860106037138697
10	0.20595475001027808	166225.60631888625	3.0477185518861267	0.02148279581579057
11	0.2207742079917807	129153.73220243535	3.3849094813120093	0.023357914881646524
12	0.25027366596623324	92172.71742247358	3.718244883069668	0.02617431492527026
13	0.2650094999698922	55999.0208627043	4.0616044816497885	0.029001072430792158
14	0.2948941249633208	20431.266293317258	4.407401061758984	0.035062014393818154
15	0.3098444169445429	-10833.428237771583	4.747703735486258	0.03581992909592471
16	0.33983829195494764	-42877.15038779787	5.101246981990047	0.043786083242814244
17	0.35429041692987084	-63482.284089793626	5.418857309457437	0.03846287815574126
18	0.3833512919372879	-105878.95302921956	5.816308013931575	0.06975106079889096
19	0.3977822089509573	-136298.1417621497	6.183420026338128	0.08875268124775926
20	0.42682275094557554	-165445.15516718634	6.5516159162296015	0.11539272201582562
30	0.6448473349446431	-350579.7620490491	10.255119173479269	0.4484000806369973
40	0.8619290429342072	-421459.89932462934	13.751078893560315	0.5363405021364189
50	1.0795747929078061	-448573.00830770266	16.54744055818531	0.5611350895584319
60	1.299905458901776	-457417.9282269649	18.290135090915875	0.5912539750493344
70	1.5214442088908982	-459470.4459948057	19.188887470249185	0.6100839545545298
80	1.742952625878388	-460622.7644120779	19.66237942800874	0.6228757147196147
90	1.96570370887639	-461033.9202236073	19.95597674254622	0.6278533795692334
100	2.183249125868315	-461035.0081064004	20.15901873388738	0.6307282619572929
110	2.3852095428446773	-460878.8018032472	20.276664468835644	0.6319195301209338
120	2.603826584847411	-461453.8558201316	20.366079087059575	0.6304125274619763
130	2.823075500840787	-461701.9301793075	20.299231752232533	0.6190931317244732
140	3.042057708837092	-459724.4349904727	20.028546190131514	0.5994253887038394
150	3.2609814588504378	-461431.7394871927	20.477384768014943	0.6268568913861632
160	3.480466125853127	-459636.7362726908	20.110468211948792	0.6061122876958349
170	3.700820416852366	-461310.7799175826	20.61672194993662	0.6340877969444876
180	3.9200790418544784	-460827.1457137977	20.49115885169407	0.6211955193818144
190	4.138896125863539	-460430.5836817255	20.48861910169121	0.6220660795247079
200	4.358730125881266	-460374.1910190719	20.389733908006885	0.6151416457221965
300	6.520067250909051	-461526.785213435	20.722066738996432	0.6326466319292527
400	8.701478250906803	-461122.33152854495	20.77244555155247	0.6336685217848903
500	10.85756058388506	-460943.1760099906	20.750432211206444	0.6346811096803061
600	13.004204917902825	-461762.20294889045	20.85524717246935	0.6435414977458178
700	15.151925625890726	-461223.41624148935	20.85736365197798	0.6419374328342908
800	17.31893937586574	-461302.19209542737	20.81312992345188	0.6415422421341048
900	19.45878462586552	-461407.2530608611	20.900440964642073	0.6444126183169897
1000	21.62027804288664	-461423.07405640424	20.695685404083534	0.6316547190755759
1100	23.778667459904682	-461589.011789816	21.745745933882823	0.685432981561493
1200	25.950686125899665	-461309.98595018405	22.14745832710594	0.7043914077060278
1300	28.100163166876882	-461697.0315436049	22.385231080543242	0.7143042483565654
1400	30.242401416879147	-461531.78368264076	22.541832116474282	0.7224424588975814
1500	32.39270162486355	-461538.10690072767	22.64645521812183	0.7269519237694728
1600	34.537707915878855	-461569.87442982866	22.724441163024373	0.7318003160519719
1700	36.727136998903006	-461477.7860838435	22.78565862055752	0.7349376500933961
1800	38.888510414923076	-461589.0032767027	22.833131526955825	0.7369447257503956
1900	41.06528037291719	-461592.8019551161	22.86720219948798	0.7369699910808241
2000	43.24792099793558	-461671.2931382996	22.89828497422631	0.7397930726479056
2100	45.42789416492451	-460840.56490629714	20.987206772671023	0.6391731232730369
2200	47.63013762392802	-461669.2950801987	20.694999058214446	0.6229691346239754
2300	49.82506574894069	-460902.85883438203	20.77260066299265	0.6357195033954708
2400	51.96890287392307	-461395.5866526011	20.833784835479896	0.6410777282575015
2500	54.16312770693912	-460757.27072984475	20.63547809429086	0.6264943252198973
2600	56.30525341493194	-461597.00777675665	20.85624413851512	0.6441020178670579
2700	58.474558122921735	-460991.38350638683	20.807674941022437	0.6398939118180478
2800	60.65462724791723	-461656.09113742976	20.862291365956057	0.6422233805785653
2900	62.82900937291561	-462117.81455525715	20.834753085953338	0.6410660575286649
3000	64.98746858190862	-461987.2694681382	20.830261031356514	0.6387421394574399
3100	67.16358774789842	-459713.5720993656	20.44446230585409	0.6106243789368002
3200	69.36328045689152	-459545.3194997052	19.90754203094443	0.581402578430159
3300	71.51838341588154	-460822.86066742876	20.484196172264596	0.6116483489071257
3400	73.65843687389861	-460850.0021878935	20.568569259001546	0.6168469857344365
3500	75.81148874890641	-460912.84538833325	20.691117526032798	0.6308260974011584
3600	77.99060841591563	-462039.64476704516	21.710284160793776	0.6836436351513866
3700	80.14434024991351	-462087.7118307385	22.124256561469693	0.7011870657420822
3800	82.29387962489272	-462030.2734790751	22.367385911276486	0.7133747254794469
3900	84.47273637491162	-462154.96238702064	22.526151412640814	0.7205986034860279
4000	86.66963637492154	-461949.4390222763	22.635592692171013	0.7261212938116686
