Merge pull request #12689 from KristofferC/patch-4

inbounds and hoist field loading for dense * sparse
JuliaLang · Aug 19, 2015 · 9c82059 · 9c82059
2 parents 3ce3f40 + 9eaa2ce
commit 9c82059
Showing 1 changed file with 4 additions and 2 deletions.
diff --git a/base/sparse/linalg.jl b/base/sparse/linalg.jl
@@ -94,8 +94,10 @@ function (*){TX,TvA,TiA}(X::StridedMatrix{TX}, A::SparseMatrixCSC{TvA,TiA})
     mX, nX = size(X)
     nX == A.m || throw(DimensionMismatch())
     Y = zeros(promote_type(TX,TvA), mX, A.n)
-    for multivec_row=1:mX, col = 1:A.n, k=A.colptr[col]:(A.colptr[col+1]-1)
-        Y[multivec_row, col] += X[multivec_row, A.rowval[k]] * A.nzval[k]
+    rowval = A.rowval
+    nzval = A.nzval
+    @inbounds for multivec_row=1:mX, col = 1:A.n, k=A.colptr[col]:(A.colptr[col+1]-1)
+        Y[multivec_row, col] += X[multivec_row, rowval[k]] * nzval[k]
     end
     Y
 end