Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “bottleneck analysis”

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 397 records · Page 22

Dynamic mode decomposition of nonequilibrium electron-phonon dynamics: accelerating the first-principles real-time Boltzmann equation

Abstract Nonequilibrium dynamics governed by electron–phonon ( e -ph) interactions plays a key role in electronic devices and spectroscopies and is central to understanding electronic excitations in materials. The real-time Boltzmann transport equation (rt-BTE) with collision processes computed from first principles can describe the coupled dynamics of electrons and atomic vibrations (phonons). Yet, a bottleneck of these simulations is the calculation of e –ph scattering integrals on dense momentum grids at each time step. Here we show a data-driven approach based on dynamic mode decomposition (DMD) that can accelerate the time propagation of the rt-BTE and identify dominant electronic processes. We apply this approach to two case studies, high-field charge transport and ultrafast excited electron relaxation. In both cases, simulating only a short time window of ~10% of the dynamics suffices to predict the dynamics from initial excitation to steady state using DMD extrapolation. Analysis of the momentum-space modes extracted from DMD sheds light on the microscopic mechanisms governing electron relaxation to a steady state or equilibrium. The combination of accuracy and efficiency makes our DMD-based method a valuable tool for investigating ultrafast dynamics in a wide range of materials.

36 MATERIALS SCIENCE↗

Quantitative insights for diagnosing performance bottlenecks in lithium–sulfur batteries

Lithium–sulfur (Li–S) batteries hold significant promise for electric vehicles and aviation due to their high energy density and cost-effectiveness. However, understanding the root causes of performance degradation remains a formidable challenge, as the interplay of multiple factors obscures key failure mechanisms. A major limitation has been the inability to quantify soluble sulfur species within practical detection limits accurately and to correlate electrochemical processes with associated physical inventory changes. Here, we introduce the high-performance liquid chromatography-ultraviolet spectroscopy and gas chromatography sequential characterization (HUGS) toolkit, capable of precisely quantifying seven distinct sulfur and polysulfide species at concentrations as low as 40 ppb. HUGS has been successfully applied to practical coin and pouch cells without requiring cell modification. Furthermore, our self-developed software, Dr HUGS, enhanced the data analysis speed by over 30 times, enabling multi-source data integration and delivering comprehensive analysis results within minutes. Using HUGS, we identify significant capacity losses from inactive lithium and sulfur during initial cycles and sulfide-rich solid–electrolyte interphase (SEI) formation on the anode during later cycles. Notably, our findings reveal that soluble polysulfides have minimal contributions to capacity loss, challenging long-standing assumptions. Moreover, HUGS demonstrates that constant-pressure setups in Li–S pouch cells improve compositional uniformity compared to constant-gap configurations. For sulfurized polyacrylonitrile (SPAN) cathodes, unique issues such as non-sulfide SEI formation and lithium pulverization are observed, which can be mitigated through localized high-concentration electrolytes to enhance lithium inventory retention. By enabling precise quantification of critical inventory components, HUGS provides transformative insights into failure mechanisms across various electrolytes and cathode chemistries, guiding rational design strategies for next-generation energy storage systems.

25 ENERGY STORAGE↗

Optimizing tertiary storage organization and access for spatio-temporal datasets

We address in this paper data management techniques for efficiently retrieving requested subsets of large datasets stored on mass storage devices. This problem represents a major bottleneck that can negate the benefits of fast networks, because the time to access a subset from a large dataset stored on a mass storage system is much greater that the time to transmit that subset over a network. This paper focuses on very large spatial and temporal datasets generated by simulation programs in the area of climate modeling, but the techniques developed can be applied to other applications that deal with large multidimensional datasets. The main requirement we have addressed is the efficient access of subsets of information contained within much larger datasets, for the purpose of analysis and interactive visualization. We have developed data partitioning techniques that partition datasets into 'clusters' based on analysis of data access patterns and storage device characteristics. The goal is to minimize the number of clusters read from mass storage systems when subsets are requested. We emphasize in this paper proposed enhancements to current storage server protocols to permit control over physical placement of data on storage devices. We also discuss in some detail the aspects of the interface between the application programs and the mass storage system, as well as a workbench to help scientists to design the best reorganization of a dataset for anticipated access patterns.

Chen, Ling Tony↗

Discovering Ni/Cu Single-Atom Alloy as a Highly Active and Selective Catalyst for Direct Methane Conversion to Ethylene: A First-Principles Kinetic Study

Direct methane conversion to liquid fuels or value-added chemicals is a promising technology to utilize natural resources without resorting to further petroleum extraction. However, discovering efficient catalysts for this reaction is challenging due to either coke formation or unfavorable C-H bond activation. Herein, we design single-atom alloy (SAA) catalysts to simultaneously eliminate the above two bottlenecks based on mechanism-guided strategies: (1) the active single atom enables favorable C-H bond breaking and (2) the less reactive host metal facilitates C-C coupling and thus avoids strong binding of carbonaceous species. Employing electronic structure theory calculations, we screened the stability of multiple SAAs with 3d-5d transition metals atomically dispersed on a copper surface in terms of avoiding dopant aggregation and segregation. We then evaluated reactivities of the stable SAAs as catalysts for direct methane conversion to C2 products, including methane dehydrogenation and C-C coupling mechanisms. Combining selectivity analysis with kinetic modeling, we predicted that nickel dispersed on copper, i.e., Ni/Cu SAA, is a highly active and selective catalyst that can efficiently transform methane to ethylene. This work designs efficient SAA catalysts for direct methane activation and provides chemical insights into engineering compositions of SAAs to tune their catalytic performances.

Kothakonda, Manish↗

Integration of RNTuple in ATLAS Athena

After using ROOT’s TTree I/O subsystem for over two decades and storing more than an exabyte of compressed High Energy Physics (HEP) data, advances in technology have motivated a complete redesign, RNTuple, which breaks backward-compatibility to take better advantage of these storage options. The RNTuple I/O subsystem has been designed to address performance bottlenecks and other shortcomings of TTree. Specifically, RNTuple comes with an updated, more compact binary data format that can be stored both in ROOT files and natively in object stores. It is designed for modern storage hardware (e.g. high-throughput low-latency NVMe SSDs), and provides robust and easy to use interfaces. The binary format of RNTuple is scheduled to become production grade in 2024, and recently has become mature enough to start exploring the integration into software used by HEP experiments. In this contribution, we discuss the developments to support the features as required by the ATLAS analysis Event Data Model (EDM) in RNTuple, which will enable its integration into the Athena software framework. With these developments in place, we evaluate the performance of the current most recent versions of RNTuple-based ATLAS data sets and compare this to that of TTree.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Experimental Systems-Biology Approaches for Clostridia-Based Bioenergy Production

This is the final project report for project "Experimental Systems-Biology Approaches for Clostridia-Based Bioenergy Production" for the funding period of 9/1/12 to 2/28/2015 (three years with a 6-month no-cost extension) OVERVIEW AND PROJECT GOALS The bottleneck of achieving higher rates and titers of toxic metabolites (such as solvents and carboxylic acids that can used as biofuels or biofuel precursors) can be overcome by engineering the stress response system. Thus, understanding and modeling the response of cells to toxic metabolites is a problem of great fundamental and practical significance. In this project, our goal is to dissect at the molecular systems level and build models (conceptual and quantitative) for the stress response of C. acetobutylicum (Cac) to its two toxic metabolites: butanol (BuOH) and butyrate (BA). Transcriptional (RNAseq and microarray based), proteomic and fluxomic data and their analysis are key requirements for this goal. Transcriptional data from mid-exponential cultures of Cac under 4 different levels of BuOH and BA stress was obtained using both microarrays (Papoutsakis group) and deep sequencing (RNAseq; Meyers and Papoutsakis groups). These two sets of data do not only serve to validate each other, but are also used for identification of stress-induced changes in transcript levels, small regulatory RNAs, & in transcriptional start sites. Quantitative proteomic data (Lee group), collected using the iTRAQ technology, are essential for understanding of protein levels and turnover under stress and the various protein-protein interactions that orchestrate the stress response. Metabolic flux changes (Antoniewicz group) of core pathways, which provide important information on the re-allocation of energy and carbon resources under metabolite stress, were examined using 13C-labelled chemicals. Omics data are integrated at different levels and scales. At the metabolic-pathway level, omics data are integrated into a 2nd generation genome-scale model (GSM) (Maranas group). Omics data are also integrated using bioinformatics (Wu and Huang group), whereby regulatory details of gene and protein expression, protein-protein interactions and metabolic flux regulation are incorporated. The PI (Papoutsakis) facilitated project integration through monthly meeting and reports, conference calls, and collaborative manuscript preparation. The five groups collaborated extensively and made a large number of presentations in national and international meetings. It has also published several papers, with several more in the preparation stage. Several PhD, MS and postdoctoral students were trained as part of this collaborative and interdisciplinary project.

09 BIOMASS FUELS↗

Temporal Genetic Dynamics of an Experimental, Biparental Field Population of Phytophthora capsici

Defining the contributions of dispersal, reproductive mode, and mating system to the population structure of a pathogenic organism is essential to estimating its evolutionary potential. After introduction of the devastating plant pathogen, Phytophthora capsici, into a grower’s field, a lack of aerial spore dispersal restricts migration. Once established, coexistence of both mating types results in formation of overwintering recombinant oospores, engendering persistent pathogen populations. To mimic these conditions, in 2008, we inoculated a field with two P. capsici isolates of opposite mating type. We analyzed pathogenic isolates collected in 2009–2013 from this experimental population, using genome-wide single-nucleotide polymorphism markers. By tracking heterozygosity across years, we show that the population underwent a generational shift; transitioning from exclusively F1 in 2009–2010, to multi-generational in 2011, and ultimately all inbred in 2012–2013. Survival of F1 oospores, characterized by heterozygosity excess, coupled with a low rate of selfing, delayed declines in heterozygosity due to inbreeding and attainment of equilibrium genotypic frequencies. Large allele and haplotype frequency changes in specific genomic regions accompanied the generational shift, representing putative signatures of selection. Finally, we identified an approximately 1.6 Mb region associated with mating type determination, constituting the first detailed genomic analysis of a mating type region (MTR) in Phytophthora. Segregation patterns in the MTR exhibited tropes of sex-linkage, where maintenance of allele frequency differences between isolates of opposite mating type was associated with elevated heterozygosity despite inbreeding. Characterizing the trajectory of this experimental system provides key insights into the processes driving persistent, sexual pathogen populations.

bottleneck↗

Development and Validation of a Lifecycle-based Prognostics Architecture with Test Bed Validation

On-line monitoring and tracking of nuclear plant system and component degradation is being investigated as a method for improving the safety, reliability, and maintainability of aging nuclear power plants. Accurate prediction of the current degradation state of system components and structures is important for accurate estimates of their remaining useful life (RUL). The correct quantification and propagation of both the measurement uncertainty and model uncertainty is necessary for quantifying the uncertainty of the RUL prediction. This research project developed and validated methods to perform RUL estimation throughout the lifecycle of plant components. Prognostic methods should seamlessly operate from beginning of component life (BOL) to end of component life (EOL). We term this "Lifecycle Prognostics." When a component is put into use, the only information available may be past failure times of similar components used in similar conditions, and the predicted failure distribution can be estimated with reliability methods such as Weibull Analysis (Type I Prognostics). As the component operates, it begins to degrade and consume its available life. This life consumption may be a function of system stresses, and the failure distribution should be updated to account for the system operational stress levels (Type II Prognostics). When degradation becomes apparent, this information can be used to again improve the RUL estimate (Type III Prognostics). This research focused on developing prognostics algorithms for the three types of prognostics, developing uncertainty quantification methods for each of the algorithms, and, most importantly, developing a framework using Bayesian methods to transition between prognostic model types and update failure distribution estimates as new information becomes available. The developed methods were then validated on a range of accelerated degradation test beds. The ultimate goal of prognostics is to provide an accurate assessment for RUL predictions, with as little uncertainty as possible. From a reliability and maintenance standpoint, there would be improved safety by avoiding all failures. Calculated risk would decrease, saving money by avoiding unnecessary maintenance. One major bottleneck for data-driven prognostics is the availability of run-to-failure degradation data. Without enough degradation data leading to failure, prognostic models can yield RUL distributions with large uncertainty or mathematically unsound predictions. To address these issues a "Lifecycle Prognostics" method was developed to create RUL distributions from Beginning of Life (BOL) to End of Life (EOL). This employs established Type I, II, and III prognostic methods, and Bayesian transitioning between each Type. Bayesian methods, as opposed to classical frequency statistics, show how an expected value, a priori, changes with new data to form a posterior distribution. For example, when you purchase a component you have a prior belief, or estimation, of how long it will operate before failing. As you operate it, you may collect information related to its condition that will allow you to update your estimated failure time. Bayesian methods are best used when limited data are available. The use of a prior also means that information is conserved when new data are available. The weightings of the prior belief and information contained in the sampled data are dependent on the variance (uncertainty) of the prior, the variance (uncertainty) of the data, and the amount of measured data (number of samples). If the variance of the prior is small compared to the uncertainty of the data, the prior will be weighed more heavily. However, as more data are collected, the data will be weighted more heavily and will eventually swamp out the prior in calculating the posterior distribution of model parameters. Fundamentally Bayesian analysis updates a prior belief with new data to get a posterior belief. The general approach to applying the Bayesian method to lifecycle prognostics consisted of identifying the prior, which is the RUL estimate and uncertainty from the previous prognostics type, and combining it with observational data related to the newer prognostics type. The resulting lifecycle prognostics algorithm uses all available information throughout the component lifecycle.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Origins of high ductility exhibited by an extruded magnesium alloy Mg-1.8Zn-0.2Ca: Experiments and crystal plasticity modeling

Low ductility and strength are major bottlenecks against Mg alloys' wide applications. In this work, we systematically design the composition and fabrication process for a low-alloyed Mg-Zn-Ca alloy, showing that it can be extruded at low temperatures (~ 250 degrees C) and high speeds (~2 mm/s). After the extrusion, this alloy exhibits a substantially weakened basal texture, relatively small grain size, very high tensile elongation (~ 30%), and good strength. The origin of the considerably improved ductility was studied using a combination of three-dimensional atom probe tomography (3D-APT), transmission electron microscopy (TEM), electron backscattered diffraction (EBSD) in conjunction with surface slip trace analysis, in-situ synchrotron X-ray diffraction, and elasto-plastic self-consistent (EPSC) modeling. Co-segregation of Zn and Ca atoms at a grain boundary is observed and associated with texture weakening and grain boundary mediated plasticity, both improving the ductility. While basal slip and prismatic slip are identified as the dominant deformation systems in the alloy, the ratio between their slip resistances is substantially reduced relative to pure Mg and most other Mg alloys, significantly contributing to the improved ductility of the alloy. This Mg-Zn-Ca alloy exhibiting excellent mechanical properties and low fabrication cost is a promising candidate for industrial productions.

36 MATERIALS SCIENCE↗

Modelling galaxy cluster triaxiality in stacked cluster weak lensing analyses

Counts of galaxy clusters offer a high-precision probe of cosmology, but control of systematic errors will determine the accuracy of this measurement. Using Buzzard simulations, we quantify one such systematic, the triaxiality distribution of clusters identified with the redMaPPer optical cluster finding algorithm, which was used in the Dark Energy Survey Year-1 (DES Y1) cluster cosmology analysis. We test whether redMaPPer selection biases the clusters’ shape and orientation and find that it only biases orientation, preferentially selecting clusters with their major axes oriented along the line of sight. Modelling the richness–mass relation as log-linear, we find that the log-richness amplitude ln (A) is boosted from the lowest to highest orientation bin with a significance of 14σ, while the orientation dependence of the richness-mass slope and intrinsic scatter is minimal. We also find that the weak lensing shear-profile ratios of cluster-associated dark haloes in different orientation bins resemble a ‘bottleneck’ shape that can be quantified with a Cauchy function. We test the correlation of orientation with two other leading systematics in cluster cosmology – miscentering and projection – and find a null correlation. The resulting mass bias predicted from our templates confirms the DES Y1 finding that triaxiality is a leading source of bias in cluster cosmology. However, the richness-dependence of the bias confirms that triaxiality does not fully resolve the tension at low-richness between DES Y1 cluster cosmology and other probes. Our model can be used for quantifying the impact of triaxiality bias on cosmological constraints for upcoming weak lensing surveys of galaxy clusters.

79 ASTRONOMY AND ASTROPHYSICS↗

Comparative Performance Analysis of Coarse Solvers for Algebraic Multigrid on Multicore and Manycore Architectures

In this paper, we study the performance of a two-level algebraic-multigrid algorithm, with a focus on the impact of the coarse-grid solver on performance. We consider two algorithms for solving the coarse-space systems: the preconditioned conjugate gradient method and a new robust HSS-embedded low-rank sparse-factorization algorithm. Our test data comes from the SPE Comparative Solution Project for oil-reservoir simulations. We contrast the performance of our code on one 12-core socket of a Cray XC30 machine with performance on a 60-core Intel Xeon Phi coprocessor. To obtain top performance, we optimized the code to take full advantage of fine-grained parallelism and made it thread-friendly for high thread count. We also developed a bounds-and-bottlenecks performance model of the solver which we used to guide us through the optimization effort, and also carried out performance tuning in the solver’s large parameter space. Finally, as a result, significant speedups were obtained on both machines.

97 MATHEMATICS AND COMPUTING↗

Accurate and complete genomes from metagenomes

Genomes are an integral component of the biological information about an organism; thus, the more complete the genome, the more informative it is. Historically, bacterial and archaeal genomes were reconstructed from pure (monoclonal) cultures, and the first reported sequences were manually curated to completion. However, the bottleneck imposed by the requirement for isolates precluded genomic insights for the vast majority of microbial life. Shotgun sequencing of microbial communities, referred to initially as community genomics and subsequently as genome-resolved metagenomics, can circumvent this limitation by obtaining metagenome-assembled genomes (MAGs); but gaps, local assembly errors, chimeras, and contamination by fragments from other genomes limit the value of these genomes. Here, we discuss genome curation to improve and, in some cases, achieve complete (circularized, no gaps) MAGs (CMAGs). To date, few CMAGs have been generated, although notably some are from very complex systems such as soil and sediment. Through analysis of about 7000 published complete bacterial isolate genomes, we verify the value of cumulative GC skew in combination with other metrics to establish bacterial genome sequence accuracy. The analysis of cumulative GC skew identified potential misassemblies in some reference genomes of isolated bacteria and the repeat sequences that likely gave rise to them. We discuss methods that could be implemented in bioinformatic approaches for curation to ensure that metabolic and evolutionary analyses can be based on very high-quality genomes.

59 BASIC BIOLOGICAL SCIENCES↗

Multilevel Analysis of Electrochemically Mediated Methanolysis of Poly(ethylene terephthalate) (PET)

Here, this study presents a multilevel analysis of electrochemically mediated methanolysis as a promising method for reducing the environmental impacts of plastic recycling, with a focus on depolymerizing poly(ethylene terephthalate) (PET) into dimethyl terephthalate (DMT). Instead of conventional chemical PET depolymerization, this electrochemical approach provides distinct technical advantages in process control and efficiency. At the process level, key operational parameters, including applied current and reaction time, were systematically investigated to optimize PET conversion and DMT selectivity. The electrochemical approach was directly compared to equivalent chemical methanolysis systems and demonstrated superior performance in terms of PET conversion and DMT selectivity. Building on these findings, a technoeconomic assessment identified the current economic bottlenecks and revealed that improvements in process design, DMT selectivity, PET conversion, and energy efficiency are key to reducing the overall process cost and enabling future implementation. While further optimization is required for market competitiveness, these results establish a performance baseline for the electrochemically mediated PET methanolysis process and underscore the importance of combining process-level innovation with systems-level evaluation in the development of sustainable recycling technologies.

chemical recycling↗

Hierarchical Data Format for Nuclear Data Sensitivities

The SCALE code system includes capabilities for sensitivity and uncertainty (S/U) analysis as part of its TSUNAMI code suite. The sensitivity of a quantity of interest (for example, an application’s $k_{eff}$) to nuclear data is stored as a profile in a text-based file, which is known as a sensitivity data file (SDF). The sensitivity profile can be used to calculate uncertainties, correlation coefficients, and similarity indices. One of the goals of the present work was to seek general performance improvements in the TSUNAMI code suite, starting with the TSUNAMI-IP code for calculating similarity indices. Through profiling, it was found that reading the text-based sensitivity files was a performance bottleneck in the TSUNAMI-IP code. In a typical TSUNAMI-IP calculation, an application might be compared to thousands of benchmarks, thus requiring the reading of thousands of SDFs. Reading of binary-based data is generally faster than reading text-based data. Hierarchical Data Format 5 (HDF5) is a binary-based format that also benefits from being portable, and it can be inspected with nonproprietary tools. This paper describes an HDF5-based file format that has been introduced for SDFs.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

X: A Comprehensive Analytic Model for Parallel Machines

To continuously comply with Moore’s Law, modern parallel machines become increasingly complex. Effectively tuning application performance for these machines therefore becomes a daunting task. Moreover, identifying performance bottlenecks at application and architecture level, as well as evaluating various optimization strategies, are becoming extremely difficult when the entanglement of numerous correlated factors is being presented. To tackle these challenges, we present a visual analytical model named “X”. It is intuitive and sufficiently flexible to track all the typical features of a parallel machine.

Modeling and Performance Analysis↗

Techno-Economic Assessment of Data Center Load Demand Powered by Small Modular Reactors and Distributed Energy Resources

The rapid increase in data center energy demand, driven by AI and large-scale data processing, poses significant challenges to global energy infrastructure. Data centers require substantial and reliable energy for continuous operations and high-performance computing. Current electrical grids face issues such as transmission bottlenecks and aging infrastructure, making it difficult to meet these demands. Integrating inverter-based-resources (IBRs) like solar and wind presents both opportunities and challenges due to their intermittent nature. Small Modular Reactors (SMRs) offer a promising solution with their enhanced safety, modularity, reliability, and scalability, providing consistent base load power ideal for data center operations. This study presents a comprehensive techno-economic assessment of powering data center load demand using a combination of SMRs and IBRs with grid-connected and islanded mode. This study utilized Idaho National Laboratory’s (INL) HPC data center hourly load profiles and Xendee microgrid optimization platform to conduct the analysis. In this configuration, SMRs serves as the primary base load power source, consistently providing a steady supply of electricity necessary to meet the minimum load demand of the data center with support from the IBRs. Key performance indicators such as Levelized Cost of Electricity (LCOE), Net Present Value (NPV) has been calculated to assess the economic feasibility. The findings from this research will underscore the strategic benefits of integrating SMR plant with DERs – particularly for critical infrastructure load such as data centers.

14 - SOLAR ENERGY↗

A GPU-based compressible combustion solver for applications exhibiting disparate space and time scales

High-speed chemically active flows pose significant computational challenges due to their disparate space and time scales, with stiff chemistry often dominating simulation time. While modern scientific computing programs achieve exascale performance by leveraging graphics processing units (GPUs), existing GPU-based compressible combustion solvers face critical limitations in memory management, load balancing, and handling the highly localized nature of chemical reactions. To this end, we present a high-performance compressible reacting flow solver built on the AMReX framework and optimized for multi-GPU settings. Here, our approach addresses three GPU performance bottlenecks: memory access patterns through column-major storage optimization, computational workload variability via a bulk-sparse integration strategy for chemical kinetics, and multi-GPU load distribution for adaptive mesh refinement applications. The solver adapts existing matrix-based chemical kinetics formulations to multi-grid contexts. Using representative combustion applications, including 2D and 3D detonations and a 3D jet-in-crossflow configuration, we demonstrate 1.4–5× performance improvements over initial implementations on an in-house cluster of NVIDIA H100 GPUs, and near-ideal weak scaling on the Frontier supercomputer (Oak Ridge Leadership Computing Facility) with up to 1024 AMD Instinct MI250X GPUs. Roofline analysis reveals substantial improvements in arithmetic intensity for both convection (∼ 10 ×) and chemistry (∼ 4 ×) routines, confirming efficient utilization of GPU memory bandwidth and computational resources.

42 ENGINEERING↗

Fast correlation function calculator: A high-performance pair-counting toolkit

A novel high-performance exact pair-counting toolkit called fast correlation function calculator (FCFC) is presented. With the rapid growth of modern cosmological datasets, the evaluation of correlation functions with observational and simulation catalogues has become a challenge. High-efficiency pair-counting codes are thus in great demand. We introduce different data structures and algorithms that can be used for pair-counting problems, and perform comprehensive benchmarks to identify the most efficient algorithms for real-world cosmological applications. We then describe the three levels of parallelisms used by FCFC, SIMD, OpenMP, and MPI, and run extensive tests to investigate the scalabilities. Finally, we compare the efficiency of FCFC with alternative pair-counting codes. The data structures and histogram update algorithms implemented in FCFC are shown to outperform alternative methods. FCFC does not benefit greatly from SIMD because the bottleneck of our histogram update algorithm is mainly cache latency. Nevertheless, the efficiency of FCFC scales well with the numbers of OpenMP threads and MPI processes, even though speedups may be degraded with over a few thousand threads in total. FCFC is found to be faster than most (if not all) other public pair-counting codes for modern cosmological pair-counting applications.

79 ASTRONOMY AND ASTROPHYSICS↗