Automatic differentiation of quadrature
نویسندگان
چکیده
We analyse the application of automatic differentiation (AD) to the quadrature (numerical integration) of a function integrand to determine the sensitivities of the integral to variation in the limits of integration. We derive an expression for the truncation errors of such ADderived sensitivities and relate them to the truncation error of the original, and a closely related, function quadrature. Our results hold provided the integrand is one degree higher continuously differentiable than that sufficient for convergence of its quadrature. Numerical results validate our analysis. However, utilising algebraic expressions for such sensitivities, instead of directly applying AD, results in an approach that proved more efficient for the tetrachoric correlation estimation example we considered using our Matlab AD framework.
منابع مشابه
Automatic Differentiation based QMOM for Population Balance Equations
The quadrature method of moments (QMOM) has emerged as a promising tool for the solution of population balance equations in the past few years. The QMOM requires solving differential algebraic equations (DAEs) consisting of ordinary differential equations related to the evolution of moments as well as nonlinear algebraic equations resulting from the quadrature approximation of moments. In this ...
متن کاملHandling accuracy in Bayesian automatic adaptive quadrature
New Bayesian inferences which significantly extend the coverage of the Bayesian automatic adaptive quadrature (BAAQ) approach to the solution of Riemann integrals are reported. The scrutiny of the possible floating-point machine number approximations of abscissa values inside an integration domain unveiled the occurrence of five classes of integration domains entering the quadrature problems: z...
متن کاملDedicated to Professor Oliviu Gherman’s 80th Anniversary GLOBAL CONTROL WITHIN THE BAYESIAN AUTOMATIC ADAPTIVE QUADRATURE
The Bayesian automatic adaptive quadrature solves numerically the Riemann integrals via the dynamic construction and handling of subrange trees the terminal nodes of which are subranges eventually characterized by integrand profiles for which the reliability of the outputs of the local quadrature rules is explicitly checked and enforced. The present paper discusses the characteristic variables ...
متن کاملQuadrature With Respect to Binomial Measures
This work is devoted to the study of integration with respect to binomial measures. We develop interpolation quadrature rules and study their properties. Applying a local error estimate based on null rules, we test two automatic integrators with local quadrature rules that generalize the five points Newton Cotes formula.
متن کاملAn Automatic Integration of Infinite Range Integrals Involving Bessel Functions
An efficient automatic quadrature procedure is developed for numerically computing the integrals 0 , where the function is smooth and nonoscillatory at infinity and is the Bessel functions of order ν =1,0 and 1/4. The procedure involves the use of an automatic integration scheme of modified FFT used for evaluating Fourier integrals and product type integration, and the modified W-transformation...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
عنوان ژورنال:
- Optimization Methods and Software
دوره 27 شماره
صفحات -
تاریخ انتشار 2012