Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Precise”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Bit Grooming: statistically accurate precision-preserving quantization with compression, evaluated in the netCDF Operators (NCO, v4.4.8+)

Abstract. Geoscientific models and measurements generate false precision (scientifically meaningless data bits) that wastes storage space. False precision can mislead (by implying noise is signal) and be scientifically pointless, especially for measurements. By contrast, lossy compression can be both economical (save space) and heuristic (clarify data limitations) without compromising the scientific integrity of data. Data quantization can thus be appropriate regardless of whether space limitations are a concern. We introduce, implement, and characterize a new lossy compression scheme suitable for IEEE floating-point data. Our new Bit Grooming algorithm alternately shaves (to zero) and sets (to one) the least significant bits of consecutive values to preserve a desired precision. This is a symmetric, two-sided variant of an algorithm sometimes called Bit Shaving that quantizes values solely by zeroing bits. Our variation eliminates the artificial low bias produced by always zeroing bits, and makes Bit Grooming more suitable for arrays and multi-dimensional fields whose mean statistics are important. Bit Grooming relies on standard lossless compression to achieve the actual reduction in storage space, so we tested Bit Grooming by applying the DEFLATE compression algorithm to bit-groomed and full-precision climate data stored in netCDF3, netCDF4, HDF4, and HDF5 formats. Bit Grooming reduces the storage space required by initially uncompressed and compressed climate data by 25–80 and 5–65 %, respectively, for single-precision values (the most common case for climate data) quantized to retain 1–5 decimal digits of precision. The potential reduction is greater for double-precision datasets. When used aggressively (i.e., preserving only 1–2 digits), Bit Grooming produces storage reductions comparable to other quantization techniques such as Linear Packing. Unlike Linear Packing, whose guaranteed precision rapidly degrades within the relatively narrow dynamic range of values that it can compress, Bit Grooming guarantees the specified precision throughout the full floating-point range. Data quantization by Bit Grooming is irreversible (i.e., lossy) yet transparent, meaning that no extra processing is required by data users/readers. Hence Bit Grooming can easily reduce data storage volume without sacrificing scientific precision or imposing extra burdens on users.

58 GEOSCIENCES↗
DOE OSTI2026

Mixed-precision numerics in scientific applications: survey and perspectives

The explosive demand for artificial intelligence (AI) workloads has led to a significant increase in silicon area dedicated to lower-precision computations on recent high-performance computing hardware designs. However, mixed-precision capabilities, which can achieve performance improvements of up to 8x compared to double-precision in extreme compute-intensive workloads, remain largely untapped in most scientific applications. A growing number of efforts have shown that mixed-precision algorithmic innovations can deliver superior performance without sacrificing accuracy. These developments should prompt computational scientists to seriously consider whether their scientific modeling and simulation applications could benefit from the acceleration offered by new hardware and mixed-precision algorithms. In this survey, we (1) review progress across diverse scientific domains—fluid dynamics, weather and climate, quantum chemistry, and computational genomics—that have begun adopting mixed-precision strategies; (2) examine state-of-the-art algorithmic techniques such as iterative refinement, splitting and emulation schemes, and adaptive precision solvers; (3) assess their implications for accuracy, performance, and resource utilization; and (4) survey the emerging software ecosystem that enables mixed-precision methods at scale. We conclude with perspectives and recommendations on cross-cutting opportunities, domain-specific challenges, and the role of co-design between application scientists, numerical analysts, and computer scientists. Collectively, this survey underscores that mixed-precision numerics can reshape computational science by aligning algorithms with the evolving landscape of hardware capabilities.

Graphics processing units↗
DOE OSTI2017

The Surface Brightness-color Relations Based on Eclipsing Binary Stars: Toward Precision Better than 1% in Angular Diameter Predictions

In this study we investigate the calibration of surface brightness–color (SBC) relations based solely on eclipsing binary stars. We selected a sample of 35 detached eclipsing binaries with trigonometric parallaxes from Gaia DR1 or Hipparcos whose absolute dimensions are known with an accuracy better than 3% and that lie within 0.3 kpc from the Sun. For the purpose of this study, we used mostly homogeneous optical and near-infrared photometry based on the Tycho-2 and 2MASS catalogs. We derived geometric angular diameters for all stars in our sample with a precision better than 10%, and for 11 of them with a precision better than 2%. The precision of individual angular diameters of the eclipsing binary components is currently limited by the precision of the geometric distances (∼5% on average). However, by using a subsample of systems with the best agreement between their geometric and photometric distances, we derived the precise SBC relations based only on eclipsing binary stars. These relations have precisions that are comparable to the best available SBC relations based on interferometric angular diameters, and they are fully consistent with them. With very precise Gaia parallaxes becoming available in the near future, angular diameters with a precision better than 1% will be abundant. At that point, the main uncertainty in the total error budget of the SBC relations will come from transformations between different photometric systems, disentangling of component magnitudes, and for hot OB stars, the main uncertainty will come from the interstellar extinction determination. We argue that all these issues can be overcome with modern high-quality data and conclude that a precision better than 1% is entirely feasible.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Statistical precision of the intensities retrieved from constrained fitting of overlapping peaks in high-resolution mass spectra

Least-squares fitting of overlapping peaks is often needed to separately quantify ions in high-resolution mass spectrometer data. A statistical simulation approach is used to assess the statistical precision of the retrieved peak intensities. The sensitivity of the fitted peak intensities to statistical noise due to ion counting is probed for synthetic data systems consisting of two overlapping ion peaks whose positions are pre-defined and fixed in the fitting procedure. The fitted intensities are sensitive to imperfections in the m/Q calibration. These propagate as a limiting precision in the fitted intensities that may greatly exceed the precision arising from counting statistics. The precision on the fitted peak intensity falls into one of three regimes. In the "counting-limited regime" (regime I), above a peak separation &chi; ~ 2 to 3 half-widths at half-maximum (HWHM), the intensity precision is similar to that due to counting error for an isolated ion. For smaller &chi; and higher ion counts (~ 1000 and higher), the intensity precision rapidly degrades as the peak separation is reduced ("calibration-limited regime", regime II). Alternatively for &chi; < 1.6 but lower ion counts (e.g. 10–100) the intensity precision is dominated by the additional ion count noise from the overlapping ion and is not affected by the imprecision in the m/Q calibration ("overlapping-limited regime", regime III). The transition between the counting and m/Q calibration-limited regimes is shown to be weakly dependent on resolving power and data spacing and can thus be approximated by a simple parameterisation based only on peak intensity ratios and separation. A simple equation can be used to find potentially problematic ion pairs when evaluating results from fitted spectra containing many ions. Longer integration times can improve the precision in regimes I and III, but a given ion pair can only be moved out of regime II through increased spectrometer resolving power. As a result, studies presenting data obtained from least-squares fitting procedures applied to mass spectral peaks should explicitly consider these limits on statistical precision.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Differentiable Neural Architecture, Mixed Precision and Accelerator Co-Search

Quantization, effective Neural Network architecture, and efficient accelerator hardware are three important design paradigms to maximize accuracy and efficiency. Mixed Precision Quantization is a process of assigning different precision to different Neural Network layers for optimized inference. Neural Architecture Search (NAS) is a process of automatically designing the neural network for a task and can also be extended to search for the precision of each weight and activation matrix. In this paper, we develop the following three methods: (i) Fast Differentiable Hardware-aware Mixed Precision Quantization Search method to find optimal precision, (ii) Joint Differentiable hardware-aware Architecture and Mixed Precision Quantization Co-search, (iii) Joint Accelerator, Architecture, and Precision triple co-search to find best possibilities in all the three worlds. We demonstrate the effectiveness of our proposed methods targeting Bitfusion accelerator by searching mixed precision models on MobilenetV2. We achieve better accuracy-latency trade-off models than the manually designed and previously proposed search methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Accuracy and Precision of Age Estimates Obtained from Three Calcified Structures for Known-Age Kokanee

Abstract Kokanee Oncorhynchus nerka growth is often density dependent; thus, proper management of kokanee populations necessitates an understanding of population dynamics using age structure data. To date, no calcified structures have been validated for kokanee. We compared the accuracy (i.e., the percentage of reconciled age estimates that matched the known ages of fish) and precision (i.e., the percentage of fish for which complete agreement was achieved on age estimates among all readers) of aging estimates for scales, sectioned otoliths, and sectioned pectoral fin rays from 455 known-age kokanee (ages 0–4) collected from five lentic waters in Idaho. Across all waters combined, mean weighted accuracy and precision were similar for scales (86% and 70%, respectively), fin rays (83% and 65%), and otoliths (82% and 65%), with no significant differences between structures. However, among water bodies, accuracy and precision of each calcified structure varied considerably. For example, scales were the most accurate and precise structure and otoliths were the least accurate and precise structure at one water body, while otoliths were the most accurate and precise structure and scales were the least accurate structure at two water bodies. Fin rays were the least precise structure at four of the five study waters, but were the least accurate structure for only one water body. Individual reader accuracy was most affected by fish age and water body, and older fish (age 3 and older) were consistently assigned incorrect ages regardless of the water body or the calcified structure. Taken collectively, all three structures produced satisfactory aging accuracy and precision for kokanee, but no structure was unequivocally best; at any individual water body, annual growth and local environmental conditions appeared to influence the readability of calcified structures.

Branigan, Philip R. (ORCID:0000000337990573)↗
DOE OSTI2013

Obtaining identical results with double precision global accuracy on different numbers of processors in parallel particle Monte Carlo simulations

We describe and compare different approaches for achieving numerical reproducibility in photon Monte Carlo simulations. Reproducibility is desirable for code verification, testing, and debugging. Parallelism creates a unique problem for achieving reproducibility in Monte Carlo simulations because it changes the order in which values are summed. This is a numerical problem because double precision arithmetic is not associative. Parallel Monte Carlo, both domain replicated and decomposed simulations, will run their particles in a different order during different runs of the same simulation because the non-reproducibility of communication between processors. In addition, runs of the same simulation using different domain decompositions will also result in particles being simulated in a different order. In [1], a way of eliminating non-associative accumulations using integer tallies was described. This approach successfully achieves reproducibility at the cost of lost accuracy by rounding double precision numbers to fewer significant digits. This integer approach, and other extended and reduced precision reproducibility techniques, are described and compared in this work. Increased precision alone is not enough to ensure reproducibility of photon Monte Carlo simulations. Non-arbitrary precision approaches require a varying degree of rounding to achieve reproducibility. For the problems investigated in this work double precision global accuracy was achievable by using 100 bits of precision or greater on all unordered sums which where subsequently rounded to double precision at the end of every time-step.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Variable Precision Computing

This report summarizes the activities and major accomplishments of the Variable Precision Computing Strategic Initiative project. The overarching goal of this project was to initiate and promote a new paradigm in High Performance Computing (HPC) that would fundamentally change how we represent and make use of representations of real numbers in finite precision and that would establish Lawrence Livermore National Laboratory (LLNL) as a leader in precision-related research for HPC. We pursued three integrated and concurrent research thrusts: new, more dynamic data representations to address data motion and capacity limitations; improved mixed precision algorithms to accelerate computation and to leverage new hardware; and new tools to help developers reason about precision and to automate code transformations. Within these thrust areas, we made significant advances in the use of compressed array data types in numerical calculations, in the theoretical justification, in improved performance, and in the prospects of hardware implementation; in the development of metrics to significantly reduce the storage needs of molecular dynamics simulations; in the development and analysis of efficient mixed-precision algorithms for time-dependent partial differential equations and for orthogonal factorization; and in the capability of tools that can analyze sensitivity of computed results to choices in precision and that can make automated code transformations based on such analysis. We established that there are definite opportunities to realize significant improvements in HPC above and beyond the 2x limits of traditional mixed precision using the technologies we have developed. The project produced numerous papers, presentations, and posters in important venues as well as open-source software contributions to the community, and many of the new capabilities have initiated or are being leveraged in ongoing projects.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Preliminary Theoretical Analysis of Mixed Precision Krylov Subspace Methods (Q3 Report)

The third quarter of the project was spent performing theoretical finite precision analysis of Krylov subspace method variants that use mixed precision. Our focus here is on the Conjugate Gradient (CG) method and the Lanczos method. We have performed an analysis of maximum attainable accuracy for the classical CG method in which 3 precisions are used: a working precision ε, a precision ε IP for the inner product computations, and a precision ε MV for the matrix-vector products. Our results show that performing inner product computations in lower precision does not affect the attainable accuracy. Further, we have performed a complete error analysis of the s-step Lanczos algorithm. In this case, we show that the numerical behavior of the algorithm can be significantly improved by using extra precision in a small part of the computation. We summarize the main theorems in the remainder of the document. Other activities include attending biweekly xSDK meetings. The subsequent quarter will be spent finalizing these results into technical reports and/or manuscripts for submission to journals, as well as identifying opportunities for future work.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

New Capability at ORNL: High Precision Uranium Titration

Destructive analytical measurements establish a nuclear facility’s nuclear material inventory and inventory differences for Nuclear Material Accountancy and Control. A nuclear laboratory’s ability to perform accurate high-precision analytical measurements is key for tracking large inventories within a facility’s material balance areas and during production to track material movement through dynamic processes. For uranium, these measurements are made by using several established high-precision measurement protocols. These include isotope dilution mass spectrometry, gravimetry, and potentiometric titrations. Whatever measurement technique a nuclear lab chooses to use, reference materials (RMs) with certified attributes and accompanying uncertainties are used to calibrate measurement systems, and are the cornerstone for accurate results. In addition to calibration, RMs provide for metrological traceability, are used for method development and validation, and thus provide critical evaluations of the appropriateness and performance of analytical processes used. Evaluations may include validation that a method is fit-for-purpose, quantification of systematic and random biases, and the evaluation of long-term and short-term performance metrics. High-precision measurement techniques require that measured attributes be certified to a high degree of precision in the RMs used—ultimately, to a higher degree than that of the measurement technique itself. The US authority on the production of special nuclear material Certified Reference Materials is the NBL Program Office (NBLPO), formally known as New Brunswick Laboratory (NBL). The NBLPO is responsible for the sales and distribution of existing NBL certified reference materials (CRMs) and for the production of the next generation of nuclear RMs. To accomplish its mission, NBLPO is establishing key base capabilities within the DOE laboratory complex that formerly existed at the NBL laboratory. The Nuclear Analytical Chemistry (NAC) section within the Chemical Sciences Division (CSD) at Oak Ridge National Laboratory (ORNL) is currently working with NBLPO to set up laboratory and measurement capabilities to provide measurements and capabilities for production and/or recertification of existing and future CRMs for uranium assay. The NBL-developed high-precision titration (HPT) method is a critically-evaluated, extremely precise and accurate primary method utilized for the determination of uranium content in a variety of uranium materials. The HPT method, combined with detailed balance weighing protocols, provides for an analytical methodology that is unsurpassed in precision and one in which all sources of error have been evaluated, a requirement of CRM certification. HPT capability within the United States was lost with the closure of the labs at NBL. The NAC has been collaborating with NBLPO to stand up and demonstrate the capability to perform uranium assay via HPT at ORNL. HPT can produce results with an expanded uncertainty of approximately 0.01% for pure uranium compounds, with typical precisions of <0.006%. The major tasks required to stand up the method at ORNL were the refurbishment of a dedicated lab and equipment setup, procedure development, analyst training, establishing method-specific quality assurance, and qualification of an analyst to perform the method. This report summarizes these tasks, outlines the documents drafted, and gives the outcome of the qualification titrations.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2022

Dynamic Differential Image Circle Diameter Measurement Precision Assessment: Application to Burning Droplets

Dynamic measurement precision assessment has been achieved for a differential circle measurement application. Differential circle diameter measurement, in image analysis, typically requires fitting a circle model that optimizes for image distortions, defects or occlusions. The differential task occurs when precise measurements of diameter change are required given object size variation with time. An automated system was designed to provide diameter measurements and associated measurement precision of images of a fuel droplet undergoing combustion in zero gravity for the FLEX-2 dataset. An image gradient-based, least-squares boundary point fitting method to a circle or ellipse model is used for diameter measurement. The presence of soot aggregates poses significant challenges for diameter measurements when it occludes part of the droplet boundary. The precision of the diameter measurements depends upon the image quality. Using synthetic image simulations that model the soot behavior, we developed a model based on image quality measures that assesses the measurement precision for each individual diameter measurement. Thus, diameter measurements with precision assessments were made available for follow-up scientific analysis. As a result, the algorithm's success rate for measurable runs was 98%. In cases of limited occlusion, a measurement precision of ±0.2 pixels for the FLEX-2 dataset was achieved.

42 ENGINEERING↗
DOE OSTI2024

The Effect of the Prior and the Experimental Design on the Inference of the Precision Matrix in Gaussian Chain Graph Models

Here, we investigate whether (and how) experimental design could aid in the estimation of the precision matrix in a Gaussian chain graph model, especially the interplay between the design, the effect of the experiment and prior knowledge about the effect. Estimation of the precision matrix is a fundamental task to infer biological graphical structures like microbial networks. We compare the marginal posterior precision of the precision matrix under four priors: flat, conjugate Normal-Wishart, Normal-MGIG and a general independent. Under the flat and conjugate priors, the Laplace-approximated posterior precision is not a function of the design matrix rendering useless any efforts to find an optimal experimental design to infer the precision matrix. In contrast, the Normal-MGIG and general independent priors do allow for the search of optimal experimental designs, yet there is a sharp upper bound on the information that can be extracted from a given experiment. We confirm our theoretical findings via a simulation study comparing (i) the KL divergence between prior and posterior and (ii) the Stein’s loss difference of MAPs between random and no experiment. Our findings provide practical advice for domain scientists conducting experiments to better infer the precision matrix as a representation of a biological network.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Accuracy and precision of ESI-Orbitrap-IRMS observations of hours to tens of hours via reservoir injection

Here, orbitrap isotope ratio mass spectrometry (Orbitrap-IRMS) has recently been applied to high-precision, natural-abundance isotope ratio measurements of a diverse range of compounds, including amino acids, oxyanions, fatty acids, and metals. These measurements can characterize many isotope ratios simultaneously at high (≈1.0‰) precision. In a successful experiment, observed precision will track the shot-noise limit and be limited by experimental time. Some isotope ratios, for example those involving 17 O in organic compounds or multiply-substituted (‘clumped’) isotopologues, require experimental times of hours to tens of hours to achieve desired precision, while current sample introduction techniques focus on observations on the order of seconds to tens of minutes. In this study, we characterize Orbitrap-IRMS performance for three long duration measurements (individual acquisitions ≥1 h and as long as 24 h) using an automated reservoir injection system coupled to a Q Exactive HF Orbitrap with an electrospray ionization (ESI) source. First, we characterize long-term intra-measurement stability through a 24-h long measurement of acetone. We report the following isotope ratios and precisions (as acquisition errors, errors on the observed ratio within this measurement (σ AE )): 13 C/ 12 C (σ AE = 0.07‰), 17 O/ 16 O (σ AE = 1.1‰), 18 O/ 16 O (σ AE = 0.3‰), and 13 C 13 C/ 12 C (σ AE = 0.65‰). The σ AE of each tracks the shot noise limit throughout and is limited by the challenging conditions (high resolution and low numbers of ions per scan) required for 17 O/ 16 O measurement in the presence of 13 C via Orbitrap. Second, we characterize inter-measurement stability via a sequence of seven 75-min analyses of perchlorate. We observe the following ratios and acquisition errors: 37 Cl/ 35 Cl (σ AE = 0.09‰); 17 O/ 16 O (σ AE = 1.6‰); 18 O/ 16 O (σ AE = 0.7‰), 37 Cl 17 O/ 35 Cl 16 O (σ AE = 2.7‰), and 37 Cl 18 O/ 35 Cl 16 O (σ AE = 1.2‰). However, we find that inter-measurement drift between acquisitions limits our accuracy and precision for standardized measurements (i.e., error on reported δ values) to ≈1‰ for the 37 Cl/ 35 Cl measurement. Hence, the benefits of low σ AE may not be fully realized. Third, we demonstrate accuracy via sample/standard comparisons of a methionine sample with 13 C enrichment of ≈20‰ relative to a known standard. Using a sequence of seven 60-min analyses, we recover the following isotope ratios and standardized precisions (i.e., error on reported δ values, denoted propagated acquisition errors, σ PAE ): 33 S/ 32 S (σ PAE = 1.0‰), 34 S/ 32 S (σ PAE = 0.7‰), 15 N/ 14 N (σ PAE = 2.1‰), 2 H/ 1 H (σ PAE = 3.2‰), 13 C/ 12 C (σ PAE = 0.4‰), 18 O/ 16 O (σ PAE = 1.6‰), & 13 C 13 C/ 12 C (σ PAE = 2.8‰) with confirmation of accurate results for the known 13 C/ 12 C and 13 C 13 C/ 12 C enrichments. Together, our results demonstrate the viability of Orbitrap-IRMS for long duration measurements of diverse sample types via an automated reservoir injection system. Inter-measurement stability remains a challenge; we expect our methods to be most applicable to extended measurements of hard-to-observe properties, such as 17 O in organics and clumped isotopologues.

58 GEOSCIENCES↗
DOE OSTI2014

Deformable Image Registration for Adaptive Radiation Therapy of Head and Neck Cancer: Accuracy and Precision in the Presence of Tumor Changes

Purpose: To compare deformable image registration (DIR) accuracy and precision for normal and tumor tissues in head and neck cancer patients during the course of radiation therapy (RT). Methods and Materials: Thirteen patients with oropharyngeal tumors, who underwent submucosal implantation of small gold markers (average 6, range 4-10) around the tumor and were treated with RT were retrospectively selected. Two observers identified 15 anatomical features (landmarks) representative of normal tissues in the planning computed tomography (pCT) scan and in weekly cone beam CTs (CBCTs). Gold markers were digitally removed after semiautomatic identification in pCTs and CBCTs. Subsequently, landmarks and gold markers on pCT were propagated to CBCTs, using a b-spline-based DIR and, for comparison, rigid registration (RR). To account for observer variability, the pair-wise difference analysis of variance method was applied. DIR accuracy (systematic error) and precision (random error) for landmarks and gold markers were quantified. Time trend of the precisions for RR and DIR over the weekly CBCTs were evaluated. Results: DIR accuracies were submillimeter and similar for normal and tumor tissue. DIR precision (1 SD) on the other hand was significantly different (P<.01), with 2.2 mm vector length in normal tissue versus 3.3 mm in tumor tissue. No significant time trend in DIR precision was found for normal tissue, whereas in tumor, DIR precision was significantly (P<.009) degraded during the course of treatment by 0.21 mm/week. Conclusions: DIR for tumor registration proved to be less precise than that for normal tissues due to limited contrast and complex non-elastic tumor response. Caution should therefore be exercised when applying DIR for tumor changes in adaptive procedures.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2012

An efficient mixed-precision, hybrid CPU-GPU implementation of a nonlinearly implicit one-dimensional particle-in-cell algorithm

Recently, a fully implicit, energy- and charge-conserving particle-in-cell method has been developed for multi-scale, full-f kinetic simulations [G. Chen, et al., J. Comput. Phys. 230, 18 (2011)]. The method employs a Jacobian-free Newton-Krylov (JFNK) solver and is capable of using very large timesteps without loss of numerical stability or accuracy. A fundamental feature of the method is the segregation of particle orbit integrations from the field solver, while remaining fully self-consistent. This provides great flexibility, and dramatically improves the solver efficiency by reducing the degrees of freedom of the associated nonlinear system. However, it requires a particle push per nonlinear residual evaluation, which makes the particle push the most time-consuming operation in the algorithm. This paper describes a very efficient mixed-precision, hybrid CPU-GPU implementation of the implicit PIC algorithm. The JFNK solver is kept on the CPU (in double precision), while the inherent data parallelism of the particle mover is exploited by implementing it in single-precision on a graphics processing unit (GPU) using CUDA. Performance-oriented optimizations, with the aid of an analytical performance model, the roofline model, are employed. Despite being highly dynamic, the adaptive, charge-conserving particle mover algorithm achieves up to 300 400 GOp/s (including single-precision floating-point, integer, and logic operations) on a Nvidia GeForce GTX580, corresponding to 20 25% absolute GPU efficiency (against the peak theoretical performance) and 50-70% intrinsic efficiency (against the algorithm s maximum operational throughput, which neglects all latencies). This is about 200-300 times faster than an equivalent serial CPU implementation. When the single-precision GPU particle mover is combined with a double-precision CPU JFNK field solver, overall performance gains 100 vs. the double-precision CPU-only serial version are obtained, with no apparent loss of robustness or accuracy when applied to a challenging long-time scale ion acoustic wave simulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Laboratory evaluation of open source and commercial electrical conductivity sensor precision and accuracy: How do they compare?

Variation in the electrical conductivity (EC) of water can reveal environmental disturbance and natural dynamics, including factors such as anthropogenic salinization. Broader application of open source (OS) EC sensors could provide an inexpensive method to measure water quality. While studies show that other water quality parameters can be robustly measured with sensors, a similar effort is needed to evaluate the performance of OS EC sensors. To address this need, we evaluated the accuracy (mean error, %) and precision (sample standard deviation) of OS EC sensors in the laboratory via comparison to EC calibration standards using three different OS and OS/commercial-hybrid (OS/C) EC sensors and data logger configurations and two commercial (C) EC sensors and data logger configurations. We also evaluated the effect of cable length (7.5 m and 30 m) and sensor calibration on OS sensor accuracy and precision. We found a significant difference between OS sensor mean accuracy (3.08%) and all other sensors combined (9.23%). Our study also found that EC sensor precision decreased across all sensor configurations with increasing calibration standard EC. There was also a significant difference between OS sensor mean precision (2.85 μS/cm) and the mean precision of all other sensors combined (9.12 μS/cm). Cable length did not affect OS sensor precision. Furthermore, our results suggest that future research should include evaluating how performance is impacted by combining OS sensors with commercial data loggers as this study found significantly decreased performance in OS/commercial-hybrid sensor configurations. To increase confidence in the reliability of OS sensor data, more studies such as ours are needed to further quantify OS sensor performance in terms of accuracy and precision across different settings and OS sensor and data collection platform configurations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Simulated Half-Precision Implementation of Blocked QR Factorization and Graph Clustering Applications

We explored half-precision implementation of blocked QR algorithms with the following motivations: 1. New GPUs perform fast half-precision arithmetic (4 to 16 times as fast as doubleprecision). 2. QR factorization is a basic linear algebra tool useful for many physics and data analysis applications. 3. Communication-avoiding, parallelizable QR algorithms already exist for talland- skinny matrices. While the standard QR algorithms are highly unstable in half-precision, our numerical simulations show that the Tall-and-Skinny QR (TSQR) algorithm can improve the backward error of QR factorization. When using subspace iteration for graph clustering applications, half-precision accuracy in forming the eigenspace is sufficient for clustering with high precision and recall for some medium-scale benchmark problems. Note that all half-precision arithmetic were simulated with conversions to single-precision floats. Therefore, the results from this work are somewhat optimistic.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Developing nanometer scale, atomically precise metallo-catalysts with Molecular Lego

We propose to create atomically precise, highly robust, nanometer scale macromolecules that organize biomimetic metal binding groups in designed three-dimensional pockets to accelerate specific polymerization reactions. These atomically precise catalysts will create atomically precise polymers with improved mechanical and environmental properties from renewable building blocks. We seek to create and study nanometer scale catalysts that mimic metallo-enzymes in their activity and selectivity because these larger catalysts will be able to make more extensive non-covalent contacts with the transition state of the growing polymer. These catalysts will have extended lifetimes because they will better protect the metals and higher activity because additional reactive groups can be brought in close contact with the active site to enhance reactivity. They will be far more tolerant to extremes of temperature and will not denature in non-aqueous solvents because they are held together by multiple, strong, covalent bonds. Our approach to catalysts could be used to develop catalysts to replace noble metal-based catalysts when targeting other reactions, because the geometry of the active site enforced by the Molecular Lego scaffolding will access new reactivity, mimicking how nature uses earth abundant metals for much of its catalysis. We will initially target Lewis acid-based catalysts that assemble aliphatic polyesters with improved stereoselectivity, tacticity and complex alternation of monomers. Aliphatic polyesters have received growing attention as attractive, environmentally benign, and sustainable alternatives to polymers developed from petroleum feedstocks. Aliphatic polyesters are environmentally friendly because they undergo facile hydrolytic degradation to benign products, and they have excellent properties and high biocompatibility. We will assemble these catalysts using unique “Spiroligomers” (aka Molecular Lego) chemistry developed in the Schafmeister group combined with the organometallic chemistry expertise of the Dobereiner group. Polymers are an attractive synthetic target of nanoscale atomically precise catalysts because their resulting structure, stereochemistry and tacticity has profound impact on their properties (crystalline vs amorphous) and provides a readout of the catalytic mechanism. Polymer synthesis also serves as an example of “atomically precise manufacturing” where atomically precise nanoscale catalysts construct polymers with atomic precision at increasing length scales.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗