# params = (ρ = 0, verbose_iter = 100, maxiter = 10000, save_results = true, save_images = true, save_images_iters = Set([5, 10000, 3000, 30, 5000, 1000, 1, 6000, 300, 2000, 2500, 7000, 3, 8000, 9000, 500, 25, 50, 7500, 2, 10, 4000, 250, 100]), pixelwise_displacement = false, dual_flow = true, prox_predict = true, handle_interrupt = true, init = :zero, plot_movement = false, stable_interval = Set([4986, 4700, 3634, 4576, 4901, 9676, 9923, 4522, 9773, 3406, 3672, 2841, 8908, 9835, 2876, 4978, 9802, 3337, 4136, 4150, 3293, 3220, 4888, 4030, 3324, 3982, 3298, 3646, 3719, 9614, 4011, 8897, 3855, 3111, 4584, 4081, 3163, 3914, 9523, 4582, 9306, 4687, 2989, 3762, 4021, 4357, 3408, 9883, 9266, 4962, 9387, 9909, 3326, 9439, 3711, 2874, 4738, 4836, 4464, 3271, 9776, 4587, 9329, 9152, 9956, 3319, 4527, 2746, 9246, 3943, 3292, 4238, 3651, 2920, 4720, 9374, 2570, 9593, 2783, 4811, 3105, 9277, 8714, 8845, 9816, 3570, 4169, 8723, 3288, 4548, 9175, 9168, 9618, 4910, 3557, 8918, 9164, 3382, 8769, 3937, 3338, 2982, 4282, 4951, 4764, 9464, 3044, 4817, 9558, 9698, 3587, 9070, 2775, 9318, 8928, 3628, 4783, 3167, 4272, 9319, 3089, 3135, 2932, 9667, 4243, 9448, 3938, 3423, 4803, 2999, 4051, 4672, 4280, 2684, 4392, 4232, 8974, 8915, 3728, 3148, 4383, 4644, 3809, 9126, 9433, 3617, 4342, 4961, 9710, 4060, 9852, 9960, 9440, 4918, 8751, 9117, 4151, 9701, 2713, 3571, 3722, 3507, 4156, 9165, 3140, 9521, 9489, 2724, 4984, 3362, 8776, 3526, 9045, 9524, 8998, 9054, 3063, 9820, 9345, 3792, 9513, 4126, 4229, 4413, 3450, 9859, 4172, 3043, 3052, 4329, 3835, 3015, 3856, 3815, 2536, 4756, 3691, 4153, 3495, 4796, 2863, 8824, 4352, 2901, 4923, 3008, 9477, 3877, 2895, 3979, 9535, 8930, 2748, 4125, 3614, 4643, 3092, 9902, 4105, 4945, 2756, 9517, 9888, 4649, 4017, 2983, 3844, 4726, 8935, 4826, 4682, 9938, 9851, 4343, 9869, 2504, 2761, 3696, 9566, 2978, 9791, 4663, 9365, 9429, 3115, 9401, 4757, 3172, 9567, 9125, 8705, 3690, 9344, 3080, 2692, 9097, 3965, 9661, 3334, 4991, 2682, 3359, 3447, 4953, 4515, 4658, 9490, 3242, 2979, 3817, 8811, 8784, 9845, 9317, 9382, 2541, 4531, 9592, 3931, 8766, 9103, 9393, 2950, 3926, 9452, 4189, 3396, 2556, 9497, 3471, 9058, 4117, 4045, 3702, 2571, 9629, 3340, 3344, 4377, 9785, 3347, 4896, 9903, 3192, 4772, 4326, 3616, 4877, 8721, 9836, 9925, 4824, 9981, 3322, 3380, 3773, 3806, 3110, 8836, 9552, 3356, 9511, 9238, 8861, 9258, 9953, 2518, 8823, 9231, 8822, 9917, 4606, 9793, 2617, 3021, 3117, 4206, 2714, 8973, 9189, 2927, 9968, 3463, 9663, 4003, 2510, 8916, 3289, 4443, 4878, 4558, 4066, 4354, 9123, 3519, 9901, 4037, 2691, 4147, 4610, 8762, 2672, 2628, 3917, 9331, 9986, 4840, 2896, 8742, 3870, 4145, 8826, 4141, 4400, 4088, 2853, 9683, 3385, 3307, 8967, 3641, 8704, 8910, 2935, 2608, 4317, 4000, 3784, 4471, 4786, 9904, 2803, 4597, 8785, 9606, 3508, 3438, 4842, 8715, 8934, 3599, 9060, 3566, 4513, 4998, 8825, 2813, 9864, 2574, 2609, 3156, 2852, 4749, 8980, 4220, 4712, 3265, 3056, 3837, 4018, 4371, 2650, 3309, 3697, 4375, 2997, 2641, 3609, 8707, 3037, 3038, 3990, 3759, 4847, 9763, 9005, 9035, 4267, 8865, 3407, 3118, 3267, 3971, 9473, 2871, 4473, 3649, 9009, 9179, 4958, 4444, 2633, 2758, 3737, 4054, 9565, 4058, 9127, 3352, 9326, 2652, 4758, 4380, 8757, 4165, 2946, 3764, 4630, 4605, 2925, 2788, 4737, 9449, 9406, 2683, 9745, 4499, 9732, 9586, 4264, 8759, 2595, 8743, 9147, 3235, 9316, 9974, 2632, 2604, 9185, 2802, 9538, 2688, 2809, 9679, 2915, 3505, 3247, 9685, 9178, 9648, 3460, 3274, 2585, 2719, 3832, 4115, 4364, 4541, 8708, 3188, 4883, 3501, 8728, 4478, 4743, 3291, 4414, 9403, 9514, 4071, 9105, 4476, 2666, 9151, 2637, 4419, 8740, 4965, 9867, 8989, 4294, 9982, 2859, 4210, 2884, 3429, 8821, 3880, 9815, 2513, 4812, 4839, 3882, 3686, 9729, 2557, 9598, 9717, 9008, 3622, 4594, 4699, 9788, 3861, 4075, 8867, 2767, 8717, 9577, 2988, 4200, 4905, 2521, 2612, 2796, 4235, 2855, 3905, 9537, 9935, 9599, 3077, 3625, 8932, 9597, 4434, 9463, 3955, 2695, 9024, 3901, 4087, 3684, 9106, 3346, 9442, 2711, 3387, 3692, 4052, 4529, 4062, 3399, 9796, 3541, 3087, 9850, 4526, 3862, 4791, 3579, 8877, 9754, 9346, 4916, 3026, 4092, 8933, 4297, 4602, 9792, 4754, 3099, 4585, 9531, 4580, 2829, 8755, 9846, 3164, 4488, 9206, 4307, 9525, 4048, 9879, 9352, 3070, 4568, 4850, 9296, 9114, 2611, 9515, 4437, 4224, 3775, 4532, 4255, 8767, 9978, 9080, 2833, 9209, 9409, 3744, 8831, 3238, 9066, 2514, 4103, 4121, 2638, 4885, 9313, 9158, 9719, 4174, 4912, 4882, 2717, 2739, 9195, 9212, 3160, 8947, 2971, 9684, 3769, 3493, 9828, 3534, 3360, 3707, 3584, 3807, 3480, 4907, 9510, 3120, 4939, 4859, 8712, 3567, 3585, 2839, 3143, 4547, 9484, 9079, 3954, 9170, 3335, 2872, 2851, 4714, 9801, 4236, 9549, 9056, 3130, 4854, 9875, 4729, 9761, 9481, 9384, 3547, 2969, 2916, 3012, 9116, 9254, 3627, 9459, 4891, 4429, 3132, 9093, 2551, 2720, 9756, 3708, 3398, 4253, 3823, 3799, 9987, 4468, 2965, 3748, 3522, 2575, 9017, 3863, 3282, 3142, 9780, 9027, 3725, 3228, 4566, 4970, 9016, 2731, 3517, 9394, 8945, 4181, 3621, 3357, 3145, 4806, 8725, 8862, 9007, 4296, 2600, 3798, 4248, 9680, 4963, 3328, 8800, 9797, 4875, 9309, 9942, 9572, 3186, 8837, 4474, 4727, 4360, 4999, 4005, 9417, 9767, 3314, 2569, 3055, 3144, 9817, 9999, 4134, 9637, 4241, 9204, 3442, 9973, 2500, 2808, 4871, 8942, 3091, 9420, 4072, 4971, 2958, 9854, 4016, 9194, 3787, 9068, 2791, 2827, 3152, 4239, 3234, 4716, 9555, 4410, 2785, 3660, 4709, 4446, 4619, 9813, 3187, 4788, 2512, 4378, 9918, 3561, 3138, 4575, 2986, 9623, 9498, 9996, 4265, 4964, 9702, 9424, 4110, 8955, 9249, 3896, 9894, 9361, 4827, 4942, 9990, 3518, 8711, 3996, 4830, 3721, 3456, 3009, 2834, 4295, 3818, 8750, 4881, 2930, 3829, 2923, 8894, 3112, 3046, 3667, 9748, 4430, 4904, 9135, 4753, 9880, 9455, 9446, 9176, 4039, 2768, 9285, 4656, 4833, 3401, 4521, 4717, 9800, 9980, 3200, 4654, 2806, 4642, 2663, 3283, 9583, 9916, 4631, 4518, 9569, 3185, 4674, 4973, 9724, 4006, 9362, 9144, 4440, 4492, 9881, 4114, 2552, 4708, 9688, 9764, 3214, 3086, 9458, 4084, 4887, 3932, 4480, 4283, 9335, 4723, 2562, 4320, 4763, 4222, 2567, 4124, 2942, 3327, 4315, 3315, 4992, 9076, 9321, 4591, 3678, 8913, 9993, 3589, 9853, 3565, 3854, 4586, 9516, 2618, 9694, 2875, 4809, 3594, 9081, 9260, 9560, 3998, 8978, 4853, 9210, 3126, 4183, 9681, 9777, 3027, 4167, 4192, 9325, 3867, 2674, 9343, 3846, 3446, 3202, 3626, 4186, 4805, 2821, 4211, 8828, 4482, 2693, 4316, 4127, 9696, 9348, 4422, 9622, 3778, 4624, 3910, 3911, 3215, 2592, 4692, 9678, 4609, 3913, 9376, 2781, 3665, 4026, 8805, 8995, 9795, 9405, 3256, 2643, 8969, 4178, 2667, 9871, 9082, 2614, 3884, 8846, 8833, 4249, 4928, 3486, 3777, 2762, 3969, 4615, 8791, 4113, 2565, 4242, 8956, 9723, 9486, 3107, 4502, 9369, 4702, 9122, 3454, 2964, 3895, 3802, 4370, 4678, 2847, 9183, 9349, 2635, 4207, 9811, 2838, 4822, 4509, 4040, 9363, 2743, 3203, 9930, 2858, 4425, 2745, 9536, 3169, 3772, 3833, 2938, 4832, 9115, 3182, 9542, 4603, 3987, 9738, 3391, 3597, 9908, 3162, 2766, 8780, 8996, 9447, 4903, 4914, 8911, 9059, 2782, 4069, 8888, 8733, 4974, 9905, 3716, 4002, 4677, 4356, 4322, 3183, 3662, 9553, 9269, 9133, 4765, 9261, 2553, 2967, 8752, 4937, 3637, 4771, 9826, 9895, 9740, 4581, 9759, 9715, 9931, 4902, 9966, 4009, 3379, 9615, 4661, 4728, 9201, 8987, 9327, 9395, 3069, 9954, 9428, 9855, 9617, 2505, 2528, 2898, 9337, 8793, 3417, 4313, 9263, 9274, 3537, 3440, 9779, 3064, 3224, 2828, 3331, 4968, 4676, 9095, 4203, 2770, 3709, 4694, 2991, 3992, 4279, 4417, 3261, 9100, 3195, 9794, 9847, 3161, 9046, 3623, 2976, 8853, 9052, 9575, 9714, 9487, 3650, 3553, 3741, 9229, 9025, 8842, 2751, 9783, 9582, 3735, 9873, 4157, 9613, 8991, 9075, 4162, 3191, 3040, 3481, 4519, 9602, 4683, 9962, 4258, 9166, 9653, 4426, 4180, 4787, 9574, 3285, 3659, 2564, 9136, 9522, 8872, 3562, 9399, 4244, 4995, 3894, 3897, 4259, 8976, 3532, 4593, 3959, 9419, 2694, 9988, 9798, 8949, 9787, 2596, 4675, 3780, 3824, 4382, 2728, 9436, 4790, 2747, 3017, 8874, 9311, 9849, 2588, 2704, 4020, 4158, 9023, 9205, 3888, 2703, 4767, 2735, 2888, 4004, 3763, 8965, 3906, 3088, 3257, 4715, 3325, 4941, 3262, 3790, 4465, 4057, 9211, 4401, 3268, 4412, 8722, 3084, 4076, 9899, 9961, 9713, 3732, 3618, 9770, 4555, 3402, 9682, 4608, 9233, 9824, 4184, 9563, 3194, 8706, 2945, 9509, 3270, 9652, 4335, 3045, 3428, 9143, 4273, 3081, 3488, 4442, 9502, 9656, 4449, 3723, 9368, 4386, 3244, 3339, 3223, 3640, 4338, 4523, 4815, 4943, 3834, 3409, 3073, 3246, 4148, 9279, 9283, 3159, 8718, 3991, 3237, 2544, 9063, 9474, 3278, 9015, 4245, 3964, 4177, 3774, 9391, 3916, 4900, 4097, 4573, 4083, 9893, 4793, 2606, 3742, 4373, 4472, 2715, 3724, 3065, 4752, 3836, 4312, 3843, 2926, 3632, 4333, 4397, 4895, 3312, 4849, 9708, 4860, 9499, 3123, 4458, 4358, 3419, 9034, 9378, 8927, 3873, 3208, 9590, 4779, 8858, 3779, 9646, 4302, 8903, 4744, 3306, 4761, 4843, 3545, 9077, 3639, 9700, 9336, 3248, 2832, 9412, 4292, 8914, 3412, 4934, 2790, 9967, 2668, 2826, 3528, 9760, 4696, 2933, 4493, 9019, 9220, 9579, 4740, 2502, 4616, 3781, 4129, 3730, 3543, 9657, 9781, 9213, 3100, 3752, 2598, 9655, 3007, 4305, 9541, 9940, 9255, 3866, 2992, 4969, 9878, 4932, 3050, 3478, 2529, 8709, 8883, 3701, 3604, 3838, 2940, 3444, 3921, 9389, 4398, 4303, 3333, 4351, 4774, 9132, 4196, 2577, 2673, 4735, 8732, 4456, 4814, 4470, 8979, 9612, 3513, 9642, 9877, 3999, 9053, 9043, 3386, 9339, 9782, 8924, 3970, 3516, 4416, 4029, 2880, 4274, 4604, 2550, 4467, 2658, 9067, 2649, 4281, 8939, 3177, 3168, 9475, 3239, 2886, 4122, 2581, 3451, 3661, 3668, 9253, 4792, 9705, 4252, 9021, 9099, 9778, 3427, 8950, 4789, 4346, 3521, 4368, 3468, 4055, 8902, 9844, 2737, 2830, 3355, 3656, 4171, 3860, 2819, 3902, 4655, 3476, 3629, 8957, 8835, 2689, 3552, 4144, 2849, 2944, 9998, 2644, 4966, 9111, 9454, 9355, 4256, 9607, 4967, 3575, 3303, 8775, 3230, 3958, 2548, 2955, 8720, 2970, 9444, 3290, 2773, 4008, 4844, 9742, 4825, 2943, 9271, 9102, 9128, 3219, 2744, 9300, 3694, 4987, 3462, 3949, 3994, 9889, 3254, 3467, 4310, 4639, 3869, 4500, 2815, 9771, 3151, 3673, 8905, 8819, 3102, 3889, 4794, 9357, 4111, 4173, 4080, 2616, 3103, 4574, 3887, 9697, 8849, 2622, 4501, 3263, 4334, 4228, 3768, 9418, 2887, 4670, 8892, 4719, 2994, 2755, 2753, 9083, 2857, 2710, 4569, 3619, 4732, 4778, 9250, 3891, 4420, 9862, 2549, 3054, 9827, 9031, 9372, 9150, 2699, 9379, 3822, 4647, 3733, 8959, 9528, 2885, 9914, 4625, 4308, 9160, 3879, 9964, 2626, 3392, 4485, 3302, 8820, 8863, 9749, 4477, 9733, 8730, 9441, 3795, 9674, 9431, 3189, 3612, 2947, 2659, 3136, 2525, 9280, 4331, 4909, 8701, 8789, 9227, 4865, 2538, 2867, 9351, 4278, 2879, 4195, 2584, 3383, 9959, 3893, 9937, 2572, 4680, 8756, 4028, 9422, 2685, 4463, 9751, 3029, 2972, 4831, 8841, 2846, 4394, 9660, 8895, 9768, 2877, 4182, 8735, 8982, 8765, 4254, 9544, 8854, 9884, 4561, 9829, 2881, 3924, 9920, 8839, 2968, 2837, 3904, 3758, 4415, 9601, 9295, 4460, 3682, 2646, 3390, 4516, 4023, 3388, 9651, 3849, 2597, 4013, 4179, 4212, 8881, 3746, 3497, 9200, 4160, 9268, 8777, 4846, 2922, 4611, 3395, 9808, 3349, 4684, 3430, 2591, 9892, 3057, 9265, 2778, 9333, 9314, 9668, 8860, 3789, 3127, 3975, 3133, 8997, 2698, 2890, 9293, 3512, 4219, 4452, 9485, 4022, 3432, 9746, 3564, 3685, 9312, 3147, 9434, 4407, 2835, 3367, 4353, 4538, 9308, 3907, 2729, 9949, 3472, 3062, 4695, 2542, 2909, 2656, 3615, 2524, 2687, 9262, 4251, 9292, 3368, 4266, 4669, 8953, 9927, 3134, 3431, 2509, 9171, 3714, 9896, 9943, 9036, 3227, 4590, 9556, 2647, 4975, 3085, 4116, 4168, 4408, 3950, 9156, 9627, 9946, 9138, 8880, 9699, 3414, 9585, 2836, 3993, 2517, 2631, 3226, 8922, 8788, 3674, 9807, 9897, 9182, 3782, 3814, 4930, 2627, 3944, 2547, 2948, 3598, 9774, 3671, 9472, 3755, 4099, 2900, 3989, 8840, 8960, 3939, 4086, 4919, 4038, 4540, 3853, 8961, 3974, 4033, 2716, 4318, 2615, 4065, 8864, 9805, 9984, 4931, 2985, 9120, 3718, 8898, 2508, 2697, 2545, 8802, 9087, 9235, 3372, 9239, 2580, 3946, 9437, 3287, 9706, 3919, 4543, 9415, 3500, 2593, 3437, 3706, 3676, 2630, 3277, 3900, 3643, 4997, 9639, 3592, 4497, 9476, 3502, 8782, 9247, 9014, 4130, 9245, 9671, 4486, 3857, 3494, 3225, 8971, 9529, 3093, 3166, 8852, 9218, 3966, 2995, 3384, 3734, 4920, 9186, 8871, 2527, 3683, 4260, 4204, 2587, 9728, 3514, 9030, 9533, 4506, 9167, 9610, 9690, 9140, 4454, 3601, 4838, 3378, 8844, 3473, 4988, 3705, 3174, 3631, 4135, 2844, 8774, 9145, 9568, 4807, 3006, 9970, 3930, 4653, 9303, 3608, 3232, 4348, 9457, 4876, 9675, 3433, 2546, 9248, 4623, 9812, 4393, 3981, 9370, 3655, 2501, 3258, 9762, 3794, 2566, 3540, 3153, 4628, 4143, 2586, 3524, 3071, 4851, 2702, 9101, 8724, 9591, 4748, 4475, 8984, 9975, 3549, 3366, 9507, 4535, 4570, 3317, 8891, 4972, 3963, 3218, 4484, 9971, 10000, 3765, 2670, 4487, 4855, 9332, 4077, 2893, 3483, 9695, 3363, 4926, 3158, 2954, 4024, 3915, 4441, 9334, 4563, 2732, 9230, 4217, 3754, 9890, 3510, 9084, 9221, 3404, 8792, 9822, 9856, 4372, 4852, 9004, 9832, 4577, 9711, 4713, 9064, 4507, 9259, 9404, 9356, 4225, 3441, 9720, 4213, 3791, 4202, 4944, 2856, 3119, 3652, 3245, 9425, 3364, 2951, 4275, 9237, 2738, 3176, 3122, 8999, 9377, 4588, 9645, 9644, 2845, 9139, 9562, 9721, 3058, 4287, 3731, 9018, 3020, 2602, 8921, 4762, 3677, 3381, 2824, 2861, 3206, 4286, 3713, 4746, 3053, 4770, 9273, 3803, 4285, 9091, 9636, 4892, 3729, 3504, 3743, 4629, 2866, 4214, 4404, 4332, 9347, 8783, 9821, 3470, 4776, 4327, 8988, 2949, 2734, 3095, 9320, 4193, 9172, 3786, 4246, 4101, 2787, 3067, 3014, 3940, 2953, 9315, 4381, 3157, 8801, 2765, 2862, 8737, 9032, 4450, 9570, 9727, 3757, 3740, 3648, 4455, 4952, 4652, 9527, 9085, 4498, 8907, 4007, 8807, 9664, 2961, 4578, 3394, 8884, 3859, 4562, 3602, 2654, 2905, 3956, 2690, 3574, 9806, 3770, 3000, 9121, 9341, 2957, 8726, 2573, 4336, 4025, 3033, 4868, 8886, 9301, 9467, 4731, 4889, 3851, 9354, 9530, 2675, 3825, 8876, 9089, 9939, 9948, 4703, 9848, 2996, 4046, 3657, 9410, 2924, 4589, 3198, 2526, 2709, 4906, 3297, 9842, 2883, 2561, 2657, 8890, 9236, 9364, 2804, 3530, 4799, 3068, 4323, 4894, 9501, 2624, 9804, 9945, 2975, 4504, 9858, 4191, 9616, 9040, 3886, 3196, 8806, 4240, 4034, 4512, 3365, 9322, 3201, 3568, 3828, 4898, 9162, 2919, 9026, 9451, 3923, 9396, 4925, 9456, 4890, 4698, 9716, 2620, 3350, 9289, 2795, 3059, 4250, 4648, 8906, 8741, 2579, 9834, 4557, 3128, 9224, 3016, 9997, 3181, 3370, 3858, 4311, 4989, 4533, 4457, 3972, 8964, 9074, 3903, 3259, 2582, 9775, 4319, 4324, 4742, 3279, 4528, 9119, 4665, 4015, 9049, 8878, 2621, 4617, 3465, 4955, 3461, 9919, 2727, 4095, 9494, 4298, 9936, 3083, 9001, 2921, 3121, 8813, 9198, 9941, 3681, 9256, 2966, 3165, 3509, 3868, 9217, 9388, 3819, 9207, 3593, 4959, 4163, 4858, 4431, 4537, 3453, 3205, 9483, 4466, 9039, 2952, 3700, 8710, 3973, 4710, 8923, 3986, 4300, 4927, 9693, 3962, 9900, 4534, 4828, 8843, 9397, 2607, 2941, 3644, 9098, 8904, 9911, 9000, 4223, 8889, 3019, 8770, 9609, 3253, 9366, 9090, 3550, 3850, 9020, 3922, 2558, 3001, 4546, 9564, 4447, 9423, 2998, 3525, 8703, 4395, 8925, 4044, 4049, 4990, 2776, 2963, 3658, 3925, 9866, 2990, 9413, 3252, 9264, 4862, 3928, 4433, 3693, 9350, 8966, 3977, 4027, 2540, 2911, 3852, 9550, 9709, 9033, 2769, 3124, 9208, 3496, 4362, 9491, 9159, 3255, 3358, 3680, 4306, 4545, 4747, 4689, 3345, 3909, 9353, 8879, 4064, 3801, 4091, 4947, 8851, 4977, 9038, 8946, 3039, 2662, 4636, 2678, 3620, 3449, 4818, 9381, 9050, 4804, 3751, 9557, 4688, 9299, 8754, 3082, 8917, 4284, 8838, 4291, 2700, 2757, 4208, 4880, 4618, 9191, 3703, 4423, 9288, 3393, 8758, 2934, 2956, 4634, 3276, 3613, 4819, 4668, 4137, 2973, 9526, 2814, 2539, 3243, 4096, 9972, 9225, 9958, 4869, 3179, 2873, 3420, 8963, 4390, 3318, 8764, 3353, 4848, 3576, 9924, 9744, 2812, 3556, 4376, 9650, 9580, 4983, 3375, 4870, 2676, 3821, 9496, 9741, 9587, 3229, 4795, 2563, 4954, 9913, 3170, 8731, 8931, 3952, 3010, 9222, 2801, 2583, 9872, 3269, 9965, 2648, 2816, 4768, 4637, 2870, 3933, 4524, 4583, 9735, 4550, 4194, 4483, 4705, 9534, 4421, 9784, 3788, 4154, 4508, 2578, 3727, 4226, 9863, 9758, 9672, 3116, 9088, 3209, 4981, 9243, 3452, 9505, 9882, 3590, 9604, 2708, 4980, 3607, 2899, 9736, 3747, 9985, 3920, 9398, 2741, 4718, 8739, 4112, 4175, 9276, 9665, 9244, 3753, 3499, 3421, 8744, 9055, 4641, 4985, 9837, 4915, 4659, 3022, 9069, 9725, 9298, 9951, 9643, 3582, 9426, 3793, 4262, 3569, 8799, 9304, 9003, 4402, 9216, 2664, 9547, 3538, 9929, 4751, 8812, 3633, 9631, 3041, 2660, 4277, 9620, 9933, 4237, 9823, 9012, 4940, 3079, 9704, 4681, 4571, 3249, 4359, 9750, 8803, 4517, 2794, 2842, 9722, 4823, 4459, 3300, 4107, 3004, 3104, 4014, 4188, 8929, 9630, 9666, 4257, 8772, 2800, 4491, 4063, 8870, 9006, 3011, 9154, 3883, 3983, 9184, 9886, 8983, 9044, 9540, 9673, 3030, 4330, 8798, 9726, 9173, 3294, 3171, 4366, 4539, 4821, 8700, 9242, 9275, 4462, 4922, 9833, 4428, 3260, 9860, 2962, 8781, 4365, 3222, 9619, 3826, 2771, 2679, 3812, 3489, 4620, 3527, 3985, 3457, 2848, 3280, 3756, 9952, 9471, 8746, 8797, 3957, 3726, 4205, 4276, 9632, 9921, 4019, 2619, 4551, 3301, 4913, 3295, 8794, 4396, 8941, 9504, 2831, 9188, 8958, 9730, 9977, 4131, 8977, 4556, 4957, 9432, 4411, 8944, 8716, 4388, 3075, 2823, 9743, 3491, 3736, 9503, 4554, 9305, 2645, 9092, 3193, 3811, 9037, 8856, 4686, 2686, 4982, 3839, 2671, 3342, 4837, 9146, 2910, 2531, 8768, 9625, 3797, 3624, 4640, 2977, 2928, 3539, 3190, 2706, 2763, 3076, 3642, 4050, 9488, 3988, 4355, 2532, 3704, 3096, 4339, 4079, 4711, 2603, 8808, 3712, 4345, 9581, 2810, 9790, 2891, 9991, 2740, 3885, 4960, 2535, 4304, 4979, 9359, 9559, 3320, 4427, 3664, 3595, 4600, 3610, 9757, 3827, 3094, 3881, 3479, 9141, 9868, 9392, 4155, 8986, 3654, 8868, 9129, 4690, 2818, 9839, 2707, 4369, 9983, 3464, 2661, 3942, 9786, 3558, 9686, 4646, 3299, 3321, 2820, 3783, 4379, 8900, 9624, 9659, 2629, 3560, 3929, 4866, 9669, 3647, 3813, 2730, 9461, 9232, 4461, 9324, 9963, 2929, 3023, 3531, 9640, 3275, 3840, 4187, 8855, 8702, 4736, 3426, 2601, 8920, 4560, 2705, 4047, 4271, 4432, 4132, 8968, 4120, 9691, 9737, 3750, 4094, 3968, 4073, 8847, 4363, 4635, 4938, 3573, 4956, 9286, 2984, 8885, 3715, 3842, 4218, 9512, 9825, 3669, 4391, 3024, 9647, 4190, 9270, 4053, 4777, 8796, 2515, 4864, 8817, 2594, 4657, 3490, 9635, 3351, 3424, 4845, 8975, 9500, 4784, 4802, 2903, 4721, 3042, 3830, 3710, 4061, 3443, 9435, 9062, 3371, 4775, 9373, 8936, 8912, 4673, 4031, 8940, 4666, 9022, 4289, 9628, 4685, 2749, 4140, 9163, 9658, 9712, 3484, 2722, 3761, 2523, 3875, 5000, 4873, 4123, 4542, 4109, 4911, 9772, 9284, 2610, 3207, 4780, 8857, 4924, 8952, 8859, 2506, 9470, 3572, 9011, 9789, 2677, 3180, 3336, 3341, 3425, 4781, 4667, 4707, 9252, 4579, 4098, 9922, 9202, 8832, 2981, 4059, 8810, 3032, 3675, 2726, 4494, 3002, 2840, 3251, 3308, 9072, 9584, 3373, 3745, 3831, 9573, 8736, 2560, 3536, 9241, 4451, 3139, 8818, 9199, 9302, 9482, 3264, 4769, 4221, 9578, 8875, 9548, 2792, 3150, 2825, 4438, 4216, 4834, 9870, 9753, 4996, 3025, 4445, 3376, 3405, 9407, 4626, 4808, 4627, 9226, 4409, 9310, 4231, 9375, 9465, 3487, 2701, 2516, 4387, 2534, 3154, 9915, 3114, 9561, 9219, 4697, 4215, 9492, 9932, 9240, 2543, 2725, 3149, 3485, 4230, 2760, 9554, 4453, 9934, 4917, 3129, 2640, 3078, 3231, 9641, 9402, 9532, 9689, 2904, 3876, 9976, 3580, 3653, 9421, 2936, 3101, 4309, 4439, 9495, 8993, 4599, 2897, 3332, 2653, 4722, 9358, 9589, 9177, 3927, 3533, 4549, 3286, 4950, 2918, 3109, 3984, 3113, 4553, 2519, 9876, 9649, 9042, 4510, 4782, 3600, 9113, 3455, 8909, 3611, 8896, 2764, 9234, 2774, 9197, 4632, 8962, 4933, 3892, 9278, 4520, 3377, 9110, 9169, 3699, 4261, 9367, 9380, 9843, 9291, 9944, 4857, 2902, 3250, 3908, 4595, 8848, 4525, 8747, 4139, 4872, 4874, 4633, 2793, 9002, 9926, 3544, 3874, 9061, 9891, 9469, 9654, 9109, 3635, 3845, 8727, 3036, 3666, 9148, 2639, 9330, 3305, 3720, 4119, 4899, 9576, 4436, 4936, 8761, 8763, 4879, 3090, 3173, 3934, 2892, 9187, 4544, 9438, 3003, 8887, 2634, 9340, 9588, 9605, 2914, 3695, 3995, 8994, 9865, 3416, 8985, 4638, 9928, 4622, 2843, 3155, 3596, 8866, 8938, 9287, 9414, 2613, 3563, 4435, 9400, 9181, 4090, 3918, 3630, 3766, 4089, 4337, 3389, 3948, 4406, 4142, 4867, 2777, 3211, 3415, 9910, 4592, 9149, 4010, 9294, 9752, 4861, 2754, 9107, 4481, 3474, 4730, 2712, 4503, 4530, 2623, 3137, 4166, 9803, 3348, 3397, 8937, 9731, 2908, 3890, 4706, 8970, 4106, 9480, 4399, 4621, 3434, 4863, 9621, 4263, 8899, 2865, 2869, 3108, 3871, 3477, 8954, 9137, 2723, 4350, 3506, 8719, 9460, 4835, 3354, 9155, 9131, 4946, 4290, 2894, 4651, 4810, 2752, 4976, 9157, 9048, 4816, 9818, 2889, 3343, 9281, 4479, 3670, 2987, 3072, 3980, 9545, 9251, 2807, 3400, 9906, 3767, 4384, 9112, 3588, 9371, 3978, 4234, 9600, 9638, 3760, 3717, 9297, 4601, 3535, 8981, 9670, 9360, 2718, 4893, 9051, 2817, 4043, 8745, 3061, 9078, 9071, 3872, 8926, 9161, 3212, 4650, 9453, 2882, 9769, 3361, 3498, 3146, 3418, 3796, 8760, 3899, 3178, 4349, 2520, 3936, 2917, 2868, 2589, 4662, 3197, 3175, 9912, 9190, 3369, 3131, 9662, 3583, 3210, 3125, 4227, 4138, 9466, 4750, 4068, 4511, 4829, 3458, 2780, 4036, 3638, 3217, 4321, 8734, 3941, 4001, 2665, 4856, 4067, 8773, 4514, 4886, 9898, 9118, 9603, 3820, 3272, 4921, 4755, 9633, 4660, 2993, 3048, 2811, 9887, 3912, 4405, 8779, 8943, 9520, 4152, 9692, 4198, 4536, 3864, 8830, 9857, 2651, 4403, 9626, 3961, 4741, 9385, 4185, 9992, 9174, 9519, 2864, 3689, 3413, 9608, 9290, 4613, 4935, 3304, 3241, 2937, 9546, 3013, 3313, 3945, 4908, 4701, 4766, 9214, 4341, 3141, 3810, 4199, 9257, 3603, 3233, 4424, 2906, 4679, 9830, 4164, 3060, 3785, 9594, 3310, 4247, 9831, 3967, 9124, 2786, 3448, 3184, 9634, 9104, 9543, 9383, 2822, 3051, 4448, 4490, 9707, 4268, 4596, 2599, 8814, 3466, 9096, 4948, 8869, 3410, 4100, 9739, 9765, 3329, 9342, 2797, 3739, 4161, 8893, 3221, 4102, 8919, 9390, 3035, 3213, 3577, 3236, 4361, 9506, 9809, 3808, 4505, 9841, 4469, 2559, 9338, 9969, 8809, 3546, 9718, 4133, 9228, 4042, 2507, 3411, 9947, 8901, 4209, 8790, 3047, 3947, 3482, 9450, 9130, 4041, 4118, 2530, 3679, 4314, 4801, 3459, 8748, 9180, 9995, 8816, 3520, 2939, 3953, 9950, 9955, 4299, 3848, 9416, 3475, 8850, 2511, 3034, 3800, 2568, 2721, 9386, 9687, 3771, 3898, 9108, 4800, 2931, 8827, 8948, 8729, 3492, 9073, 3311, 4104, 4785, 8951, 2680, 2912, 3018, 4820, 2554, 8873, 4085, 3636, 4489, 9571, 9747, 4035, 3439, 4559, 2681, 3296, 4176, 3976, 4993, 4949, 4760, 9907, 3688, 3316, 4093, 4367, 3445, 4813, 8778, 3605, 8990, 9134, 9799, 3551, 9028, 3559, 9677, 3663, 4645, 3031, 8786, 9328, 4607, 4344, 4385, 9445, 4773, 2537, 3542, 9065, 3841, 2742, 9307, 4797, 2503, 2733, 3878, 4552, 4170, 3284, 3435, 9734, 4733, 3997, 3529, 2878, 2576, 4671, 4745, 4288, 9323, 9838, 3199, 8804, 2772, 4614, 3503, 9595, 3049, 9010, 3586, 4128, 9057, 9013, 8882, 9267, 9468, 3865, 4078, 4565, 9223, 4841, 2779, 4725, 3266, 3554, 4374, 8815, 2907, 2590, 2913, 9203, 9029, 4693, 4347, 9193, 4056, 9427, 8829, 3951, 2605, 2850, 4340, 4572, 9703, 4994, 3005, 9282, 9041, 2959, 9979, 9508, 8738, 2789, 2696, 8771, 9766, 3804, 2860, 4664, 2555, 4032, 9611, 4012, 4108, 4325, 4389, 2533, 3074, 3581, 3749, 3776, 4495, 4567, 4691, 4884, 4897, 9192, 2974, 3935, 2798, 3330, 4082, 2522, 9814, 9874, 4739, 2805, 4328, 4293, 9596, 9086, 9994, 8787, 8753, 9411, 9493, 4201, 9478, 3240, 4734, 9819, 4197, 3066, 3204, 3738, 9957, 9153, 9408, 2759, 2736, 3469, 4564, 3645, 3436, 3523, 4070, 4269, 8795, 2636, 3511, 9810, 3097, 8713, 9430, 9989, 4149, 3606, 4146, 9196, 9142, 3281, 9443, 9861, 2642, 3578, 9551, 3098, 2854, 2655, 4798, 3323, 3687, 4612, 3403, 9272, 8972, 3816, 2980, 4074, 4929, 8749, 3548, 8834, 3106, 2750, 4704, 9518, 9462, 2669, 9047, 4270, 2960, 4598, 3374, 3805, 4418, 9094, 9479, 3273, 3216, 9840, 9539, 3960, 2799, 2625, 2784, 3698, 4759, 9885, 3515, 3422, 9215, 3028, 3555, 4233, 3591, 9755, 3847, 4724, 4496, 4159, 8992, 4301]), noise_level = 0.5, shake_noise_level = 0.025, shake = 2.0, α = 0.25, ρ̃₀ = 1.0, σ̃₀ = 1.0, δ = 0.9, σ₀ = 1.0, τ₀ = 0.01, save_prefix = "img/lighthouse200x300_pdps_known_rotation_25_10000_100", dynrange = 1.0, Λ = 1.0)
iter	time	function_value	psnr	ssim
1	1.3750017387792468e-5	17036.45803843617	4.935377882355947	0.010713389784152977
2	0.004920750012388453	17005.162552740978	4.998210440499573	0.01989122696950344
3	0.010721042024670169	16846.818461004546	5.055130160748168	0.03054886195733225
4	0.0170157500251662	16768.734385888558	5.134264470140812	0.043244460047165756
5	0.024443750036880374	16614.4643946921	5.212944293661926	0.052190240222083864
6	0.032168875040952116	16380.3147716353	5.299194970533515	0.06308969273832403
7	0.03861062502255663	16053.23262867989	5.400566474280241	0.07261375093333282
8	0.04666866603656672	16076.400605852617	5.499302418067636	0.08055837990327151
9	0.05124358303146437	15952.504053366345	5.573432631246778	0.0964266833683668
10	0.05832174903480336	15597.274901445502	5.653632755915028	0.10069933517519858
11	0.06309908305411227	15550.985648362292	5.7490286005168825	0.11622998390510217
12	0.07047241705004126	15337.440166346032	5.860228898311723	0.12187901167820091
13	0.07679025005199946	15167.553657159804	5.922774201379464	0.13012303898893607
14	0.08428600005572662	14999.174013444996	6.031521599107197	0.1429323740219062
15	0.09067116706864908	14762.305701093555	6.152355813717552	0.14961300825231866
16	0.09697050109389238	14741.21007108452	6.215864853497651	0.16610520967183795
17	0.1033847090729978	14628.440095849526	6.312156587416757	0.1679701689035668
18	0.11017887605703436	14482.811673135338	6.416641379429033	0.17516199089166742
19	0.1156167090521194	14280.342620340105	6.4787715172627305	0.18708308756574782
20	0.12173920904751867	14118.326829413492	6.543632756135588	0.20293973746121416
30	0.1931887500686571	12978.708933978507	7.459517435633619	0.2773971809900377
40	0.2721750830824021	12173.662082118079	8.257517669084411	0.36956329087544526
50	0.34751362408860587	11232.036046791727	9.096177393511633	0.42014998639690276
60	0.41832374909427017	10673.975102352151	9.843508583982425	0.4287012368091597
70	0.49229529008152895	10189.430591554255	10.666649521031857	0.5077569758111443
80	0.5699706650630105	9743.007965258097	11.440070951420147	0.5013094861477418
90	0.6472419570782222	9286.218809155098	12.244106624497233	0.5586846831250571
100	0.7194712480704766	9008.379919356885	12.975524014667432	0.6044460865014308
110	0.7920764140726533	8922.248232006757	13.655256211487254	0.5970357240913172
120	0.8677855800779071	8771.007230863823	14.411777668677829	0.6065822738669273
130	0.9461131220741663	8546.534625048092	15.104769001165842	0.6085303497056861
140	1.0204639550647698	8452.17869313293	15.742960452596863	0.5986425766619422
150	1.0965052470564842	8273.334177702494	16.336860539730196	0.6122987015057609
160	1.1710666220460553	8227.22360412647	16.96433047740871	0.6155380015768278
170	1.2457705810375046	8249.433730702849	17.31601285661983	0.5986001220548078
180	1.3228841640229803	8153.361299876968	17.735422638130615	0.5869171051945036
190	1.4000980390119366	8115.390308026218	18.000810070987363	0.5628028754902441
200	1.4810212480369955	8008.953703709957	18.48296652166439	0.5780319997085015
300	2.2694766650383826	7864.716227283669	22.404435033000098	0.6678863165876615
400	3.0329098740185145	7805.167254259783	22.58350355195061	0.6754551755233579
500	3.7822160820360295	7743.450583002042	22.966030125422186	0.6904744183536446
600	4.5412534980569035	7865.561095728671	22.294250873040102	0.6330837334179605
700	5.289570082037244	7717.566527550956	23.68544784005303	0.701862934997791
800	6.027116790035507	7761.019636381813	22.777283447519952	0.7418297553556962
900	6.7586812900262885	7766.930368349117	22.276495178511894	0.7074732388306058
1000	7.482037332025357	7779.748186711526	22.416894130567066	0.7767105196624577
1100	8.207273707026616	7791.890065200269	23.886977567711046	0.7727509901293816
1200	8.924417248024838	7769.67673036067	23.425727801853718	0.7304467701671189
1300	9.657945165003184	7890.5112901295015	19.442100180990096	0.6950244358724238
1400	10.41649820600287	7871.501487643188	20.651684842336685	0.6957533637791328
1500	11.148728956002742	7615.646383161898	25.75430023499205	0.7528913606114618
1600	11.862543747993186	7720.792063225015	23.53610459938785	0.6721291644397175
1700	12.59042812298867	7805.576315982072	24.221699011551774	0.723825530581459
1800	13.310956581000937	7649.40745807506	25.230318552887194	0.778776096125585
1900	14.028655205998803	7778.663616140591	21.696275973460544	0.6816480033477612
2000	14.77224420598941	7705.580953736973	24.225004823502005	0.7452701594624009
2100	15.500101497978903	7752.430891346984	23.35226450060185	0.7507944308710588
2200	16.21086574799847	7635.825929571018	25.190361582815175	0.775070269956368
2300	16.94130933098495	7771.121316449181	23.658661543672146	0.7163326103177566
2400	17.660661288973643	7647.860030135753	26.43437532040743	0.7707311407507962
2500	18.38265545599279	7713.945414995487	24.174801988709284	0.7607487045087408
2600	19.089490246988134	7579.975105373484	27.052625020180102	0.7975841965892861
2700	19.788468454993563	7568.012162536432	27.768663687253778	0.8022057504672713
2800	20.499084164010128	7622.611800080488	27.969613875466088	0.8032922037738249
2900	21.201529414014658	7570.927329892415	28.017874950881403	0.8034305602010582
3000	21.90538537199609	7607.069311815061	28.019304938770556	0.8034136304664135
3100	22.608946912980173	7586.512462578789	28.077938657850684	0.8049267482971748
3200	23.311048996984027	7643.084929893055	28.059664023475378	0.804604898573818
3300	24.015238204970956	7609.953214924527	28.06527796589579	0.8037830211659431
3400	24.717693829967175	7624.419962092668	28.060429374058273	0.8038102161983628
3500	25.419979829981457	7577.882517782774	28.043648518760182	0.8039929120394936
3600	26.121845496993046	7646.712629692958	28.068345886392212	0.8040518249938594
3700	26.82603891400504	7656.054469269335	28.095440002925752	0.8047721540775693
3800	27.523926955997013	7510.406721564839	28.060989258701877	0.804677278797235
3900	28.226154080970446	7601.828344557941	28.027506082591692	0.8032508480004357
4000	28.92533833096968	7581.912298292904	28.05379622297165	0.8035450143228106
4100	29.629142871970544	7530.283998533004	28.096061192335018	0.8044567191240999
4200	30.332633996964432	7656.934378647605	28.049968522723866	0.8031879007919659
4300	31.038346996967448	7653.212086871685	28.02582289758309	0.8029905166281951
4400	31.74718257994391	7540.288319809373	28.01759616113773	0.8038210857644369
4500	32.44692287192447	7659.183264813418	28.05616701618598	0.8040212836013508
4600	33.14698032994056	7610.805598058812	28.042263763017058	0.8043519256133137
4700	33.84526212193305	7575.006939356938	28.036746282488494	0.8040909801591346
4800	34.543415621941676	7670.283781295864	28.05256772017376	0.8041269400904157
4900	35.24569545593113	7565.443945573589	28.06733421004376	0.8045377979901996
5000	35.94865695593762	7585.598798337533	28.041976279159698	0.8039938508170368
5100	36.64385833093547	7709.127375234199	24.50960143019089	0.7094510703338114
5200	37.35519220592687	7658.149774432704	26.48177158373968	0.7598527084896053
5300	38.06190274792607	7742.067632993954	23.24956241729561	0.681931709876919
5400	38.78211224792176	7716.1781422267895	22.382770513673652	0.632886481669512
5500	39.51845195592614	7913.742073286904	20.705599037960965	0.5681054863557181
5600	40.27400016391766	7720.22217704841	23.127478161232958	0.6894610011026423
5700	41.012527622893685	7801.6852901223465	21.194890987200328	0.716480869301299
5800	41.76354908090434	7828.406828045541	21.702068382977796	0.7544527935229437
5900	42.49736987191136	7777.060502936925	21.734264584144434	0.736214081508787
6000	43.222460746910656	7677.765966439724	24.514012780551994	0.7826458895727514
6100	43.95602057990618	7791.387267199786	23.376090024170555	0.7541308000311356
6200	44.667465579899726	7680.664974355961	27.134107717256878	0.8396760491301912
6300	45.3718635799014	7717.061319060069	25.997154317221426	0.8310954912305509
6400	46.037683871923946	7621.048080955539	27.894672360439706	0.8491646983228088
6500	46.7179864969512	7629.091665758529	26.42799298614296	0.8370169458191061
6600	47.40297162195202	7731.708871737757	24.320007209256353	0.8003130742724529
6700	48.11471095494926	7710.354205764435	24.60846701932997	0.7978521258659961
6800	48.80738828793983	7641.074446132592	27.874753065026866	0.8291985371938344
6900	49.5028448299272	7678.023866103773	27.025541270677955	0.8447787476317818
7000	50.16905378791853	7590.283913172578	30.73596939152722	0.8761603709897844
7100	50.865957745932974	7647.6513126922055	24.175295697632954	0.8177428460490632
7200	51.54533387091942	7635.0508448486735	25.639086644278912	0.8248456482593033
7300	52.21913837091415	7589.136568252853	29.841611816612556	0.8414954446661531
7400	52.856322787934914	7605.811992347608	29.272441879919814	0.8563161569226613
7500	53.52465466293506	7687.8681688042325	28.20332260938512	0.844571966775588
7600	54.17140687193023	7609.435428080351	33.73519039666935	0.859984226689566
7700	54.840452079923125	7616.219546333647	31.062428065600088	0.8373021110528851
7800	55.51693766290555	7628.682965882036	28.334729318753148	0.8134221799471503
7900	56.20739562189556	7830.13359032119	21.255518114432434	0.7773517558355888
8000	56.90076666290406	7802.571223070503	26.338202723183453	0.8233714602159099
8100	57.60952812191681	7753.006599990971	24.26475523421105	0.8242574600045239
8200	58.31158153890283	7635.708333720564	26.748417837748022	0.8343308393487306
8300	58.99051137189963	7624.280142354264	27.670576074143007	0.8438721207599575
8400	59.67515299687511	7686.32360371171	25.832098899967264	0.8038625187541055
8500	60.36724691386917	7660.509355813461	27.95216414346488	0.851386435507331
8600	61.05988349686959	7633.120634491248	27.932110538900794	0.830098736230559
8700	61.85974878887646	7590.301771287001	27.691104956011905	0.8384011971074047
8800	62.65195307985414	7490.1785362291475	30.37862920310014	0.8786085539962765
8900	63.388242204877315	7508.880763883343	31.127140080390404	0.8817201729422717
9000	64.11541057989234	7588.365995961253	31.342185468511445	0.8818393015165621
9100	64.8263862048916	7597.410618367154	31.429867056208707	0.883017462672304
9200	65.54513499588938	7640.4266480526985	31.43250401941845	0.8825890063322326
9300	66.25543587087304	7586.031233503669	31.461872880542	0.8830534105516887
9400	66.96558774585719	7562.0066573993345	31.478672693449422	0.8828081286790433
9500	67.67042362084612	7586.509470272601	31.42167833449961	0.8822838970771084
9600	68.37141157884616	7557.039272817189	31.457856408010663	0.8833273532445818
9700	69.13009249581955	7510.636997546954	31.482019311329957	0.8834098592515948
9800	69.96108212083345	7527.542152379646	31.50324196022173	0.8836778103649127
9900	70.67019437084673	7617.082410684228	31.49403024403778	0.8829037353457182
10000	71.37938737083459	7644.247425320034	31.49353815051584	0.883210938613134
