Fix diagonal operator spawning

joachimbrand · Aug 4, 2023 · 8ac9b08 · 8ac9b08
1 parent e3730f8
commit 8ac9b08
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/StochasticStyles/spawning.jl b/src/StochasticStyles/spawning.jl
@@ -173,7 +173,7 @@ end
 
 @inline function spawn!(s::Exact, w, offdiags::AbstractVector, add, val, boost=1)
     T = valtype(w)
-    spawns = sum(offdiags) do (new_add, mat_elem)
+    spawns = sum(offdiags; init=zero(T)) do (new_add, mat_elem)
         abs(projected_deposit!(
             w, new_add, val * mat_elem, add => val, s.threshold
         ))