Fix preconditioning with ldiv=true in BiLQ and QMR

JuliaSmoothOptimizers · May 20, 2024 · b7c6774 · b7c6774
1 parent 9f25aec
commit b7c6774
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/src/bilq.jl b/src/bilq.jl
@@ -172,7 +172,7 @@ kwargs_bilq = (:c, :transfer_to_bicg, :M, :N, :ldiv, :atol, :rtol, :itmax, :time
       @kaxpby!(n, one(FC), b, -one(FC), r₀)
     end
     if !MisI
-      mul!(p, M, r₀)
+      mulorldiv!(p, M, r₀, ldiv)
       r₀ = p
     end
 

diff --git a/src/qmr.jl b/src/qmr.jl
@@ -178,7 +178,7 @@ kwargs_qmr = (:c, :M, :N, :ldiv, :atol, :rtol, :itmax, :timemax, :verbose, :hist
       @kaxpby!(n, one(FC), b, -one(FC), r₀)
     end
     if !MisI
-      mul!(p, M, r₀)
+      mulorldiv!(p, M, r₀, ldiv)
       r₀ = p
     end