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.

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

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

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

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

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

Accelerating Gravitational Microlensing Simulations Using the Xeon Phi Coprocessor

GPU Accelerated Anomaly Detection of Large Scale Light Curves

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

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

GPU Parallel Acceleration of Beam Tracing for Propagation Prediction in Urban Environments

View Dependent Layer Sampling: an Approach to Hardware Implementation of Volume Ray Casting

VBMicroLensing: three algorithms for multiple lensing with contour integration

Efficient Modeling and Real-Time Rendering of Massive Urban Pipelines Based on GPU

An Improved Back-Projection Algorithm for GNSS-R BSAR Imaging Based on CPU and GPU Platform

Microlensing with advanced contour integration algorithm: Green's theorem to third order, error control, optimal sampling and limb darkening

GAMER: a GPU-Accelerated Adaptive Mesh Refinement Code for Astrophysics

Gravitational Microlensing of Gamma-Ray Bursts

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

An efficiency-improved GPU algorithm for the 2 + 2 + 1 method in nonlinear beamforming

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

Parallel GPU Optimization of the Shooting and Bouncing Ray Tracing Methodology for Propagation Modeling