Reference Lanczos doesn't overwrite v

SunnySuite · Jan 6, 2025 · d04ece1 · d04ece1
1 parent c475263
commit d04ece1
Showing 1 changed file with 1 addition and 0 deletions.
diff --git a/src/KPM/Lanczos.jl b/src/KPM/Lanczos.jl
@@ -93,6 +93,7 @@ function lanczos(mulA!, mulS!, v; min_iters, resolution=Inf, lhs=zeros(length(v)
     βs = Float64[]
     lhs_adj_Q = Vector{ComplexF64}[]
 
+    v = copy(v)
     vp  = zero(v)
     Sv  = zero(v)
     Svp = zero(v)