Proposed numerical methods address the cardiac electrophysiology model, which involves the electrical activity of the heart described by a nonlinear reaction-diffusion PDE coupled with ODEs for electrochemical reactions in cells.
The methods integrate an operator splitting technique for the PDE with primal hybrid methods for spatial discretization, allowing for flexible approximations of the Lagrange multiplier.
Results from convergence studies demonstrate optimal rates of convergence, and numerical experiments indicate that these new methods are often more efficient than traditional numerical techniques when using preconditioned iterative approaches for solving linear systems.