Efficient sum-of-exponentials approximations for the heat kernel and their applications

Shidong Jiang, Leslie Greengard, Shaobo Wang

Research output: Contribution to journalArticlepeer-review

Abstract

In this paper, we show that efficient separated sum-of-exponentials approximations can be constructed for the heat kernel in any dimension. In one space dimension, the heat kernel admits an approximation involving a number of terms that is of the order (Formula Presented.) for any x∈ℝR and δ≤t≤T, where (Formula Presented.) is the desired precision. In all higher dimensions, the corresponding heat kernel admits an approximation involving only (Formula Presented.) terms for fixed accuracy (Formula Presented.). These approximations can be used to accelerate integral equation-based methods for boundary value problems governed by the heat equation in complex geometry. The resulting algorithms are nearly optimal. For NS points in the spatial discretization and NT time steps, the cost is (Formula Presented.) in terms of both memory and CPU time for fixed accuracy (Formula Presented.). The algorithms can be parallelized in a straightforward manner. Several numerical examples are presented to illustrate the accuracy and stability of these approximations.

Original languageEnglish (US)
Pages (from-to)529-551
Number of pages23
JournalAdvances in Computational Mathematics
Volume41
Issue number3
DOIs
StatePublished - Jun 22 2015

Keywords

  • Heat kernels
  • Heat potentials
  • Inverse laplace transform
  • Sum-of-exponentials approximation

ASJC Scopus subject areas

  • Computational Mathematics
  • Applied Mathematics

Fingerprint Dive into the research topics of 'Efficient sum-of-exponentials approximations for the heat kernel and their applications'. Together they form a unique fingerprint.

Cite this