Cosmological Calculations on the GPU

Abstract

Cosmological measurements require the calculation of nontrivial quantities over large datasets. The next generation of survey telescopes (such as DES, PanSTARRS, and LSST) will yield measurements of billions of galaxies. The scale of these datasets, and the nature of the calculations involved, make cosmological calculations ideal models for implementation on graphics processing units (GPUs). We consider two cosmological calculations, the two-point angular correlation function and the aperture mass statistic, and aim to improve the calculation time by constructing code for calculating them on the GPU. Using CUDA, we implement the two algorithms on the GPU and compare the calculation speeds to comparable code run on the CPU. We obtain a code speed-up of between 10 - 180x faster, compared to performing the same calculation on the CPU. The code has been made publicly available.

10 pages, 1 figure; accepted by Astronomy and Computing; code available at https://github.com/djbard/ccogs

References in corpus (14)

• Weak Gravitational Lensing
• The 2dF Galaxy Redshift Survey: Power-spectrum analysis of the final dataset and cosmological implications
• LSST Science Book, Version 2.0
• Detection of (dark) matter concentrations via weak gravitational lensing
• Accelerating NBODY6 with Graphics Processing Units
• A sparse octree gravitational N-body code that runs entirely on the GPU processor
• Cosmology with the shear-peak statistics
• Probing Cosmology with Weak Lensing Peak Counts
• Phantom-GRAPE: numerical software library to accelerate collisionless $N$-body simulation with SIMD instruction set on x86 architecture
• Baryon Acoustic Oscillations
• Accelerating Radio Astronomy Cross-Correlation with Graphics Processing Units
• Application of GPUs for the Calculation of Two Point Correlation Functions in Cosmology
• Challenges for LSST scale data sets
• Massive Datasets in Astronomy