On computing modified moments for half-range Hermite weights (Q6200847)
From MaRDI portal
scientific article; zbMATH DE number 7807009
Language | Label | Description | Also known as |
---|---|---|---|
English | On computing modified moments for half-range Hermite weights |
scientific article; zbMATH DE number 7807009 |
Statements
On computing modified moments for half-range Hermite weights (English)
0 references
20 February 2024
0 references
This interesting paper studies the computation of modified moments for the half-range Hermite weight. Integrals defined with the Hermite weight on the real semiaxis appear in many applications for example lattice Bolzman models, harmonically oscillation rarefied gas flows and many others. These integrals are typically computed using modified moments via product rules. It is known, however that there is a loss of accuracy in computing such modified moments, having only three correct decimal digits for the size of the problem equal to 6, and Gautschi has described a way to overcome this problem by computing the modified moments with extended precision arithmetic. In the paper under review, the authors study the problem of computing the modified moments for the system of orthonormal Laguerre polynomials on the real semiaxis with the Hermite weight, also called half-range Hermite weights. The orthonomal Laguerre polynomials are orthonormal polynomials on the interval \([0,\infty)\) satisfying the three term recurrence relation: \[ \left\{ \begin{array}{ll} p_{0}=1, \\ \beta_{1}p_1(x)=(x-\alpha_{1})p_{0}(x) \\ \beta_{l}p_{l}(x)=(x-\alpha_{l})p_{l-1}(x)-\beta_{l-1} p_{l-2}(x),\, l\geq 2 \end{array} \right. \] with \(\beta_l=-l\) and \(\alpha_{l}=2l-1, l=0.1,...\) They are orthonormal with respect to the weight function \(\exp(-x)\), i.e. \[ \int_{0}^{\infty}p_{i}(x)p_{j}(x)\exp(-x)dx=\delta_{i,j}. \] The authors propose a new method to compute these modified moments, based on the construction of the null-space of a rectangular matrix derived from the three-term recurrence relation of the system of orthonormal Laguerre polynomials. They show that that their algorithm computes the modified moments with high relative accuracy in floating point arithmetic, thereby avoiding the use of extended precision arithmetic. Numerical examples show the effectiveness of the proposed approach. The paper is well written with a good set of references.
0 references
null-space
0 references
Gaussian quadrature rule
0 references
modified moments
0 references
product rule
0 references
0 references