Prev

Next

Index-> contents reference index search external

Up-> CppAD library RombergMul

CppAD-> Install Introduction AD ADFun multi_thread library cppad_ipopt_nlp Example preprocessor Appendix

library-> ErrorHandler NearEqual speed_test SpeedTest time_test NumericType CheckNumericType SimpleVector CheckSimpleVector nan pow_int Poly LuDetAndSolve RombergOne RombergMul Runge45 Rosen34 OdeErrControl OdeGear OdeGearControl BenderQuad opt_val_hes LuRatio CppAD_vector thread_alloc memory_leak

RombergMul-> RombergMul.cpp

Headings-> Syntax Description Include m r F a b n p e Float FloatVector Example Source Code

Multi-dimensional Romberg Integration
Syntax
# include <cppad/romberg_mul.hpp>
RombergMul<Fun, SizeVector, FloatVector, m> R
r = R(F, a, b, n, p, e)

Description
Returns the Romberg integration estimate $\mathit{r}$ for the multi-dimensional integral $$\mathit{r}={\int }_{\mathit{a}[0]}^{\mathit{b}[0]}\cdots {\int }_{\mathit{a}[\mathit{m}-1]}^{\mathit{b}[\mathit{m}-1]}\mathit{F}(\mathit{x})\mathbf{d}{\mathit{x}}_0\cdots \mathbf{d}{\mathit{x}}_{\mathit{m}-1}+\sum _{\mathit{i}=0}^{\mathit{m}-1}\mathit{O}{\left[(\mathit{b}[\mathit{i}]-\mathit{a}[\mathit{i}])/2^{\mathit{n}[\mathit{i}]-1}\right]}^{2(\mathit{p}[\mathit{i}]+1)}$$
Include
The file cppad/romberg_mul.hpp is included by cppad/cppad.hpp but it can also be included separately with out the rest of the CppAD routines.

m
The template parameter m must be convertible to a size_t object with a value that can be determined at compile time; for example 2. It determines the dimension of the domain space for the integration.

r
The return value r has prototype
     Float r
It is the estimate computed by RombergMul for the integral above (see description of Float below).

F
The object F has the prototype
     Fun &F
It must support the operation
     F(x)
The argument x to F has prototype
     const Float &x
The return value of F is a Float object

a
The argument a has prototype
     const FloatVector &a
It specifies the lower limit for the integration (see description of FloatVector below).

b
The argument b has prototype
     const FloatVector &b
It specifies the upper limit for the integration.

n
The argument n has prototype
     const SizeVector &n
A total number of $2^{\mathit{n}[\mathit{i}]-1}+1$ evaluations of F(x) are used to estimate the integral with respect to $\mathbf{d}{\mathit{x}}_{\mathit{i}}$ .

p
The argument p has prototype
     const SizeVector &p
For $\mathit{i}=0,\dots ,\mathit{m}-1$ , $\mathit{n}[\mathit{i}]$ determines the accuracy order in the approximation for the integral that is returned by RombergMul. The values in p must be less than or equal n; i.e., p[i] <= n[i].

e
The argument e has prototype
     Float &e
The input value of e does not matter and its output value is an approximation for the absolute error in the integral estimate.

Float
The type Float is defined as the type of the elements of FloatVector . The type Float must satisfy the conditions for a NumericType type. The routine CheckNumericType will generate an error message if this is not the case. In addition, if x and y are Float objects,
     x < y
returns the bool value true if x is less than y and false otherwise.

FloatVector
The type FloatVector must be a SimpleVector class. The routine CheckSimpleVector will generate an error message if this is not the case.

Example
The file RombergMul.cpp contains an example and test a test of using this routine. It returns true if it succeeds and false otherwise.

Source Code
The source code for this routine is in the file cppad/romberg_mul.hpp.

---

Input File: cppad/romberg_mul.hpp