Update interpolation.py
fixes lowerbound interpolation for non-monotonic functions
@@ -1041,7 +1041,7 @@ class CubicInterp(HARKinterpolator1D):
- self.coeffs[self.n, 2] * np.exp(alpha * self.coeffs[self.n, 3])
)
- y[x == self.x_list.min()] = self.y_list.min()
+ y[x == self.x_list[0]] = self.y_list[0]
return y
