Search NASA⌕ Search

SEARCH · Search NASA

Results for “concurrent data structures”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

Data-flow parallelism for high-energy and nuclear physics computing frameworks

The processing tasks of a scientific workflow in high-energy and nuclear physics (HENP) can typically be represented as a directed acyclic graph formed according to the data flow—i.e. the data dependencies among algorithms executed as part of the workflow. With this representation, an HENP computing framework can optimally execute a workflow, exploiting the parallelism inherent among independent tasks. Despite such a natural description of a workflow, most HENP frameworks do not make use of technologies that provide concurrent execution of graph-based tasking structures. In this session, we describe Fermilab efforts to adopt a graph-based technology (specifically Intel’s oneTBB flow graph) for meeting the framework needs of its experiments, notably DUNE. After introducing the physics DUNE intends to explore, we will show that all common processing idioms supported by current HENP frameworks can naturally be supported by oneTBB’s data-flow technology, optimally leveraging the concurrent capabilities of the machine. In addition, we discuss collaborative efforts between Fermilab and the Intel oneTBB development team, who is considering improvements to the flow-graph technology to better support HENP use cases.

43 PARTICLE ACCELERATORS↗

Three-dimensional lattice modulations in the charge density wave system Lu 2 Ir 3 Si 5

Using total and resonant x-ray scattering coupled to large-scale computer modeling, we study the lattice modulations in the complex charge density wave (CDW) material Lu 2 ⁢Ir 3 ⁢Si 5 . Here, we find that it is a unique quantum system where periodic lattice modulations related to emergent CDW order occur in three orthogonal atomic planes of the crystal lattice, leading to the emergence of an unusual three-dimensional (3D) pattern of short and long Ir-Ir and Lu-Lu bonds. The 3D character of observed lattice modulations explains the largely isotropic character of the changes in the electronic properties occurring when the CDW order sets in, demonstrating the strong electron-lattice coupling in Lu 2 ⁢Ir 3 ⁢Si 5 . The result is supported by DFT calculations based on the experimental structure data. Altogether, our work provides strong evidence for the presence of a relationship between the dimensionality of emergent lattice distortions and that of concurrent changes in the electronic properties of CDW materials. The relationship may need to be accounted for when these materials are explored for practical applications.

36 MATERIALS SCIENCE↗

Ultrafast x-ray scattering of photodissociation dynamics in 2-iodothiophene

Here, a time-resolved x-ray scattering (TRXS) investigation of the photodissociation dynamics of gas-phase 2-iodothiophene molecules following 252 nm excitation is presented. Structural evolution of the molecule and dynamical information on the resulting photofragments were captured using femtosecond x-ray free-electron laser pulses. Two dissociation pathways were identified, arising via excitation to ππ* and (n/π)σ* states, respectively, yielding distinct interfragment recoil velocities of ∼6.4 Åps −1 and 17.0 Åps −1 . A comparison of asymptotic scattering data with simulated patterns indicates that the thiophene ring remains closed following dissociation at this wavelength. Modeling the experimental data yields a branching ratio of ∼3:1 in favor of the high velocity channel. These findings demonstrate the capability of TRXS to disentangle concurrent ultrafast reaction pathways and provide detailed structural insight into energy redistribution during photoinduced bond fission in complex molecular systems.

Razmus, Weronika O. [Univ. of Southampton (United ↗

PSelInv—A Distributed Memory Parallel Algorithm for Selected Inversion: The Symmetric Case

We describe an efficient parallel implementation of the selected inversion algorithm for distributed memory computer systems, which we call PSelInv. The PSelInv method computes selected elements of a general sparse matrix Athat can be decomposed as A= LU, where L is lower triangular and U is upper triangular. The implementation described in this article focuses on the case of sparse symmetric matrices. It contains an interface that is compatible with the distributed memory parallel sparse direct factorization SuperLU-DIST. Yet, the underlying data structure and design of PSelInv allows it to be easily combined with other factorization routines, such as PARDISO. We discuss general parallelization strategies such as data and task distribution schemes. Specifically, we describe how to exploit the concurrency exposed by the elimination tree associated with the LU factorization of A. We demonstrate the efficiency and accuracy of PSelInv by presenting several numerical experiments. In particular, we show that PSelInv can run efficiently on more than 4,000 cores for a modestly sized matrix. We also show how PSelInv can be used to accelerate large-scale electronic structure calculations.

97 MATHEMATICS AND COMPUTING↗

Scalable Incremental Checkpointing using GPU-Accelerated De-Duplication

Writing large amounts of data concurrently to stable storage is a typical I/O pattern of many HPC workflows. This pattern introduces high I/O overheads and results in increased storage space utilization especially for workflows that need to capture the evolution of data structures with high frequency as checkpoints. In this context, many applications, such as graph pattern matching, perform sparse updates to large data structures between checkpoints. For these applications, incremental checkpointing techniques that save only the differences from one checkpoint to another can dramatically reduce the checkpoint sizes, I/O bottlenecks, and storage space utilization. However, such techniques are not without challenges: it is non-trivial to transparently determine what data has changed since a previous checkpoint and assemble the differences in a compact fashion that does not result in excessive metadata. State-of-art data reduction techniques (e.g., compression and de-duplication) have significant limitations when applied to modern HPC applications that leverage GPUs: slow at detecting the differences, generate a large amount of metadata to keep track of the differences, and ignore crucial spatiotemporal checkpoint data redundancy. This paper addresses these challenges by proposing a Merkle tree-based incremental checkpointing method to exploit GPUs' high memory bandwidth and massive parallelism. Experimental results at scale show a significant reduction of the I/O overhead and space utilization of checkpointing compared with state-of-the-art incremental checkpointing and compression techniques.

Tan, Nigel↗

Accelerating nuclear configuration interaction calculations through a preconditioned block iterative eigensolver

We describe a number of recently developed techniques for improving the performance of large-scalenuclear configuration interaction calculations on high performance parallel computers. We show thebenefit of using a preconditioned block iterative method to replace the Lanczos algorithm that hastraditionally been used to perform this type of computation. The rapid convergence of the block iterativemethod is achieved by a proper choice of starting guesses of the eigenvectors and the construction of aneffective preconditioner. These acceleration techniques take advantage of special structure of the nuclearconfiguration interaction problem which we discuss in detail. The use of a block method also allows usto improve the concurrency of the computation, and take advantage of the memory hierarchy of modernmicroprocessors to increase the arithmetic intensity of the computation relative to data movement. Wealso discuss the implementation details that are critical to achieving high performance on massivelyparallel multi-core supercomputers, and demonstrate that the new block iterative solver is two to threetimes faster than the Lanczos based algorithm for problems of moderate sizes on a Cray XC30 system.

Shao, Meiyue↗

Combining the Best of Two Chelating Titans: A Hydroxypyridinone-Decorated Macrocyclic Ligand for Efficient and Concomitant Complexation and Sensitized Luminescence of f-Elements

An ideal chelator for f-elements features rapid kinetics of complexation, high thermodynamic stability, and slow kinetics of dissociation. Here we present the facile synthesis of a macrocyclic ligand bearing four 1-hydroxy-2-pyridinone units linked to a cyclen scaffold that rapidly forms thermodynamically stable complexes with lanthanides (Sm 3+ , Eu 3+ , Tb 3+ , Dy 3+ ) and a representative late actinide (Cm 3+ ) in aqueous media and concurrently sensitizes them. Extended X-ray absorption fine structure (EXAFS) spectroscopy revealed an increase in the Ln/An-O bond lengths following the trend Cm>Eu>Tb and EXAFS data were compatible with time-resolved luminescence studies, which indicated one to two water molecules in the inner metal coordination sphere of Eu(III) and two water molecules for the Cm(III) complex. Spectrofluorimetric ligand competition titrations against DTPA confirmed the high thermodynamic stability of DOTHOPO complexes, with pM values between 19.9(1) and 21.9(2).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

In Situ Insights into Cathode Calcination for Predictive Synthesis: Kinetic Crystallization of LiNiO 2 from Hydroxides

Abstract Calcination is a solid‐state synthesis process widely deployed in battery cathode manufacturing. However, its inherent complexity associated with elusive intermediates hinders the predictive synthesis of high‐performance cathode materials. Here, correlative in situ X‐ray absorption/scattering spectroscopy is used to investigate the calcination of nickel‐based cathodes, focusing specifically on the archetypal LiNiO 2 from Ni(OH) 2 . Combining in situ observation with data‐driven analysis reveals concurrent lithiation and dehydration of Ni(OH) 2 and consequently, the low‐temperature crystallization of layered LiNiO 2 alongside lithiated rocksalts. Following early nucleation, LiNiO 2 undergoes sluggish crystallization and structural ordering while depleting rocksalts; ultimately, it turns into a structurally‐ordered layered phase upon full lithiation but remains small in size. Subsequent high‐temperature sintering induces rapid crystal growth, accompanied by undesired delithiation and structural degradation. These observations are further corroborated by mesoscale modeling, emphasizing that, even though calcination is thermally driven and favors transformation towards thermodynamically equilibrium phases, the actual phase propagation and crystallization can be kinetically tuned via lithiation, providing freedom for structural and morphological control during cathode calcination.

36 MATERIALS SCIENCE↗

National Cancer Institute Think-Tank Meeting Report on Proteomic Cartography and Biomarkers at the Single-Cell Level: Interrogation of Premalignant Lesions

A Think-Tank Meeting was convened by the National Cancer Institute (NCI) to solicit experts’ opinion on the development and application of multi-omic single-cell analyses, and especially single-cell proteomics, to improve the development of a new generation of biomarkers for cancer risk, early detection, diagnosis, prognosis as well as to discuss the discovery of new targets for prevention and therapy. It is anticipated that such markers and targets will be based on cellular, subcellular, molecular and functional aberrations within the lesion and within individual cells. Additionally, single-cell proteomic data will be essential for the establishment of new tools that with searchable and scalable features that include spatial and temporal cartographies of premalignant and malignant lesions. Challenges and potential solutions that were discussed included: • The best way/s to analyze single-cells from fresh and preserved tissue • Detection and analysis of secreted molecules and from single cells, especially from a tissue slice • Detection of new, previously undocumented cell type/s in the premalignant and early stage cancer tissue microenvironment • Multi-omic integration of data to support and inform proteomic measurements • Subcellular organelles – identifying abnormal structure, function, distribution and location within individual premalignant and malignant cells • How to improve the dynamic range of single-cell proteomic measurements for discovery of differentially expressed proteins and their post translational modifications (PTM) • The depth of coverage measured concurrently using single-cell techniques • Quantitation - absolute or semiquantitative? • Single methodology or multiplexed combinations? • Application of analytical methods for identification of biologically significant subsets • Data visualization of N-dimensional datasets • How to construct intercellular signaling networks in individual cells within premalignant tumor microenvironments (TME) • Associations between intrinsic cellular processes and extrinsic stimuli • How to predict cellular responses to stress inducing stimuli • Identification of new markers for prediction of progression from precursor, benign and localized lesions to invasive cancer, based on spatial and temporal changes within individual cells • Identification of new targets for immunoprevention or immunotherapy – identification of neoantigens and surfactome of individual cells within a lesion

59 BASIC BIOLOGICAL SCIENCES↗

Glass-forming ability of TeO{sub 2} and temperature induced changes on the structure of the glassy, supercooled, and molten states

Polarized (VV) and depolarized (VH) Raman spectra are obtained for glassy, supercooled, and molten TeO{sub 2} at temperatures up to 1000 K in order to elucidate the temperature evolution of the pertinent structural and vibrational properties. The intrinsic tendency of the system for crystallization is avoided by means of a newly applied protocol, thereby enabling the recording of Raman spectra of pure TeO{sub 2} on going from the molten to the supercooled liquid and to the room temperature glass states. Following an appropriate fitting procedure, the revealed bands are assigned to specific modes of structural polymorphs. A weak polarised band at ∼880 cm{sup −1} is assigned to Te=O terminal stretching in agreement with the literature ab initio molecular orbital calculations. Subtle changes to the relative band intensities within the 550-900 cm{sup −1} stretching region are caused by temperature increase. The network-like structure of the glass/melt is composed by TeO{sub 4} trigonal bipyramid and TeO{sub 3} trigonal pyramid units. With increasing temperature, TeO{sub 4} units convert to TeO{sub 3} units with a concurrent increase in the number of Te=O sites resulting from cleavages within the network structure. The fraction of such terminal oxygen atoms has been directly estimated from the spectroscopic data. The relative populations of the basic building blocks and the average number of O atoms around Te have been estimated for a wide temperature range directly from the Raman spectra, implying a gradual transformation of TeO{sub 4/2} to TeO{sub 2/2}(= O) trigonal pyramid units. The results are discussed in the context of the current phenomenological and theoretical status of the field.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

GENFIRE: A generalized Fourier iterative reconstruction algorithm for high-resolution 3D imaging

Tomography has made a radical impact on diverse fields ranging from the study of 3D atomic arrangements in matter to the study of human health in medicine. Despite its very diverse applications, the core of tomography remains the same, that is, a mathematical method must be implemented to reconstruct the 3D structure of an object from a number of 2D projections. Here, we present the mathematical implementation of a tomographic algorithm, termed GENeralized Fourier Iterative REconstruction (GENFIRE), for high-resolution 3D reconstruction from a limited number of 2D projections. GENFIRE first assembles a 3D Fourier grid with oversampling and then iterates between real and reciprocal space to search for a global solution that is concurrently consistent with the measured data and general physical constraints. The algorithm requires minimal human intervention and also incorporates angular refinement to reduce the tilt angle error. We demonstrate that GENFIRE can produce superior results relative to several other popular tomographic reconstruction techniques through numerical simulations and by experimentally reconstructing the 3D structure of a porous material and a frozen-hydrated marine cyanobacterium. As a result, equipped with a graphical user interface, GENFIRE is freely available from our website and is expected to find broad applications across different disciplines.

36 MATERIALS SCIENCE↗

Vectorization of Random Number Generation and Reproducibility of Concurrent Particle Transport Simulation

Efficient random number generation with high quality statistical properties and exact reproducibility of Monte Carlo simulations are important requirements in many areas of computational science. VecRNG is a package providing pseudo-random number generation (pRNG) in the context of a new library VecMath. This library bundles up several general-purpose mathematical utilities, data structures, and algorithms having both SIMD and SIMT (GPUs) support based on VecCore.~Several state-of-the-art RNG algorithms are implemented as kernels supporting parallel generation of random numbers in scalar, vector, and Cuda workflows. In this report, we will present design considerations, implementation details, and computing performance of parallel pRNG engines on both CPU and GPU. Reproducibility of propagating multiple particles in parallel for HEP event simulation is demonstrated, using GeantV-based examples, for both sequential and fine-grain track-level concurrent simulation wor kflows. Strategies for efficient uses of vectorized pRNG and non-overlapping streams of random number sequences in concurrent computing environments is discussed as well.

Jun, S. Y.↗

The Influence of Environment on Post-Detonation Chemistry and Debris Formation (Abbreviated Final Report: 20-SI-006)

Predicting, responding to, or interpreting the chemical record preserved in debris derived from nuclear events can be challenging due to chemical fractionation. Chemical fractionation is where different species of the evolving radionuclide inventory segregate and/or are lost from the system over the timescales of debris formation. Both historic data and recent research suggest that the interaction and character of the local environment may exert controls on chemical fractionation by influencing the cooling and evolution of the associated fireball as well as the composition of the vapor term and resultant speciation. Prior to this work, an integrated platform permitting dynamic and concurrent consideration of physical and chemical evolution of early time post-detonation event environments did not exist. Our work merged historic data and experimental approaches to support development of a computational framework able to simulate fundamental processes (e.g., entrainment of local environment, oxidation chemistry, and cooling time scales) that may perturb the radionuclide inventory captured in post-detonation debris. Work with historic debris confirmed that entrained environmental material affect debris composition, structure, and radionuclide incorporation. Complementary work utilizing a readily controllable and tunable benchtop setup (a plasma flow reactor) simulated the late cooling of a nuclear fireball (e.g., T < 6000 K) and bounded the sensitivity of actinide speciation and particle size distribution to variations in oxygen concentration and cooling rates. Concurrent laser ablation and laser heating experiments were used to investigate the chemistry and physics of processes occurring in vaporized and/or rapidly heated actinides and other elements in the presence of oxygen. A more computationally efficient microphysical model was developed for predicting and evolving size distributions of particles forming from mixed vapor terms and simulating particle formation processes under a variety of extreme conditions. Continued study of historic nuclear event film confirmed that shockwave data and physics codes agree to within the uncertainty of the data. Good agreement was achieved for thermal emission from an airburst, however the paucity of low-temperature molecular opacity data for mixtures of air, bomb debris, entrained dirt, and water vapor complicate agreement for more elaborate scenarios. A multiphysics code (ALE3D) was modified to bring the necessary physics and chemistry, including these new data and insights, onto a single platform. Code development included improved initialization of large physical systems, modernization of chemistry capabilities, and modifications to enable inclusion of particle transport.

07 ISOTOPE AND RADIATION SOURCES↗

Enabling Scalable and Extensible Memory-mapped Datastores in Userspace

Exascale workloads are expected to incorporate data-intensive processing in close coordination with traditional physics simulations. These emerging scientific, data-analytics and machine learning applications need to access a wide variety of datastores in flat files and structured databases. Programmer productivity is greatly enhanced by mapping datastores into the application process's virtual memory space to provide a unified “in-memory” interface. Currently, memory mapping is provided by system software primarily designed for generality and reliability. However, scalability at high concurrency is a formidable challenge on exascale systems. Also, there is a need for extensibility to support new datastores potentially requiring HPC data transfer services. In this article, we present UMap , a scalable and extensible userspace service for memory-mapping datastores. Furthermore, through decoupled queue management, concurrency aware adaptation, and dynamic load balancing, UMap enables application performance to scale even at high concurrency. We evaluate UMap in data-intensive applications, including sorting, graph traversal, database operations, and metagenomic analytics. Our results show that UMap as a userspace service outperforms an optimized kernel-based service across a wide range of intra-node concurrency by 1.22-1.9 × . We performed two case studies to demonstrate UMap 's extensibility. First, a new datastore residing in remote memory is incorporated into UMap as an application-specific plugin. Second, we present a persistent memory allocator Metall built atop UMap for unified storage/memory.

97 MATHEMATICS AND COMPUTING↗

MassiveGNN: Efficient Training via Prefetching for Massively Connected Distributed Graphs

Graph Neural Networks (GNN) are indispensable in learning from graph-structured data, yet their rising computational costs, especially on massively connected graphs, pose significant challenges in terms of execution performance. To tackle this, distributed-memory solutions such as partitioning the graph to concurrently train multiple replicas of GNNs are in practice. However, approaches requiring a partitioned graph usually suffer from communication overhead and load imbalance, even under optimal partitioning and communication strategies due to irregularities in the neighborhood minibatch sampling. This paper proposes practical trade-offs for improving the sampling and communication overheads for representation learn- ing on distributed graphs (using popular GraphSAGE architecture) by developing a parameterized prefetch and eviction scheme on top of the state-of-the-art Amazon DistDGL distributed GNN framework, demonstrating about 15–40% improvement in end-to-end training performance on the NERSC Perlmutter supercomputer for various OGB datasets.

Machine Leanring, high performance comptuing, grap↗

Hardware accelerated dynamic work creation on a graphics processing unit

A processor core is configured to execute a parent task that is described by a data structure stored in a memory. A coprocessor is configured to dispatch a child task to the at least one processor core in response to the coprocessor receiving a request from the parent task concurrently with the parent task executing on the at least one processor core. In some cases, the parent task registers the child task in a task pool and the child task is a future task that is configured to monitor a completion object and enqueue another task associated with the future task in response to detecting the completion object. The future task is configured to self-enqueue by adding a continuation future task to a continuation queue for subsequent execution in response to the future task failing to detect the completion object.

Gutierrez, Anthony↗

Hardware accelerated dynamic work creation on a graphics processing unit

A processor core is configured to execute a parent task that is described by a data structure stored in a memory. A coprocessor is configured to dispatch a child task to the at least one processor core in response to the coprocessor receiving a request from the parent task concurrently with the parent task executing on the at least one processor core. In some cases, the parent task registers the child task in a task pool and the child task is a future task that is configured to monitor a completion object and enqueue another task associated with the future task in response to detecting the completion object. The future task is configured to self-enqueue by adding a continuation future task to a continuation queue for subsequent execution in response to the future task failing to detect the completion object.

Gutierrez, Anthony T.↗

Hardware-Accelerated Ray Tracing of CAD-Based Geometry for Monte Carlo Radiation Transport

Monte Carlo radiation transport (MCRT) methods have been used to simulate radiation environments for many decades by tracking individual particles through a model to accumulate statistical information. MCRT geometry is historically formed using the constructive solid geometry (CSG). Recently, significant work has been performed to support simulations using computer-aided design (CAD)-based tessellated surfaces to support highly complex geometries. Ray tracing acceleration data structures from the rendering and visualization community are applied to accelerate particle tracking in CAD-based models. Despite these efforts, CSG representations provide the superior performance in surface intersection operations during particle flight. Concurrently, pseudo Monte Carlo methods have become prevalent in rendering applications to support more realistic models for scattering media, motivating innovations that are advantageous for MCRT simulations. Finally, the authors’ work extends these innovations by employing Intel’s Embree ray tracing kernel within a geometry toolkit for Monte Carlo to improve the simulation performance using CAD-based models by factors of 1.5 to 2.

42 ENGINEERING↗