Skip deepcopy in addmul! in case of aliasing

thofma · Oct 25, 2024 · 9ddf26e · 9ddf26e
1 parent ea21c55
commit 9ddf26e
Showing 1 changed file with 6 additions and 0 deletions.
diff --git a/src/Sparse/Row.jl b/src/Sparse/Row.jl
@@ -839,10 +839,16 @@ function mul!(z::SRow{T}, c, y::SRow{T}) where T
 end
 
 function addmul!(z::SRow{T}, x::SRow{T}, y) where T
+  if z === x
+    return scale_row_right!(x, y+1)
+  end
   return add_right_scaled_row!(x, z, y)
 end
 
 function addmul!(z::SRow{T}, x, y::SRow{T}) where T
+  if z === x
+    return scale_row_left!(y, x+1)
+  end
   return add_left_scaled_row!(y, z, x)
 end