Stimulus Protocol Determines the Most Computationally Efficient Preconditioner for the Bidomain Equations

被引：8

作者：

Bernabeu, Miguel O. ^{[1
]}

Pathmanathan, Pras ^{[1
]}

Pitt-Francis, Joe ^{[1
]}

Kay, David ^{[1
]}

机构：

[1] Univ Oxford, Comp Lab, Oxford OX1 3QD, England

来源：

IEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING | 2010年 / 57卷 / 12期

基金：

英国工程与自然科学研究理事会;

关键词：

Bidomain equations; computational efficiency; preconditioning; whole heart geometries; REACTION-DIFFUSION SYSTEMS; MULTIGRID PRECONDITIONER; CARDIAC VULNERABILITY; ELECTRICAL BEHAVIOR; MONODOMAIN; TISSUE;

D O I：

10.1109/TBME.2010.2078817

中图分类号：

R318 [生物医学工程];

学科分类号：

0831 ;

摘要：

The efficient solution of the bidomain equations is a fundamental tool in the field of cardiac electrophysiology. When choosing a finite element discretization of the coupled system, one has to deal with the solution of a large, highly sparse system of linear equations. The conjugate gradient algorithm, along with suitable preconditioning, is the natural choice in this scenario. In this study, we identify the optimal preconditioners with respect to both stimulus protocol and mesh geometry. The results are supported by a comprehensive study of the mesh-dependence properties of several preconditioning techniques found in the literature. Our results show that when only intracellular stimulus is considered, incomplete LU factorization remains a valid choice for current cardiac geometries. However, when extracellular shocks are delivered to tissue, preconditioners that take into account the structure of the system minimize execution time and ensure mesh-independent convergence.

引用

页码：2806 / 2815

页数：10

共 28 条

[1]

[Anonymous], 2003, ITERATIVE METHODS SP, DOI DOI 10.1137/1.9780898718003

[2] Development of an anatomically detailed MRI-derived rabbit ventricular model and assessment of its impact on simulations of electrophysiological function [J].