test_nint
test_nint,
a C++ code which
defines a set of test problems
for the approximate computation of integrals over multidimensional
regions.
Routines are available to evaluate the integrand, return the exact
value of the integral, report the name of the problem,
report the integration limits, get, set or modify a base point.
The integrands is assigned an index. The integrands can be invoked
by index. Most integrands may be defined for any value of the
spatial dimension, which we denote here by m. Most integrands
are defined on the unit mdimensional hypercube. Some
integrands include one or more parameters. These generally have
default values, which can be altered by the user.
For each problem, a set of routines are available with a standard
interface, for manipulating and evaluating the problem. For a
problem with index "87", for instance, we might have the following
set of routines. The most important is P87_F which evaluates
the integrand. We probably also need P87_LIM to determine
the limits of integration, and P87_EXACT to get the exact
value of the integral (if known). A number of routines are available
to set, get, or randomize parameters associated with the problem.

P87_DEFAULT sets default values for problem 87.

P87_EXACT returns the exact integral for problem 87.

P87_F evaluates the integrand for problem 87.

P87_I4 sets or gets integer scalar parameters for problem 87.

P87_I4VEC sets or gets integer vector parameters for problem 87.

P87_LIM returns the integration limits for problem 87.

P87_NAME returns the name of problem 87.

P87_R8 sets or gets real scalar parameters for problem 87.

P87_R8VEC sets or gets real vector parameters for problem 87.

P87_REGION returns the name of the integration region
for problem 87.

P87_TITLE prints a title for problem 87.
The list of integrand functions includes:

f(x) = ( sum ( x(1:m) ) )^2;

f(x) = ( sum ( 2 * x(1:m)  1 ) )^4;

f(x) = ( sum ( x(1:m) ) )^5;

f(x) = ( sum ( 2 * x(1:m)  1 ) )^6;

f(x) = 1 / ( 1 + sum ( 2 * x(1:m) ) );

f(x) = product ( 2 * abs ( 2 * x(1:m)  1 ) );

f(x) = product ( pi / 2 ) * sin ( pi * x(1:m) );

f(x) = ( sin ( (pi/4) * sum ( x(1:m) ) ) )**2;

f(x) = exp ( sum ( c(1:m) * x(1:m) ) );

f(x) = sum ( abs ( x(1:m)  0.5 ) );

f(x) = exp ( sum ( abs ( 2 * x(1:m)  1 ) ) );

f(x) = product ( 1 <= i <= m ) ( i * cos ( i * x(i) ) );

f(x) = product ( 1 <= i <= m ) t(n(i))(x(i)), t(n(i))
is a Chebyshev polynomial;

f(x) = sum ( 1 <= i <= m ) (1)^i * product ( 1 <= j <= i ) x(j);

f(x) = product ( 1 <= i <= order ) x(mod(i1,m)+1);

f(x) = sum ( abs ( x(1:m)  x0(1:m) ) );

f(x) = sum ( ( x(1:m)  x0(1:m) )^2 );

f(x) = 1 inside an mdimensional sphere around x0(1:m), 0 outside;

f(x) = product ( sqrt ( abs ( x(1:m)  x0(1:m) ) ) );

f(x) = ( sum ( x(1:m) ) )^power;

f(x) = c * product ( x(1:m)^e(1:m) ) on the surface of
an mdimensional unit sphere;

f(x) = c * product ( x(1:m)^e(1:m) ) in an mdimensional ball;

f(x) = c * product ( x(1:m)^e(1:m) ) in the unit mdimensional simplex;

f(x) = product ( abs ( 4 * x(1:m)  2 ) + c(1:m) )
/ ( 1 + c(1:m) ) );

f(x) = exp ( c * product ( x(1:m) ) );

f(x) = product ( c(1:m) * exp (  c(1:m) * x(1:m) ) );

f(x) = cos ( 2 * pi * r + sum ( c(1:m) * x(1:m) ) ),
Genz "Oscillatory";

f(x) = 1 / product ( c(1:m)**2 + (x(1:m)  x0(1:m))^2),
Genz "Product Peak";

f(x) = 1 / ( 1 + sum ( c(1:m) * x(1:m) ) )^(m+r),
Genz "Corner Peak";

f(x) = exp(sum(c(1:m)^2 * ( x(1:m)  x0(1:m))^2 ) ),
Genz "Gaussian";

f(x) = exp (  sum ( c(1:m) * abs ( x(1:m)  x0(1:m) ) ) ),
Genz "Continuous";

f(x) = exp(sum(c(1:m)*x(1:m)) for x(1:m) <= x0(1:m), 0 otherwise,
Genz "Discontinuous";
An Important Quote:
"When good results are obtained in integrating a highdimensional
function, we should conclude first of all that an especially tractable
integrand was tried and not that a generally successful method has
been found. A secondary conclusion is that we might have made a
very good choice in selecting an integration method to exploit whatever
features of f made it tractable."
Art Owen,
Latin Supercube Sampling for Very High Dimensional Simulation,
ACM Transactions on Modeling and Computer Simulations,
Volume 8, Number 1, January 1998, pages 71102.
Licensing:
The computer code and data files described and made available on this web page
are distributed under
the GNU LGPL license.
Languages:
test_nint is available in
a C++ version and
a FORTRAN90 version and
a MATLAB version.
Related Data and Programs:
CLENSHAW_CURTIS_RULE,
a C++ code which
sets a Clenshaw Curtis quadrature grid in multiple dimensions.
NINTLIB,
a C++ code which
numerically estimates integrals in multiple dimensions.
PRODUCT_RULE,
a C++ code which
creates a multidimensional quadrature rule as a product of
one dimensional rules.
STROUD,
a C++ code which
defines quadrature rules for a variety of unusual areas, surfaces and volumes in 2D,
3D and Ndimensions.
test_nint_test
TESTPACK,
a C++ code which
defines a set of integrands used to test multidimensional quadrature.
Reference:

JD Beasley, SG Springer,
Algorithm AS 111:
The Percentage Points of the Normal Distribution,
Applied Statistics,
Volume 26, 1977, pages 118121.

Paul Bratley, Bennett Fox, Harald Niederreiter,
Implementation and Tests of LowDiscrepancy Sequences,
ACM Transactions on Modeling and Computer Simulation,
Volume 2, Number 3, July 1992, pages 195213.

Roger Broucke,
Algorithm 446:
Ten Subroutines for the Manipulation of Chebyshev Series,
Communications of the ACM,
Volume 16, 1973, pages 254256.

William Cody, Kenneth Hillstrom,
Chebyshev Approximations for the Natural Logarithm of the
Gamma Function,
Mathematics of Computation,
Volume 21, Number 98, April 1967, pages 198203.

Richard Crandall,
Projects in Scientific Computing,
Springer, 2005,
ISBN: 0387950095,
LC: Q183.9.C733.

Philip Davis, Philip Rabinowitz,
Methods of Numerical Integration,
Second Edition,
Dover, 2007,
ISBN: 0486453391,
LC: QA299.3.D28.

Gerald Folland,
How to Integrate a Polynomial Over a Sphere,
American Mathematical Monthly,
Volume 108, Number 5, May 2001, pages 446448.

Leslie Fox, Ian Parker,
Chebyshev Polynomials in Numerical Analysis,
Oxford Press, 1968,
LC: QA297.F65.

Alan Genz,
Testing Multidimensional Integration Routines,
in Tools, Methods, and Languages for Scientific and
Engineering Computation,
edited by B Ford, JC Rault, F Thomasset,
NorthHolland, 1984, pages 8194,
ISBN: 0444875700,
LC: Q183.9.I53.

Alan Genz,
A Package for Testing Multiple Integration Subroutines,
in Numerical Integration:
Recent Developments, Software and Applications,
edited by Patrick Keast, Graeme Fairweather,
Reidel, 1987, pages 337340,
ISBN: 9027725144,
LC: QA299.3.N38.

Kenneth Hanson,
QuasiMonte Carlo: halftoning in high dimensions?
in Computatinal Imaging,
Edited by CA Bouman, RL Stevenson,
Proceedings SPIE,
Volume 5016, 2003, pages 161172.

John Hart, Ward Cheney, Charles Lawson, Hans Maehly,
Charles Mesztenyi, John Rice, Henry Thatcher,
Christoph Witzgall,
Computer Approximations,
Wiley, 1968,
LC: QA297.C64.

Claude Itzykson, JeanMichel Drouffe,
Statistical Field Theory,
Volume 1,
Cambridge, 1991,
ISBN: 0521408067,
LC: QC174.8.I89.

Stephen Joe, Frances Kuo
Remark on Algorithm 659:
Implementing Sobol's Quasirandom Sequence Generator,
ACM Transactions on Mathematical Software,
Volume 29, Number 1, March 2003, pages 4957.

David Kahaner, Cleve Moler, Steven Nash,
Numerical Methods and Software,
Prentice Hall, 1989,
ISBN: 0136272584,
LC: TA345.K34.

Bradley Keister,
Multidimensional Quadrature Algorithms,
Computers in Physics,
Volume 10, Number 2, March/April, 1996, pages 119122.

Arnold Krommer, Christoph Ueberhuber,
Numerical Integration on Advanced Compuer Systems,
Springer, 1994,
ISBN: 3540584102,
LC: QA299.3.K76.

Anargyros Papageorgiou, Joseph Traub,
Faster Evaluation of Multidimensional Integrals,
Computers in Physics,
Volume 11, Number 6, November/December 1997, pages 574578.

Thomas Patterson,
On the Construction of a Practical ErmakovZolotukhin
Multiple Integrator,
in Numerical Integration:
Recent Developments, Software and Applications,
edited by Patrick Keast and Graeme Fairweather,
D. Reidel, 1987, pages 269290.

Arthur Stroud,
Approximate Calculation of Multiple Integrals,
Prentice Hall, 1971,
ISBN: 0130438936,
LC: QA311.S85.

Arthur Stroud, Don Secrest,
Gaussian Quadrature Formulas,
Prentice Hall, 1966,
LC: QA299.4G3S7.

Xiaoqun Wang, KaiTai Fang,
The Effective Dimension and quasiMonte Carlo Integration,
Journal of Complexity,
Volume 19, pages 101124, 2003.
Source Code:
Last revised on 23 April 2020.