I forgot to check earlier, np.percentile is available in our minimum numpy version.
np.percentile is vectorized
`np.percentile(10*np.arange(6), [0.1, 0.5, 0.75])`
which is faster because it doesn't do repeated sorting (starting in 1.8 it will do only partial sorting which is even faster)