| 240 // This is much easier with K_μ = 0, which is the only reason why are enforcing it. |
240 // This is much easier with K_μ = 0, which is the only reason why are enforcing it. |
| 241 // TODO: Write a version of initial_transport that can deal with K_μ ≠ 0. |
241 // TODO: Write a version of initial_transport that can deal with K_μ ≠ 0. |
| 242 |
242 |
| 243 //dbg!(&μ); |
243 //dbg!(&μ); |
| 244 |
244 |
| 245 γ.initial_transport(&μ, τ, &mut θ_or_adaptive, v); |
245 γ.initial_transport(&μ, τ, &mut θ_or_adaptive, v, &config.transport); |
| 246 |
246 |
| 247 let mut attempts = 0; |
247 let mut attempts = 0; |
| 248 |
248 |
| 249 // Solve finite-dimensional subproblem several times until the dual variable for the |
249 // Solve finite-dimensional subproblem several times until the dual variable for the |
| 250 // regularisation term conforms to the assumptions made for the transport above. |
250 // regularisation term conforms to the assumptions made for the transport above. |