[rrss]

> if you change their source code you can have them be faster, that's an obvious to me? The question is, what kind of performance can you expect as a user using them as a library, and this benchmark is revealing of that in this case. What is wrong with it?

You don't need to change the numpy or cupy source code, just write your own version of the function in a few lines of code using other primitives that numpy/cupy provides.

Note that is exactly what is done for the neanderthal version.

[didibus]

How am I supposed to know that corrcoef is slow and I should rewrite it? I didn't know before, now thanks to this article I do.