Structured inversion of the Bernstein mass matrix

by   Larray 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. Previously, we gave block-structured algorithms for inverting the Bernstein mass matrix on simplicial cells, but did not study fast alorithms for the univariate case. Here, we give several approaches to inverting the univariate mass matrix based on exact formulae for the inverse; decompositions of the inverse in terms of Hankel, Toeplitz, and diagonal matrices; and a spectral decomposition. In particular, the eigendecomposition can be explicitly constructed in O(n^2) operations, while its accuracy for solving linear systems is comparable to that of the Cholesky decomposition. Moreover, we study conditioning and accuracy of these methods from the standpoint of the effect of roundoff error in the L^2 norm on polynomials, showing that the conditioning in this case is far less extreme than in the standard 2-norm.



There are no comments yet.


page 1

page 2

page 3

page 4


Structured inversion of the Bernstein-Vandermonde Matrix

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

A High-Order Lower-Triangular Pseudo-Mass Matrix for Explicit Time Advancement of hp Triangular Finite Element Methods

Explicit time advancement for continuous finite elements requires the in...

Easy and Efficient preconditioning of the Isogeometric Mass Matrix

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

Sublinear-Time Quadratic Minimization via Spectral Decomposition of Matrices

We design a sublinear-time approximation algorithm for quadratic functio...

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...

Spectral analysis of continuous FEM for hyperbolic PDEs: influence of approximation, stabilization, and time-stepping

We study continuous finite element dicretizations for one dimensional hy...

Mixed precision recursive block diagonalization for bivariate functions of matrices

Various numerical linear algebra problems can be formulated as evaluatin...
This week in AI

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