# HG changeset patch
# User Tuomo Valkonen <tuomov@iki.fi>
# Date 1772128538 18000
# Node ID ec5f160c413ba69d29b3b7845c26172999f7d6af
# Parent  7a8a55fd41c039bc633e9e6ad53b5558fcc2999d
Oops, a τ had dropped in forward_pdps.

diff -r 7a8a55fd41c0 -r ec5f160c413b src/forward_pdps.rs
--- a/src/forward_pdps.rs	Thu Feb 26 11:36:22 2026 -0500
+++ b/src/forward_pdps.rs	Thu Feb 26 12:55:38 2026 -0500
@@ -199,7 +199,7 @@
     // Run the algorithm
     for state in iterator.iter_init(|| full_stats(&μ, &z, ε, stats.clone())) {
         // Calculate initial transport
-        let Pair(mut τv, τz) = f.differential(Pair(&μ, &z));
+        let Pair(mut τv, τz) = f.differential(Pair(&μ, &z)) * τ;
         let μ_base = μ.clone();
 
         // Construct μ^{k+1} by solving finite-dimensional subproblems and insert new spikes.