Search NASA⌕ Search

SEARCH · Search NASA

Results for “Arithmetic”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Performance Characterization of De Novo Genome Assembly on Leading Parallel Systems.

De novo genome assembly is one of the most important and challenging computational problems in modern genomics; further, it shares algorithms and communication patterns important to other graph analytic and irregular applications. Unlike simulations, it has no floating point arithmetic and is dominated by small memory transactions within and between computing nodes. In this work, we focus on the highly scalable HipMer assembler and identify the dominant algorithms and communication patterns, also using microbenchmarks to capture the workload. We evaluate HipMer on a variety of platforms from the latest HPC systems to ethernet clusters. HipMer performs well on all single node systems, including the Xeon Phi manycore architecture. Given large enough problems, it also demonstrates excellent scaling across nodes in an HPC system, but requires a high speed network with low overhead and high injection rates. Our results shed light on the architectural features that are most important for achieving good parallel efficiency on this and related problems.

Ellis, Marquita↗
DOE OSTI2024

Two-loop master integrals for leading-color $$ pp\to t\overline{t}H $$ amplitudes with a light-quark loop

Abstract We compute the two-loop master integrals for leading-color QCD scattering amplitudes including a closed light-quark loop in$$ t\overline{t}H $$ t t ¯ H production at hadron colliders. Exploiting numerical evaluations in modular arithmetic, we construct a basis of master integrals satisfying a system of differential equations inϵ-factorized form. We present the analytic form of the differential equations in terms of a minimal set of differential one-forms. We explore properties of the function space of analytic solutions to the differential equations in terms of iterative integrals which can be exploited for studying the analytic form of related scattering amplitudes. Finally, we solve the differential equations using generalized series expansions to numerically evaluate the master integrals in physical phase space. As the first computation of a set of two-loop seven-scale master integrals, our results provide valuable input for analytic studies of scattering amplitudes in processes involving massive particles and a large number of kinematic scales.

Physics↗
DOE OSTI2024

Measurement of the effective leptonic weak mixing angle

Using pp collision data at $\sqrt{s}$ = 13 TeV, recorded by the LHCb experiment between 2016 and 2018 and corresponding to an integrated luminosity of 5.4 fb −1 , the forward-backward asymmetry in the pp → Z/γ* → μ + μ − process is measured. The measurement is carried out in ten intervals of the difference between the muon pseudorapidities, within a fiducial region covering dimuon masses between 66 and 116 GeV, muon pseudorapidities between 2.0 and 4.5 and muon transverse momenta above 20 GeV. These forward-backward asymmetries are compared with predictions, at next-to-leading order in the strong and electroweak couplings. The measured effective leptonic weak mixing angle is ${\sin}^2{\theta}_{\textrm{eff}}^{\ell }=0.23147\pm 0.00044\pm 0.00005\pm 0.00023,$ where the first uncertainty is statistical, the second arises from systematic uncertainties associated with the asymmetry measurement, and the third arises from uncertainties in the fit model used to extract ${\sin}^2{\theta}_{\textrm{eff}}^{\ell }$ from the asymmetry measurement. This result is based on an arithmetic average of results using the CT18, MSHT20, and NNPDF31 parameterisations of the proton internal structure, and is consistent with previous measurements and with predictions from the global electroweak fit.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Computational Errors in the Calculation of Long Radioactive Decay Chains

For safety security of objects utilizing atomic energy, the effect of numerical errors on the calculation of the activity and integrals of the activity in time when using Bateman’s formula and Siewers’ algorithm to calculate radioactive decay chains is investigated. The smallest errors obtain when using Siewers’ algorithm and Bateman’s formula together with long arithmetic. It is shown that the standard calculations of radioactive decay chains according to Bateman’s formula with double precision lead to significant errors. However, these errors have almost no effect on the dose characteristics of mixtures of radionuclides. The analysis is based on data contained in ICRP Publication No. 107.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2019

Application of nonlinear photonic crystal ring resonators in realizing all optical OR/NOT/AND gates

Arithmetic logic unit (ALU) is the core of any digital processing systems. For creating an all optical ALU one needs basic logic gates such as optical NOT, OR and AND gates along with an optical full adder and optical multiplexer. In this paper we aim to design and propose these three basic logic gates using nonlinear photonic crystal ring resonators. The simulations will be done using plane wave expansion and finite difference time domain methods. The proposed structures work based on controlling the optical behavior of the nonlinear resonators using optical power. The simulation results show that the switching threshold for this resonant ring is about 15 mW/μm{sup 2}. Reduced delay times and switching optical power are the main advantages of the proposed structures. Also by adding two resonant rings in input ports we reduced the cross reflection between the input ports.

42 ENGINEERING↗
DOE OSTI2022

Enhancing Surface Finish of Additively Manufactured 316L Stainless Steel with Pulse/Pulse Reverse Electropolishing

To improve the reliable integration of additively manufactured (AM) metal components, there is an urgent need for surface finishing methods capable of processing rough, “as-printed” surfaces. Here, we present a new and scalable electropolishing (EP) technique for use on AM 316L stainless steel (SS) that is capable of achieving a surface roughness equivalent to a machined finish, with arithmetic mean heights (S a ) as low as 1.0 µm. This technique utilizes the sequential application of a new, NaCl/H 2 O pulse/pulse reverse (P/PR) EP step along with a conventional DC (direct current) EP step, where P/PR expands the capability of EP by removing roughness features > 1 µm, which has been previously unattainable with DC EP alone without substantial material removal. This is achieved by independently varying pulse parameters, such as pulse width and height, allowing more control over electrodynamic processes at the part-electrolyte interface. Additionally, P/PR EP resulted in higher smoothing efficiency, where a higher degree of smoothing was achieved with less material removal when compared with conventional DC EP alone. This minimizes unwanted change in geometry from excessive metal dissolution.

36 MATERIALS SCIENCE↗
DOE OSTI2018

A cooperative conjugate gradient method for linear systems permitting efficient multi-thread implementation

This paper revisits, in a multi-thread context, the so-called multi-parameter or block conjugate gradient (B-CG) methods, first proposed as sequential algorithms by O’Leary and Brezinski, for the solution of the linear system Ax=b, for an n-dimensional symmetric positive definite matrix A. Instead of the scalar parameters of the classical CG algorithm, which minimizes a scalar functional at each iteration, multiple descent and conjugate directions are updated simultaneously. Implementation involves the use of multiple threads and the algorithm is referred to as cooperative CG (CCG) to emphasize that each thread now uses information that comes from the other threads. It is shown that for a sufficiently large matrix dimension n, the use of an optimal number of threads results in a worst case flop count of O(n{sup 7/3}) in exact arithmetic. Numerical experiments on a multi-core, multi-thread computer, for synthetic and real matrices, illustrate the theoretical results.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A structural equation modeling approach to leveraging the power of extant sentiment analysis tools

Machine-derived sentiment analysis has become a pervasive and useful tool to address a wide array of issues in natural language processing. Leading technology companies such as Google now provide sentiment analysis tools (SATs) as readily accessible online products. Academic researchers develop and make available SATs to support the research enterprise. One of the major challenges with SATs is the inconsistencies in results among the various SATs. Consequently, the selection of a SAT for a specific purpose may significantly impact the application. This study addresses the foregoing problem by utilizing structural equation modeling to merge the outputs of SATs to develop a combined sentiment metric without the need for a labeled training dataset. This method is applicable to a wide range of text-based problems, is data-driven, and replicable. It was tested using three publicly available datasets and compared against seven different SATs. The results indicate that as a continous measure, the proposed method outperformed other SATs in the movie reviews and SemEval datasets, and achieved a tie for first place with IBM Watson on the Sentiment 140 dataset. Also, compared to the published major alternatives, the arithmetic mean solution, this approach performed better across these three datasets.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Parametric model-order reduction for radiation transport using multi-resolution proper orthogonal decomposition

For parametric high-fidelity simulations, it is often desirable to utilize a reduced-order model (ROM) to emulate, at a reduced computational cost, parametric solutions of the governing partial differential equations (PDEs) for unseen parameter values. One commonly employed option is to utilize a data-driven, projection-based ROM supplemented with subspace identification via proper orthogonal decomposition (POD). POD discovers the ROM subspace by computing the singular value decomposition (SVD) of a set of training data from the full-order model (FOM). In streaming-dominated radiation transport simulations with localized sources, solutions often greatly vary over the spatial domain by many orders of magnitude. In such cases, machine-precision arithmetic can be insufficient to obtain an accurate SVD, resulting in a poorly performing ROM. We present a method called multiresolution POD (mrPOD) that mitigates these inaccuracies. The mrPOD method works by decomposing the spatial domain into regions and performing proper orthogonal decomposition on the training dataset separately in each region. In conclusion, mrPOD is tested on single energy group and multigroup atmospheric shielding transport problems and is shown to outperform classic POD.

42 ENGINEERING↗
DOE OSTI2019

Generalized Bose–Hubbard Hamiltonians exhibiting a complete non-Hermitian degeneracy

Highlights: • Hot topic: physics with non-Hermitian Hamiltonians. • Computer algebra used to construct new, next-to-BH models. • Degeneracy verification via variable-precision arithmetics. • Unfolding patterns at N>5 (recent problem solved). - Abstract: The method of construction of the tridiagonal and symmetric complex-matrix Hamiltonians H{sup (N)}(z) exhibiting an exceptional-point (EP) degeneracy of the Nth (i.e., maximal) order at a preselected parameter z=z{sup (EPN)}=1 is proposed and tested. In general, the implementation of the method requires the use of computer-assisted symbolic manipulations, especially at the larger matrix dimensions N. The well known PT−symmetric N-by-N-matrix Bose–Hubbard Hamiltonians as well as their recent N≤5 non-Bose–Hubbard alternatives are obtained as special solutions. Several other N≥6 non-Bose–Hubbard models H{sup (N)}(z) exhibiting the maximal EPN degeneracy are also constructed and analyzed in detail. In particular, their z−dependent, real or complex spectra of energies are displayed and discussed near as well as far from z{sup (EPN)}.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

A simulation-based efficiency comparison of AC and DC power distribution networks in commercial buildings

Direct current (DC) power distribution has recently gained traction in buildings research due to the proliferation of on-site electricity generation and battery storage, and an increasing prevalence of internal DC loads. The research discussed in this paper uses Modelica-based simulation to compare the efficiency of DC building power distribution with an equivalent alternating current (AC) distribution. The buildings are all modeled with solar generation, battery storage, and loads that are representative of the most efficient building technology. A variety of paramet ric simulations determine how and when DC distribution proves advantageous. These simulations also validate previous studies that use simpler approaches and arithmetic efficiency models. This work shows that using DC distribution can be considerably more efficient: a medium sized office building using DC distribution has an expected baseline of 12% savings, but may also save up to 18%. In these results, the baseline simulation parameters are for a zero net energy (ZNE) building that can island as a microgrid. DC is most advantageous in buildings with large solar capacity, large battery capacity, and high voltage distribution.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Robust trajectory-constrained frequency control for microgrids considering model linearization error

Grid supportive modes integrated within inverter-based resources can improve the frequency response of renewable-rich microgrids. The synthesis of grid supportive modes to guarantee frequency trajectory constraints under a predefined disturbance set is challenging but essential. To tackle this challenge, a numerical optimal control (NOC)-based control synthesis methodology is proposed. Without loss of generality, a wind-diesel fed microgrid is studied, where we aim to design grid supportive functions in the wind turbine. In the control design, linearized models are used, and the linearization-induced errors are quantitatively analyzed by reachability and interval arithmetics and represented in the form of interval uncertainties. Then, the NOC problem can be formulated into a robust mixed-integer linear program. The control structure is strategically configured into two levels to realize online deployment. The proposed control is verified on the modified 33-node microgrid with a full-order three-phase nonlinear model in Simulink. In conclusion, the simulation results show the effectiveness of the proposed control paradigm and the necessity of considering linearization-induced uncertainty.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2017

Parallel accelerated cyclic reduction preconditioner for three-dimensional elliptic PDEs with variable coefficients

We present a robust and scalable preconditioner for the solution of large-scale linear systems that arise from the discretization of elliptic PDEs amenable to rank compression. The preconditioner is based on hierarchical low-rank approximations and the cyclic reduction method. The setup and application phases of the preconditioner achieve log-linear complexity in memory footprint and number of operations, and numerical experiments exhibit good weak and strong scalability at large processor counts in a distributed memory environment. Numerical experiments with linear systems that feature symmetry and nonsymmetry, definiteness and indefiniteness, constant and variable coefficients demonstrate the preconditioner applicability and robustness. Furthermore, it is possible to control the number of iterations via the accuracy threshold of the hierarchical matrix approximations and their arithmetic operations, and the tuning of the admissibility condition parameter. Together, these parameters allow for optimization of the memory requirements and performance of the preconditioner.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Lightweight runtime checking of C programs with RTC

The C Programming Language is known for being an efficient language that can be compiled on almost any architecture and operating system. However the absence of dynamic safety checks and a relatively weak type system allows programmer oversights that are hard to spot. In this paper, we present RTC, a runtime monitoring tool that instruments unsafe code and monitors the program execution. RTC is built on top of the ROSE compiler infrastructure. RTC finds memory bugs and arithmetic overflows and underflows, and run-time type violations. Most of the instrumentations are directly added to the source file and only require a minimal runtime system. As a result, the instrumented code remains portable. In tests against known error detection benchmarks, RTC found 98% of all memory related bugs and had zero false positives. Finally, in performance tests conducted with well known algorithms, such as binary search and MD5, we determined that our tool has an average run-time overhead rate of 9.7× and memory overhead rate of 3.5×.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

An efficient and portable SIMD algorithm for charge/current deposition in Particle-In-Cell codes

In current computer architectures, data movement (from die to network) is by far the most energy consuming part of an algorithm (≈20pJ/word on-die to ≈10,000 pJ/word on the network). To increase memory locality at the hardware level and reduce energy consumption related to data movement, future exascale computers tend to use many-core processors on each compute nodes that will have a reduced clock speed to allow for efficient cooling. To compensate for frequency decrease, machine vendors are making use of long SIMD instruction registers that are able to process multiple data with one arithmetic operator in one clock cycle. SIMD register length is expected to double every four years. As a consequence, Particle-In-Cell (PIC) codes will have to achieve good vectorization to fully take advantage of these upcoming architectures. In this paper, we present a new algorithm that allows for efficient and portable SIMD vectorization of current/charge deposition routines that are, along with the field gathering routines, among the most time consuming parts of the PIC algorithm. Our new algorithm uses a particular data structure that takes into account memory alignment constraints and avoids gather/scat;ter instructions that can significantly affect vectorization performances on current CPUs. The new algorithm was successfully implemented in the 3D skeleton PIC code PICSAR and tested on Haswell Xeon processors (AVX2-256 bits wide data registers). Results show a factor of ×2 to ×2.5 speed-up in double precision for particle shape factor of orders 1–3. The new algorithm can be applied as is on future KNL (Knights Landing) architectures that will include AVX-512 instruction sets with 512 bits register lengths (8 doubles/16 singles). Program summary Program Title: vec_deposition Program Files doi:http://dx.doi.org/10.17632/nh77fv9k8c.1 Licensing provisions: BSD 3-Clause Programming language: Fortran 90 External routines/libraries: OpenMP > 4.0 Nature of problem: Exascale architectures will have many-core processors per node with long vector data registers capable of performing one single instruction on multiple data during one clock cycle. Data register lengths are expected to double every four years and this pushes for new portable solutions for efficiently vectorizing Particle-In-Cell codes on these future many-core architectures. One of the main hotspot routines of the PIC algorithm is the current/charge deposition for which there is no efficient and portable vector algorithm. Solution method: Here we provide an efficient and portable vector algorithm of current/charge deposition routines that uses a new data structure, which significantly reduces gather/scatter operations. Vectorization is controlled using OpenMP 4.0 compiler directives for vectorization which ensures portability across different architectures. Restrictions: Here we do not provide the full PIC algorithm with an executable but only vector routines for current/charge deposition. These scalar/vector routines can be used as library routines in your 3D Particle-In-Cell code. However, to get the best performances out of vector routines you have to satisfy the two following requirements: (1) Your code should implement particle tiling (as explained in the manuscript) to allow for maximized cache reuse and reduce memory accesses that can hinder vector performances. The routines can be used directly on each particle tile. (2) You should compile your code with a Fortran 90 compiler (e.g Intel, gnu or cray) and provide proper alignment flags and compiler alignment directives (more details in README file).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

GPU-accelerated Tersoff potentials for massively parallel Molecular Dynamics simulations

The Tersoff potential is one of the empirical many-body potentials that has been widely used in simulation studies at atomic scales. Unlike pair-wise potentials, the Tersoff potential involves three-body terms, which require much more arithmetic operations and data dependency. In this contribution, we have implemented the GPU-accelerated version of several variants of the Tersoff potential for LAMMPS, an open-source massively parallel Molecular Dynamics code. Compared to the existing MPI implementation in LAMMPS, the GPU implementation exhibits a better scalability and offers a speedup of 2.2X when run on 1000 compute nodes on the Titan supercomputer. On a single node, the speedup ranges from 2.0 to 8.0 times, depending on the number of atoms per GPU and hardware configurations. The most notable features of our GPU-accelerated version include its design for MPI/accelerator heterogeneous parallelism, its compatibility with other functionalities in LAMMPS, its ability to give deterministic results and to support both NVIDIA CUDA- and OpenCL-enabled accelerators. Our implementation is now part of the GPU package in LAMMPS and accessible for public use.

Computer Science↗
DOE OSTI2021

Caravel: A C++ framework for the computation of multi-loop amplitudes with numerical unitarity

We present the first public version of Caravel, a C++17 framework for the computation of multi-loop scattering amplitudes in quantum field theory, based on the numerical unitarity method. Caravel is composed of modules for the D-dimensional decomposition of integrands of scattering amplitudes into master and surface terms, the computation of tree-level amplitudes in floating point or finite-field arithmetic, the numerical computation of one- and two-loop amplitudes in QCD and Einstein gravity, and functional reconstruction tools. Here, we provide programs that showcase Caravel's main functionalities and allow to compute selected one- and two-loop amplitudes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Thyroid antagonists and thyroid indicators in U.S. pregnant women in the Vanguard Study of the National Children's Study

The sodium iodide-symporter (NIS) mediates uptake of iodide into thyroid follicular cells. This key step in thyroid hormone synthesis is inhibited by perchlorate, thiocyanate (SCN) and nitrate (NO{sub 3}) anions. When these exposures occur during pregnancy the resulting decreases in thyroid hormones may adversely affect neurodevelopment of the human fetus. Our objectives were to describe and examine the relationship of these anions to the serum thyroid indicators, thyroid stimulating hormone (TSH) and free thyroxine (FT4), in third trimester women from the initial Vanguard Study of the National Children's Study (NCS); and to compare urine perchlorate results with those in pregnant women from the National Health and Nutritional Examination Survey (NHANES). Urinary perchlorate, SCN, NO{sub 3}, and iodine, serum TSH, FT4, and cotinine were measured and a food frequency questionnaire (FFQ) was administered to pregnant women enrolled in the initial Vanguard Study. We used multiple regression models of FT4 and TSH that included perchlorate equivalent concentration (PEC, which estimates combined inhibitory effects of the anions perchlorate, SCN, and NO{sub 3} on the NIS). We used multiple regression to model predictors of each urinary anion, using FFQ results, drinking water source, season of year, smoking status, and demographic characteristics. Descriptive statistics were calculated for pregnant women in NHANES 2001–2012. The geometric mean (GM) for urinary perchlorate was 4.04 µg/L, for TSH 1.46 mIU/L, and the arithmetic mean for FT4 1.11 ng/dL in 359 NCS women. In 330 women with completed FFQs, consumption of leafy greens, winter season, and Hispanic ethnicity were significant predictors of higher urinary perchlorate, which differed significantly by study site and primary drinking water source, and bottled water was associated with higher urinary perchlorate compared to filtered tap water. Leafy greens consumption was associated with higher urinary NO{sub 3} and higher urinary SCN. There was no association between urinary perchlorate or PEC and TSH or FT4, even for women with urinary iodine <100 µg/L. GM urinary perchlorate concentrations in the full sample (n=494) of third trimester NCS women (4.03 µg/L) were similar to pregnant women in NHANES (3.58 µg/L). - Highlights: • Limited biomonitoring data are available in pregnant women. • We found no effect of perchlorate and related anions on thyroid indicators in third trimester women. • Leafy green vegetables were a dietary source for perchlorate, nitrate, and thiocyanate.

54 ENVIRONMENTAL SCIENCES↗