Bug in binned.py

[floern1]

Line 147 in binned.py uses inside = (knot_y[..., 0] < y) & (y <= knot_x[..., -1]) but should be inside = (knot_y[..., 0] < y) & (y <= knot_y[..., -1]) instead. This results sometimes in crashes in the line 173 when rev=True. The reason is that np.searchsorted cannot find the value and index=upper is out of bounds.