High performance BLAS formulation of the multipole-to-local operator in the Fast Multipole Method - Inria - Institut national de recherche en sciences et technologies du numérique Accéder directement au contenu
Article Dans Une Revue Journal of Computational Physics Année : 2008

High performance BLAS formulation of the multipole-to-local operator in the Fast Multipole Method

Résumé

The multipole-to-local (M2L) operator is the most time-consuming part of the far field computation in the Fast Multipole Method for Laplace equation. Its natural expression, though commonly used, does not respect a sharp error bound: we here first prove the correctness of a second expression. We then propose a matrix formulation implemented with BLAS (Basic Linear Algebra Subprograms) routines in order to speed up its computation for these two expressions. We also introduce special data storages in memory to gain greater computational efficiency. This BLAS scheme is finally compared, for uniform distributions, to other M2L improvements such as block FFT, rotations and plane wave expansions. When considering runtime, extra memory storage, numerical stability and common precisions for Laplace equation, the BLAS version appears as the best one.
Fichier principal
Vignette du fichier
CoulaudFortinRoman05_V2.pdf (405.66 Ko) Télécharger le fichier
Origine : Fichiers produits par l'(les) auteur(s)
Loading...

Dates et versions

inria-00000957 , version 1 (19-12-2005)
inria-00000957 , version 2 (04-01-2007)

Identifiants

Citer

Olivier Coulaud, Pierre Fortin, Jean Roman. High performance BLAS formulation of the multipole-to-local operator in the Fast Multipole Method. Journal of Computational Physics, 2008, 227 (3), pp.1836-1862. ⟨10.1016/j.jcp.2007.09.027⟩. ⟨inria-00000957v2⟩
214 Consultations
256 Téléchargements

Altmetric

Partager

Gmail Facebook X LinkedIn More