You can loop on i in the inner loop, it helps quite a bit: @inbounds for j = 1:size(b,2)
for k = 1:size(b,1)
@simd for i =1:size(a,1)
c[i,j] += a[i,k]*b[k,j]
end
end
end
You can loop on i in the inner loop, it helps quite a bit: @inbounds for j = 1:size(b,2)
for k = 1:size(b,1)
@simd for i =1:size(a,1)
c[i,j] += a[i,k]*b[k,j]
end
end
end