Also,
The getindex method for the InterpGrid object can accept more than one point. There is no need to evaluate the points using a comprehension: julia> y2 = ig[Xq - 2.0]; julia> y=[ig[x - 2.0] for x in Xq]; julia> maximum(y - y2) 0.0 julia> minimum(y - y2) 0.0