Structured inversion of the Bernstein-Vandermonde Matrix

by   Larry Allen, et al.

Bernstein polynomials, long a staple of approximation theory and computational geometry, have also increasingly become of interest in finite element methods. Many fundamental problems in interpolation and approximation give rise to interesting linear algebra questions. When attempting to find a polynomial approximation of boundary or initial data, one encounters the Bernstein-Vandermonde matrix, which is found to be highly ill-conditioned. Previously, we used the relationship between monomial Bezout matrices and the inverse of Hankel matrices to obtain a decomposition of the inverse of the Bernstein mass matrix in terms of Hankel, Toeplitz, and diagonal matrices. In this paper, we use properties of the Bernstein-Bezout matrix to factor the inverse of the Bernstein-Vandermonde matrix into a difference of products of Hankel, Toeplitz, and diagonal matrices. We also use a nonstandard matrix norm to study the conditioning of the Bernstein-Vandermonde matrix, showing that the conditioning in this case is better than in the standard 2-norm. Additionally, we use properties of multivariate Bernstein polynomials to derive a block LU decomposition of the Bernstein-Vandermonde matrix corresponding to equispaced nodes on the d-simplex.



There are no comments yet.


page 1

page 2

page 3

page 4


Structured inversion of the Bernstein mass matrix

Bernstein polynomials, long a staple of approximation theory and computa...

On the Closed Form Expression of Elementary Symmetric Polynomials and the Inverse of Vandermonde Matrix

Inverse Vandermonde matrix calculation is a long-standing problem to sol...

Graph-Theoretical Based Algorithms for Structural Optimization

Five new algorithms were proposed in order to optimize well conditioning...

Mixed precision recursive block diagonalization for bivariate functions of matrices

Various numerical linear algebra problems can be formulated as evaluatin...

On the minimal ranks of matrix pencils and the existence of a best approximate block-term tensor decomposition

Under the action of the general linear group, the ranks of matrices A an...

Easy and Efficient preconditioning of the Isogeometric Mass Matrix

This paper deals with the fast solution of linear systems associated wit...

Modifying AMG coarse spaces with weak approximation property to exhibit approximation in energy norm

Algebraic multigrid (AMG) coarse spaces are commonly constructed so that...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.