Optimized M2L Kernels for the Chebyshev Interpolation based Fast Multipole Method

Abstract : A fast multipole method (FMM) for asymptotically smooth kernel functions (1/r, 1/r^4, Gauss and Stokes kernels, radial basis functions, etc.) based on a Chebyshev interpolation scheme has been introduced in [Fong et al., 2009]. The method has been extended to oscillatory kernels (e.g., Helmholtz kernel) in [Messner et al., 2012]. Beside its generality this FMM turns out to be favorable due to its easy implementation and its high performance based on intensive use of highly optimized BLAS libraries. However, one of its bottlenecks is the precomputation of the multiple-to-local (M2L) operator, and its higher number of floating point operations (flops) compared to other FMM formulations. Here, we present several optimizations for that operator, which is known to be the costliest FMM operator. The most efficient ones do not only reduce the precomputation time by a factor up to 340 but they also speed up the matrix-vector product. We conclude with comparisons and numerical validations of all presented optimizations.


https://hal.inria.fr/hal-00746089
Contributeur : Matthias Messner <>
Soumis le : mardi 20 novembre 2012 - 01:44:16
Dernière modification le : samedi 17 septembre 2016 - 01:36:54
Document(s) archivé(s) le : jeudi 21 février 2013 - 11:50:42

Fichiers

main.pdf
Fichiers produits par l'(les) auteur(s)

Identifiants

  • HAL Id : hal-00746089, version 2
  • ARXIV : 1210.7292

Collections

Citation

Matthias Messner, Bérenger Bramas, Olivier Coulaud, Eric Darve. Optimized M2L Kernels for the Chebyshev Interpolation based Fast Multipole Method. [Research Report] 2012, pp.22. <hal-00746089v2>

Exporter

Partager

Métriques

Consultations de
la notice

227

Téléchargements du document

224