Dear users, I would like to speed up the following code.
function essai(n,s1,s2) a = zeros(Int64,n) ss1::Float64 = 0. ss2::Float64 = 0. @inbounds begin for k=1:n for i=1:length(s1) ss1 = s1[i] for j=1:length(s2) ss2 = s2[j] if ss1>ss2 a[k] += 1 end end end end end end a = rand(1,5000) b = rand(1,5000) essai(1,a,b) @time essai(200,a,b) On my computer, if I remove the *if*, it runs in 2ms and 13s otherwise. Could it be there is a way to speed it up without change the "algorithm"? Thank you for your help