3 ms·

So following the numpy example the author gives one using scipy. This should be the same if not faster than a pure numpy implementation