Wavemoth -- Fast spherical harmonic transforms by butterfly matrix compression

We present Wavemoth, an experimental open source code for computing scalar spherical harmonic transforms (SHTs). Such transforms are ubiquitous in astronomical data analysis. Our code performs substantially better than existing publicly available codes due to improvements on two fronts. First, the computational core is made more efficient by using small amounts of precomputed data, as well as paying attention to CPU instruction pipelining and cache usage. Second, Wavemoth makes use of a fast and numerically stable algorithm based on compressing a set of linear operators in a precomputation step. The resulting SHT scales as O(L^2 (log L)^2) for the resolution range of practical interest, where L denotes the spherical harmonic truncation degree. For low and medium-range resolutions, Wavemoth tends to be twice as fast as libpsht, which is the current state of the art implementation for the HEALPix grid. At the resolution of the Planck experiment, L ~ 4000, Wavemoth is between three and six times faster than libpsht, depending on the computer architecture and the required precision. Due to the experimental nature of the project, only spherical harmonic synthesis is currently supported, although adding support or spherical harmonic analysis should be trivial.Comment: 13 pages, 6 figures, accepted by ApJ

Scale-discretised ridgelet transform on the sphere

We revisit the spherical Radon transform, also called the Funk-Radon transform, viewing it as an axisymmetric convolution on the sphere. Viewing the spherical Radon transform in this manner leads to a straightforward derivation of its spherical harmonic representation, from which we show the spherical Radon transform can be inverted exactly for signals exhibiting antipodal symmetry. We then construct a spherical ridgelet transform by composing the spherical Radon and scale-discretised wavelet transforms on the sphere. The resulting spherical ridgelet transform also admits exact inversion for antipodal signals. The restriction to antipodal signals is expected since the spherical Radon and ridgelet transforms themselves result in signals that exhibit antipodal symmetry. Our ridgelet transform is defined natively on the sphere, probes signal content globally along great circles, does not exhibit blocking artefacts, supports spin signals and exhibits an exact and explicit inverse transform. No alternative ridgelet construction on the sphere satisfies all of these properties. Our implementation of the spherical Radon and ridgelet transforms is made publicly available. Finally, we illustrate the effectiveness of spherical ridgelets for diffusion magnetic resonance imaging of white matter fibers in the brain.Comment: 5 pages, 4 figures, matches version accepted by EUSIPCO, code available at http://www.s2let.or

A Fast and Accurate Algorithm for Spherical Harmonic Analysis on HEALPix Grids with Applications to the Cosmic Microwave Background Radiation

The Hierarchical Equal Area isoLatitude Pixelation (HEALPix) scheme is used extensively in astrophysics for data collection and analysis on the sphere. The scheme was originally designed for studying the Cosmic Microwave Background (CMB) radiation, which represents the first light to travel during the early stages of the universe's development and gives the strongest evidence for the Big Bang theory to date. Refined analysis of the CMB angular power spectrum can lead to revolutionary developments in understanding the nature of dark matter and dark energy. In this paper, we present a new method for performing spherical harmonic analysis for HEALPix data, which is a central component to computing and analyzing the angular power spectrum of the massive CMB data sets. The method uses a novel combination of a non-uniform fast Fourier transform, the double Fourier sphere method, and Slevinsky's fast spherical harmonic transform (Slevinsky, 2019). For a HEALPix grid with $N$ pixels (points), the computational complexity of the method is $\mathcal{O}(N\log^2 N)$, with an initial set-up cost of $\mathcal{O}(N^{3/2}\log N)$. This compares favorably with $\mathcal{O}(N^{3/2})$ runtime complexity of the current methods available in the HEALPix software when multiple maps need to be analyzed at the same time. Using numerical experiments, we demonstrate that the new method also appears to provide better accuracy over the entire angular power spectrum of synthetic data when compared to the current methods, with a convergence rate at least two times higher

Libpsht - algorithms for efficient spherical harmonic transforms

Libpsht (or "library for Performant Spherical Harmonic Transforms") is a collection of algorithms for efficient conversion between spatial-domain and spectral-domain representations of data defined on the sphere. The package supports transforms of scalars as well as spin-1 and spin-2 quantities, and can be used for a wide range of pixelisations (including HEALPix, GLESP and ECP). It will take advantage of hardware features like multiple processor cores and floating-point vector operations, if available. Even without this additional acceleration, the employed algorithms are among the most efficient (in terms of CPU time as well as memory consumption) currently being used in the astronomical community. The library is written in strictly standard-conforming C90, ensuring portability to many different hard- and software platforms, and allowing straightforward integration with codes written in various programming languages like C, C++, Fortran, Python etc. Libpsht is distributed under the terms of the GNU General Public License (GPL) version 2 and can be downloaded from http://sourceforge.net/projects/libpsht.Comment: 9 pages, 8 figures, accepted by A&

On the computation of directional scale-discretized wavelet transforms on the sphere

We review scale-discretized wavelets on the sphere, which are directional and allow one to probe oriented structure in data defined on the sphere. Furthermore, scale-discretized wavelets allow in practice the exact synthesis of a signal from its wavelet coefficients. We present exact and efficient algorithms to compute the scale-discretized wavelet transform of band-limited signals on the sphere. These algorithms are implemented in the publicly available S2DW code. We release a new version of S2DW that is parallelized and contains additional code optimizations. Note that scale-discretized wavelets can be viewed as a directional generalization of needlets. Finally, we outline future improvements to the algorithms presented, which can be achieved by exploiting a new sampling theorem on the sphere developed recently by some of the authors.Comment: 13 pages, 3 figures, Proceedings of Wavelets and Sparsity XV, SPIE Optics and Photonics 2013, Code is publicly available at http://www.s2dw.org

Magnetic field topology of the unique chemically peculiar star CU Virginis

The late-B magnetic chemically peculiar star CU Vir is one of the fastest rotators among the intermediate-mass stars with strong fossil magnetic fields. It shows a prominent rotational modulation of the spectral energy distribution and absorption line profiles due to chemical spots and exhibits a unique strongly beamed variable radio emission. Little is known about the magnetic field topology of CU Vir. In this study we aim to derive, for the first time, detailed maps of the magnetic field distribution over the surface of this star. We use high-resolution spectropolarimetric observations covering the entire rotational period. These data are interpreted using a multi-line technique of least-squares deconvolution (LSD) and a new Zeeman Doppler imaging code based on detailed polarised radiative transfer modelling of the Stokes I and V LSD profiles. This new magnetic inversion approach relies on the spectrum synthesis calculations over the full wavelength range covered by observations and does not assume that the LSD profiles behave as a single spectral line with mean parameters. We present magnetic and chemical abundance maps derived from the Si and Fe lines. Mean polarisation profiles of both elements reveal a significant departure of the magnetic field topology of CU Vir from the commonly assumed axisymmetric dipolar configuration. The field of CU Vir is dipolar-like, but clearly non-axisymmetric, showing a large difference of the field strength between the regions of opposite polarity. The main relative abundance depletion features in both Si and Fe maps coincide with the weak-field region in the magnetic map. Detailed information on the distorted dipolar magnetic field topology of CU Vir provided by our study is essential for understanding chemical spot formation, radio emission, and rotational period variation of this star.Comment: 14 pages, 14 figures; accepted for publication in A&

A Large Sky Simulation of the Gravitational Lensing of the Cosmic Microwave Background

Large scale structure deflects cosmic microwave background (CMB) photons. Since large angular scales in the large scale structure contribute significantly to the gravitational lensing effect, a realistic simulation of CMB lensing requires a sufficiently large sky area. We describe simulations that include these effects, and present both effective and multiple plane ray-tracing versions of the algorithm, which employs spherical harmonic space and does not use the flat sky approximation. We simulate lensed CMB maps with an angular resolution of ~0.9 arcmin. The angular power spectrum of the simulated sky agrees well with analytical predictions. Maps generated in this manner are a useful tool for the analysis and interpretation of upcoming CMB experiments such as PLANCK and ACT.Comment: 14 pages, 12 figures, replaced with version accepted for publication by the AP

A simple model potential for hollow nanospheres

A new model potential is introduced to describe the hollow nanospheres such as fullerene and molecular structures and to obtain their electronic properties. A closed analytical solution of the corresponding treatment is given within the framework of supersymmetric perturbation theory.Comment: 7 pages, 3 figure