Mercurial > repos > pointsource_algs / changeset

--- a/src/sliding_pdps.rs	Tue Dec 31 09:25:45 2024 -0500
+++ b/src/sliding_pdps.rs	Mon Jan 06 11:32:57 2025 -0500
@@ -214,7 +214,7 @@
     let β = σ_p * config.σd0 * nKz / a; // σ_p * σ_d * (nKz * nK_z) / a;
     assert!(β < 1.0);
     // Now we need κ‖K_μ(π_♯^1 - π_♯^0)γ‖^2 ≤ (1/θ - τ[ℓ_v + ℓ]) ∫ c_2 dγ for κ defined as:
-    let κ = σ_d * ψ / ((1.0 - β) * ψ - τ * σ_d * bigM);
+    let κ = τ * σ_d * ψ / ((1.0 - β) * ψ - τ * σ_d * bigM);
     //  The factor two in the manuscript disappears due to the definition of 𝚹 being
     // for ‖x-y‖₂² instead of c_2(x, y)=‖x-y‖₂²/2.
     let calculate_θ = |ℓ_v, max_transport| {