# 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/brainphantom256x256_pdps_known_rotation_20_10000_9000", dynrange = 1.0, Λ = 1.0)
iter	time	function_value	psnr	ssim
1	9.66701190918684e-6	563529.9047468598	-0.02636674332749045	0.0041734697881560985
2	0.02877075001015328	478860.2653100332	0.4729792052599031	0.007393427206780912
3	0.04311683401465416	446626.2876367058	0.8541619936646793	0.010070943327301872
4	0.0722115840180777	400243.37978844263	1.1819917687548782	0.01207211383280447
5	0.08647975101484917	362903.75240876776	1.4927527272327812	0.014045577004022168
6	0.11542437601019628	323517.27880956477	1.7768935562065704	0.015131728686981212
7	0.1297638760006521	284278.7388255165	2.084840612954693	0.016023772194017228
8	0.1585740430164151	246266.23350080894	2.3992492144808355	0.01777398624054255
9	0.17293237600824796	205496.1747394166	2.7261910521095767	0.019900948428165973
10	0.20242170902201906	166722.51245009553	3.0487589816465	0.02154424058935638
11	0.21695545903639868	129595.94543380497	3.385893078486955	0.023407573158945073
12	0.24588008402497508	92589.05286840507	3.7193978096084193	0.026250425132621193
13	0.26070950002758764	56390.57057934398	4.062803280783469	0.029096309789118112
14	0.2903478340303991	20805.631665172947	4.408741438532354	0.03517006620894968
15	0.3051171670376789	-10472.99776275068	4.748968791332911	0.035925703155609244
16	0.3348170420504175	-42532.320262397945	5.102500432878775	0.04389283941473977
17	0.34966662505757995	-63145.15661342107	5.420055379279055	0.038570550470243284
18	0.37952020805096254	-105581.3321454509	5.817590366720566	0.0699281291347063
19	0.3943494990526233	-136013.77302614387	6.184812400389773	0.08905387566479556
20	0.4239988740591798	-165173.91865567074	6.5532405370099225	0.11585897025230812
30	0.6470629570831079	-350410.3092912984	10.257725273273355	0.44985329359900694
40	0.8692494570859708	-421323.05182432005	13.75695137502648	0.5367817455531115
50	1.0912682490888983	-448462.47855350334	16.557305969414603	0.5615964673132685
60	1.3135556660708971	-457329.7274138537	18.299840086810896	0.5918948102927029
70	1.5357775830780156	-459392.49313135236	19.196551218168445	0.6105889882148734
80	1.757238083082484	-460551.21970442025	19.66725132732669	0.6230537796755524
90	1.978470583067974	-460963.92788536975	19.95815120963067	0.6275886510593224
100	2.199967500055209	-460968.5136157653	20.158816726108316	0.6303481318114454
110	2.406070500059286	-460815.0629428459	20.273378620119086	0.6315258766460485
120	2.6272950830461923	-461391.0023669823	20.361316525602497	0.6299580032555294
130	2.8441925410588738	-461640.39797228074	20.29331353840849	0.6187626891320379
140	3.0578791240695864	-459662.8562124126	20.02285313904919	0.5990111512646492
150	3.270875666057691	-461370.35583897616	20.46773211035933	0.6262588424897644
160	3.484006332058925	-459576.05869975314	20.102864960529942	0.6056310111676781
170	3.7020841650664806	-461248.63484423445	20.603684600068245	0.6333167666452466
180	3.9180504980613478	-460765.0076530978	20.47711976856618	0.6204972887990944
190	4.132013540074695	-460368.46346053574	20.4734589290177	0.6212314667185186
200	4.345908081071684	-460313.5973109162	20.37684210688802	0.6144396614448289
300	6.515198498091195	-461463.3258037997	20.70079018185767	0.6313648991441537
400	8.674415331101045	-461055.71188931796	20.744151767916932	0.6322798743157789
500	10.809096914104884	-460878.2453239834	20.72149423121758	0.6333073627426274
600	12.9502694561088	-461695.6020332558	20.82718945356902	0.6423045950762352
700	15.075433039106429	-461156.3582565605	20.83044564694633	0.640576703831315
800	17.23009541409556	-461236.9267132035	20.788354681123767	0.6404501123864409
900	19.384372956119478	-461340.8723486018	20.875470490310303	0.6432821078461084
1000	21.53408408112591	-461358.67638304026	20.67569407347818	0.6306102009986503
1100	23.731297456135508	-461481.1798425294	21.58477855156007	0.6796034121746438
1200	25.93844637312577	-461195.35640629195	21.888883089594948	0.6948712563042606
1300	28.142789956124034	-461581.4433667799	22.055022593892648	0.7019982403113245
1400	30.35619770613266	-461415.64552988095	22.14989706517471	0.7073923196451856
1500	32.5616323731374	-461419.70460595266	22.205496368837508	0.7099419292457212
1600	34.781294956133934	-461451.33927862835	22.24298252407411	0.7129063330327251
1700	37.013551790150814	-461358.65101609955	22.269585673477646	0.7145236707879179
1800	39.23503729014192	-461467.9945638083	22.290460608771888	0.7152698571700562
1900	41.47022087313235	-461471.4919548442	22.302591681800166	0.7142836224093747
2000	43.68976458214456	-461551.3574866272	22.314072776195385	0.7160223160518155
2100	45.92007637413917	-460768.7764849436	20.92667388700574	0.6362685557474019
2200	48.14280408312334	-461604.14464767574	20.6642175219841	0.6214101178266135
2300	50.36941704212222	-460837.8507565543	20.746081630153025	0.6343319435963962
2400	52.5926038751204	-461329.4369833698	20.807704620648636	0.6398955555873252
2500	54.81101500010118	-460692.66350547015	20.612767421013633	0.6255728642417271
2600	57.04147754111909	-461529.55400986254	20.826545746024557	0.6427887784086174
2700	59.25178045811481	-460926.41354568733	20.780211563323174	0.6386339330146128
2800	61.46456374911941	-461589.13178382907	20.835574460037808	0.641021392507188
2900	63.68609691612073	-462050.4510823144	20.807188809408345	0.6399229006264731
3000	65.89733320812229	-461920.2802446627	20.802676800846037	0.6375449286622252
3100	68.11656291713007	-459649.8329814513	20.419635332777666	0.6093226284572139
3200	70.31011962614139	-459483.942951215	19.890985331381728	0.580683784125292
3300	72.50147446015035	-460761.88374322385	20.461488858229913	0.61048027836114
3400	74.7120566261583	-460789.333355378	20.54488986265257	0.6156325162898197
3500	76.88998116817675	-460849.63754284545	20.666478495916927	0.6297139766029789
3600	79.08756770918262	-461930.4543876718	21.551416725550652	0.6780864912552933
3700	81.21363162618945	-461973.01148512	21.876340408887742	0.6924491947011493
3800	83.34992416817113	-461915.365298515	22.052824827253797	0.70207661919722
3900	85.47732883517165	-462037.79113190103	22.154002577818247	0.7069076855544503
4000	87.61745541816344	-461832.4036635608	22.2171888164634	0.7104784823402773
