# 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.15, ρ̃₀ = 1.0, σ̃₀ = 1.0, δ = 0.9, σ₀ = 1.0, τ₀ = 0.01, save_prefix = "img/lighthouse200x300_pdps_known_rotation_15_10000_100", dynrange = 1.0, Λ = 1.0)
iter	time	function_value	psnr	ssim
1	1.158300437964499e-5	16984.761585634194	4.935377882355947	0.010713389784152977
2	0.004544499999610707	16913.18898930401	4.998340859594887	0.02003995164343634
3	0.008661250001750886	16750.193803920756	5.055288302008346	0.03079688152893856
4	0.014498083997750655	16671.121063485978	5.134428792845256	0.043574713333592
5	0.019589208997786045	16516.463176692465	5.2131053745546465	0.05259420146766833
6	0.02791233401512727	16271.937825809135	5.299457076455073	0.06366235779386183
7	0.03420216700760648	15945.210944968887	5.4008452669668605	0.07328524208133452
8	0.04249612599960528	15974.181367505127	5.499620661364012	0.08125474959430541
9	0.04726033398765139	15858.445837265735	5.573792168708679	0.09722364809866409
10	0.05276854199473746	15503.39117557451	5.654036886849081	0.10149296859783202
11	0.05727383401244879	15454.61012967366	5.7495501901365795	0.11720193284467302
12	0.06327795900870115	15239.939731625193	5.860857413649669	0.12294234651272244
13	0.07100850000279024	15059.040094573	5.9234178858606175	0.13136065627873253
14	0.07828570899437182	14891.562142176792	6.03211256308824	0.14425750978442975
15	0.08526483399327844	14664.861022472149	6.152844068588522	0.15093438857769398
16	0.09174675098620355	14641.178398699278	6.21646557885688	0.1676228676829198
17	0.09678504199837334	14529.909939062802	6.31286379644095	0.1694373998033091
18	0.10385145799955353	14379.705098056555	6.4172091469064805	0.17670629832741128
19	0.1087840830150526	14179.453368008217	6.479520197657715	0.18872684684363394
20	0.11608445801539347	14020.848105462535	6.5445567863701415	0.20474623122689758
30	0.18559258300228976	12871.947746376674	7.461218626935155	0.28007571114924296
40	0.25807166600134224	12060.288474931844	8.259252278318392	0.37338219579177834
50	0.3265702080097981	11122.276492601344	9.099572090485927	0.42488574827384507
60	0.39536050002789125	10550.270315586282	9.849158830871863	0.43426822276460514
70	0.4645764590532053	10080.060984796544	10.673770295353753	0.5131519326020315
80	0.5362553340673912	9626.534298702694	11.45027937859924	0.5081658075098451
90	0.6114502920827363	9183.233269141487	12.255799110936499	0.5644123621891661
100	0.6808032090775669	8898.236219351487	12.989862807095253	0.6103747461746691
110	0.7479969590785913	8809.951001932704	13.67371323918665	0.6031451615463254
120	0.8172319590521511	8655.101365187198	14.435264239369562	0.613849601581059
130	0.8885022090689745	8422.620311449538	15.13092155241641	0.616054959486306
140	0.9546161250618752	8336.41548129948	15.773278297029155	0.6063176961436858
150	1.025674291071482	8153.405878967763	16.373128127399685	0.6203923581372525
160	1.0964637490687892	8112.835874755481	17.00323299214496	0.6230752602336734
170	1.1670600400539115	8125.618209675765	17.35866400319437	0.6061958342907529
180	1.2390473310369998	8029.364254689228	17.78630083944302	0.5949506209377361
190	1.3132644560246263	7994.885627673553	18.06344125308478	0.5718875725311316
200	1.3959701220155694	7899.7964941151595	18.546335188073314	0.585502759924342
300	2.1401939970091917	7740.584377754145	22.52005713644982	0.675358561272889
400	2.862409455003217	7695.582884385696	22.6816696892138	0.6831752813596924
500	3.5613134970190004	7625.610711925934	23.077632175339083	0.6991805129064713
600	4.269134539004881	7754.982198379753	22.403300052815393	0.6418902787118214
700	4.973046997009078	7611.83610185616	23.801131359836205	0.7113256628240129
800	5.684840663016075	7650.488589214438	22.8698319727053	0.7504247246884851
900	6.3655843710002955	7646.509402356019	22.359241698611658	0.716466590880605
1000	7.0389545380021445	7668.886121640976	22.499936169658188	0.7853196626102886
1100	7.710475829982897	7673.451893024856	24.0368183386424	0.7828333033288555
1200	8.618623746966477	7655.715964911601	23.568287590306376	0.7390551983573022
1300	9.456737246946432	7784.653621026143	19.456528654944872	0.7030213859686036
1400	10.175538413925096	7753.379311440898	20.702903460303673	0.7036526505404448
1500	10.864930871903198	7507.315937184592	25.929679171191363	0.7616761375968653
1600	11.56571253889706	7603.154306937039	23.66237538886828	0.6815647175140157
1700	12.26827358090668	7683.057115332622	24.389110422338604	0.7348239981336342
1800	12.953262372902827	7546.5661948313245	25.390888822762697	0.788160811522912
1900	13.650335039885249	7662.961502436569	21.760714631784843	0.6883832313453777
2000	14.332668331888271	7595.559557308641	24.361751322668983	0.7557768556634176
2100	15.0158604978933	7642.641510804748	23.44259446268145	0.759704765065882
2200	15.812293080904055	7529.946423293651	25.337402283282948	0.7837707411232193
2300	16.529103621898685	7650.293092387081	23.782945135291826	0.7267884490047118
2400	17.24880120489979	7531.350293034386	26.627333993596952	0.7789931463732542
2500	18.236935371911386	7604.937350810198	24.254093954731324	0.7688797243041007
2600	18.98574237190769	7506.979159490371	28.163287861835187	0.8213692548912471
2700	19.692724830907537	7495.132521588308	29.338972573578445	0.8303435550424318
2800	20.418023830919992	7550.861495358559	29.688668681401357	0.8325261142712462
2900	21.151178455911577	7497.681378049184	29.77113674490053	0.833510219676084
3000	22.469975080894073	7534.864279839007	29.760176356969055	0.8325596838218031
3100	23.242205705901142	7512.8016621416955	29.849760190752402	0.8343977637627975
3200	23.96700045588659	7568.977120466229	29.82377919540966	0.8339655492645061
3300	24.685648372891592	7536.47575250437	29.82779250651126	0.8326487245505965
3400	25.410894205910154	7550.12613902099	29.847425664276194	0.8339694125858834
3500	26.134183455898892	7506.795761958375	29.79776179508788	0.8336850913317675
3600	26.84944883090793	7572.2284505058515	29.822090778266116	0.8334390425971695
3700	27.579378205904504	7582.1704330241855	29.871070442641503	0.8346460044610935
3800	28.292889121890767	7439.094626072915	29.82233457823815	0.8342533887885921
3900	29.005965621909127	7528.041189369876	29.78374802783843	0.8325943398371105
4000	29.705760246899445	7509.590722262952	29.811820790718734	0.8332053765410026
4100	30.4276719969057	7458.171668837847	29.864071116486485	0.8345310229759666
4200	31.129338496917626	7584.3237834324555	29.794273535196947	0.831842501391206
4300	31.8562886629079	7582.441780888883	29.783391509370336	0.8321317870898365
4400	32.55630512192147	7467.954865776653	29.763698262445182	0.8331304531614087
4500	33.26248699691496	7588.595344826527	29.815481580634405	0.8332284369903453
4600	33.98749416292412	7538.761649853208	29.806839312729867	0.8339655553542691
4700	34.71784066292457	7501.114274613001	29.782025957865628	0.8334443996982707
4800	35.445598745922325	7596.040010093486	29.81471283899427	0.8333600795048546
4900	36.15405020391336	7491.107641665366	29.82440375617628	0.8335393489258388
5000	36.88177916192217	7512.010735555482	29.795031549003607	0.8336169240254815
5100	37.553861577936914	7601.731823796803	24.639734206497295	0.718731217377206
5200	38.212563535926165	7550.010610191385	26.61879096064414	0.7692474935887758
5300	38.883532618929166	7628.710669576085	23.338469364047487	0.6916891649268602
5400	39.57633686895133	7628.12346561555	22.47477835923926	0.638533817330561
5500	40.263060326949926	7801.965719204189	20.799948537751334	0.5759658265441671
5600	40.97751503495965	7590.406281269364	23.272389689544376	0.7004891061299592
5700	41.69070515997009	7686.640032197544	21.286703262117506	0.7249615495895064
5800	42.390688867948484	7715.908051437104	21.772581186282732	0.7622136205232116
5900	43.07864670193521	7677.099181847955	21.80427298856596	0.743074175589619
6000	43.76831682692864	7573.898277157426	24.68497084974471	0.7894670390122641
6100	44.466636034921976	7670.000686942537	23.509744642920406	0.7633259721080018
6200	45.143922784918686	7577.16436369657	27.32180112366382	0.8482051710105608
6300	45.79350349292508	7610.5301635939595	26.156708159483948	0.8396642143201595
6400	46.41413886792725	7516.608587698038	28.109309748526748	0.8578501417257888
6500	47.0422783679096	7526.628236196506	26.61345359182814	0.8451422101286203
6600	47.66718915992533	7617.665484831259	24.412500471573534	0.8095191043759433
6700	48.32594545194297	7592.001954983449	24.67670711495551	0.8065243710810321
6800	48.98214024392655	7531.233189425531	28.13211188152673	0.839074880063195
6900	49.63292824392556	7571.139238522656	27.163933757288206	0.8543925739075514
7000	50.25030528593925	7496.563887401963	30.956149582928155	0.884712322939132
7100	50.896915077930316	7551.366553567551	24.244664321093147	0.8260577388043638
7200	51.51143224394764	7534.571709550444	25.756620870381255	0.8344890876806853
7300	52.131177827948704	7485.960519148183	30.076868751489837	0.8511779879818084
7400	52.716625452943845	7509.252296110658	29.365195117396894	0.865193489458899
7500	53.32932041093591	7590.377196099804	28.292229004600685	0.853547915409306
7600	53.93155649493565	7511.920612134252	33.998559464748375	0.8692588979560306
7700	54.55736841194448	7520.4289745423375	31.21873851121314	0.8460522801868691
7800	55.173953369958326	7515.01843521902	28.48016880924737	0.8232343078425448
7900	55.82728507797583	7721.380953219343	21.31383065319813	0.785493461480182
8000	56.45445670298068	7695.563799071023	26.512129585646598	0.8321527447930043
8100	57.099798118986655	7649.981373505653	24.37150784924491	0.8330535929121899
8200	57.76061307798955	7532.126977190399	26.948906460515303	0.843123783252039
8300	58.38335895299679	7520.173879660432	27.94082916917048	0.8526875900285286
8400	59.037121661007404	7568.537250784265	26.000420865404866	0.8132389288311548
8500	59.690731120004784	7556.984142310944	28.182215189603905	0.8599580091419193
8600	60.33046382802422	7519.835801920714	28.20112401935014	0.8397698933960218
8700	60.949555327999406	7484.455131969999	27.94604407057809	0.8474665442397825
8800	61.65961128601339	7424.784368236039	31.580729098513117	0.8821582644915212
8900	62.36635357703199	7443.628080317521	32.75742362999206	0.8869630853595883
9000	63.06816161904135	7524.912001889027	33.09153311089601	0.8875289565822205
9100	63.75655032804934	7530.055202891102	33.21476256733631	0.8884281788172025
9200	64.45943362003891	7574.220896121868	33.252934067227194	0.8889270675816986
9300	65.14593112003058	7519.54053286266	33.27710940188125	0.8890492259587879
9400	65.84583432803629	7497.476373314503	33.30475288461811	0.8887149741862514
9500	66.55636366203544	7520.959882347643	33.24703235288527	0.8881605973020652
9600	67.2616211620334	7489.835206310586	33.29580354370719	0.8894888876893868
9700	67.94827216205886	7444.416495291221	33.3172274862889	0.889383063864621
9800	68.635197412048	7461.26930177474	33.33262109262642	0.889888290143119
9900	69.34838907804806	7550.116901443209	33.321838888452206	0.8890961509296528
10000	70.06482066205353	7576.1456499203905	33.310095426055014	0.8886832740202457
