In Speeding up isotonic regression in scikit-learn, we dropped down into Cython to improve the performance of a regression algorithm. I thought it would be interesting to compare the performance of this (optimized) code in Python against the naive Julia implementation. This article continues on from the previous one, so it may be worth reading that before continuing here to obtain the necessary ba