Abstract
The numerical solution of the bidomain equations is a fundamental tool in the field of computational cardiac electrophysiology. The multi-scale nature of the processes involved and the increasing complexity of the tissue and ionic models of interest make the use of High Performance Computing techniques mandatory. The solution of the linear system arising from the Finite Element discretisation of the bidomain equations is the main bottleneck in terms of performance and parallel scalability. In this work we present and evaluate parallel implementations of two problemspecific block preconditioners, namely BD and LDU, in the framework of the open source cardiac electrophysiology simulation package Chaste. Both rely upon providing a reliable and practical method for calculating the action of the inverse of certain matrix subblocks. Our results show that: a) total execution time can be reduced by using different approximation methods for different matrix subblocks, b) there exists a good correlation between the coefficient Δt=h2 and the cases where BD and LDU outperform generic Algebraic Multigrid (AMG) preconditioning, and c) when applied to simulations with current state-of-the-art cardiac geometries, BD and LDU are over 3 times faster than AMG preconditioning and show similar scaling.
Published Version (
Free)
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a callDisclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.
