Published February 15, 2015
| Submitted
Journal Article
Open
Short note on the mass matrix for Gauss–Lobatto grid points
- Creators
- Teukolsky, Saul A.
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(N^2).
Additional Information
© 2014 Elsevier Inc. Received 16 September 2014, Revised 20 November 2014, Accepted 6 December 2014, Available online 12 December 2014. I thank Jan Hesthaven for several helpful comments, including pointing out that interpolation from a coarse grid to a finer one introduces no aliasing. This work was supported in part by NSF Grants PHY-1306125 and AST-1333129 at Cornell University, and by a grant from the Sherman Fairchild Foundation.Attached Files
Submitted - 1412.2276.pdf
Files
1412.2276.pdf
Files
(74.8 kB)
Name | Size | Download all |
---|---|---|
md5:71f2952abafebe0021b1f0b1cc19923b
|
74.8 kB | Preview Download |
Additional details
- Eprint ID
- 86801
- DOI
- 10.1016/j.jcp.2014.12.012
- Resolver ID
- CaltechAUTHORS:20180605-142355801
- PHY-1306125
- NSF
- AST-1333129
- NSF
- Sherman Fairchild Foundation
- Created
-
2018-06-05Created from EPrint's datestamp field
- Updated
-
2021-11-15Created from EPrint's last_modified field