[shele]

The Julia-code has some performance flaws in the hot loop (like making a superfluous copy of a matrix before just reading it.)