Modeling excitable cells with the EMI equations: spectral analysis and iterative solution strategy

by   Pietro Benedusi, et al.

In this work, we are interested in solving large linear systems stemming from the Extra-Membrane-Intra (EMI) model, which is employed for simulating excitable tissues at a cellular scale. After setting the related systems of partial differential equations (PDEs) equipped with proper boundary conditions, we provide numerical approximation schemes for the EMI PDEs 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. We show, through numerical experiments, that the presented solution strategy is robust w.r.t. problem and discretization parameters, efficient and scalable.


page 1

page 2

page 3

page 4


A systematic approach to reduced GLT

This paper concerns the spectral analysis of matrix-sequences that are g...

Spectral Analysis and Preconditioned Iterative Solvers for Large Structured Linear Systems

In this thesis, the numerical solution of three different classes of pro...

Learning to correct spectral methods for simulating turbulent flows

Despite their ubiquity throughout science and engineering, only a handfu...

Recursive blocked algorithms for linear systems with Kronecker product structure

Recursive blocked algorithms have proven to be highly efficient at the n...

Evaluating Accuracy and Efficiency of HPC Solvers for Sparse Linear Systems with Applications to PDEs

Partial Differential Equations (PDEs) describe several problems relevant...

Please sign up or login with your details

Forgot password? Click here to reset