Skip to main content Accessibility help

A Fast Direct Solver for a Class of 3-D Elliptic Partial Differential Equation with Variable Coefficient

  • Beibei Huang (a1) (a2), Bin Tu (a1) and Benzhuo Lu (a1)


We propose a direct solver for the three-dimensional Poisson equation with a variable coefficient, and an algorithm to directly solve the associated sparse linear systems that exploits the sparsity pattern of the coefficient matrix. Introducing some appropriate finite difference operators, we derive a second-order scheme for the solver, and then two suitable high-order compact schemes are also discussed. For a cube containing N nodes, the solver requires arithmetic operations and memory to store the necessary information. Its efficiency is illustrated with examples, and the numerical results are analysed.


Corresponding author



Hide All
[1]Martinsson, P.‐G., A fast direct solver for a class of elliptic partial differential equations, J. Sci. Comput. 38 (2009), pp. 316–330.
[2]Engquist, B. and Ying, L., Sweeping preconditioner for the Helmholtz equation: Moving perfectly matched layers, Multiscale Model. Sim. 9(2) (2011), 686–710.
[3]Chandrasekaran, S., Gu, M., and Lyons, W., A fast adaptive solver for hierarchically semisep-arable representations, Calcolo 42(3-4) (2005), pp. 171–185.
[4]Chandrasekaran, S., Gu, M., and Lyons, W., Superfast multifrontal method for structured linear systems of equations, Private Communication (2007).
[5]Grasedyck, L., Kriemann, R., and Le Borne, S., Domain-decomposition based H-matrix pre-conditioners, Proceedings of DD16. LNSCE, Vol. 55, Springer, Berlin (2006), pp. 661–668, pp. 471–476.
[6]Le Tallec, P., Domain decomposition methods in computational mechanical, Computational Mechanics Advances (1986), pp. 121–220.
[7]Pack, G.R., Garrett, G.A., Wong, L., and Lamm, G., The Effect of a variable dielectric coefficient and finite Ion size on Poisson-Boltzmann calculations of DNA-electrolyte systems, Biophysical J. 65 (1993), pp. 1363–1370.
[8]Warwicker, J. and Watson, H.C., Calculation of the electric potential in the active site cleft due to α-helical dipoles, J. Mol. Biol. 157 (1982), pp. 671–679.
[9]Lu, B.Z. and McCammon, J.A., Molecular surface-free continuum model for electrodiffusion processes, Chem. Phys. Lett. 451(4-6) 2008, pp. 282–286.
[10]Warming, R.F. and Hyett, B.J., The modified equation approach to the stability and accuracy analysis of finite-difference methods, J. Comput. Phys. 14(2) (1974), pp. 159–179.
[11]Hackbusch, W., Elliptic Differential Equations Theory and Numerical Treatment, Springer-Verlag Berlin Heidelberg, Berlin, 1992, pp. 40–72. (Electronic).
[12]Iserles, A., A First Course in the Numerical Analysis of Differential Equations, Tinghua University Press, 2005, pp. 101–127. (Chinese version).
[13]Yu, D.H. and Tang, H.Z., Numerical Solution of Differential Equations, Science Press, 2007, pp. 300–307. (Chinese version).
[14]Liberty, E., Woolfe, F., Martinsson, P., Rokhlin, V., and Tygert, M., Randomized algorithms for the low-rank approximation of matrices, Proc. Natl. Acad. Sci. USA 104 (2007), pp. 20167– 0172.
[15]Woolfe, F., Liberty, E., Rokhlin, V., and Tygert, M., A fast randomized algorithm for the approximation of matrices, Appl. Comput. Harmon. Anal. 25 (2008), pp. 335–366.
[16]Lin, L., Lu, J., and Ying, L., Fast construction of hierarchical matrix representation from matrix-vector multiplication, J. Comput. Phys. 230(10) (2010), pp. 4071–4087.
[17]Spotz, W.F. and Carey, G.F., A high-order compact formulation for the 3-D Poisson equation, Numer. Methods Part. Diff. Eqs. 12 (1996), pp. 235–243.
[18]Spotz, W. F. and Carey, G. F., High-order compact scheme for the steady stream-function vorticity equations, Int. J. Numer. Methods Eng. 38 (1995), pp. 3497–3512.
[19]Gupta, M.M. and Kouatchou, J., Symbolic derivation of finite difference approximations for three dimensional Poisson equation, Int. J. Numer. Methods Part. Diff. Eqs. 14(5) (1998), pp. 593–606.
[20]Ananthakrishnaiah, U., Manohar, R., and Stephenson, J. W., Fourth-order finite difference for three-dimensional general elliptic problems with variable coefficients, Numer. Methods Part. Diff. Eqs. 3 (1987), pp. 229–240.
[21]Zhang, J., An explicit fourth-order compact finite difference scheme for three dimensional convection-diffusion equation, Commun. Numer. Methods Eng. 14 (1997), pp. 263–280.



Full text views

Total number of HTML views: 0
Total number of PDF views: 0 *
Loading metrics...

Abstract views

Total abstract views: 0 *
Loading metrics...

* Views captured on Cambridge Core between <date>. This data will be updated every 24 hours.

Usage data cannot currently be displayed