N-D

usage:

std::vector<NumericVector> grid = {x_pts, y_pts};
NumericVector values(4);  // Fortran order
values << 0, 1, 1, 2;     // z(x,y) = x + y
 
janus::Interpolator interp(grid, values);
 
NumericVector query(2);
query << 0.5, 0.5;
double result = interp(query);  // Query at (0.5, 0.5)

Note

Interpolator caches numeric table data. For symbolic table values that must remain optimization variables, use the free interpn() overloads instead.