Mercurial > repos > pointsource_algs / file comparison
comparison: src/sliding_pdps.rs
src/sliding_pdps.rs
- branch
- dev
- changeset 36
- fb911f72e698
- parent 35
- b087e3eab191
- child 37
- c5d8bd1a7728
equal
deleted
inserted
replaced
| 35:b087e3eab191 | 36:fb911f72e698 |
|---|---|
| 212 let τ = config.τ0 * φ / ( σ_d * bigM * a + φ * l ); | 212 let τ = config.τ0 * φ / ( σ_d * bigM * a + φ * l ); |
| 213 let ψ = 1.0 - τ * l; | 213 let ψ = 1.0 - τ * l; |
| 214 let β = σ_p * config.σd0 * nKz / a; // σ_p * σ_d * (nKz * nK_z) / a; | 214 let β = σ_p * config.σd0 * nKz / a; // σ_p * σ_d * (nKz * nK_z) / a; |
| 215 assert!(β < 1.0); | 215 assert!(β < 1.0); |
| 216 // Now we need κ‖K_μ(π_♯^1 - π_♯^0)γ‖^2 ≤ (1/θ - τ[ℓ_v + ℓ]) ∫ c_2 dγ for κ defined as: | 216 // Now we need κ‖K_μ(π_♯^1 - π_♯^0)γ‖^2 ≤ (1/θ - τ[ℓ_v + ℓ]) ∫ c_2 dγ for κ defined as: |
| 217 let κ = σ_d * ψ / ((1.0 - β) * ψ - τ * σ_d * bigM); | 217 let κ = τ * σ_d * ψ / ((1.0 - β) * ψ - τ * σ_d * bigM); |
| 218 // The factor two in the manuscript disappears due to the definition of 𝚹 being | 218 // The factor two in the manuscript disappears due to the definition of 𝚹 being |
| 219 // for ‖x-y‖₂² instead of c_2(x, y)=‖x-y‖₂²/2. | 219 // for ‖x-y‖₂² instead of c_2(x, y)=‖x-y‖₂²/2. |
| 220 let calculate_θ = |ℓ_v, max_transport| { | 220 let calculate_θ = |ℓ_v, max_transport| { |
| 221 config.transport.θ0 / (τ*(ℓ + ℓ_v) + κ * bigθ * max_transport) | 221 config.transport.θ0 / (τ*(ℓ + ℓ_v) + κ * bigθ * max_transport) |
| 222 }; | 222 }; |
