would be nice to have a "mapper" function

[doutriaux1]

given a list of intervals and a value,return the index of the interval where the value falls.

I'm sure it already exsits somewhere in numpy.