AKRA 3.0: A matrix-free Inversion Framework for Weak Lensing Mass Mapping and Its Application to DES Y3 Data
Abstract
Weak gravitational lensing mass mapping offers a direct probe of the matter distribution. Accurate reconstruction of mass maps from masked shear catalogs remains challenging due to survey boundaries and spatially varying noise. In AKRA 2.0, we addressed the mask problem on the curved sky by constructing and inverting the normal-equation matrix H ATN-1 A explicitly, necessitating a split-scale strategy that reconstructed different angular scales independently to reach high resolution. Here we present AKRA 3.0, in which H is treated as a linear operator and the normal equations are solved by the conjugate gradient (CG) method. This reformulation reduces the memory requirement from O(N2) to O(N) and the inversion cost from O(N3) to O(N iterN3/2), N max2 for full-sky (SHT-based) operations. Such optimizations render high-resolution full-sky reconstruction tractable for Stage~III and Stage~IV surveys. Applying AKRA 3.0 to the DES Y3 METACALIBRATION catalog, we produce the highest-resolution convergence map of this dataset to date at HEALPix Nnside= 2048 without imposing any prior assumptions. We extract the convergence power spectrum directly from the reconstructed map and demonstrate that unbiased two-point measurements can be obtained directly from the reconstructed map. The reconstructed E-mode convergence map will be publicly released as data products to enable future studies of non-Gaussian statistics, higher-order moments, and cross-correlations with external datasets.
Turn this paper into a lesson
ArcXiv compiles a structured reading guide from this paper's metadata: plain-English importance, contributions, prerequisite concepts, which sections to read first, flashcards, and a quiz. Grounded in the abstract, never invented.