tmatlablike.cc - pism - [fork] customized build of PISM, the parallel ice sheet model (tillflux branch)
HTML git clone git://src.adamsgaard.dk/pism
DIR Log
DIR Files
DIR Refs
DIR LICENSE
---
tmatlablike.cc (1250B)
---
1 /* Copyright (C) 2004-2009, 2014, 2015, 2017, 2018 Ed Bueler and Constantine Khroulev
2
3 This file is part of PISM.
4
5 PISM is free software; you can redistribute it and/or modify it under the
6 terms of the GNU General Public License as published by the Free Software
7 Foundation; either version 3 of the License, or (at your option) any later
8 version.
9
10 PISM is distributed in the hope that it will be useful, but WITHOUT ANY
11 WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
12 FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
13 details.
14
15 You should have received a copy of the GNU General Public License
16 along with PISM; if not, write to the Free Software
17 Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
18 */
19
20 #include "matlablike.hh"
21
22 double dblquad_cubature(integrand f, double ax, double bx, double ay, double by,
23 double reqRelError, void *fdata) {
24
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;
29
30 /* see cubature.h: */
31 adapt_integrate(f, fdata, 2, xmin, xmax,
32 maxEval, 0.0, reqRelError, &result, &estimated_error);
33 return result;
34 }
35