PISM, A Parallel Ice Sheet Model  stable v2.1-1-g6902d5502 committed by Ed Bueler on 2023-12-20 08:38:27 -0800

◆ dblquad_cubature()

double dblquad_cubature ( integrand  f,
double  ax,
double  bx,
double  ay,
double  by,
double  reqRelError,
void *  fdata 
)

Functions like Matlab's dblquad() for 2D integrals. Wrapper for cubature by Steven Johnson.

The cubature code is used for double integration. It replaces iterated use of a single-variable adaptive rule as in Matlab's dblquad. cubature.{h|c} by Johnson is based on HIntLib by Rudolf Schuerer as well on GSL itself. cubature is GPLed.

Definition at line 22 of file matlablike.cc.

References adapt_integrate().

Referenced by pism::bed::LingleClarkSerial::compute_load_response_matrix().