I converted the authors c version of the code to Java and ran it through JMH, on my machine it runs about 30% slower. (This is before he added the metaprogramming that essentially precalculates some of the matrices at compile-time).

That's also without making any attempt at optimizing, basically a straight conversion that took ten minutes.

It's not the speed of C but it's pretty close. And it has memory safety.