This seems to be one of those “poster child” projects for Julia. As in there are likely lots of custom primitives or numerical routine that’d be impossible to have perform well in Python/Numpy but still saving development effort over coding a custom C++ case (much less C++ expertise needed to do so).