Merge d814a0d into 566dc79

src/linear_algebra.jl

@@ -22,21 +22,20 @@ function mutable_operate_to!(C::Vector, ::typeof(*), A::AbstractMatrix, B::Abstr
 
     Astride = size(A, 1)
     @inbounds begin
-    for i = 1:mA
-        C[i] = zero!(C[i])
-    end
+        for i = 1:mA
+            C[i] = zero!(C[i])
+        end
 
-    # We need a buffer to hold the intermediate multiplication.
-    mul_buffer = zero(promote_operation(*, eltype(A), eltype(B)))
-    @inbounds for k = Base.OneTo(mB)
-        aoffs = (k-1) * Astride
-        b = B[k]
-        for i = Base.OneTo(mA)
-            # `C[i] = muladd_buf!(mul_buffer, C[i], A[aoffs + i], b)`
-            mutable_buffered_operate!(mul_buffer, add_mul, C[i], A[aoffs + i], b)
+        # We need a buffer to hold the intermediate multiplication.
+        mul_buffer = zero(promote_operation(*, eltype(A), eltype(B)))
+        for k = Base.OneTo(mB)
+            aoffs = (k-1) * Astride
+            b = B[k]
+            for i = Base.OneTo(mA)
+                mutable_buffered_operate!(mul_buffer, add_mul, C[i], A[aoffs + i], b)
+            end
         end
-    end
-    end
+    end # @inbounds
     return C
 end