Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “test vector”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
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

Deploy Nalu/Kokkos algorithmic infrastructure with performance benchmarking.

The former Nalu interior heterogeneous algorithm design, which was originally designed to manage matrix assembly operations over all elemental topology types, has been modified to operate over homogeneous collections of mesh entities. This newly templated kernel design allows for removal of workset variable resize operations that were formerly required at each loop over a Sierra ToolKit (STK) bucket (nominally, 512 entities in size). Extensive usage of the Standard Template Library (STL) std::vector has been removed in favor of intrinsic Kokkos memory views. In this milestone effort, the transition to Kokkos as the underlying infrastructure to support performance and portability on many-core architectures has been deployed for key matrix algorithmic kernels. A unit-test driven design effort has developed a homogeneous entity algorithm that employs a team-based thread parallelism construct. The STK Single Instruction Multiple Data (SIMD) infrastructure is used to interleave data for improved vectorization. The collective algorithm design, which allows for concurrent threading and SIMD management, has been deployed for the core low-Mach element- based algorithm. Several tests to ascertain SIMD performance on Intel KNL and Haswell architectures have been carried out. The performance test matrix includes evaluation of both low- and higher-order methods. The higher-order low-Mach methodology builds on polynomial promotion of the core low-order control volume nite element method (CVFEM). Performance testing of the Kokkos-view/SIMD design indicates low-order matrix assembly kernel speed-up ranging between two and four times depending on mesh loading and node count. Better speedups are observed for higher-order meshes (currently only P=2 has been tested) especially on KNL. The increased workload per element on higher-order meshes bene ts from the wide SIMD width on KNL machines. Combining multiple threads with SIMD on KNL achieves a 4.6x speedup over the baseline, with assembly timings faster than that observed on Haswell architecture. The computational workload of higher-order meshes, therefore, seems ideally suited for the many-core architecture and justi es further exploration of higher-order on NGP platforms. A Trilinos/Tpetra-based multi-threaded GMRES preconditioned by symmetric Gauss Seidel (SGS) represents the core solver infrastructure for the low-Mach advection/diffusion implicit solves. The threaded solver stack has been tested on small problems on NREL's Peregrine system using the newly developed and deployed Kokkos-view/SIMD kernels. fforts are underway to deploy the Tpetra-based solver stack on NERSC Cori system to benchmark its performance at scale on KNL machines.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

A heterologous prime-boosting strategy with replicating Vaccinia virus vectors and plant-produced HIV-1 Gag/dgp41 virus-like particles

Showing modest efficacy, the RV144 HIV-1 vaccine clinical trial utilized a non-replicating canarypox viral vector and a soluble gp120 protein boost. Here we built upon the RV144 strategy by developing a novel combination of a replicating, but highly-attenuated Vaccinia virus vector, NYVAC-KC, and plant-produced HIV-1 virus-like particles (VLPs). Both components contained the full-length Gag and a membrane anchored truncated gp41 presenting the membrane proximal external region with its conserved broadly neutralizing epitopes in the pre-fusion conformation. We tested different prime/boost combinations of these components in mice and showed that the group primed with NYVAC-KC and boosted with both the viral vectors and plant-produced VLPs have the most robust Gag-specific CD8 T cell responses, at 12.7% of CD8 T cells expressing IFN-γ in response to stimulation with five Gag epitopes. The same immunization group elicited the best systemic and mucosal antibody responses to Gag and dgp41 with a bias towards IgG1. - Highlights: • We devised a prime/boost anti HIV-1 vaccination strategy modeled after RV144. • We used plant-derived virus-like particles (VLPs) consisting of Gag and dgp41. • We used attenuated, replicating vaccinia virus vectors expressing the same antigens. • The immunogens elicited strong cellular and humoral immune responses.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

Tidal Effects in Some Regular Black Holes

This paper is aimed to study the tidal forces produced by a class of regular black holes. We consider the radial infall of test particle and find radial as well as angular components of tidal forces by taking geodesic deviation equations. We also compute geodesic deviation vector by solving geodesic deviation equation numerically. It is concluded that a particle undergos either compression or stretching in radial or angular direction due to tidal forces.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Testing Lorentz invariance of dark matter with satellite galaxies

We develop the framework for testing Lorentz invariance in the dark matter sector using galactic dynamics. We consider a Lorentz violating (LV) vector field acting on the dark matter component of a satellite galaxy orbiting in a host halo. We introduce a numerical model for the dynamics of satellites in a galactic halo and for a galaxy in a rich cluster to explore observational consequences of such an LV field. The orbital motion of a satellite excites a time dependent LV force which greatly affects its internal dynamics. Our analysis points out key observational signatures which serve as probes of LV forces. These include modifications to the line of sight velocity dispersion, mass profiles and shapes of satellites. With future data and a more detailed modeling these signatures can be exploited to constrain a new region of the parameter space describing the LV in the dark matter sector.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Transverse force induced by a magnetized wake

The force on a test charge moving through a strongly magnetized plasma is calculated using linear response theory. Strong magnetization is found to generate a component of the force perpendicular to the velocity of the particle in the plane formed by the velocity and magnetic field vectors. Here, this transverse force is generated by an asymmetry with respect to the velocity vector in the induced electrostatic wake potential that is associated with the action of the Lorentz force on the background plasma. The direction depends on the speed of the test charge. If it is faster than a critical speed characteristic of the sound speed, it acts to reduce the component of velocity parallel to the magnetic field and to increase the gyroradius. In contrast, if the speed is below this critical speed, it acts to increase the velocity parallel to the magnetic field and to decrease the gyroradius. Because the transverse force is perpendicular to the velocity, it does not directly influence the total energy of the test charge. Nevertheless, it significantly alters the trajectory on a timescale associated with the Coulomb collision time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

SU-G-TeP4-12: Individual Beam QA for a Robotic Radiosurgery System Using a Scintillator Cone

Purpose: The targeting accuracy of the Cyberknife system is measured by end-to-end tests delivering multiple isocentric beams to a point in space. While the targeting accuracy of two representative beams can be determined by a Winston-Lutz-type test, no test is available today to determine the targeting accuracy of each clinical beam. We used a scintillator cone to measure the accuracy of each individual beam. Methods: The XRV-124 from Logos Systems Int’l is a scintillator cone with an imaging system that is able to measure individual beam vectors and a resulting error between planned and measured beam coordinates. We measured the targeting accuracy of isocentric and non-isocentric beams for a number of test cases using the Iris and the fixed collimator. The average difference between plan and measured beam position was 0.8–1.2mm across the collimator sizes and plans considered here. The max error for a single beam was 2.5mm for the isocentric plans, and 1.67mm for the non-isocentric plans. The standard deviation of the differences was 0.5mm or less. Conclusion: The CyberKnife System is specified to have an overall targeting accuracy for static targets of less than 0.95mm. In E2E tests using the XRV124 system we measure average beam accuracy between 0.8 to 1.23mm, with maximum of 2.5mm. We plan to investigate correlations between beam position error and robot position, and to quantify the effect of beam position errors on patient specific plans. Martina Descovich has received research support and speaker honoraria from Accuray.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Ectopic expression of anti-HIV-1 shRNAs protects CD8{sup +} T cells modified with CD4ζ CAR from HIV-1 infection and alleviates impairment of cell proliferation

Chimeric antigen receptors (CARs) are artificially engineered receptors that confer a desired specificity to immune effector T cells. As an HIV-1-specific CAR, CD4ζ CAR has been extensively tested in vitro as well as in clinical trials. T cells modified with this CAR mediated highly potent anti-HIV-1 activities in vitro and were well-tolerated in vivo, but exerted limited effects on viral load and reservoir size due to poor survival and/or functionality of the transduced cells in patients. We hypothesize that ectopic expression of CD4ζ on CD8{sup +} T cells renders them susceptible to HIV-1 infection, resulting in poor survival of those cells. To test this possibility, highly purified CD8{sup +} T cells were genetically modified with a CD4ζ-encoding lentiviral vector and infected with HIV-1. CD8{sup +} T cells were vulnerable to HIV-1 infection upon expression of CD4ζ as evidenced by elevated levels of p24{sup Gag} in cells and culture supernatants. Concurrently, the number of CD4ζ-modified CD8{sup +} T cells was reduced relative to control cells upon HIV-1 infection. To protect these cells from HIV-1 infection, we co-expressed two anti-HIV-1 shRNAs previously developed by our group together with CD4ζ. This combination vector was able to suppress HIV-1 infection without impairing HIV-1-dependent effector activities of CD4ζ. In addition, the number of CD4ζ-modified CD8{sup +} T cells maintained similar levels to that of the control even under HIV-1 infection. These results suggest that protecting CD4ζ-modified CD8{sup +} T cells from HIV-1 infection is required for prolonged HIV-1-specific immune surveillance. - Highlights: • Ectopic expression of CD4ζ CAR in CD8{sup +} T cells renders them susceptible to HIV-1 infection. • Co-expression of two anti-HIV-1 shRNAs protects CD4ζ CAR-modified CD8{sup +} T cells from HIV-1 infection. • Protecting CD4ζ CAR-modified CD8{sup +} T cells from HIV-1 infection suppresses its cytopathic effect.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Intrafraction Variability and Deformation Quantification in the Breast

Purpose: To evaluate intrafraction variability and deformation of the lumpectomy cavity (LC), breast, and nearby organs. Methods and Materials: Sixteen left-sided postlumpectomy and 1 bilateral breast cancer cases underwent free-breathing CT (FBCT) and 10-phase 4-dimensional CT (4DCT). Deformable image registration was used for deformation analysis and contour propagation of breast, heart, lungs, and LC between end-exhale and end-inhale 4DCT phases. Respiration-induced motion was calculated via centroid analysis. Two planning target volumes (PTVs) were compared: PTV{sub FBCT} from the FBCT volume with an isotropic 10 mm expansion (5 mm excursion and 5 mm setup error) and PTV{sub 4DCT} generated from the union of 4DCT contours with isotropic 5 mm margin for setup error. Volume and geometry were evaluated via percent difference and bounding box analysis, respectively. Deformation correlations between breast/cavity, breast/lung, and breast/heart were evaluated. Associations were tested between cavity deformation and proximity to chest wall and breast surface. Results: Population-based 3-dimensional vector excursions were 2.5 ± 1.0 mm (range, 0.8-3.8 mm) for the cavity and 2.0 ± 0.8 mm (range, 0.7-3.0 mm) for the ipsilateral breast. Cavity excursion was predominantly in the anterior and superior directions (1.0 ± 0.8 mm and −1.8 ± 1.2 mm, respectively). Similarly, for all cases, LCs and ipsilateral breasts yielded median deformation values in the superior direction. For 14 of 17 patients, the LCs and breast interquartile ranges tended toward the anterior direction. The PTV{sub FBCT} was 51.5% ± 10.8% larger (P<.01) than PTV{sub 4DCT}. Bounding box analysis revealed that PTV{sub FBCT} was 9.8 ± 1.2 (lateral), 9.0 ± 2.2 (anterior–posterior), and 3.9 ± 1.8 (superior–inferior) mm larger than PTV{sub 4DCT}. Significant associations between breast and cavity deformation were found for 6 of 9 axes. No dependency was found between cavity deformation and proximity to chest wall or breast surface. Conclusions: Lumpectomy cavity and breast deformation and motion demonstrated large variability. A PTV{sub 4DCT} approach showed value in patient-specific margins, particularly if robust interfraction setup analysis can be performed.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2021

EpitopeVec: linear epitope prediction using deep protein sequence embeddings

Motivation: B-cell epitopes (BCEs) play a pivotal role in the development of peptide vaccines, immuno-diagnostic reagents and antibody production, and thus in infectious disease prevention and diagnostics in general. Experimental methods used to determine BCEs are costly and time-consuming. Therefore, it is essential to develop computational methods for the rapid identification of BCEs. Although several computational methods have been developed for this task, generalizability is still a major concern, where cross-testing of the classifiers trained and tested on different datasets has revealed accuracies of 51–53%. Results: We describe a new method called EpitopeVec, which uses a combination of residue properties, modified antigenicity scales, and protein language model-based representations (protein vectors) as features of peptides for linear BCE predictions. Extensive benchmarking of EpitopeVec and other state-of-the-art methods for linear BCE prediction on several large and small datasets, as well as cross-testing, demonstrated an improvement in the performance of EpitopeVec over other methods in terms of accuracy and area under the curve. As the predictive performance depended on the species origin of the respective antigens (viral, bacterial and eukaryotic), we also trained our method on a large viral dataset to create a dedicated linear viral BCE predictor with improved cross-testing performance.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

CLMM : a LSST-DESC cluster weak lensing mass modeling library for cosmology

ABSTRACT We present the v1.0 release of CLMM, an open source python library for the estimation of the weak lensing masses of clusters of galaxies. CLMM is designed as a stand-alone toolkit of building blocks to enable end-to-end analysis pipeline validation for upcoming cluster cosmology analyses such as the ones that will be performed by the Vera C. Rubin Legacy Survey of Space and Time-Dark Energy Science Collaboration (LSST-DESC). Its purpose is to serve as a flexible, easy-to-install, and easy-to-use interface for both weak lensing simulators and observers and can be applied to real and mock data to study the systematics affecting weak lensing mass reconstruction. At the core of CLMM are routines to model the weak lensing shear signal given the underlying mass distribution of galaxy clusters and a set of data operations to prepare the corresponding data vectors. The theoretical predictions rely on existing software, used as backends in the code, that have been thoroughly tested and cross-checked. Combined theoretical predictions and data can be used to constrain the mass distribution of galaxy clusters as demonstrated in a suite of example Jupyter Notebooks shipped with the software and also available in the extensive online documentation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Confinement, quark mass functions, and spontaneous chiral symmetry breaking in Minkowski space

We formulate the covariant equations for quark-antiquark bound states in Minkowski space in the framework of the Covariant Spectator Theory. The quark propagators are dressed with the same kernel that describes the interaction between different quarks. We show that these equations are charge conjugation invariant, and that in the chiral limit of vanishing bare quark mass, a massless pseudoscalar bound state is produced in a Nambu--Jona-Lasinio (NJL) mechanism, which is associated with the Goldstone boson of spontaneous chiral symmetry breaking. In this introductory paper we test the formalism by using a simplified kernel consisting of a momentum-space $\delta$-function with a vector Lorentz structure, to which one adds a mixed scalar and vector confining interaction. The scalar part of the confining interaction is not chirally invariant by itself, but decouples from the equations in the chiral limit and therefore allows the NJL mechanism to work. With this model we calculate the quark mass function, and we compare our Minkowski-space results to LQCD data obtained in Euclidean space. In a companion paper we apply this formalism to a calculation of the pion form factor.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Performance of an Astrophysical Radiation Hydrodynamics Code under Scalable Vector Extension Optimization

We present results of a performance study of an astrophysical radiation hydrodynamics code, V2D, on the Arm-based A64FX processor developed by Fujitsu. The code solves sparse linear systems, a task for which the A64FX architecture should be well suited. Here, we performed the performance analysis study on Ookami, an Apollo 80 platform utilizing the A64FX processor. We explored several compilers and performance anal-ysis packages and found the code did not perform as expected under scalable vector extension optimization, suggesting that a “deeper dive” into analyzing the code is worthwhile. However, a simple driver program that exercised basic sparse linear algebra routines used by V2D did show significant speedup with the use of the scalable vector extension optimization. We present the initial results from the study which used V2D on a relatively simple test problem that emphasized the repeated solution of sparse linear systems.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

TU-H-CAMPUS-JeP2-03: Machine-Learning-Based Delineation Framework of GTV Regions of Solid and Ground Glass Opacity Lung Tumors at Datasets of Planning CT and PET/CT Images

Purpose: In radiation treatment planning, delineation of gross tumor volume (GTV) is very important, because the GTVs affect the accuracies of radiation therapy procedure. To assist radiation oncologists in the delineation of GTV regions while treatment planning for lung cancer, we have proposed a machine-learning-based delineation framework of GTV regions of solid and ground glass opacity (GGO) lung tumors following by optimum contour selection (OCS) method. Methods: Our basic idea was to feed voxel-based image features around GTV contours determined by radiation oncologists into a machine learning classifier in the training step, after which the classifier produced the degree of GTV for each voxel in the testing step. Ten data sets of planning CT and PET/CT images were selected for this study. The support vector machine (SVM), which learned voxel-based features which include voxel value and magnitudes of image gradient vector that obtained from each voxel in the planning CT and PET/CT images, extracted initial GTV regions. The final GTV regions were determined using the OCS method that was able to select a global optimum object contour based on multiple active delineations with a level set method around the GTV. To evaluate the results of proposed framework for ten cases (solid:6, GGO:4), we used the three-dimensional Dice similarity coefficient (DSC), which denoted the degree of region similarity between the GTVs delineated by radiation oncologists and the proposed framework. Results: The proposed method achieved an average three-dimensional DSC of 0.81 for ten lung cancer patients, while a standardized uptake value-based method segmented GTV regions with the DSC of 0.43. The average DSCs for solid and GGO were 0.84 and 0.76, respectively, obtained by the proposed framework. Conclusion: The proposed framework with the support vector machine may be useful for assisting radiation oncologists in delineating solid and GGO lung tumors.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

EKAT v.1.0

E3SM Kokkos Application Toolkit (EKAT) is a collection of C++, Fortran, and CMake utilities for providing a single implementation of common kernels based on the Kokkos programming model. The library contains utilities for vectorization, tridiagonal linear system solvers, and linear interpolation as well as some general-purpose utilities such as testing utilities, parameter lists, representation of physical units, and additional interfaces. The goal is to provide a centralized implementation for high-performance computing structures and common utilities that reduce code duplication and streamline maintenance efforts. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-NA0003525. SAND2022-1327 O

Bertagna, Luca↗
DOE OSTI2021

Overview of Algorithms for Using Particle Morphology in Pre-Detonation Nuclear Forensics

A major goal in pre-detonation nuclear forensics is to infer the processing conditions and/or facility type that produced radiological material. This review paper focuses on analyses of particle size, shape, texture (“morphology”) signatures that could provide information on the provenance of interdicted materials. For example, uranium ore concentrates (UOC or yellowcake) include ammonium diuranate (ADU), ammonium uranyl carbonate (AUC), sodium diuranate (SDU), magnesium diuranate (MDU), and others, each prepared using different salts to precipitate U from solution. Once precipitated, UOCs are often dried and calcined to remove adsorbed water. The products can be allowed to react further, forming uranium oxides UO3, U3O8, or UO2 powders, whose surface morphology can be indicative of precipitation and/or calcination conditions used in their production. This review paper describes statistical issues and approaches in using quantitative analyses of measurements such as particle size and shape to infer production conditions. Statistical topics include multivariate t tests (Hotelling’s T 2 ), design of experiments, and several machine learning (ML) options including decision trees, learning vector quantization neural networks, mixture discriminant analysis, and approximate Bayesian computation (ABC). ABC is emphasized as an attractive option to include the effects of model uncertainty in the selected and fitted forward model used for inferring processing conditions.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2017

Black holes in vector-tensor theories

We study static and spherically symmetric black hole (BH) solutions in second-order generalized Proca theories with nonminimal vector field derivative couplings to the Ricci scalar, the Einstein tensor, and the double dual Riemann tensor. We find concrete Lagrangians which give rise to exact BH solutions by imposing two conditions of the two identical metric components and the constant norm of the vector field. These exact solutions are described by either Reissner-Nordström (RN), stealth Schwarzschild, or extremal RN solutions with a non-trivial longitudinal mode of the vector field. We then numerically construct BH solutions without imposing these conditions. For cubic and quartic Lagrangians with power-law couplings which encompass vector Galileons as the specific cases, we show the existence of BH solutions with the difference between two non-trivial metric components. The quintic-order power-law couplings do not give rise to non-trivial BH solutions regular throughout the horizon exterior. The sixth-order and intrinsic vector-mode couplings can lead to BH solutions with a secondary hair. For all the solutions, the vector field is regular at least at the future or past horizon. The deviation from General Relativity induced by the Proca hair can be potentially tested by future measurements of gravitational waves in the nonlinear regime of gravity.

79 ASTRONOMY AND ASTROPHYSICS↗