Accurate recovery of recursion coefficients from Gaussian quadrature formulas
We present an algorithm for the accurate recovery of recursion coefficients from quadrature formulas with positive weights, based on the differential form of the quotient-difference algorithm. The process is slightly faster than the Gragg–Harrod algorithm and is forward stable in the sense of compon...
Saved in:
Published in | Journal of computational and applied mathematics Vol. 112; no. 1; pp. 165 - 180 |
---|---|
Main Author | |
Format | Journal Article |
Language | English |
Published |
Elsevier B.V
30.11.1999
|
Online Access | Get full text |
Cover
Loading…
Summary: | We present an algorithm for the accurate recovery of recursion coefficients from quadrature formulas with positive weights, based on the differential form of the quotient-difference algorithm. The process is slightly faster than the Gragg–Harrod algorithm and is forward stable in the sense of componentwise relative error when the gaps between the nodes are available as data. This result shows constructively that the converse problem is well posed when the data are required to be accurate floating point numbers. This is not a contradiction of the examples given by Gragg and Harrod which are ill posed in a vector norm, because in that case very small numbers are not required to have any significant digits. |
---|---|
ISSN: | 0377-0427 1879-1778 |
DOI: | 10.1016/S0377-0427(99)00228-9 |