Chebyshev multidomain pseudospectral method to solve cardiac wave equations with rotational anisotropy

The implementation of numerical methods to solve and study equations for cardiac wave propagation in realistic geometries is very costly, in terms of computational resources. The aim of this work is to show the improvement that can be obtained with Chebyshev polynomials-based methods over the classical finite difference schemes to obtain numerical solutions of cardiac models. To this end, we present a Chebyshev multidomain (CMD) Pseudospectral method to solve a simple two variable cardiac models on three-dimensional anisotropic media and we show the usefulness of the method over the traditional finite differences scheme widely used in the literature.