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

02/12/2018
by   Fredrik Johansson, et al.
0

We describe a strategy for rigorous arbitrary-precision evaluation of Legendre polynomials on the unit interval and its application in the generation of Gauss-Legendre quadrature rules. Our focus is on making the evaluation practical for a wide range of realistic parameters, corresponding to the requirements of numerical integration to an accuracy of about 100 to 100 000 bits. Our algorithm combines the summation by rectangular splitting of several types of expansions in terms of hypergeometric series with a fixed-point implementation of Bonnet's three-term recurrence relation. We then compute rigorous enclosures of the Gauss-Legendre nodes and weights using the interval Newton method. We provide rigorous error bounds for all steps of the algorithm. The approach is validated by an implementation in the Arb library, which achieves order-of-magnitude speedups over previous code for computing Gauss-Legendre rules with simultaneous high degree and precision.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
02/22/2018

Numerical integration in arbitrary-precision ball arithmetic

We present an implementation of arbitrary-precision numerical integratio...
research
05/09/2017

Computing the Lambert W function in arbitrary-precision complex interval arithmetic

We describe an algorithm to evaluate all the complex branches of the Lam...
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
06/12/2019

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

Methods for the computation of classical Gaussian quadrature rules are d...
research
04/16/2021

Recent advances in rigorous computation of Poincaré maps

In this article we present recent advances on interval methods for rigor...
research
07/07/2016

Rigorous Multiple-Precision Evaluation of D-Finite Functions in SageMath

We present a new open source implementation in the SageMath computer alg...
research
02/12/2023

Numerical methods and arbitrary-precision computation of the Lerch transcendent

We examine the use of the Euler-Maclaurin formula and new derived unifor...

Please sign up or login with your details

Forgot password? Click here to reset