An order optimal solver for the discretized bidomain equations

The electrical activity in the heart is governed by the bidomain equations. In this paper, we analyse an order optimal method for the algebraic equations arising from the discretization of this model. Our scheme is defined in terms of block Jacobi or block symmetric Gauss–Seidel preconditioners. Furthermore, each block in these methods is based on standard preconditioners for scalar elliptic or parabolic partial differential equations (PDEs). Such preconditioners can be realized in terms of multigrid or domain decomposition schemes, and are thus readily available by applying ‘off-the-shelves’ software. Finally, our theoretical findings are illuminated by a series of numerical experiments. Copyright © 2006 John Wiley & Sons, Ltd.

[1]  Kent-André Mardal,et al.  Uniform preconditioners for the time dependent Stokes problem , 2004, Numerische Mathematik.

[2]  T. Chan,et al.  Domain decomposition algorithms , 1994, Acta Numerica.

[3]  A. Tveito,et al.  Multigrid Block Preconditioning for a Coupled System of Partial Differential Equations Modeling the Electrical Activity in the Heart , 2002, Computer methods in biomechanics and biomedical engineering.

[4]  O. Widlund.  Domain Decomposition Algorithms , 1993 .

[5]  Randolph E. Bank,et al.  An optimal order process for solving finite element equations , 1981 .

[6]  A. Tveito,et al.  Electrical Activity in the Human Heart , 2003 .

[7]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[8]  W. Hackbusch.  Iterative Solution of Large Sparse Systems of Equations , 1993 .

[9]  D. Rose,et al.  Linear algebraic transformations of the bidomain equations: implications for numerical methods. , 1994, Mathematical biosciences.

[10]  James P. Keener,et al.  Mathematical physiology , 1998 .

[11]  Ragnar Winther,et al.  A Preconditioned Iterative Method for Saddlepoint Problems , 1992, SIAM J. Matrix Anal. Appl..

[12]  Denis Noble,et al.  Simulating cardiac sinus and atrial network dynamics on the Connection Machine , 1993 .

[13]  Jinchao Xu,et al.  Iterative Methods by Space Decomposition and Subspace Correction , 1992, SIAM Rev..

[14]  A. Bruaset.  A survey of preconditioned iterative methods , 1995 .

[15]  D. Arnold,et al.  Preconditioning discrete approximations of the Reissner-Mindlin plate model , 1997 .

[16]  Maxim A. Olshanskii,et al.  On the Convergence of a Multigrid Method for Linear Reaction-Diffusion Problems , 2000, Computing.

[17]  Zhong-Zhi Bai.  A class of modified block SSOR preconditioners for symmetric positive definite systems of linear equations , 1999, Adv. Comput. Math..

[18]  H. Elman,et al.  Block preconditioners for the discrete incompressible Navier–Stokes equations , 2002 .

[19]  A. Garfinkel,et al.  An advanced algorithm for solving partial differential equation in cardiac conduction , 1999, IEEE Transactions on Biomedical Engineering.

[20]  G. Meurant.  Computer Solution of Large Linear Systems , 1999 .

[21]  Xiao-Chuan Cai,et al.  Additive Schwarz algorithms for parabolic convection-diffusion equations , 1991 .