# 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.25, ρ̃₀ = 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_25_10000_9000", dynrange = 1.0, Λ = 1.0)
iter	time	function_value	psnr	ssim
1	1.0707997716963291e-5	565295.3803130721	-0.02636674332749045	0.0041734697881560985
2	0.028988374979235232	479936.6940496358	0.47303816102183643	0.007393645149152622
3	0.043390999955590814	447337.29465607984	0.854226057935519	0.010071411555931864
4	0.0721233339572791	400782.8347744507	1.18213124589612	0.012074005506735772
5	0.08648137596901506	363325.6211593988	1.4928857970847578	0.014046910434111803
6	0.11522816796787083	323918.3063501865	1.7772320220108926	0.015142105452964183
7	0.1294931679731235	284613.91669050674	2.0852821537068147	0.01604272054989946
8	0.15891237696632743	246556.16302408997	2.3996586133703444	0.017791738245918693
9	0.17371837695827708	205754.27803143955	2.726575250612619	0.01991963038910463
10	0.20407216894091107	166966.02838565764	3.0492558541473107	0.021570786476998133
11	0.2188391689560376	129814.00315071123	3.3863627763101163	0.023429474952561685
12	0.24864716894808225	92792.84521357977	3.7199365131118367	0.02628347653179347
13	0.263487293937942	56582.91155478681	4.063364569210158	0.029137144336865067
14	0.29324075192562304	20988.567789692377	4.409386636798599	0.0352167010386286
15	0.30805529293138534	-10295.26284233826	4.74956646748376	0.035968334540495285
16	0.3378625009499956	-42361.61350861602	5.103090713525757	0.04393899951957969
17	0.35265508395968936	-62975.750637252146	5.420604745988592	0.03861925498906002
18	0.38227987597929314	-105433.89629825394	5.818198893069195	0.07001252443904628
19	0.39705070995842107	-135873.24882035606	6.185481103643483	0.08918035545676753
20	0.42667195995454676	-165041.67073430505	6.554008333369941	0.11606388338067497
30	0.6487572929472663	-350329.12585618434	10.25892748309177	0.45051988347641203
40	0.8703106269531418	-421257.56846832554	13.759697799257946	0.5369905547189335
50	1.0910189599380828	-448409.6850226374	16.561888222150902	0.5618204857355307
60	1.311915167956613	-457287.31000737386	18.30424182807881	0.5921628166553292
70	1.5334339599648956	-459354.6405065188	19.19989020173965	0.6107655567814587
80	1.7532776269654278	-460516.4579955916	19.66928519780336	0.6230980035009128
90	1.9744140849797986	-460929.99868724524	19.95896735480432	0.6274534537436341
100	2.195762876974186	-460936.0032709116	20.158559352195283	0.6301436186768921
110	2.402005168987671	-460783.9362610121	20.272062278269246	0.63135545930282
120	2.6233021689695306	-461360.1949246803	20.35918100523153	0.6297748047702523
130	2.844342960976064	-461610.14580466604	20.29046849102887	0.6186381399891705
140	3.057389544963371	-459632.46062880853	20.020121611450442	0.5988465647798757
150	3.270347752957605	-461340.3505556192	20.46326349816168	0.6260325189681373
160	3.4883546699711587	-459546.2758162923	20.099110619501104	0.6054418538570915
170	3.709388503979426	-461218.4538160261	20.59787271038905	0.6329780642308959
180	3.9320040459861048	-460734.6943387905	20.470688923867325	0.6201743838546125
190	4.154303753981367	-460338.07616257825	20.466438994650158	0.6208264349748359
200	4.374596586974803	-460283.77306294645	20.370670854944866	0.6141050527619258
300	6.500672378984746	-461432.23007193283	20.69048484786439	0.6307340247157932
400	8.656305420998251	-461023.5199178825	20.73159050842538	0.631642938062736
500	10.805220463022124	-460846.7197019425	20.708273943783404	0.632650913786067
600	12.939912880043266	-461663.29846633866	20.814703316201324	0.6417158248320463
700	15.105015671055298	-461123.6259968327	20.81817921046001	0.6399316739720301
800	17.27866017105407	-461205.0380922281	20.776946567530402	0.6399586271099702
900	19.40851621303591	-461308.68257504	20.86445143583836	0.6427741187143806
1000	21.559861588058993	-461327.08990290936	20.666444093330387	0.6301150855901796
1100	23.70320417106268	-461433.7143383608	21.510821130362135	0.6764269050732626
1200	25.835706171084894	-461144.6388166647	21.77400835736371	0.6899647970249609
1300	27.98031025409	-461530.0409253144	21.91222094443246	0.6958856169016354
1400	30.167393129086122	-461363.87551534275	21.987995213620447	0.7002111032517541
1500	32.319458087091334	-461367.2535378648	22.030290138091843	0.7020665634450887
1600	34.49852171208477	-461398.669522093	22.058277156581234	0.704431698123562
1700	36.665211254090536	-461305.7005398814	22.07830437140478	0.7056477572908388
1800	38.835525087080896	-461414.1410352026	22.095690347338014	0.7061579182720408
1900	41.03871287908987	-461417.82311967056	22.103590902828223	0.704949959443096
2000	43.240473629091866	-461498.15629946435	22.110776223279657	0.7064599196032031
2100	45.432963504077634	-460733.8197462368	20.895466690283705	0.6349044884150639
2200	47.6270135450759	-461572.39957110415	20.648874187635677	0.62064623032515
2300	49.7922158370784	-460806.46629412984	20.73444678788794	0.633691292241924
2400	51.99507387907943	-461297.1195402388	20.795833169076037	0.6393525357171586
2500	54.14621908808476	-460661.1095265538	20.60243646360756	0.6251195477566365
2600	56.300207796070026	-461496.8557794344	20.812971685825232	0.642148012456536
2700	58.4340616290865	-460894.8052721762	20.767953596914978	0.6380449865468143
2800	60.56881929509109	-461556.56929934415	20.823589076899868	0.6404609974890327
2900	62.762007337092655	-462017.73445431684	20.79475546518493	0.6393959026714731
3000	64.89748317009071	-461887.81847307173	20.790415051174463	0.6370245991265441
3100	67.06207783607533	-459618.59477879986	20.407548253791102	0.608709929122851
3200	69.19232391906553	-459453.69076501305	19.882424966557693	0.5802867688684055
3300	71.35106079405523	-460732.2678339048	20.45128690878247	0.6100237598179833
3400	73.5085932520451	-460759.9102170203	20.53433018917123	0.6151239379826552
3500	75.65015529404627	-460818.8676992728	20.65554863860025	0.6292105119175335
3600	77.86505233607022	-461882.57213369926	21.477584827838008	0.6750485073236959
3700	80.02096908606472	-461922.4051090157	21.7638070933217	0.6879424660128474
3800	82.1604639190773	-461864.01715284266	21.914493874493665	0.6962527742787488
3900	84.3396060440864	-461985.57766614074	21.996430609233528	0.7001251679692586
4000	86.48145687708165	-461780.1805682422	22.04599090529063	0.703083117774591
