At this point, all three methods are CONSISTENT
This does not necessarily mean that there are no bugs left, but we can
now directly compare results between the Python, vectorized, and Cython
versions of the code [ci skip]