Abstract
We describe a fast and robust method for solving the large sparse linear systems that arise upon the discretization of elliptic partial differential equations such as Laplace’s equation and the Helmholtz equation at low frequencies. While most existing fast schemes for this task rely on so called “iterative” solvers, the method described here solves the linear system directly (to within an arbitrary predefined accuracy). The method is described for the particular case of an operator defined on a square uniform grid, but can be generalized other geometries. For a grid containing N points, a single solve requires O(Nlog 2 N) arithmetic operations and \(O(\sqrt{N}\log N)\) storage. Storing the information required to perform additional solves rapidly requires O(Nlog N) storage. The scheme is particularly efficient in situations involving domains that are loaded on the boundary only and where the solution is sought only on the boundary. In this environment, subsequent solves (after the first) can be performed in \(O(\sqrt{N}\log N)\) operations. The efficiency of the scheme is illustrated with numerical examples. For instance, a system of size 106×106 is directly solved to seven digits accuracy in four minutes on a 2.8 GHz P4 desktop PC.
Similar content being viewed by others
References
Börm, S.: H 2-matrix arithmetics in linear complexity. Computing 77(1), 1–28 (2006)
Börm, S.: Approximation of solution operators of elliptic partial differential equations by H and H 2-matrices. Tech. Report 85/2007, Max Planck Institute for Mathematics in the Sciences (2007)
Chandrasekaran, S., Gu, M., Li, X.S., Xia, J.: Some fast algorithms for hierarchically semiseparable matrices. Private Communication (2007)
Chandrasekaran, S., Gu, M., Li, X.S., Xia, J.: Superfast multifrontal method for structured linear systems of equations. Private Communication (2007)
Chandrasekaran, S., Gu, M., Lyons, W.: A fast adaptive solver for hierarchically semiseparable representations. Calcolo 42(3–4), 171–185 (2005)
Chandrasekaran, S., Gu, M., Pals, T.: A fast ULV decomposition solver for hierarchically semiseparable representations. SIAM J. Matrix Anal. Appl. 28(3), 603–622 (2006). (Electronic)
George, A.: Nested dissection of a regular finite element mesh. SIAM J. Numer. Anal. 10, 345–363 (1973)
Grasedyck, L., Kriemann, R., Le Borne, S.: Domain-decomposition based H-matrix preconditioners. In: Proceedings of DD16. LNSCE, vol. 55, pp. 661–668. Springer, Berlin (2006)
Hackbusch, W.: A sparse matrix arithmetic based on H-matrices. I. Introduction to H-matrices. Computing 62(2), 89–108 (1999)
Hackbusch, W., Khoromskij, B., Sauter, S.A.: On H 2-matrices. In: Lectures on Applied Mathematics (Munich, 1999), pp. 9–29. Springer, Berlin (2000)
Hoffman, A.J., Martin, M.S., Rose, D.J.: Complexity bounds for regular finite difference and finite element grids. SIAM J. Numer. Anal. 10, 364–369 (1973)
Martinsson, P.G., Rokhlin, V.: A fast direct solver for boundary integral equations in two dimensions. J. Comput. Phys. 205(1), 1–23 (2005)
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Martinsson, PG. A Fast Direct Solver for a Class of Elliptic Partial Differential Equations. J Sci Comput 38, 316–330 (2009). https://doi.org/10.1007/s10915-008-9240-6
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10915-008-9240-6