Replace mul! by broadcasting in exp! (JuliaLang#54493)

lazarusA · Jul 12, 2024 · e43989d · e43989d
1 parent 45fd21a
commit e43989d
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/stdlib/LinearAlgebra/src/dense.jl b/stdlib/LinearAlgebra/src/dense.jl
@@ -705,8 +705,8 @@ function exp!(A::StridedMatrix{T}) where T<:BlasFloat
         V = mul!(C[3]*P, true, C[1]*I, true, true) #V = C[1]*I + C[3]*P
         for k in 2:(div(length(C), 2) - 1)
             P *= A2
-            mul!(U, C[2k + 2], P, true, true) # U += C[2k+2]*P
-            mul!(V, C[2k + 1], P, true, true) # V += C[2k+1]*P
+            U .+= C[2k + 2] .* P
+            V .+= C[2k + 1] .* P
         end
 
         U = A * U