Disentangling complexity in Bayesian automatic adaptive quadrature

Jul 5, 2017, 9:00 AM
30m
LIT Conference Hall

LIT Conference Hall

Speaker

Prof. Gheorghe Adam (LIT-JINR and IFIN-HH, Bucharest, Romania)

Description

The present report discusses a Bayesian automatic adaptive quadrature (BAAQ) solution for numerical integration which is simultaneously robust, reliable, and efficient, hence able to yield guaranteed output in numerical experiments involving sudden unexpected modifications of the behavior of the integrand function. An essential ingredient of the solution is the multiscale approach [1]. Within it, for integration ranges of macroscopic length, which are of primary practical interest, an early decision path for the integrand profile (IP) scrutiny is defined which enables fast solution of four basic problems: (i) identification of simple integrals; (ii) check of the need to relax the user requested accuracy parameters; (iii) end of computation diagnostic for simple integrals; (iv) hints on manifestly ill-conditioning IP features. For integrals which are neither trivial, nor manifestly ill-conditioned, the Clenshaw-Curtis quadrature is activated within the approach discussed in [2]. This enables further identification of unresolved ill-conditioning features. We are thus left either with a hopefully well-conditioned integral, for which the standard automatic adaptive quadrature [3] is expected to yield reliable output, or with a manifestly ill-conditioned problem for which an improved version of the full BAAQ machinery [4] is activated. References [1] Gh. Adam, S. Adam, “Length Scales in Bayesian Automatic Adaptive Quadrature”, in EPJ Web of Conferences, vol. 108, 2016, 02002, 1-6; DOI: 10.1051/epjconf/201610802002. [2] S. Adam, Gh. Adam, “Summation Paths in Clenshaw-Curtis Quadrature”, in EPJ Web of Conferences, vol. 108, 2016, 02003, 1-6; DOI: 10.1051/epjconf/201610802003. [3] A.R. Krommer and C.W. Ueberhuber, Computational Integration SIAM, Philadelphia, 1998. [4] Gh. Adam and S. Adam, in Mathematical Modeling and Computational Science (MMCP2011), Gh. Adam, J. Buša, and M. Hnatič (Eds.), Springer, Heidelberg, LNCS 7125, (2012) pp.1–16.

Primary authors

Prof. Gheorghe Adam (LIT-JINR and IFIN-HH, Bucharest, Romania) Dr Sanda Adam (LIT-JINR and IFIN-HH, Bucharest, Romania)

Presentation materials