# 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_greedy_15_10000_100", dynrange = 1.0, Λ = 1.0)
iter	time	function_value	psnr	ssim
1	9.459006832912564e-6	16984.761585634194	4.935377882355947	0.010713389784152977
2	0.0027665839879773557	16913.30637293789	4.998339471715698	0.020038460127813636
3	0.004400917998282239	16751.874548588195	5.055282677025294	0.030773503136381827
4	0.0074785009783227	16675.59318989853	5.134405200055879	0.04348267882173725
5	0.009071666980162263	16519.627023682413	5.213092574877857	0.05250707337965823
6	0.011895957955857739	16283.880989107562	5.299151092604632	0.06313438010905992
7	0.013695832953089848	15956.7559086914	5.400513910462024	0.07276459073603986
8	0.016906540957279503	15984.286601574102	5.499275613913268	0.08077218823403695
9	0.018733498960500583	15862.034115549892	5.5734809701428585	0.09697670465303475
10	0.02270516596036032	15505.973513050476	5.653702898047572	0.1013115025623551
11	0.024344290955923498	15457.665117831717	5.7491675725363836	0.1169407551820293
12	0.02728391595883295	15244.052573853878	5.860435185270816	0.12255516472196427
13	0.02985308296047151	15067.288572432637	5.922866449282072	0.13061598788737439
14	0.03276779197040014	14902.02122627321	6.031645949725149	0.1433379387176869
15	0.03455779197975062	14669.508552244712	6.152522947001659	0.15063014383534631
16	0.03803112599416636	14644.807368261618	6.21603819836686	0.16736886395539902
17	0.03989208399434574	14533.760754238756	6.312268536142362	0.16907540018380526
18	0.04292820897535421	14384.678516758857	6.416757464309095	0.1762789165922429
19	0.044509708968689665	14182.897167955603	6.47894409921748	0.18836096699180307
20	0.04720804197131656	14023.868552695294	6.543847427858075	0.20448784195041242
30	0.07033574994420633	12880.595060249443	7.459508950299395	0.2781634584727047
40	0.09238099993672222	12068.489015562742	8.256153778039806	0.37002593930122496
50	0.11472004192182794	11128.81971841798	9.094205226938325	0.4201305009805389
60	0.1372836669033859	10558.363284464049	9.841755340051355	0.4286222935145585
70	0.15999600090435706	10084.916766784721	10.664921806472535	0.5085967277125533
80	0.1824061258812435	9632.188346162278	11.438587529886286	0.5021162328825586
90	0.20564962588832714	9185.87847980603	12.241266182626946	0.5596539213224592
100	0.22782233389443718	8902.535295382342	12.97424502922588	0.6061899419412976
110	0.2499515418894589	8812.924290869876	13.655475091728793	0.5992575123150529
120	0.2739960418839473	8661.616220113547	14.41142472921532	0.6077387165768228
130	0.2977702918869909	8430.662622999736	15.10399522447688	0.6093523779938992
140	0.3208194168691989	8340.36694964283	15.745784568180035	0.6006538027357621
150	0.3441482498892583	8159.236554504958	16.338233739474813	0.6139047787612464
160	0.36803204187890515	8116.41787607613	16.968007272086652	0.6179461678249274
170	0.3921203338832129	8131.50472796069	17.324559371508972	0.600737677141733
180	0.41673616686603054	8032.228967796908	17.742890976021734	0.5883617683545839
190	0.44063199986703694	7999.167683045084	18.015894447404705	0.5656991910872208
200	0.46547587486566044	7901.567434435487	18.493142674052464	0.5806969243242441
300	0.6945271658769343	7744.621460586899	22.404752505262316	0.6694631402686483
400	0.9234754998760764	7696.81597259264	22.60408950014969	0.6790517873710666
500	1.1520704168651719	7630.2008052740985	22.980102228934246	0.6925100185875607
600	1.3927676258608699	7756.69258828954	22.32609699883373	0.6369383533210603
700	1.622385292837862	7614.017033178141	23.698982937813543	0.7061012270203996
800	1.8614029598538764	7654.437758992583	22.755057612627386	0.7450735649327236
900	2.1005915848654695	7651.495737251431	22.27639095540639	0.7088548187668683
1000	2.3409649598470423	7672.829010181706	22.427315388898627	0.7810282343560052
1100	2.571770001843106	7679.3724531220805	23.9358912928825	0.7770259834542077
1200	2.807847876829328	7659.447697142759	23.466717385335258	0.7347375741583264
1300	3.0473836678429507	7784.862069192328	19.449758045706083	0.6996499872601374
1400	3.2807555428298656	7758.770230472157	20.65231111770588	0.6977938287876079
1500	3.5184316678496543	7511.422315503292	25.757579246820022	0.7558377846789588
1600	3.7457244598772377	7608.196035937882	23.562523248434903	0.6756810875418341
1700	3.977651709865313	7697.670157743667	24.24426209811344	0.723212156276774
1800	4.213088626856916	7547.637793041002	25.279817533654704	0.7844018484872726
1900	4.451237751869485	7667.920506487475	21.72051921460519	0.6828336218227026
2000	4.689798876846908	7598.746227277767	24.241558266295268	0.749841753394096
2100	4.922756168845808	7645.811341523108	23.36983983920483	0.7543507847196569
2200	5.150594709848519	7533.228792593347	25.231301822335865	0.7796706764201016
2300	5.38284004284651	7654.749020668592	23.686393087174295	0.7203110763089673
2400	5.615964250842808	7541.466072073888	26.465311631930962	0.7710626209228894
2500	5.848620083823334	7608.358721589208	24.196907814250626	0.7645283809720801
2600	6.075843750819331	7507.225162855352	28.127467782523894	0.8207545711576181
2700	6.30493600081536	7495.117097665414	29.324728348553137	0.8301583874190869
2800	6.536273667821661	7550.855809042114	29.683734586072937	0.8324660773963289
2900	6.759037000825629	7497.685231385834	29.769377215163825	0.8334900625709835
3000	6.972873875842197	7534.864430960992	29.759554753554895	0.8325528306710228
3100	7.197159583825851	7512.801082337991	29.849542422120578	0.8343954699788185
3200	7.423133958829567	7568.976972948631	29.823703354094597	0.8339647583523062
3300	7.652460166835226	7536.475714561033	29.827765786478686	0.8326484526386115
3400	7.871725166827673	7550.126132618042	29.84741588148506	0.8339693162678066
3500	8.090837041818304	7506.795761201511	29.7977583465735	0.8336850579828082
3600	8.312458582833642	7572.2284465970415	29.822089586913094	0.8334390307888344
3700	8.526169582823059	7582.170432515587	29.871070017584263	0.8346460003409135
3800	8.743197165837046	7439.094625706412	29.82233442097835	0.8342533872936406
3900	8.966374123847345	7528.041189336874	29.783747970150003	0.8325943393028667
4000	9.198021790856728	7509.590722264806	29.81182076994608	0.8332053763533972
4100	9.419515956862597	7458.17166883916	29.864071109316384	0.8345310229110278
4200	9.642372290865751	7584.323783433058	29.794273532559682	0.8318425013677145
4300	9.863351332867751	7582.441780888057	29.783391508421445	0.8321317870813338
4400	10.085081249853829	7467.954865775342	29.76369826211031	0.8331304531583429
4500	10.300521957862657	7588.595344826685	29.815481580511694	0.8332284369892616
4600	10.521645541855833	7538.7616498531115	29.806839312685028	0.8339655553538801
4700	10.744349125860026	7501.114274612981	29.782025957849072	0.8334443996981304
4800	10.972111875860719	7596.040010093481	29.81471283898813	0.8333600795048038
4900	11.20129612585879	7491.107641665356	29.82440375617413	0.8335393489258202
5000	11.424383292847779	7512.010735555474	29.7950315490028	0.8336169240254758
5100	11.659154209832195	7606.346505955395	24.547897757737758	0.7126314879405296
5200	11.887422251835233	7554.757629823769	26.490087323391343	0.763200425997098
5300	12.124965626833728	7634.839195777455	23.251826692348565	0.6840233743595976
5400	12.359079917834606	7629.750243265272	22.38585060968374	0.6328488100805307
5500	12.582586500851903	7805.177821390203	20.712711549312544	0.5705742044424992
5600	12.818468542856863	7603.905952370702	23.143413289079085	0.6884330594749877
5700	13.05466633386095	7688.578233813721	21.213717009785185	0.721636141504465
5800	13.290357458870858	7719.635211776798	21.71195955129489	0.7570478347851294
5900	13.526611833862262	7676.6816135402105	21.750619127706408	0.7403546375004602
6000	13.758223708864534	7575.080335613914	24.573381519916417	0.7875561199842938
6100	13.997269041836262	7676.284675923737	23.439093458133975	0.757009339674707
6200	14.232671874837251	7578.364923457294	27.231637038386012	0.8455471305468374
6300	14.481987041828688	7613.979370346609	26.07013462848363	0.8359441011382875
6400	14.70771487482125	7519.046861264674	27.997874286985255	0.8545032122601275
6500	14.944584833836416	7528.3952922768185	26.50449677595438	0.8424162952163371
6600	15.177699666848639	7624.93080874274	24.350779383899273	0.803258899485134
6700	15.410235249844845	7613.077355013846	24.621546037280556	0.7919224088874013
6800	15.649307665851666	7531.9505519229115	28.000964755702395	0.836617921369054
6900	15.883560040849261	7576.760593454588	27.05046723720748	0.8486065004053568
7000	16.12480108285672	7497.54148329205	30.78922398830416	0.8826064995701333
7100	16.350105957855703	7553.837092613434	24.16853115179478	0.8219914838713879
7200	16.577384874835843	7538.355087947535	25.648969528966887	0.830198149976945
7300	16.816331999842077	7487.136842246847	29.873573183838023	0.8488374521166463
7400	17.055520124849863	7513.217166823415	29.300587411867088	0.8616637913384613
7500	17.288585708854953	7592.57548697497	28.216167493337903	0.8507529113053007
7600	17.528063500853023	7515.16736000597	33.77506940517368	0.8654520584919614
7700	17.759960917843273	7524.851477228005	31.073310800919533	0.8415409435023122
7800	17.992375834845006	7520.021629120392	28.32135266766251	0.8169240420477784
7900	18.230295167857548	7729.597545694648	21.26367207111414	0.7793687715483483
8000	18.45639720984036	7704.974258865148	26.378426936477663	0.8253986411804386
8100	18.68567287584301	7651.753488931303	24.30611111149455	0.82978453686767
8200	18.92381733385264	7534.056869442974	26.8491896478354	0.8406205162350359
8300	19.157327583874576	7523.38132614155	27.7902514166489	0.8489746192796985
8400	19.395289333886467	7575.363440646889	25.906706169408498	0.808142994842427
8500	19.63425841688877	7558.923700348188	28.063546182832262	0.8577812486485874
8600	19.872082958900137	7526.868528131662	28.0533371061733	0.8338569528062263
8700	20.104079666896723	7486.049333603622	27.826949037583827	0.8450839363743374
8800	20.32895895888214	7424.8996538189	31.519319399570023	0.8818385116225501
8900	20.54868645890383	7443.628293714584	32.73464134397366	0.8868612227294423
9000	20.7744442919211	7524.891511050345	33.083683954022504	0.8874964247053455
9100	21.001210333924973	7530.054099647805	33.212080535295605	0.8884177095532856
9200	21.226242750912206	7574.218966637169	33.2519638533287	0.8889233763902179
9300	21.450743792927824	7519.540326693112	33.2767698414882	0.8890479553394782
9400	21.67771950090537	7497.476207127806	33.30463126513708	0.8887145333590175
9500	21.906040000903886	7520.959849162117	33.24698797583543	0.8881604407095081
9600	22.1302481679013	7489.83517071967	33.2957869569424	0.8894888318274736
9700	22.353870833903784	7444.416504670087	33.31722158513732	0.889383044141604
9800	22.583873624913394	7461.2693002844435	33.33261897821747	0.8898882831980676
9900	22.81429399989429	7550.116901014406	33.32183813719376	0.889096148465819
10000	23.031792416906683	7576.14565024196	33.31009514882632	0.888683273133075
