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.

Gravity.jl: Fast and accurate gravitational lens modeling in Julia. I. Point-like and linearized extended sources

Gravity.jl: fast and accurate gravitational lens modeling in Julia

Gravitational lens modelling in a citizen science context

Caustics: A Python Package for Accelerated Strong Gravitational Lensing Simulations

A model for galaxy-galaxy strong lensing statistics in surveys

Lensing of gravitational waves: efficient wave-optics methods and validation with symmetric lenses

Extrapolating the projected potential of gravitational lens models: property-preserving degeneracies

A model for galaxy–galaxy strong lensing statistics in surveys

Determination of Lens Mass Density Profile from Strongly-Lensed Gravitational-Wave Signals

DeepGraviLens: a Multi-Modal Architecture for Classifying Gravitational Lensing Data

Comparison of Strong Gravitational Lens Model Software III. Direct and indirect semi-independent lens model comparisons of COSMOS J095930+023427, SDSS J1320+1644, SDSSJ1430+4105 and J1000+0021

Effective Models for Statistical Studies of Galaxy-Scale Gravitational Lensing

PyAutoLens: Open-Source Strong Gravitational Lensing

Unveiling lens light complexity with a novel multi-Gaussian expansion approach for strong gravitational lensing

On the simulation of gravitational lensing

Along the Primary Curve: Simultaneous Source and Lens Reconstruction of Bright Arcs in Cluster Lenses

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

Strong Gravitational Lensing in Horndeski theory

Reconstructions of the Strong Gravitational Lenses MS 2137 and MS 1455 Using a Two-Stage Inversion Algorithm

Through the lens of Sgr : Identifying and resolving strongly lensed continuous gravitational waves beyond the Einstein radius