The ResponseSurface and PolyharmonicSpline essentially do the exact same thing with different basis. I suggest we refactor them into a single type that is passed the desired basis. This can include other bases like a Gaussian radial basis.

We can then use these in the IntervalPredictorModel.