How to implement a fast linear interpolator in Python for this case?

Question

I have a precomputed array of x_val and its corresponding y_val. These are 1-d array of length about 5000. I have another multi-dimens