Except Julia vectorised is not as fast as Julia looped.

so this :

    for i = 1:size(a)[1]
        a[i] = 2a[i] + 3b[i]
    end

is faster than

    a = 2a .+ 3b

(although I haven't benchmarked it)

Here's someone that did : http://aflyax.github.io/de_vectorization-runtime/

but in 0.5 there are speedups, apparently.

Here's a stack overflow discussion on it

http://stackoverflow.com/questions/27928502/julia-vectorized...