Skip to Main content Skip to Navigation
Journal articles

Efficient Spherical Harmonic Transforms aimed at pseudo-spectral numerical simulations

Nathanaël Schaeffer 1
1 Géodynamo
ISTerre - Institut des Sciences de la Terre
Abstract : In this paper, we report on very efficient algorithms for the spherical harmonic transform (SHT). Explicitly vectorized variations of the algorithm based on the Gauss-Legendre quadrature are discussed and implemented in the SHTns library which includes scalar and vector transforms. The main breakthrough is to achieve very efficient on-the-fly computations of the Legendre associated functions, even for very high resolutions, by taking advantage of the specific properties of the SHT and the advanced capabilities of current and future computers. This allows us to simultaneously and significantly reduce memory usage and computation time of the SHT. We measure the performance and accuracy of our algorithms. Even though the complexity of the algorithms implemented in SHTns are in $O(N^3)$ (where N is the maximum harmonic degree of the transform), they perform much better than any third party implementation, including lower complexity algorithms, even for truncations as high as N=1023. SHTns is available at as open source software.
Complete list of metadata

Cited literature [19 references]  Display  Hide  Download
Contributor : Nathanaël Schaeffer <>
Submitted on : Tuesday, November 4, 2014 - 10:47:35 PM
Last modification on : Monday, July 19, 2021 - 11:22:02 AM
Long-term archiving on: : Thursday, February 5, 2015 - 11:36:06 AM


Files produced by the author(s)


Distributed under a Creative Commons Attribution - NonCommercial - ShareAlike 4.0 International License




Nathanaël Schaeffer. Efficient Spherical Harmonic Transforms aimed at pseudo-spectral numerical simulations. Geochemistry, Geophysics, Geosystems, AGU and the Geochemical Society, 2013, 14 (3), pp.751-758. ⟨10.1002/ggge.20071⟩. ⟨insu-00675145v4⟩



Record views


Files downloads