Weights and indices for linear interpolation/extrapolation

[stevengj]

@fredrikpaues

Implement binary search as
        #     `index = searchsortedlast(gridpoints, value)`
        # but with support for tuple gridpoints

I’m confused, why can’t you use searchsortedlast? You almost certainly don’t want to use a tuple to store a long array of grid points.