for me personally that is not a good enough reason to do it the same in numpy, especially as implemented in this PR the multidimensional result is inconsistent.

by just skipping (or implementing differently) the limit parameter it scipys scoreatpercentile can still be reduced to two lines.