Abstract:Strong gravitational lensing is a powerful probe of cosmology and the dark matter distribution. Efficient lensing software is already a necessity to fully use its potential and the performance demands will only increase with the upcoming generation of telescopes. In this paper, we present a proof-of-concept study on the impact of High Performance Computing techniques on a performance-critical part of the widely used lens modeling software LENSTOOL. We implement the algorithm once as a highly optimized CPU version and once with graphics card acceleration for a simple parametric lens model. In addition, we study the impact of finite machine precision on the lensing algorithm. While double precision is the default choice for scientific applications, we find that single precision can be sufficiently accurate for our purposes and lead to a big speedup. Therefore we develop and present a mixed precision algorithm which only uses double precision when necessary. We measure the performance of the different implementations and find that the use of High Performance Computing Techniques dramatically improves the code performance both on CPUs and GPUs. Compared to the current LENSTOOL implementation on 12 CPU cores, we obtain speedup factors of up to 170. We achieve this optimal performance by using our mixed precision algorithm on a high-end GPU which is common in modern supercomputers. We also show that these techniques reduce the energy consumption by up to 98%. Furthermore, we demonstrate that a highly competitive speedup can be reached with consumer GPUs. While they are an order of magnitude cheaper than the high-end graphics cards, they are rarely used for scientific computations due to their low double precision performance. However, our mixed precision algorithm unlocks their full potential. Consequently, the consumer GPU delivers a speedup which is only a factor of four lower than the best speedup achieved by a high-end GPU.

Computational advances in gravitational microlensing: a comparison of CPU, GPU, and parallel, large data codes

An Improved GPU-based Ray-shooting Code for Gravitational Microlensing

Teraflop per second gravitational lensing ray-shooting using graphics processing units

High Performance Computing for gravitational lens modeling: Single vs double precision on GPUs and CPUs

HyGrid: A CPU-GPU Hybrid Convolution-Based Gridding Algorithm in Radio Astronomy.

Accelerating Gravitational Microlensing Simulations Using the Xeon Phi Coprocessor

Parallelized Inference for Gravitational-Wave Astronomy

Simulation Of The Microlensing Effect Near The Critical Curve Of The Galaxy Cluster

Accelerating Astrophysical Particle Simulations with Programmable Hardware (FPGA and GPU)

Wave effect of gravitational waves intersected with a microlens field II: an adaptive hierarchical tree algorithm and population study

Performance Analysis of Parallel Gravitationaln-Body Codes on Large GPU Clusters

A GPU-Enabled, High-Resolution Cosmological Microlensing Parameter Survey

Study of A Gpu-Based Parallel Computing Method for the Monte Carlo Program

High Performance Direct Gravitational N-body Simulations on Graphics Processing Unit I: An implementation in Cg

High Performance Direct Gravitational N-body Simulations on Graphics Processing Units

High Performance Direct Gravitational N-body Simulations on Graphics Processing Units -- II: An implementation in CUDA

Astrophysical applications of gravitational microlensing

Scalable matched-filtering pipeline for gravitational-wave searches of compact binary mergers

Accelerating the 3D Elastic Wave Forward Modeling on GPU and MIC

VBMicroLensing: three algorithms for multiple lensing with contour integration

Supermassive Black Hole Binaries in High Performance Massively Parallel Direct N-body Simulations on Large GPU Clusters