ELUCID-DESI I: A Parallel MPI Implementation of the Initial Condition Solver for Large-Scale Reconstruction Simulations

Clicks: 2
ID: 318755
2026
Free PDF
Article Quality & Performance Metrics
Overall Quality
Not rated
Combines reader engagement with the AI quality analysis. This article has not been analysed, so there is no overall score — reader engagement is measured and shown alongside.
AI Quality Assessment
Not analyzed
Readership in this journal
Emerging

Ranked #842 of 913 articles by views in monthly notices of the royal astronomical society

Most read Least read

Bar heights use a square-root scale. Only the 120 most-read articles are drawn; the journal has 913 in total.

Mint this article as an NFT
Not yet minted

Create a permanent, verifiable on-chain record of this article on the Scimatic Network. The NFT is held in your Journament account, and you can withdraw it to your own wallet at any time.

5 SUSD one-off · no wallet required
Sign in to mint
Abstract
Abstract We present a highly scalable, MPI-parallelized framework for reconstructing the initial cosmic density field, designed to meet the computational demands of next-generation cosmological simulations, particularly the upcoming ELUCID-DESI simulation based on DESI BGS data. Building upon the Hamiltonian Monte Carlo approach and the FastPM solver, our code employs domain decomposition to efficiently distribute memory between nodes. Although communication overhead increases the per-step runtime of the MPI version by roughly a factor of eight relative to the shared-memory implementation, our scaling tests-spanning different particle numbers, core counts, and node layouts-show nearly linear scaling with respect to both the number of particles and the number of CPU cores. Furthermore, to significantly reduce computational costs during the initial burn-in phase, we introduce a novel ‘guess’ module that rapidly generates a high-quality initial density field. The results of the simulation test confirm substantial efficiency gains: for 2563 particles, 53 steps (∼ 54 core hours) are saved, accelerating convergence by a factor of ∼ 18; for 10243, 106 steps (∼7500 core hours), achieving a speedup factor of ∼ 3. The total core hour gain grows with the number of particles, rendering large-volume reconstructions computationally practical for upcoming surveys, including our planned ELUCID-DESI reconstruction simulation with 40963 particles. We estimate that achieving convergence for this scale (targeting DESI-BGS data) requires about 800 HMCMC steps (∼ 5 million core hours). Our initial guess module will save approximately 360 steps (∼2.3 million core hours), reducing the total computational time by about 45 %.
Reference Key
openalex_W7165857413 Use this key to autocite in the manuscript while using SciMatic Manuscript Manager or Thesis Manager
Authors Wensheng Hong, X L Yang, J Li, H G Wang, Zhao Chen, Hong-Ming Zhu, Qingyang Li, Yizhou Gu, Youcai Zhang, Feng Shi, Jiaxin Han, Yu Yu, Zhongxu Zhai
Journal monthly notices of the royal astronomical society
Year 2026
DOI
10.1093/mnras/stag1211
URL
https://doi.org/10.1093/mnras/stag1211
Keywords Keywords not found

Citations

No citations found. To add a citation, contact the admin at info@scimatic.org

Comments

No comments yet. Be the first to comment on this article.