Richardson method and totally nonnegative linear systems (Q603132)

From MaRDI portal
scientific article
Language Label Description Also known as
English
Richardson method and totally nonnegative linear systems
scientific article

    Statements

    Richardson method and totally nonnegative linear systems (English)
    0 references
    0 references
    0 references
    0 references
    5 November 2010
    0 references
    Totally nonnegative stochastic matrices arise as collocation matrices in computer-aided design interpolating curves in a shape preserving manner [see e. g. \textit{J. M. Peña}. Shape preserving representations in computer aided geometric design. Huntington, NY: Nova Science. (1999; Zbl 1005.68163)] by polynomial and spline functions. The corresponding linear systems \(A x = b\) for the computation of the control points can be solved by the Richardson iterative method which converges if and only if the spectral radius of \(\rho(I-A)\), \(I\) identity matrix, is less than \(1\). In order to accelerate the convergence the authors use a relaxation of the method. They prove the convergence of this modified Richardson method for nonsingular totally nonnegative stochastic matrices and derive the parameter for the optimal convergence speed. Additionally, a modified Richardson method is presented, that converges for any nonsingular totally nonnegative matrix, including a procedure for the estimation of the optimal convergence parameter. The results are verified by numerical examples arising in the interpolation by spline and polynomial functions. Discussed are Vandermonde matrices, Bernstein-Vandermonde matrices, and quadratic and cubic B-spline collocation matrices. The choice of the optimal parameter reduces the number of iterations to a factor greater than a half. The convergence is low in the case of bad conditioning using higher degrees in the polynomial interpolation but, the convergence is ensured for any dimension.
    0 references
    Richardson method
    0 references
    nonsingular totally nonnegative matrices
    0 references
    nonsingular totally nonnegative stochastic matrices
    0 references
    shape preserving modeling of curves
    0 references
    convergence acceleration
    0 references
    collocation matrices
    0 references
    computer-aided design
    0 references
    numerical examples
    0 references
    interpolation
    0 references
    spline
    0 references
    Vandermonde matrices
    0 references
    Bernstein-Vandermonde matrices
    0 references
    cubic B-spline
    0 references

    Identifiers

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