Modeling Excitable Cells with the EMI Equations: Spectral Analysis and Iterative Solution Strategy

In this work, we are interested in solving large linear systems stemming from the extra–membrane–intra model, which is employed for simulating excitable tissues at a cellular scale. After setting the related systems of partial differential equations equipped with proper boundary conditions, we provide its finite element discretization and focus on the resulting large linear systems. We first give a relatively complete spectral analysis using tools from the theory of Generalized Locally Toeplitz matrix sequences. The obtained spectral information is used for designing appropriate preconditioned Krylov solvers. Through numerical experiments, we show that the presented solution strategy is robust w.r.t. problem and discretization parameters, efficient and scalable.