A sparse Laplacian in tensor product wavelet coordinates

作者：Dijkema Tammo Jan; Stevenson Rob^*

来源：Numerische Mathematik, 2010, 115(3): 433-449.

DOI：10.1007/s00211-010-0288-5

• Sparse representations
• Tensor product approximation
• Adaptive wavelet scheme
• Riesz bases
• Cubic Hermite splines

We construct a wavelet basis on the unit interval with respect to which both the (infinite) mass and stiffness matrix corresponding to the one-dimensional Laplacian are (truly) sparse and boundedly invertible. As a consequence, the (infinite) stiffness matrix corresponding to the Laplacian on the n-dimensional unit box with respect to the n-fold tensor product wavelet basis is also sparse and boundedly invertible. This greatly simplifies the implementation and improves the quantitative properties of an adaptive wavelet scheme to solve the multi-dimensional Poisson equation. The results extend to any second order partial differential operator with constant coefficients that defines a boundedly invertible operator.

• 出版日期2010-5