add max_points_used to evaluate expression from a context using only a subset of all points contained in the context
work on error estimation from eim approximation