Document Type
Article
Publication Date
2014
Publication Title
SIAM Journal on Scientific Computing
Department
Department of Mathematics
Abstract
A numerical method for solving elliptic PDEs with variable coefficients on two-dimensional domains is presented. The method is based on high-order composite spectral approximations and is designed for problems with smooth solutions. The resulting system of linear equations is solved using a direct (as opposed to iterative) solver that has optimal $O(N)$ complexity for all stages of the computation when applied to problems with nonoscillatory solutions such as the Laplace and the Stokes equations. Numerical examples demonstrate that the scheme is capable of computing solutions with a relative accuracy of $10^{-10}$ or better for challenging problems such as highly oscillatory Helmholtz problems and convection-dominated convection-diffusion equations. In terms of speed, it is demonstrated that a problem with a nonoscillatory solution that was discretized using $10^{8}$ nodes can be solved in 115 minutes on a personal workstation with two quad-core 3.3 GHz CPUs. Since the solver is direct, and the “solution operator” fits in RAM, any solves beyond the first are very fast. In the example with $10^{8}$ unknowns, solves require only 30 seconds.
DOI
10.1137/130918988
Dartmouth Digital Commons Citation
Gillman, A. and Martinsson, P. G., "A Direct Solver with O(N) Complexity for Variable Coefficient Elliptic PDEs Discretized via a High-Order Composite Spectral Collocation Method" (2014). Dartmouth Scholarship. 2049.
https://digitalcommons.dartmouth.edu/facoa/2049
