Related papers: THOR: a GPU-accelerated and MPI-parallel radiative…
The numerical simulations of massive collisional stellar systems, such as globular clusters (GCs), are very time-consuming. Until now, only a few realistic million-body simulations of GCs with a small fraction of binaries (5%) have been…
We present the newly developed code, GAMER (GPU-accelerated Adaptive MEsh Refinement code), which has adopted a novel approach to improve the performance of adaptive mesh refinement (AMR) astrophysical simulations by a large factor with the…
The Lyman alpha (LyA) line of neutral hydrogen plays a central role in observations of star-forming galaxies. However, resonant scattering makes it difficult to directly interpret LyA signatures. Monte Carlo radiative transfer (MCRT)…
The emergence of three-dimensional magneto-hydrodynamic (MHD) simulations of stellar atmospheres has sparked a need for efficient radiative transfer codes to calculate detailed synthetic spectra. We present RH 1.5D, a massively parallel…
The Event Horizon Telescope (EHT) has recently released high-resolution images of accretion flows onto two supermassive black holes. Our physical understanding of these images depends on accuracy and precision of numerical models of plasma…
Ray-tracing is a central tool for constructing mock observations of compact object emission and for comparing physical emission models with observations. We present Arcmancer, a publicly available general ray-tracing and tensor algebra…
Aims. We present MCFOST-art, a new non-local thermodynamic equilibrium radiative transfer solver for multilevel atomic systems. The code is embedded in the 3D radiative transfer code MCFOST and is compatible with most of the MCFOST modules.…
The study of planetary atmospheres is crucial for understanding the origin, evolution, and processes that shape celestial bodies like planets, moons and comets. The interpretation of planetary spectra requires a detailed understanding of…
MPI derived datatypes are an abstraction that simplifies handling of non-contiguous data in MPI applications. These datatypes are recursively constructed at runtime from primitive Named Types defined in the MPI standard. More recently, the…
The multi-resolution approximation (MRA) of Gaussian processes was recently proposed to conduct likelihood-based inference for massive spatial data sets. An advantage of the methodology is that it can be parallelized. We implemented the MRA…
Using the radiative transfer code Torus, we produce spectral-line cubes of the predicted HI profile from global SPH simulations of spiral galaxies. Torus grids the SPH galaxy using Adaptive Mesh Refinement, then applies a ray-tracing method…
We develop a new three-dimensional time-dependent radiative transfer code, TRINITY (Time-dependent Radiative transfer In Near-Infrared TomographY), for in-vivo diffuse optical tomography (DOT). The simulation code is based on the design of…
Higher resolution telescopes as well as 3D numerical simulations will require the development of detailed 3D radiative transfer calculations. Building upon our previous work we extend our method to include both continuum and line transfer.…
We present an implementation of an adaptive ray tracing (ART) module in the Athena hydrodynamics code that accurately and efficiently handles the radiative transfer involving multiple point sources on a three-dimensional Cartesian grid. We…
In this paper we focus on the integration of high-performance numerical libraries in ab initio codes and the portability of performance and scalability. The target of our work is FLEUR, a software for electronic structure calculations…
The development of radiation hydrodynamical methods that are able to follow gas dynamics and radiative transfer self-consistently is key to the solution of many problems in numerical astrophysics. Such fluid flows are highly complex, rarely…
We introduce Monte Carlo-based non-LTE line radiative transfer calculations in the 3D dust radiative transfer code SKIRT, which was originally set up as a dust radiative transfer code. By doing so, we develop a generic and powerful 3D…
We present EOS, a procedure for determining the Outgoing Longwave Radiation (OLR) and top-of-atmosphere (TOA) albedo for a wide range of conditions expected to be present in the atmospheres of rocky planets with temperate conditions. EOS is…
Star-forming regions in galaxies are surrounded by vast reservoirs of gas capable of both emitting and absorbing Lyman-alpha (Lya) radiation. Observations of Lya emitters and spatially extended Lya haloes indeed provide insights into the…
We describe a new implementation of a parallel Tree-SPH code with the aim to simulate Galaxy Formation and Evolution. The code has been parallelized using SHMEM, a Cray proprietary library to handle communications between the 256 processors…