PISM, A Parallel Ice Sheet Model  stable v2.1-1-g6902d5502 committed by Ed Bueler on 2023-12-20 08:38:27 -0800
matlablike.cc
Go to the documentation of this file.
1 /* Copyright (C) 2004-2009, 2014, 2015, 2017, 2018, 2023 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 "pism/earth/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 
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)
Definition: cubature.c:701
double(* integrand)(unsigned ndim, const double *x, void *)
Definition: cubature.h:60
double dblquad_cubature(integrand f, double ax, double bx, double ay, double by, double reqRelError, void *fdata)
Definition: matlablike.cc:22