Mercurial > repos > pointsource_algs / file diff

--- a/src/sliding_pdps.rs	Thu Jan 23 23:48:52 2025 +0100
+++ b/src/sliding_pdps.rs	Sun Jan 26 11:29:57 2025 -0500
@@ -162,7 +162,7 @@
 
     // Set up parameters
     // TODO: maybe this PairNorm doesn't make sense here?
-    let opAnorm = opA.opnorm_bound(PairNorm(Radon, L2, L2), L2);
+    // let opAnorm = opA.opnorm_bound(PairNorm(Radon, L2, L2), L2);
     let bigθ = 0.0; //opKμ.transport_lipschitz_factor(L2Squared);
     let bigM = 0.0; //opKμ.adjoint_product_bound(&op𝒟).unwrap().sqrt();
     let nKz = opKz.opnorm_bound(L2, L2);
@@ -245,9 +245,7 @@
         // TODO: Write a version of initial_transport that can deal with K_μ ≠ 0.
  
         let (μ_base_masses, mut μ_base_minus_γ0) = initial_transport(
-            &mut γ1, &mut μ, |ν| opA.apply(Pair(ν, &z)),
-            ε, τ, &mut θ_or_adaptive, opAnorm,
-            v, &config.transport,
+            &mut γ1, &mut μ, τ, &mut θ_or_adaptive, v,
         );
 
         // Solve finite-dimensional subproblem several times until the dual variable for the