Skip to Main content Skip to Navigation
New interface
Journal articles

ARKCoS: artifact-suppressed accelerated radial kernel convolution on the sphere

Abstract : We describe a hybrid Fourier/direct space convolution algorithm for compact radial (azimuthally symmetric) kernels on the sphere. For high resolution maps covering a large fraction of the sky, our implementation takes advantage of the inexpensive massive parallelism afforded by consumer graphics processing units (GPUs). Its applications include modeling of instrumental beam shapes in terms of compact kernels, computation of fine-scale wavelet transformations, and optimal filtering for the detection of point sources. Our algorithm works for any pixelization where pixels are grouped into isolatitude rings. Even for kernels that are not bandwidth-limited, ringing features are completely absent on an ECP grid. We demonstrate that they can be highly suppressed on the popular HEALPix pixelization, for which we develop a freely available implementation of the algorithm. As an example application, we show that running on a high-end consumer graphics card our method speeds up beam convolution for simulations of a characteristic Planck high frequency instrument channel by two orders of magnitude compared to the commonly used HEALPix implementation on one CPU core, while typically maintaining a fractional RMS accuracy of about 1 part in 105.
Complete list of metadata
Contributor : Hal Sorbonne Université Gestionnaire Connect in order to contact the contributor
Submitted on : Friday, April 22, 2022 - 11:47:22 AM
Last modification on : Friday, August 5, 2022 - 3:37:49 PM
Long-term archiving on: : Saturday, July 23, 2022 - 6:29:31 PM


Publisher files allowed on an open archive



F. Elsner, B. D. Wandelt. ARKCoS: artifact-suppressed accelerated radial kernel convolution on the sphere. Astronomy and Astrophysics - A&A, 2011, 532, ⟨10.1051/0004-6361/201116963⟩. ⟨insu-03645943⟩



Record views


Files downloads