Stability of rootfinding for barycentric Lagrange interpolants (Q2454385): Difference between revisions

From MaRDI portal
Changed an Item
ReferenceBot (talk | contribs)
Changed an Item
(2 intermediate revisions by 2 users not shown)
Property / MaRDI profile type
 
Property / MaRDI profile type: MaRDI publication profile / rank
 
Normal rank
Property / full work available at URL
 
Property / full work available at URL: https://doi.org/10.1007/s11075-013-9770-3 / rank
 
Normal rank
Property / OpenAlex ID
 
Property / OpenAlex ID: W1969652671 / rank
 
Normal rank
Property / cites work
 
Property / cites work: LAPACK Users' Guide / rank
 
Normal rank
Property / cites work
 
Property / cites work: ON MATRICES DEPENDING ON PARAMETERS / rank
 
Normal rank
Property / cites work
 
Property / cites work: Rational functions for guaranteed and experimentally well-conditioned global interpolation / rank
 
Normal rank
Property / cites work
 
Property / cites work: Barycentric Lagrange Interpolation / rank
 
Normal rank
Property / cites work
 
Property / cites work: Q3447166 / rank
 
Normal rank
Property / cites work
 
Property / cites work: Roots of Polynomials Expressed in Terms of Orthogonal Polynomials / rank
 
Normal rank
Property / cites work
 
Property / cites work: Polynomial Roots from Companion Matrix Eigenvalues / rank
 
Normal rank
Property / cites work
 
Property / cites work: On eigenvectors and adjoints of modified matrices / rank
 
Normal rank
Property / cites work
 
Property / cites work: The Derivative of a Determinant / rank
 
Normal rank
Property / cites work
 
Property / cites work: THE COLLEAGUE MATRIX, A CHEBYSHEV ANALOGUE OF THE COMPANION MATRIX / rank
 
Normal rank
Property / cites work
 
Property / cites work: Q4320535 / rank
 
Normal rank
Property / cites work
 
Property / cites work: Backward Error of Polynomial Eigenproblems Solved by Linearization / rank
 
Normal rank
Property / cites work
 
Property / cites work: Fast Reduction of Generalized Companion Matrix Pairs for Barycentric Lagrange Interpolants / rank
 
Normal rank
Property / cites work
 
Property / cites work: Numerical stability of barycentric Hermite root-finding / rank
 
Normal rank
Property / cites work
 
Property / cites work: Balancing Regular Matrix Pencils / rank
 
Normal rank
Property / cites work
 
Property / cites work: Balancing a matrix for calculation of eigenvalues and eigenvectors / rank
 
Normal rank
Property / cites work
 
Property / cites work: Die Lage der Nullstellen eines Polynoms / rank
 
Normal rank
Property / cites work
 
Property / cites work: Pseudozeros of polynomials and pseudospectra of companion matrices / rank
 
Normal rank
Property / cites work
 
Property / cites work: Q4904857 / rank
 
Normal rank
Property / cites work
 
Property / cites work: Balancing the Generalized Eigenvalue Problem / rank
 
Normal rank
Property / cites work
 
Property / cites work: The Matrix Eigenvalue Problem / rank
 
Normal rank
Property / cites work
 
Property / cites work: The evaluation of the zeros of ill-conditioned polynomials. I, II / rank
 
Normal rank

Revision as of 14:41, 8 July 2024

scientific article
Language Label Description Also known as
English
Stability of rootfinding for barycentric Lagrange interpolants
scientific article

    Statements

    Stability of rootfinding for barycentric Lagrange interpolants (English)
    0 references
    0 references
    0 references
    13 June 2014
    0 references
    Given a univariate polynomial defined by \(n+1\) values \([f_0,\dots,f_n]\) at distinct nodes \([x_0,\dots,x_n]\), define the nodal polynomial \(\ell(z)\) by \(\ell(z) = \prod_{i=0}^n(z-x_i)\), and the barycentric weights \(w_j\) by \(w_j = \prod_{k \neq j}(x_j-x_k)^{-1}\), for \(0 \leq j \leq n\). The barycentric interpolation formula is then: \[ p(z) = \ell(z) \sum_{j=0}^{n} \frac{w_jf_j}{(z-x_j)}. \] The roots of polynomials expressed in this basis can be found via the eigenvalues of a generalized companion matrix pair \((\mathbf{A},\mathbf{B})\), and it is often advantageous to do so for reasons of numerical stability. While the stability of this method previously has been verified empirically, the present paper provides a detailed analysis of the stability of this approach. The authors show, in particular, that this method is normwise backward stable, and they propose a bound on the backward error which is ``typically only an order of magnitude larger than the actual error''. The analysis is rigorous as well as practical. For example, the computed eigenvalues of \((\mathbf{A},\mathbf{B})\) are assumed to be the exact eigenvalues of the matrix pair \((\mathbf{A+\Delta A}, \mathbf{B+\Delta B})\), where the perturbations satisfy \[ ||(\mathbf{\Delta A},\mathbf{\Delta B})||_F \leq \sigma(n) \varepsilon_M||(\mathbf{A}, \mathbf{B})||_F \] ``(as guaranteed by LAPACK...)'', where \(\varepsilon_M\) is machine precision, and \(||(\mathbf{A}, \mathbf{B})||_F\) is the Frobenius norm of a matrix pair \((\mathbf{A}, \mathbf{B})\). Further details as to how this is calculated can be found in the LAPACK Users' Guide. This leads them to an explicit bound on, essentially, the ratio of the backward error in the coefficients of the polynomial \(p(z)\) to the backward error of the generalized eighenvalue problem. The authors note that for their experiments, they cho0se \(\sigma(n) = \sqrt{n}\) because it seemed to be a ``good experimental fit, but [they] admit that [they] have no theoretical reason for choosing this.'' The authors give a full discussion of balancing and scaling of the companion matrix to bring it closer to a ``normal pair'', which allows them to produce roots with small backward error. They conclude with numerous experiments, comparing their bound to the test problems given in the paper by \textit{A. Edelman} and \textit{H. Murakami} [Math. Comput. 64, No. 210, 763--776 (1995; Zbl 0833.65041)], Chebyshev polynomials of the first and second kinds, the scaled Wilkinson polynomial and a Wilkinson filter example introduced on Page 178 of the article by \textit{J. H. Wilkinson} [Numer. Math. 1, 150--166, 167--180 (1959; Zbl 0202.43701)].
    0 references
    stability
    0 references
    backward error
    0 references
    Lagrange interpolation
    0 references
    barycentric formula
    0 references
    generalized companion matrices
    0 references
    polynomial roots
    0 references
    eigenvalue problem
    0 references
    numerical examples
    0 references

    Identifiers

    0 references
    0 references
    0 references
    0 references
    0 references
    0 references