Automatic differentiation of quadrature

نویسندگان

  • Marina Menshikova
  • Shaun A. Forth
چکیده

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.

برای دانلود رایگان متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

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