Fast parallel transient electromagnetic modeling using a uniform-in-time approximation to the exponential

Abstract

A new approach for the parallel forward modeling of transient electromagnetic (TEM) fields is presented. It is based on a family of uniform-in-time rational approximants to the matrix exponential that share a common denominator independent of the evaluation time points. The partial fraction decomposition of this family is exploited to devise a fast solver with high parallel efficiency. The number of shifted linear systems that need to be solved in parallel does not depend on the number of required time channels nor the spatial discretization. We also argue that similar parallel efficiency gains can be expected when solving the inverse TEM problem.