Efficient Spherical Harmonic Transforms aimed at pseudo-spectral numerical simulations - INSU - Institut national des sciences de l'Univers Accéder directement au contenu
Pré-Publication, Document De Travail Année : 2012

Efficient Spherical Harmonic Transforms aimed at pseudo-spectral numerical simulations

Nathanaël Schaeffer

Résumé

In this paper, we report on very efficient algorithms for the spherical harmonic transform (SHT) that can be used in numerical simulations of partial differential equations. Explicitly vectorized variations of the Gauss-Legendre algorithm are discussed and implemented in the open-source library SHTns which includes scalar and vector transforms. This library is especially suitable for direct numerical simulations of non-linear partial differential equations in spherical geometry, like the Navier-Stokes equation. The performance of our algorithms is compared to third party SHT implementations, including fast algorithms. Even though the complexity of the algorithms implemented in SHTns are of order O(N^3) (where N is the maximum harmonic degree of the transform), they perform much better than the available implementations of asymptotically fast algorithms, even for a truncation as high as N=1023. In our performance tests, the best performance for SHT on the x86 platform is delivered by SHTns, which is available at https://bitbucket.org/nschaeff/shtns as open source software.
Fichier principal
Vignette du fichier
sht.pdf (185.73 Ko) Télécharger le fichier
Origine : Fichiers produits par l'(les) auteur(s)

Dates et versions

insu-00675145 , version 1 (29-02-2012)
insu-00675145 , version 2 (08-12-2012)
insu-00675145 , version 3 (29-01-2013)
insu-00675145 , version 4 (04-11-2014)

Identifiants

Citer

Nathanaël Schaeffer. Efficient Spherical Harmonic Transforms aimed at pseudo-spectral numerical simulations. 2012. ⟨insu-00675145v2⟩
606 Consultations
1186 Téléchargements

Altmetric

Partager

Gmail Facebook X LinkedIn More