Eva Sciacca, Valentina Cesare, Nicola Tuccari, Fabio Vitello, Iacopo Colonnelli, Carmelita Carbone, Emiliano Tramontana, Ugo Becciani
Toward a Portable and Reproducible High-performance Visualization for Astronomy & Cosmology
16th International Workshop on Science Gateways (IWSG2024), Toulouse, France.doi.org/10.5281/zenodo.13858498
Modern Astronomy and Cosmology (A&C) generate petabyte-scale data volumes requiring the development of a new eneration of software tools to access, store, and process them. The Visualization Interface for the Virtual Observatory (VisIVO) tool performs multi-dimensional data analysis and knowledge discovery in multivariate astrophysical datasets. Thanks to containerization and virtualization technologies, VisIVO has already been exploited on top of several distributed computing infrastructures including the European Open Science Cloud (EOSC). In this work, we outline the execution strategies designed to enhance the portability and reproducibility of the VisIVO modular applications for high performance visualization of data generated on the (pre-)exascale systems by HPC A&C simulations performed with the GADGET (GAlaxies with Dark matter and Gas) code.
Download the article.
Giuliano Taffoni, Andrea Mignone, Luca Tornatore, Eva Sciacca, Massimiliano Guarrasi, Giovanni Lapenta, Lubomir Riha, Radim Vavrik, Ondrej Vysocky, Kristian Kadlubiak, Petr Strakos, Milan Jaros, Klaus Dolag, Benoit Commercon, Luciano Rezzolla, Khalil Pierre, Georgios Doulis, Sijing Shen, Manolis Marazakis, Daniele Gregori, Elisabetta Boella, Gino Perna, Marisa Zanotti, Erwan Raffin, Kai Polsterer, Sebastian Trujillo Gomez, Guillermo Marin
Redesign of astrophysical codes for exascale computing: the SPACE experience
Proceedings Volume 13101, Software and Cyberinfrastructure for Astronomy VIII; 131010X (2024) doi.org/10.1117/12.3020361
Event: SPIE Astronomical Telescopes + Instrumentation, 2024, Yokohama, Japan
High Performance Computing based simulations are crucial in Astrophysics and Cosmology, helping scientists investigate and understand complex astrophysical phenomena. Taking advantage of Exascale computing capabilities is essential for these efforts. However, the unprecedented architectural complexity of exascale systems impacts simulation codes. The SPACE Center of Excellence aims to re-engineer key astrophysical codes to adapt to these new computational challenges by adopting innovative programming paradigms and software solutions. Through co-design activities, SPACE brings together scientists, code developers, HPC experts, hardware manufacturers, and software developers. This collaboration enhances exascale astrophysics and cosmology applications, promoting the use of exascale and post-exascale computing capabilities. Additionally, SPACE addresses high-performance data analysis for the massive data outputs from exascale simulations, using machine learning and visualization tools. The project facilitates application deployment across platforms by focusing on code repositories and data sharing, integrating European astrophysical communities around exascale computing with standardized software and data protocols. In this paper, we present the SPACE Center of Excellence and the preliminary results achieved by the project.
Find out more here.
CINECA HPC Team
HPC Annual Report 2023/2024
CINECA HPC is the main Centre for scientific computing in Italy.
It runs a large computing infrastructure and makes it available to Italian and European Researchers, as well as to important Italian Companies within a programme for supporting national industrial competitiveness.
This report presents data about usage and users, projects, events and classes offered in the field of High Performance Computing.
SPACE CoE is one of the EU projects CINECA is involved in and such engagement and information about our Centre of Excellence is reported in this document (pags. 88-91).
Download the full report (pdf format) or find out more here.
V. Berta, A. Mignone, M. Bugli, G. Mattia
“A 4th-order accurate finite volume method for ideal classical and special relativistic MHD based on pointwise reconstructions”,
Journal of Computational Physics 499 (2024) 112701
We present a novel implementation of a genuinely 4th-order accurate finite volume scheme
for multidimensional classical and special relativistic magnetohydrodynamics (MHD) based on
the constrained transport (CT) formalism.
The scheme introduces several novel aspects when
compared to its predecessors yielding a more efficient computational tool.
Among the most relevant ones, our scheme exploits pointwise to pointwise reconstructions (rather than onedimensional
finite volume ones), employs the generic upwind constrained transport averaging
and sophisticated limiting strategies that include both a discontinuity detector and an order
reduction procedure. Selected numerical benchmarks demonstrate the accuracy and robustness
of the method.
Download the article (pdf format).
Bernhard Baumschlager, Sijing Shen, and James W. Wadsley
“Spectral reconstruction for radiation hydrodynamic simulations of galaxy evolution”,
Astronomy & Astrophysics manuscript no. 10BandRT - October 27, 2023
Radiation from stars and active galactic nuclei (AGN) plays an important role in galaxy formation and evolution, and profoundly
transforms the intergalactic, circumgalactic and interstellar medium (IGM, CGM & ISM). On-the-fly radiative transfer (RT) has
started being incorporated in cosmological simulations, but the complex, evolving radiation spectra are often crudely approximated
with a small number of broad bands with piece-wise constant intensity and a fixed photo-ionisation cross-section.
Such a treatment is unable to capture the changes to the spectrum as light is absorbed while it propagates through a medium with non-zero opacity. This can lead to large errors in photo-ionisation and heating rates. In this work, we present a novel approach of discretising the radiation field in narrow bands, located at the edges of the typically used bands, in order to capture the power-law slope of the radiation field.
In combination with power-law approximations for the photo-ionisation cross-sections, this model allows us to self-consistently combine
radiation from sources with different spectra and accurately follow the ionisation states of primordial and metal species through time.
The method is implemented in Gasoline2 in connection with Trevr2, a fast reverse RT algorithm with O(Nactive log2 N) scaling.
We compare our new piece-wise power-law reconstruction to the piece-wise constant method in calculating the primordial chemistry
photo-ionisation and heating rates under an evolving UV-background (UVB) and stellar spectrum, and find that our method reduces
errors significantly, up to two orders of magnitude in the case of HeII ionisation.
We apply our new spectral reconstruction method in RT post-processing of a cosmological zoom-in simulation, MUGS2 g1536 (Keller et al. 2016), including radiation from stars and a
live UVB, and find a significant increase in total neutral hydrogen (HI) mass in the ISM and the CGM due to shielding of the UVB
and a low escape fraction of the stellar radiation. This demonstrates the importance of RT and an accurate spectral approximation in
simulating the CGM-galaxy ecosystem.
Download the article (pdf format).
M. Zanotti, A. Mignone, and M. Marazakis
“SPACE: Space Parallel Astrophysical Codes for Exascale”,
ERCIM News No. 134, July 2023, pp. 43-44
In astrophysics and cosmology (A&C), high-performance computing (HPC)-based numerical simulations are invaluable
instruments to support scientific discovery.
The efficient and effective exploitation of exascale (and beyond)
computing capabilities will be key to paving the way for scientific discovery in this scientific domain, but requires a
coordinated effort towards adapting A&C codes for space applications on exascale HPC systems.
Go to the article webpage or download the magazine (pdf format).
