--- a/src/pdps.rs	Thu Feb 26 11:38:43 2026 -0500
+++ b/src/pdps.rs	Thu Feb 26 11:36:22 2026 -0500
@@ -147,7 +147,7 @@
     pdpsconfig: &PDPSConfig<F>,
     iterator: I,
     mut plotter: Plot,
-    μ0 : Option<RNDM<N, F>>,
+    μ0: Option<RNDM<N, F>>,
 ) -> DynResult<RNDM<N, F>>
 where
     F: Float + ToNalgebraRealField,
@@ -206,15 +206,15 @@
         let μ_base = μ.clone();
 
         // Insert and reweigh
-        let (maybe_d, _within_tolerances) = prox_penalty.insert_and_reweigh(
-            &mut μ, &mut τv, &μ_base, None, τ, ε, config, &reg, &state, &mut stats,
-        )?;
+        let (maybe_d, _within_tolerances) = prox_penalty
+            .insert_and_reweigh(&mut μ, &mut τv, τ, ε, config, &reg, &state, &mut stats)?;
 
         // Prune and possibly merge spikes
         if config.merge_now(&state) {
-            stats.merged += prox_penalty
-                .merge_spikes_no_fitness(&mut μ, &mut τv, &μ_base, None, τ, ε, config, &reg);
+            stats.merged +=
+                prox_penalty.merge_spikes_no_fitness(&mut μ, &mut τv, &μ_base, τ, ε, config, &reg);
         }
+        stats.inserted += μ.len() - μ_base.len();
         stats.pruned += prune_with_stats(&mut μ);
 
         // Update step length parameters