This function implements the interpolation search algorithm, that starts from a given feasible point and looks in the given direction for the root of the violation curve.
gs.interpolation(evaluations, dir, iter.nbr, hp)
The next point to evaluate.
Matrix of evaluated test statistics and critical values.
Search direction.
Iteration number.
List of hyperparameters.