Correct MulAddMul in triangular * Diagonal __muldiag_nonzeroalpha!

JuliaLang · Oct 28, 2024 · a6557e2 · a6557e2
1 parent 14c01fa
commit a6557e2
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/stdlib/LinearAlgebra/src/diagonal.jl b/stdlib/LinearAlgebra/src/diagonal.jl
@@ -489,7 +489,7 @@ function __muldiag_nonzeroalpha!(out, A::UpperOrLowerTriangular, D::Diagonal, _a
         if !_has_matching_zeros(out, A)
             rowrange = _rowrange_tri_zeros(A, j)
             @inbounds @simd for i in rowrange
-                _modify!(_add, A[i,j] * dja, out, (i,j))
+                _modify!(_add_aisone, A[i,j] * dja, out, (i,j))
             end
         end
     end