Search NASA⌕ Search

SEARCH · Search NASA

Results for “Tensor compression”

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 55 records · Page 3

Towards Compact Neural Networks via End-to-End Training: A Bayesian Tensor Approach with Automatic Rank Determination

Post-training model compression can reduce the inference costs of deep neural networks, but uncompressed training still consumes enormous hardware resources and energy. To enable low-energy training on edge devices, it is highly desirable to directly train a compact neural network from scratch with a low memory cost. Low-rank tensor decomposition is an effective approach to reduce the memory and computing costs of large neural networks. However, directly training low-rank tensorized neural networks is a very challenging task because it is hard to determine a proper tensor rank a priori, and the tensor rank controls both model complexity and accuracy. Here, this paper presents a novel end-to-end framework for low-rank tensorized training. We first develop a Bayesian model that supports various low-rank tensor formats (e.g., CANDECOMP/PARAFAC, Tucker, tensor-train, and tensor-train matrix) and reduces neural network parameters with automatic rank determination during training. Then we develop a customized Bayesian solver to train large-scale tensorized neural networks. Our training methods shows orders-of-magnitude parameter reduction and little accuracy loss (or even better accuracy) in the experiments. On a very large deep learning recommendation system with over 4.2 ×10 9 model parameters, our method can reduce the parameter number to 1.6 ×10 5 automatically in the training process (i.e., by 2.6 ×10 4 times) while achieving almost the same accuracy. Code is available at https://github.com/colehawkins/bayesian-tensor-rank-determination.

compact neural networks↗

Kinetic study of strong shock waves in fully ionized plasmas

Fully kinetic simulations are used in this work to study strong shock waves in fully ionized plasmas for charge state 𝑍 ranging from 1 to 18. The ion stress tensor amplitude agrees well with classical predictions within the compression layer; however, its spatial structure is noticeably shaped by nonlocal transport effects. The electron stress tensor remains significantly smaller than that of the ions for all values of 𝑍 considered. Ion overheating in the compression layer is observed to scale approximately linearly with 𝑍. Both findings diverge from hydrodynamic models based on classical electron stress formulations. In conclusion, the electrons exhibit strong nonlocal behavior in the compression layer, with Knudsen numbers approaching unity for 𝑍 ≥ 6.

Physics - Plasma physics↗

Generic, Sparse Tensor Core for Neural Networks

Sparse neural network attracts more attention for model compression, fast execution, and power reduction. The state-of-the-art designed sparse tensor core for structured and static sparsity, which did not support well for generic or dynamic sparsity. We design a sparse tensor core to support generic sparsity pruning with a novel hybrid and blocked sparse matrix storage format, HB-ELL, which saves computation and storage while keeping the most significant elements, as well as supporting dynamic sparsity for data flow in neural networks. We achieve better performance with preliminary results than the state-of- the-art on an NVIDIA GPU simulator.

Wu, Xiaolong↗

Coronal heating by the resonant absorption of Alfven waves: The effect of viscous stress tensor

The time-dependent linearized magnetohydrodynamics (MHD) equations for a fully compressible, low-beta, viscoresistive plasma are solved numerically using an implicit integration scheme. The full viscosity stress tensor (Braginskii 1965) is included with the five parameters eta(sub i) i = 0 to 4. In agreement with previous studies, the numerical simulations demonstrate that the dissipation on inhomogeneities in the background Alfven speed occurs in a narrow resonant layer. For an active region in the solar corona the values of eta(sub i) are eta(sub o) = 0.65 g/cm/s, eta(sub 1) = 3.7 x 10(exp -12) g/cm/s, eta(sub 2) = 4 eta(sub 1), eta(sub 3) = 1.4 x 10(exp -6) g/cm/s, eta(sub 4) = 2 eta(sub 3), with n = 10(exp 10)/cu cm, T = 2 x 10(exp 6) K, and B = 100 G. When the Lundquist number S = 10(exp 4) and R(sub 1) much greater than S (where R(sub 1) is the dimensionless shear viscous number) the width of the resistive dissipation layer d(sub r) is 0.22a (where a is the density gradient length scale) and d(sub r) approximately S(exp -1/3). When S much greater than R(sub 1) the shear viscous dissipation layer width d(sub r) scales as R(sub 1)(exp -1/3). The shear viscous and the resistive dissipation occurs in an overlapping narrow region, and the total heating rate is independent of the value of the dissipation parameters in agreement with previous studies. Consequently, the maximum values of the perpendicular velocity and perpendicular magnetic field scale as R(sub 1)(exp -1/3). It is evident from the simulations that for solar parameters the heating due to the compressive viscosity (R(sub 0) = 560) is negligible compared to the resistive and the shear viscous (R(sub 1)) dissipation and it occurs in a broad layer of order a in width. In the solar corona with S approximately equals 10(exp 4) and R(sub 1) approximately equals 10(exp 14) (as calculated from the Braginskii expressions), the shear viscous resonant heating is of comparable magnitude to the resistive resonant heating.

Ofman, L.↗

Analysis of kink band formation under compression

The kink band formation in unidirectional composites under compression is analyzed in the present paper. The kinematics of kink band formation is described in terms of a deformation tensor. Equilibrium conditions are then applied to relate the compression load to the deformation of fibers. Since the in situ shear behavior of the matrix resin is not known, an analysis-experiment correlation is used to find the shear failure strain in the kink band. The present analysis thus elucidates the mechanisms and identifies the controlling parameters, of compression failure.

Hahn, H. Thomas↗

Characterization of anisotropic mechanical properties of polymer composites from a hybrid additive manufacturing-compression molding process using x-ray computer tomography

In this study, anisotropic stiffness tensors were reconstructed based on fiber orientation distributions obtained from X-ray computer tomography (xCT). A preform was manufactured via a big area additive manufacturing (BAAM) system with carbon fiber (CF) filled acrylonitrile butadiene styrene (ABS). The tailored preform from additive manufacturing (AM) was used in the compression molding (CM) process to produce a low-void high-performance thermoplastic composite panel. An xCT technique was employed to detect the fiber orientations in CF/ABS composites manufactured via three different methods: AM from BAAM, extrusion compression molding (ECM), and AM-CM. The anisotropic stiffness tensor was obtained from the composite panel manufactured via the three manufacturing methods (AM, ECM, and AMCM). A micromechanics theory was used to obtain the orthotropic stiffness tensors of the composite panels and compared with the experimental values. The predicted stiffness tensors of AM and AM-CM composite panels were used to study the deformation characteristics of a steering wheel during airbag deployment by performing finite element analysis (FEA). The approach developed in this study can be utilized for evaluating high-performance composites.

Pokkalla, Deepak↗

Local well-posedness and singularity formation in non-Newtonian compressible fluids

We investigate the initial value problem of a very general class of 3 + 1 non-Newtonian compressible fluids in which the viscous stress tensor with shear and bulk viscosity relaxes to its Navier–Stokes values. These fluids correspond to the non-relativistic limit of well-known Israel–Stewart-like theories used in the relativistic fluid dynamic simulations of high-energy nuclear and astrophysical systems. After establishing the local well-posedness of the Cauchy problem, we show for the first time in the literature that there exists a large class of initial data for which the corresponding evolution breaks down in finite time due to the formation of singularities. This implies that a large class of non-Newtonian fluids do not have finite solutions defined at all times.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Neutron (and other Particle) Transport at LANL: An Overview [Presentation]

For decades, Los Alamos National Laboratory has been at the forefront of neutron transport methods research and code development. One such code is PARTISN, the LANL parallel time-dependent discrete ordinate neutron transport code. In this presentation, we describe the various research efforts currently underway by the PARTISN and other code teams. Some examples of current research are a block automated mesh refinement scheme, the application of tensor trains to the discretized neutron transport equation, and GPU code porting. The block automated mesh refinement scheme uses cross section information to refine and coarsen the solution mesh to improve time to solution and reduce memory. The tensor train approach expresses discretized transport operators as tensor products of vectors and matrices to compress the size of linear systems being solved by transport codes. Rather than relying on matrix-free methods such as the transport sweep, we have access to an operator that can be inverted, reshaped, or manipulated algebraically. Finally, we describe how PARTISN is used, what problems we are looking to solve, and what the future holds for neutron transport at LANL. In addition to this, we briefly describe the various research efforts in other particle transport teams using both deterministic and Monte Carlo methods. In the presentation, we list possible opportunities for collaboration between the laboratory and faculty and students.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Efficient 3-D finite element failure analysis of compression loaded angle-ply plates with holes

Finite element stress analysis and the tensor polynomial failure criterion predict that failure always initiates at the interface between layers on the hole edge for notched angle-ply laminates loaded in compression. The angular location of initial failure is a function of the fiber orientation in the laminate. The dominant stress components initiating failure are shear. It is shown that approximate symmetry can be used to reduce the computer resources required for the case of unaxial loading.

Burns, S. W.↗

Fast increased fidelity samplers for approximate Bayesian Gaussian process regression

Gaussian processes (GPs) are common components in Bayesian non-parametric models having a rich methodological literature and strong theoretical grounding. The use of exact GPs in Bayesian models is limited to problems containing several thousand observations due to their prohibitive computational demands. We develop a posterior sampling algorithm using H-matrix approximations that scales at O(n log 2 n). We show that this approximation’s Kullback-Leibler divergence to the true posterior can be made arbitrarily small. Though multidimensional GPs could be used with our algorithm, d-dimensional surfaces are modeled as tensor products of univariate GPs to minimize the cost of matrix construction and maximize computational efficiency. We illustrate the performance of this fast increased fidelity approximate GP, FIFA-GP, using both simulated and non-synthetic data sets

97 MATHEMATICS AND COMPUTING↗

TensorID v1.0

This Python software package includes new and efficient algorithms for satellite and core interpolative decomposition of tensor data. In general, these algorithms target high-dimensional data reduction and compression. The software is purely numerical and can be applied by others to many important sources of tensor data generated by computation or experiment.

Zhang, Yifan [Lawrence Berkeley National Laborator↗

Failure analysis of composite laminates including biaxial compression

This report describes a continued effort on the development and application of the tensor polynomial failure criterion for composite laminate analysis. In particular, emphasis is given to the design, construction and testing of a cross-beam laminate configuration to obtain "pure' biaxial compression failure. The purpose of this test case was to provide to permit "closure' of the cubic form of the failure surface in the 1-2 compression-compression quadrant. This resulted in a revised set of interaction strength parameters and the construction of a failure surface which can be used with confidence for strength predictions, assuming a plane stress state exists. Furthermore, the problem of complex conjugate roots which can occur in some failure regions is addressed and an "engineering' interpretation is provided. Results are presented illustrating this behavior and the methodology for overcoming this problem is discussed.

Tennyson, R. C.↗

SymProp: Scaling Sparse Symmetric Tucker Decomposition via Symmetry Propagation

Sparse symmetric tensors are an important class of tensors, and their decompositions serve as powerful tools for revealing low-rank structures. This paper introduces SymProp, a novel approach for scaling sparse symmetric Tucker decomposition by propagating symmetry through intermediate computations. SymProp optimizes two key computational kernels: Sparse Symmetric Tensor Times Same Matrix chain (S3 TTMc) for Higher-Order Orthogonal Iteration (HOOI) and Sparse Symmetric Tensor Times Same Matrix chain Times Core (S3 TTMcTC) for Higher-Order QR Iteration (HOQRI). Our method employs a metaprogramming-based index iteration approach to efficiently handle the upper triangular parts of intermediate dense symmetric tensors. SymProp achieves up to 50.9× speedup over SPLATT and up to 360.8× over Compressed Sparse Symmetric (CSS) format on the S3 TTMc operation. Moreover, our S3 TTMc and S3 TTMcTC implementations support tensor orders four levels higher than state-of-the-art methods. Our HOQRI demonstrates superior scalability and up to a 33.6× speedup over optimized HOOI. By enabling more scalable Tucker decompositions for higher orders, decomposition ranks, and dimension sizes, SymProp opens new possibilities for analyzing complex hypergraph structures in fields such as network science, data mining, and machine learning.

Li, Zecheng [North Carolina State University]↗

Quantitative kinetic rules for plastic strain-induced α - ω phase transformation in Zr under high pressure

Plastic strain-induced phase transformations (PTs) and chemical reactions under high pressure are broadly spread in modern technologies, friction and wear, geophysics, and astrogeology. However, because of very heterogeneous fields of plastic strain $E$ p and stress σ tensors and volume fraction c of phases in a sample compressed in a diamond anvil cell (DAC) and impossibility of measurements of σ and $E$ p , there are no strict kinetic equations for them. Here, we develop a kinetic model, finite element method (FEM) approach, and combined FEM-experimental approaches to determine all fields in strongly plastically predeformed Zr compressed in DAC, and specific kinetic equation for α-ω PT consistent with experimental data for the entire sample. Since all fields in the sample are very heterogeneous, data are obtained for numerous complex 7D paths in the space of 3 components of the plastic strain tensor and 4 components of the stress tensor. Kinetic equation depends on accumulated plastic strain (instead of time) and pressure and is independent of plastic strain and deviatoric stress tensors, i.e., it can be applied for various above processes. Our results initiate kinetic studies of strain-induced PTs and provide efforts toward more comprehensive understanding of material behavior in extreme conditions.

36 MATERIALS SCIENCE↗

Direct numerical simulation of buoyantly driven turbulence

Numerical simulations of homogeneous turbulence subject to buoyant forcing were performed. The presence of a mean temperature gradient combined with a gravitational field results in a forcing term in the momentum equations. The development of the turbulence was studied and compared to the decay of similar fields in the absence of gravity. In the buoyantly driven field, the vorticity is preferentially aligned with the intermediate eigenvector of the strain-rate tensor and the local temperature gradient is more likely to be aligned with the most compressive eigenvector. These relationships are qualitatively similar to those observed in previous shear flow results studied by Ashurst (1987). A tensor diffusivity model for passive scalar transport developed from shear flow results in Rogers, Moin, and Reynolds (1986) also predicts this buoyant scalar transport, indicating that the relationship between the scalar flux and the Reynolds stress is similar in both flows.

Ashurst, William T.↗

Direct interpolative construction of the discrete Fourier transform as a matrix product operator

The quantum Fourier transform (QFT), which can be viewed as a reindexing of the discrete Fourier transform (DFT), has been shown to be compressible as a low-rank matrix product operator (MPO) or quantized tensor train (QTT) operator. However, the original proof of this fact does not furnish a construction of the MPO with a guaranteed error bound. Meanwhile, the existing practical construction of this MPO, based on the compression of a quantum circuit, is not as efficient as possible. We present a simple closed-form construction of the QFT MPO using the interpolative decomposition, with guaranteed near-optimal compression error for a given rank. This construction can speed up the application of the QFT and the DFT, respectively, in quantum circuit simulations and QTT applications. We also connect our interpolative construction to the approximate quantum Fourier transform (AQFT) by demonstrating that the AQFT can be viewed as an MPO constructed using a different interpolation scheme.

97 MATHEMATICS AND COMPUTING↗

Residual strain orientation in rolled titanium determined with synchrotron X-ray Laue microdiffraction

Previously, synchrotron X-ray Laue microdiffraction has been used to measure the magnitudes of residual strain in materials. Recently the method was advanced to determine the orientation of the strain ellipsoid and applied to naturally deformed quartzites; however, the deformation history of these quartzites is ambiguous due to their natural origin. Here, in this study, synchrotron X-ray Laue microdiffraction (µXRD) is used to measure the residual strain for the first time in a sample with known stress history, rolled titanium. A deviatoric strain tensor is calculated from each Laue diffraction image collected with two µXRD scans of a rolled titanium sheet in different sample orientations. The principal strain axes are calculated using an eigen decomposition of the deviatoric strain tensors. The results show that the principal axis of compression is aligned with the normal direction of the titanium sheet, and the principal axis of extension is aligned with the rolling direction. Pole figures are used to represent the 3D distribution of residual strain axes.

36 MATERIALS SCIENCE↗