Article ID | Journal | Published Year | Pages | File Type |
---|---|---|---|---|
6931833 | Journal of Computational Physics | 2015 | 6 Pages |
Abstract
The mass matrix for Gauss-Lobatto grid points is usually approximated by Gauss-Lobatto quadrature because this leads to a diagonal matrix that is easy to invert. The exact mass matrix and its inverse are full. We show that the exact mass matrix and its inverse differ from the approximate diagonal ones by a simple rank-1 update (outer product). They can thus be applied to an arbitrary vector in O(N) operations instead of O(N2).
Related Topics
Physical Sciences and Engineering
Computer Science
Computer Science Applications
Authors
Saul A. Teukolsky,