I see that c is a constant array of Ints, and its elements multiply ux, uy and uz in a loop, where ux, uy and uz are arrays of floats, so there's a type stability problem.
On Sunday, May 8, 2016 at 9:18:09 AM UTC+2, feza wrote: > > https://gist.github.com/musmo/27436a340b41c01d51d557a655276783 > > On Sunday, May 8, 2016 at 3:17:39 AM UTC-4, feza wrote: >> >> I have read the performance section and believe I have followed all the >> suggested guidelines >> >> The same matlab script takes less than 3 seconds, julia 0.45 9.7 seconds >> (julia 0.5 is even worse...) >> >> <script src=" >> https://gist.github.com/musmo/27436a340b41c01d51d557a655276783.js >> "></script> >> >>