Someone posted a benchmark for Numba and Fortran https://www.reddit.com/r/Python/comments/431tsm/numba_applied_to_high_intensity_computations_a/
As a Friday afternoon project I ported the Python version to Julia https://github.com/abieler/LJSim.jl (please go to the reddit link above to also download the python script to generate the necessary input data if you want to replicate the benchmark) My Julia implementation has about the same performance as the Numba version. (It is slightly slower though) Both are well within a factor of 2 compared to the fortran implementation, which is very good in my mind! If you have suggestions for speed up please share :) I quickly profiled and @code_warntyped and didn't see anything worrying with my code. Have a good weekend everybody!