Publications (42)
The Mira-Titan Universe II: Matter Power Spectrum Emulation
Earl Lawrence, Katrin Heitmann, Juliana Kwan +7
We introduce a new cosmic emulator for the matter power spectrum covering eight cosmological parameters. Targeted at optical surveys, the emulator provides accurate predictions out…
ClangJIT: Enhancing C++ with Just-in-Time Compilation
Hal Finkel, David Poliakoff, David F. Richards
The C++ programming language is not only a keystone of the high-performance-computing ecosystem but has proven to be a successful base for portable parallel-programming frameworks.…
A Proposal for Loop-Transformation Pragmas
Michael Kruse, Hal Finkel
Pragmas for loop transformations, such as unrolling, are implemented in most mainstream compilers. They are used by application programmers because of their ease of use compared to…
Simulations of the Pairwise Kinematic Sunyaev-Zel'dovich Signal
Samuel Flender, Lindsey Bleem, Hal Finkel +3
The pairwise kinematic Sunyaev-Zel'dovich (kSZ) signal from galaxy clusters is a probe of their line-of-sight momenta, and thus a potentially valuable source of cosmological inform…
Really Embedding Domain-Specific Languages into C++
Hal Finkel, Alexander McCaskey, Tobi Popoola +2
Domain-specific languages (DSLs) are both pervasive and powerful, but remain difficult to integrate into large projects. As a result, while DSLs can bring distinct advantages in pe…
The Q Continuum Simulation: Harnessing the Power of GPU Accelerated Supercomputers
Katrin Heitmann, Nicholas Frontiere, Chris Sewell +6
Modeling large-scale sky survey observations is a key driver for the continuing development of high resolution, large-volume, cosmological simulations. We report the first results…
Large-Scale Structure Formation with Massive Neutrinos and Dynamical Dark Energy
Amol Upadhye, Rahul Biswas, Adrian Pope +4
Over the next decade, cosmological measurements of the large-scale structure of the Universe will be sensitive to the combined effects of dynamical dark energy and massive neutrino…
Autotuning PolyBench Benchmarks with LLVM Clang/Polly Loop Optimization Pragmas Using Bayesian Optimization (extended version)
Xingfu Wu, Michael Kruse, Prasanna Balaprakash +4
In this paper, we develop a ytopt autotuning framework that leverages Bayesian optimization to explore the parameter space search and compare four different supervised learning met…
CosmoDC2: A Synthetic Sky Catalog for Dark Energy Science with LSST
Danila Korytov, Andrew Hearin, Eve Kovacs +27
This paper introduces cosmoDC2, a large synthetic galaxy catalog designed to support precision dark energy science with the Large Synoptic Survey Telescope (LSST). CosmoDC2 is the…
Halo Profiles and the Concentration-Mass Relation for a ÎCDM Universe
Hillary L. Child, Salman Habib, Katrin Heitmann +4
Profiles of dark matter-dominated halos at the group and cluster scales play an important role in modern cosmology. Using results from two very large cosmological -body simulati…
HACC: Simulating Sky Surveys on State-of-the-Art Supercomputing Architectures
Salman Habib, Adrian Pope, Hal Finkel +11
Current and future surveys of large-scale cosmic structure are associated with a massive and complex datastream to study, characterize, and ultimately understand the physics behind…
The Outer Rim Simulation: A Path to Many-Core Supercomputers
Katrin Heitmann, Hal Finkel, Adrian Pope +10
We describe the Outer Rim cosmological simulation, one of the largest high-resolution N-body simulations performed to date, aimed at promoting science to be carried out with large-…
Design and Use of Loop-Transformation Pragmas
Michael Kruse, Hal Finkel
Adding a pragma directive into the source code is arguably easier than rewriting it, for instance for loop unrolling. Moreover, if the application is maintained for multiple platfo…
Quantum Sensing for High Energy Physics
Zeeshan Ahmed, Yuri Alexeev, Giorgio Apollinari +113
Report of the first workshop to identify approaches and techniques in the domain of quantum sensing that can be utilized by future High Energy Physics applications to further the s…
HACC Cosmological Simulations: First Data Release
Katrin Heitmann, Thomas D. Uram, Hal Finkel +10
We describe the first major public data release from cosmological simulations carried out with Argonne's HACC code. This initial release covers a range of datasets from large gravi…
Oscillons After Inflation
Mustafa A. Amin, Richard Easther, Hal Finkel +2
Oscillons are massive, long-lived, localized excitations of a scalar field. We show that in a large class of well-motivated single-field models, inflation is followed by self-reson…
The Mira-Titan Universe. III. Emulation of the Halo Mass Function
Sebastian Bocquet, Katrin Heitmann, Salman Habib +5
We construct an emulator for the halo mass function over group and cluster mass scales for a range of cosmologies, including the effects of dynamical dark energy and massive neutri…
GRChombo : Numerical Relativity with Adaptive Mesh Refinement
Katy Clough, Pau Figueras, Hal Finkel +3
In this work, we introduce GRChombo: a new numerical relativity code which incorporates full adaptive mesh refinement (AMR) using block structured Berger-Rigoutsos grid generation.…
High Energy Physics Forum for Computational Excellence: Working Group Reports (I. Applications Software II. Software Libraries and Tools III. Systems)
Salman Habib, Robert Roser, Tom LeCompte +22
Computing plays an essential role in all aspects of high energy physics. As computational technology evolves rapidly in new directions, and data throughput and volume continue to f…
User-Directed Loop-Transformations in Clang
Michael Kruse, Hal Finkel
Directives for the compiler such as pragmas can help programmers to separate an algorithm's semantics from its optimization. This keeps the code understandable and easier to optimi…
GRChombo: An adaptable numerical relativity code for fundamental physics
Tomas Andrade, Llibert Areste Salo, Josu C. Aurrekoetxea +29
GRChombo is an open-source code for performing Numerical Relativity time evolutions, built on top of the publicly available Chombo software for the solution of PDEs. Whilst GRChomb…
Autotuning Search Space for Loop Transformations
Michael Kruse, Hal Finkel, Xingfu Wu
One of the challenges for optimizing compilers is to predict whether applying an optimization will improve its execution speed. Programmers may override the compiler's profitabilit…
Amplitude-Aware Lossy Compression for Quantum Circuit Simulation
Xin-Chuan Wu, Sheng Di, Franck Cappello +3
Classical simulation of quantum circuits is crucial for evaluating and validating the design of new quantum algorithms. However, the number of quantum state amplitudes increases ex…
Loop Optimization Framework
Michael Kruse, Hal Finkel
The LLVM compiler framework supports a selection of loop transformations such as vectorization, distribution and unrolling. Each transformation is carried-out by specialized passes…
Autotuning PolyBench Benchmarks with LLVM Clang/Polly Loop Optimization Pragmas Using Bayesian Optimization
Xingfu Wu, Michael Kruse, Prasanna Balaprakash +4
An autotuning is an approach that explores a search space of possible implementations/configurations of a kernel or an application by selecting and evaluating a subset of implement…
The Last Journey. I. An Extreme-Scale Simulation on the Mira Supercomputer
Katrin Heitmann, Nicholas Frontiere, Esteban Rangel +10
The Last Journey is a large-volume, gravity-only, cosmological N-body simulation evolving more than 1.24 trillion particles in a periodic box with a side-length of 5.025Gpc. It was…
Full-State Quantum Circuit Simulation by Using Data Compression
Xin-Chuan Wu, Sheng Di, Emma Maitreyee Dasgupta +4
Quantum circuit simulations are critical for evaluating quantum algorithms and machines. However, the number of state amplitudes required for full simulation increases exponentiall…
Stochastic Evolution of Graphs using Local Moves
Hal Finkel
Inspired by theories such as Loop Quantum Gravity, a class of stochastic graph dynamics was studied in an attempt to gain a better understanding of discrete relational systems unde…
Cosmic Emulation: Fast Predictions for the Galaxy Power Spectrum
Juliana Kwan, Katrin Heitmann, Salman Habib +4
The halo occupation distribution (HOD) approach has proven to be an effective method for modeling galaxy clustering and bias. In this approach, galaxies of a given type are probabi…
The Universe at Extreme Scale: Multi-Petaflop Sky Simulation on the BG/Q
Salman Habib, Vitali Morozov, Hal Finkel +9
Remarkable observational advances have established a compelling cross-validated model of the Universe. Yet, two key pillars of this model -- dark matter and dark energy -- remain m…
The differential transformation method and Miller's recurrence
Hal Finkel
The differential transformation method (DTM) enables the easy construction of a power-series solution to a nonlinear differential equation. The exponentiation operation has not bee…
An Iterated, Multipoint Differential Transform Method for Numerically Evolving PDE IVPs
Hal Finkel
Traditional numerical techniques for solving time-dependent partial-differential-equation (PDE) initial-value problems (IVPs) store a truncated representation of the function value…
The Mira-Titan Universe: Precision Predictions for Dark Energy Surveys
Katrin Heitmann, Derek Bingham, Earl Lawrence +8
Ground and space-based sky surveys enable powerful cosmological probes based on measurements of galaxy properties and the distribution of galaxies in the Universe. These probes inc…
Gravitational Waves from Oscillon Preheating
Shuang-Yong Zhou, Edmund J. Copeland, Richard Easther +3
Oscillons are long-lived, localized excitations of nonlinear scalar fields which may be copiously produced during preheating after inflation, leading to a possible oscillon-dominat…
Memory-Efficient Quantum Circuit Simulation by Using Lossy Data Compression
Xin-Chuan Wu, Sheng Di, Franck Cappello +3
In order to evaluate, validate, and refine the design of new quantum algorithms or quantum computers, researchers and developers need methods to assess their correctness and fideli…
Report of the Workshop on Program Synthesis for Scientific Computing
Hal Finkel, Ignacio Laguna
Program synthesis is an active research field in academia, national labs, and industry. Yet, work directly applicable to scientific computing, while having some impressive successe…
Inflaton Fragmentation and Oscillon Formation in Three Dimensions
Mustafa A. Amin, Richard Easther, Hal Finkel
Analytical arguments suggest that a large class of scalar field potentials permit the existence of oscillons -- pseudo-stable, non-topological solitons -- in three spatial dimensio…
Extending C++ for Heterogeneous Quantum-Classical Computing
Thien Nguyen, Anthony Santana, Tyler Kharazi +3
We present qcor - a language extension to C++ and compiler implementation that enables heterogeneous quantum-classical programming, compilation, and execution in a single-source co…
The Borg Cube Simulation: Cosmological Hydrodynamics with CRK-SPH
J. D. Emberson, Nicholas Frontiere, Salman Habib +4
A challenging requirement posed by next-generation observations is a firm theoretical grasp of the impact of baryons on structure formation. Cosmological hydrodynamic simulations m…
PSpectRe: A Pseudo-Spectral Code for (P)reheating
Richard Easther, Hal Finkel, Nathaniel Roth
PSpectRe is a C++ program that uses Fourier-space pseudo-spectral methods to evolve interacting scalar fields in an expanding universe. PSpectRe is optimized for the analysis of pa…
Redshift-space distortions in massive neutrino and evolving dark energy cosmologies
Amol Upadhye, Juliana Kwan, Adrian Pope +4
Large-scale structure surveys in the coming years will measure the redshift-space power spectrum to unprecedented accuracy, allowing for powerful new tests of the LambdaCDM picture…
Doing Moore with Less -- Leapfrogging Moore's Law with Inexactness for Supercomputing
Sven Leyffer, Stefan M. Wild, Mike Fagan +4
Energy and power consumption are major limitations to continued scaling of computing systems. Inexactness, where the quality of the solution can be traded for energy savings, has b…