From d89c432838b6504e42a64e4013ddb42a0cdfe648 Mon Sep 17 00:00:00 2001
From: Alexis Montoison <alexis.montoison@polymtl.ca>
Date: Tue, 12 Sep 2023 14:52:51 -0500
Subject: [PATCH] Improve computation

---
 src/block_krylov_processes.jl | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/src/block_krylov_processes.jl b/src/block_krylov_processes.jl
index 9263ab552..7351c4899 100644
--- a/src/block_krylov_processes.jl
+++ b/src/block_krylov_processes.jl
@@ -22,7 +22,7 @@ function gs(A::AbstractMatrix{FC}) where FC <: FloatOrComplex
       @kaxpy!(n, -R[i,j], qᵢ, qⱼ)  # qⱼ = qⱼ - rᵢⱼqᵢ
     end
     R[j,j] = @knrm2(n, qⱼ)  # rⱼⱼ = ‖qⱼ‖
-    qⱼ .= qⱼ ./ R[j,j]      # qⱼ = qⱼ / rⱼⱼ
+    qⱼ ./= R[j,j]           # qⱼ = qⱼ / rⱼⱼ
   end
   return Q, R
 end
@@ -44,7 +44,7 @@ function mgs(A::AbstractMatrix{FC}) where FC <: FloatOrComplex
   for i = 1:k
     qᵢ = view(Q,:,i)
     R[i,i] = @knrm2(n, qᵢ)  # rᵢᵢ = ‖qᵢ‖
-    qᵢ .= qᵢ ./ R[i,i]      # qᵢ = qᵢ / rᵢᵢ
+    qᵢ ./= R[i,i]           # qᵢ = qᵢ / rᵢᵢ
     for j = i+1:k
       qⱼ = view(Q,:,j)
       R[i,j] = @kdot(n, qᵢ, qⱼ)    # rᵢⱼ = ⟨qᵢ , qⱼ⟩