SIAM Journal on Scientific and Statistical Computing
A fast method for solving the heat equation by layer potentials
Journal of Computational Physics
Nyström discretization of parabolic boundary integral equations
Applied Numerical Mathematics
Multidimensional Fast Gauss Transforms by Chebyshev Expansions
SIAM Journal on Scientific Computing
High Order Accurate Methods for the Evaluation of Layer Heat Potentials
SIAM Journal on Scientific Computing
Hi-index | 31.45 |
An efficient scheme for solving boundary integral equations of the heat equation based on the Galerkin method is introduced. The parabolic fast multipole method (pFMM) is applied to accelerate the evaluation of the thermal layer potentials. In order to remain attractive for a wide range of applications, a key issue is to ensure efficiency for a big variety of temporal to spatial mesh ratios. Within the parabolic Galerkin FMM (pGFMM) it turns out that the temporal nearfield can become very costly. To that end, a modified fast Gauss transform (FGT) is developed. The complexity and convergence behavior of the method are analyzed and numerically investigated on a range of model problems. The results demonstrate that the complexity is nearly optimal in the number of discretization parameters while the convergence rate of the Galerkin method is preserved.