CMACIONIZE 2.0: a novel task-based approach to Monte Carlo radiation transfer

Author:

Vandenbroucke B.ORCID,Camps P.ORCID

Abstract

Context. Monte Carlo radiative transfer (MCRT) is a widely used technique to model the interaction between radiation and a medium. It plays an important role in astrophysical modelling and when these models are compared with observations. Aims. We present a novel approach to MCRT that addresses the challenging memory-access patterns of traditional MCRT algorithms, which prevent an optimal performance of MCRT simulations on modern hardware with a complex memory architecture. Methods. We reformulated the MCRT photon-packet life cycle as a task-based algorithm, whereby the computation is broken down into small tasks that are executed concurrently. Photon packets are stored in intermediate buffers, and tasks propagate photon packets through small parts of the computational domain, moving them from one buffer to another in the process. Results. Using the implementation of the new algorithm in the photoionization MCRT code CMACIONIZE 2.0, we show that the decomposition of the MCRT grid into small parts leads to a significant performance gain during the photon-packet propagation phase, which constitutes the bulk of an MCRT algorithm because memory caches are used more efficiently. Our new algorithm is faster by a factor 2 to 4 than an equivalent traditional algorithm and shows good strong scaling up to 30 threads. We briefly discuss adjustments to our new algorithm and extensions to other astrophysical MCRT applications. Conclusions. We show that optimising the memory access patterns of a memory-bound algorithm such as MCRT can yield significant performance gains.

Funder

Belgian Federal Science Policy Office

Science and Technology Facilities Council

Publisher

EDP Sciences

Subject

Space and Planetary Science,Astronomy and Astrophysics

Reference32 articles.

1. Agostinelli S., Allison J., Amako K., et al. 2003, Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectorsand Associated Equipment, 506, 250

2. Random forcing of three-dimensional homogeneous turbulence

3. The impact of Lyman-α radiative transfer on large-scale clustering in the Illustris simulation

4. starbench: the D-type expansion of an H ii region

5. Bordner J., & Norman M. L. 2012, Proceedings of the Extreme Scaling Workshop, BW-XSEDE’12 (USA: University of Illinois at Urbana-Champaign)

Cited by 3 articles. 订阅此论文施引文献 订阅此论文施引文献,注册后可以免费订阅5篇论文的施引文献,订阅后可以查看论文全部施引文献

1. The persistence of high altitude non-equilibrium diffuse ionized gas in simulations of star-forming galaxies;Monthly Notices of the Royal Astronomical Society;2024-04-10

2. AREPO-MCRT: Monte Carlo Radiation Hydrodynamics on a Moving Mesh;The Astrophysical Journal;2020-12-09

3. Photoionization feedback in turbulent molecular clouds;Monthly Notices of the Royal Astronomical Society;2020-10-30

同舟云学术

1.学者识别学者识别

2.学术分析学术分析

3.人才评估人才评估

"同舟云学术"是以全球学者为主线,采集、加工和组织学术论文而形成的新型学术文献查询和分析系统,可以对全球学者进行文献检索和人才价值评估。用户可以通过关注某些学科领域的顶尖人物而持续追踪该领域的学科进展和研究前沿。经过近期的数据扩容,当前同舟云学术共收录了国内外主流学术期刊6万余种,收集的期刊论文及会议论文总量共计约1.5亿篇,并以每天添加12000余篇中外论文的速度递增。我们也可以为用户提供个性化、定制化的学者数据。欢迎来电咨询!咨询电话:010-8811{复制后删除}0370

www.globalauthorid.com

TOP

Copyright © 2019-2024 北京同舟云网络信息技术有限公司
京公网安备11010802033243号  京ICP备18003416号-3