Fast, reliable and unrestricted iterative computation of Gauss--Hermite and Gauss--Laguerre quadratures

06/12/2019
by   A. Gil, et al.
0

Methods for the computation of classical Gaussian quadrature rules are described which are effective both for small and large degree. These methods are reliable because the iterative computation of the nodes has guaranteed convergence, and they are fast due to their fourth-order convergence and its asymptotic exactness for an appropriate selection of the variables. For Gauss--Hermite and Gauss--Laguerre quadratures, local Taylor series can be used for computing efficiently the orthogonal polynomials involved, with exact initial values for the Hermite case and first values computed with a continued fraction for the Laguerre case. The resulting algorithms have almost unrestricted validity with respect to the parameters. Full relative precision is reached for the Hermite nodes, without any accuracy loss and for any degree, and a mild accuracy loss occurs for the Hermite and Laguerre weights as well as for the Laguerre nodes. These fast methods are exclusively based on convergent processes, which, together with the high order of convergence of the underlying iterative method, makes them particularly useful for high accuracy computations. We show examples of very high accuracy computations (of up to 1000 digits of accuracy).

READ FULL TEXT
research
08/19/2020

Fast and reliable high accuracy computation of Gauss–Jacobi quadrature

Iterative methods with certified convergence for the computation of Gaus...
research
02/12/2018

Fast and rigorous arbitrary-precision computation of Gauss-Legendre quadrature nodes and weights

We describe a strategy for rigorous arbitrary-precision evaluation of Le...
research
12/28/2021

Gaussian quadrature rules for composite highly oscillatory integrals

Highly oscillatory integrals of composite type arise in electronic engin...
research
11/10/2021

Fast Computation of Hahn Polynomials for High Order Moments

Discrete Hahn polynomials (DHPs) and their moments are considered to be ...
research
10/08/2019

Efficient implementation of adaptive order reconstructions

Including polynomials with small degree and stencil when designing very ...
research
08/11/2018

Fast RodFIter for Attitude Reconstruction from Inertial Measurements

Attitude computation is of vital importance for a variety of application...
research
02/11/2019

Computing Extremely Accurate Quantiles Using t-Digests

We present on-line algorithms for computing approximations of rank-based...

Please sign up or login with your details

Forgot password? Click here to reset