Twinkle: A GPU-based binary-lens microlensing code with contour integration method
arXiv:2501.03322 · doi:10.3847/1538-4365/ad9b8d
Abstract
With the rapidly increasing rate of microlensing planet detections, microlensing modeling software faces significant challenges in computation efficiency. Here, we develop the Twinkle code, an efficient and robust binary-lens modeling software suite optimized for heterogeneous computing devices, especially GPUs. Existing microlensing codes have the issue of catastrophic cancellation that undermines the numerical stability and precision, and Twinkle resolves them by refining the coefficients of the binary-lens equation. We also devise an improved method for robustly identifying ghost images, thereby enhancing computational reliability. We have advanced the state of the art by optimizing Twinkle specifically for heterogeneous computing devices by taking into account the unique task and cache memory dispatching patterns of GPUs, while the compatibility with the traditional computing architectures of CPUs is still maintained. Twinkle has demonstrated an acceleration of approximately 2 orders of magnitude (>~100 times) on contemporary GPUs. The enhancement in computational speed of Twinkle will translate to the delivery of accurate and highly efficient data analysis for ongoing and upcoming microlensing projects. Both GPU and CPU versions of Twinkle are open-source and publicly available.
Accepted by ApJS, GitHub link: https://github.com/AsterLight0626/Twinkle
References in corpus (13)
- Real-Time Difference Imaging Analysis of MOA Galactic Bulge Observations During 2000
- Predictions of the WFIRST Microlensing Survey I: Bound Planet Detection Rates
- Planetary Detection Efficiency of the Magnification 3000 Microlensing Event OGLE-2004-BLG-343
- Microlensing Event MOA-2007-BLG-400: Exhuming the Buried Signature of a Cool, Jovian-Mass Planet
- Properties of Planetary Caustics in Gravitational Microlensing
- Extended-Source Effect and Chromaticity in Two-Point-Mass Microlensing
- Microlensing with advanced contour integration algorithm: Green's theorem to third order, error control, optimal sampling and limb darkening
- Hexadecapole Approximation in Planetary Microlensing
- An Efficient Method for Modeling High Magnification Planetary Microlensing Events
- VBBinaryLensing: a public package for microlensing light curve computation
- Adaptive Contouring -- an efficient way to calculate microlensing light curves of extended sources
- Light Curve Calculations for Triple Microlensing Systems
- Fast computation of quadrupole and hexadecapole approximations in microlensing with a single point-source evaluation