20 #include "pism/earth/matlablike.hh"
23 double reqRelError,
void *fdata) {
25 double xmin[2] = {ax, ay};
26 double xmax[2] = {bx, by};
27 unsigned maxEval = 5000;
28 double result = 0.0, estimated_error = 0.0;
32 maxEval, 0.0, reqRelError, &result, &estimated_error);
int adapt_integrate(integrand f, void *fdata, unsigned dim, const double *xmin, const double *xmax, unsigned maxEval, double reqAbsError, double reqRelError, double *val, double *estimated_error)
double(* integrand)(unsigned ndim, const double *x, void *)
double dblquad_cubature(integrand f, double ax, double bx, double ay, double by, double reqRelError, void *fdata)