SEARCH · Search NASA

Results for “Tensor factorizations”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 397 records · Page 22

All sources NASA NTRS DOE OSTI

DOE OSTI2018

High-resolution sub-Doppler infrared spectroscopy of atmospherically relevant Criegee precursor CH 2 I radicals: CH 2 stretch vibrations and “charge-sloshing” dynamics

The combination of a pulsed supersonic slit-discharge source and single-mode difference frequency direct absorption infrared spectroscopy permit first high resolution infrared study of the iodomethyl (CH 2 I) radical, with the CH 2 I radical species generated in a slit jet Ne/He discharge and cooled to 16 K in the supersonic expansion. Dual laser beam detection and collisional collimation in the slit expansion yield sub-Doppler linewidths (60 MHz), an absolute frequency calibration of 13 MHz, and absorbance sensitivities within a factor of two of the shot-noise limit. Fully rovibrationally resolved direct absorption spectra of the CH 2 symmetric stretch mode (ν 2 ) are obtained and fitted to a Watson asymmetric top Hamiltonian with electron spin-rotation coupling, providing precision rotational constants and spin-rotation tensor elements for the vibrationally excited state. Assessment of the asymmetric top rotational constants confirms a vibrationally averaged planar geometry in both the ground- and first-excited vibrational levels. Sub-Doppler resolution permits additional nuclear spin hyperfine structures to be observed, with splittings in excellent agreement with microwave measurements on the ground state. Spectroscopic data on CH 2 I facilitate systematic comparison with previous studies of halogen-substituted methyl radicals, with the periodic trends strongly correlated with the electronegativity of the halogen atom. Surprisingly, we do not observe any asymmetric CH 2 stretch transitions, despite S/N ≈ 25:1 on strongest lines in the corresponding symmetric CH 2 stretch manifold. This dramatic reversal of the more typical 3:1 antisymmetric/symmetric CH 2 stretch intensity ratio signals a vibrational transition moment poorly described by simple “bond-dipole” models. Instead, the data suggest that this anomalous intensity ratio arises from “charge sloshing” dynamics in the highly polar carbon-iodine bond, as supported by ab initio electron differential density plots and indeed consistent with observations in other halomethyl radicals and protonated cluster ions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

High-performance strategies for the recent MRSF-TDDFT in GAMESS

Multiple ERI (Electron Repulsion Integral) tensor contractions (METC) with several matrices are ubiquitous in quantum chemistry. In response theories, the contraction operation, rather than ERI computations, can be the major bottleneck, as its computational demands are proportional to the multiplicatively combined contributions of the number of excited states and the kernel pre-factors. Here, this paper presents several high-performance strategies for METC. Optimal approaches involve either the data layout reformations of interim density and Fock matrices, the introduction of intermediate ERI quartet buffer, and loop-reordering optimization for a higher cache hit rate. The combined strategies remarkably improve the performance of the MRSF (mixed reference spin flip)-TDDFT (time-dependent density functional theory) by nearly 300%. The results of this study are not limited to the MRSF-TDDFT method and can be applied to other METC scenarios.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Real-time scattering in Ising field theory using matrix product states

We study scattering in Ising field theory (IFT) using matrix product states and the time-dependent variational principle. IFT is a one-parameter family of strongly coupled nonintegrable quantum field theories in 1+1 dimensions, interpolating between massive free fermion theory and Zamolodchikov's integrable massive 𝐸 8 theory. Particles in IFT may scatter either elastically or inelastically. In the postcollision wave function, particle tracks from all final-state channels occur in superposition; processes of interest can be isolated by projecting the wave function onto definite particle sectors, or by evaluating energy density correlation functions. Using numerical simulations we determine the time delay of elastic scattering and the probability of inelastic particle production as a function of collision energy. We also study the mass and width of the lightest resonance near the 𝐸 8 point in detail. Close to both the free fermion and 𝐸 8 theories, our results for both elastic and inelastic scattering are in good agreement with expectations from form-factor perturbation theory. Using numerical computations to go beyond the regime accessible by perturbation theory, we find that the high-energy behavior of the two-to-two particle scattering probability in IFT is consistent with a conjecture of Zamolodchikov. Our results demonstrate the efficacy of tensor-network methods for simulating the real-time dynamics of strongly coupled quantum field theories in 1+1 dimensions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Reexamining the relation between the binding energy of finite nuclei and the equation of state of infinite nuclear matter

The energy density is calculated in coordinate space for 12 C, 40 Ca, 48 Ca, and 208 Pb using a dispersive optical model constrained by all relevant data including the corresponding energy of the ground state. The energy density of 8 Be is also calculated using the Green’s function Monte-Carlo method employing the Argonne/Urbana two and three-body interactions. The nuclear interior minimally contributes to the total binding energy due to the 4πr 2 phase space factor. Thus, the volume contribution to the energy in the interior is not well constrained. The dispersive-optical-model energy densities are in good agreement with ab initio self-consistent Green’s function calculations of infinite nuclear matter restricted to treat only short-range and tensor correlations. These results call into question the degree to which the equation of state for nuclear matter is constrained by the empirical mass formula. In particular, the results in this work indicate that saturated nuclear matter does not require the canonical value of 16 MeV binding per particle but only about 13-14 MeV when the interior of 208 Pb is considered.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2011

Embedding DBI inflation in scalar-tensor theory

The Dirac-Born-Infeld (DBI) action has been widely studied as an interesting example of a model of k-inflation in which the sound speed of the cosmological perturbations differs from unity. In this article we consider a scalar-tensor theory in which the matter component is a field with a DBI action. Transforming to the Einstein frame, we explore the effect of the resulting coupling on the background dynamics of the fields and the first-order perturbations. We find that the coupling forces the scalar field into the minimum of its effective potential. While the additional scalar field contributes significantly to the energy density during inflation, the dynamics are determined by the DBI field, which has the interesting effect of increasing the number of efolds of inflation and decreasing the boost factor of the DBI field. Focusing on this case, we show, with the benefit of numerical examples, that the power spectrum of the primordial perturbations is determined by the behaviour of the perturbations of the modified DBI field.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

M5-brane and D-brane scattering amplitudes

We present tree-level $\mathcal{n}$-particle on-shell scattering amplitudes of various brane theories with 16 conserved supercharges. These include the world-volume theory of a probe D3-brane or D5-brane in 10D Minkowski spacetime as well as a probe M5-brane in 11D Minkowski spacetime, which describes self interactions of an abelian tensor supermultiplet with 6D (2, 0) supersymmetry. Twistor-string-like formulas are proposed for tree-level scattering amplitudes of all multiplicities for each of these theories. The R symmetry of the D3-brane theory is shown to be SU(4) × U(1), and the U(1) factor implies that its amplitudes are helicity conserving. Each of 6D theories (D5-brane and M5-brane) reduces to the D3-brane theory by dimensional reduction. Conclusively, as special cases of the general M5-brane amplitudes, we present compact formulas for examples involving only the self-dual B field with $\mathcal{n}$ = 4, 6, 8.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2017

Tensor-structured coupled cluster theory

We derive and implement a new way of solving coupled cluster equations with lower computational scaling. Our method is based on the decomposition of both amplitudes and two electron integrals, using a combination of tensor hypercontraction and canonical polyadic decomposition. While the original theory scales as O ( N 6 ) with respect to the number of basis functions, we demonstrate numerically that we achieve sub-millihartree difference from the original theory with O ( N 4 ) scaling. This is accomplished by solving directly for the factors that decompose the cluster operator. Here, the proposed scheme is quite general and can be easily extended to other many-body methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Faster Johnson–Lindenstrauss transforms via Kronecker products

The Kronecker product is an important matrix operation with a wide range of applications in signal processing, graph theory, quantum computing and deep learning. In this work, we introduce a generalization of the fast Johnson–Lindenstrauss projection for embedding vectors with Kronecker product structure, the Kronecker fast Johnson–Lindenstrauss transform (KFJLT). The KFJLT reduces the embedding cost by an exponential factor of the standard fast Johnson–Lindenstrauss transform’s cost when applied to vectors with Kronecker structure, by avoiding explicitly forming the full Kronecker products. Here, we prove that this computational gain comes with only a small price in embedding power: consider a finite set of $p$ points in a tensor product of $d$ constituent Euclidean spaces $\bigotimes _{k=d}^{1}{\mathbb{R}}^{n_k}$, and let $N = \prod _{k=1}^{d}n_k$. With high probability, a random KFJLT matrix of dimension $m \times N$ embeds the set of points up to multiplicative distortion $(1\pm \varepsilon )$ provided $m \gtrsim \varepsilon ^{-2} \, \log ^{2d - 1} (p) \, \log N$. We conclude by describing a direct application of the KFJLT to the efficient solution of large-scale Kronecker-structured least squares problems for fitting the CP tensor decomposition.

Kronecker structure↗

NASA NTRSDate not supplied

Assessing the Relationships Between Sensorimotor Biomarkers and Post-Landing Functional Task Performance

Spaceflight drives adaptive changes in healthy individuals appropriate for sensorimotor function in a microgravity environment. These changes are maladaptive for return to earth's gravity. The inter-individual variability of sensorimotor decrements is striking, although poorly understood. The goal of this study was to identify a set of behavioral, neuroimaging and genetic measures that can be used to predict early post-flight performance on a set of sensorimotor tasks. Astronauts were recruited who previously participated in sensorimotor field tests and/or posturography soon after long-duration spaceflight. Behavioral tests included assessments of sensory dependency and adaptability. Visual dependency involved treadmill walking while viewing a moving virtual visual scene. Vestibular perceptual thresholds were measured while seated during lateral translations. Proprioception dependency was measured during one-legged stance on a horizontal air-bearing surface. Ground assessment of adaptability was performed (1) during treadmill walking with a virtual linear hallway and a moving walking surface, and (2) during multiple trials of navigating an obstacle course while wearing reversing prisms (adaptive Functional Mobility Test, aFMT). The neuroimaging tests characterized individual differences in regional brain volumes (using Structural MRI) and white matter microstructure (using Diffusion Tensor Imaging) to serve as potential predictors of adaptive capacity. The genetic tests utilized saliva samples to examine variations in four genes chosen because of their ability to differentiate sensorimotor adaptation ability in a normative population, including Catechol-O-methyltransferase (COMT), Dopamine Receptor D2 (DRD2), Brain-derived neurotrophic factor (BDNF) and the α2-adrenergic receptor. Thirty ISS crewmembers were recruited, including 12 first-time fliers, 6F, and mission durations lasting 185.5 ± 45.5 days, mean ± std. We utilized a combination of three post-flight functional task outcomes: tandem walk, recovery from fall and dynamic posturography. There was considerable variability among all post-flight performance outcomes. Based on an ordinal scale survey, 72% indicated their ability to perform functional tasks were more impacted postflight relative to inflight with 50% indicating they needed to restrict movements for a longer period postflight relative to inflight. While there is a strong association within tests obtained at different R+0 timepoints, by R+24 hr performance on one post-flight test does not necessarily correlate with performance on other post-flight tests. Preliminary statistical analysis indicates combining biomarkers will increase predictive power and this will be explored with future analyses. Our preliminary findings underscore the importance of a comprehensive post-flight test battery including different types of tasks with varying sensory feedback. Understanding the relationships between these sensorimotor biomarkers and post-flight functional task performance improve both our understanding of the individual variability and our strategy to optimize sensorimotor countermeasures.

biomarkers↗

DOE OSTI2015

Similarity and differences between the radion and Higgs boson production and decay processes involving off-shell fermions

The radion is a scalar particle that occurs in brane world models and interacts with the trace of the energy–momentum tensor of the Standard Model (SM). The radion–SM fermion interaction Lagrangian differs from the Higgs boson–fermion interaction Lagrangian for off-shell fermions. It is shown that all additional, as compared to the Higgs boson, contributions to the amplitudes of radion production and decay processes involving off-shell fermions are canceled out for both massless and massive fermions. Thus, additional terms in the interaction Lagrangian do not change properties of these processes for the radion and the Higgs boson, except for the general normalization factors. This similarity is a consequence of gauge invariance for the processes with production of gauge bosons. When an additional scalar particle is produced, there are no apparent reasons for the above cancellation, as confirmed, for example, by the process with production of two scalar particles, which features an additional contribution of the radion in comparison with the Higgs boson.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2019

Determining Gyration Tensor of Orienting Macromolecules through Their Scattering Signature

A method is presented for quantitatively evaluating the shape and size of deformed particles in dispersion from their two-dimensional anisotropic spectra by small-angle scattering. By means of real spherical harmonic expansion, we derive analytical expressions of the gyration tensor R in terms of experimentally measured anisotropic scattering functions, yielding a tensorial extension of the Guinier law. We demonstrate the usefulness of this approach by a model study of an affinely deformed Gaussian chain. We further show that radius of gyration R g is the source term of intraparticle structure factor at the mean-field limit, and from this perspective, we address the connection between R and conformation asphericity. As a result, the developed method not only facilitates quantitative scattering studies of deforming materials, but also provides insightful information regarding their deformation behavior at the molecular level based on the symmetric properties of real spherical harmonics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1974

Geodesic synchrotron radiation in the Kerr geometry by the method of asymptotically factorized Green's functions

The scalar, electromagnetic, and gravitational geodesic-synchrotron-radiation (GSR) spectra are determined for the case of a test particle moving on a highly relativistic circular orbit about a rotating (Kerr) black hole. It is found that the spectral shape depends only weakly on the value of the angular-momentum parameter (a/M) of the black hole, but the total radiated power drops unexpectedly for a value of at least 0.95 and vanishes as the value approaches unity. A spin-dependent factor (involving the inner product of the polarization of a radiated quantum with the source) is isolated to explain the dependence of the spectral shape on the spin of the radiated field. Although the scalar wave equation is solved by separation of variables, this procedure is avoided for the vector and tensor cases by postulating a sum-over-states expansion for the Green's function similar to that found to hold in the scalar case. The terms in this sum, significant for GSR, can then be evaluated in the geometric-optics approximation without requiring the use of vector or tensor spherical harmonics.

Chrzanowski, P. L.↗

NASA NTRSDate not supplied

Assessing the Relationships Between Sensorimotor Biomarkers and Post-Landing Functional Task Performance

Spaceflight drives adaptive changes in healthy individuals appropriate for sensorimotor function in a microgravity environment. These changes are maladaptive for return to earth's gravity. The inter-individual variability of sensorimotor decrements is striking, although poorly understood. The goal of this study is to identify a set of behavioral, neuroimaging and genetic measures that can be used to predict early post-flight performance on a set of sensorimotor tasks. Astronauts are recruited who previously participated in sensorimotor field tests and/or posturography soon after long-duration spaceflight. Behavioral tests include assessments of sensory dependency and adaptability. Visual dependency involves treadmill walking while viewing a moving virtual visual scene. Vestibular thresholds are measured while seated during lateral translations. Proprioception dependency is measured during one-legged stance on a horizontal air-bearing surface. Ground assessment of adaptability is performed (1) during treadmill walking with a virtual linear hallway and a moving walking surface, and (2) during multiple trials of navigating an obstacle course while wearing reversing prisms. The neuroimaging tests will characterize individual differences in regional brain volumes (using Structural MRI) and white matter microstructure (using Diffusion Tensor Imaging) to serve as potential predictors of adaptive capacity. The genetic tests will utilize saliva samples to examine variations in four genes chosen because of their ability to differentiate sensorimotor adaptation ability in a normative population, including Catechol-O-methyltransferase (COMT), Dopamine Receptor D2 (DRD2), Brain-derived neurotrophic factor (BDNF) and the α2-adrenergic receptor. Twenty-one ISS crewmembers have been tested to date, including 6 from this past year after testing resumed post-COVID. This cohort includes 9 first-time fliers, 4F, and mission durations lasting 182 ± 32 days, mean ± std. We are utilizing a combination of three post-flight functional task outcomes: tandem walk, recovery from fall and dynamic posturography. There is considerable variability among the post-flight performance outcomes for the 21 participants to date. Based on a partial sample using an ordinal scale survey, 80% indicated their ability to perform functional tasks were more impacted postflight relative to inflight with 50% indicating they needed to restrict movements for a longer period postflight relative to inflight. While there is a strong association within tests obtained at different R+0 timepoints, by R+24 hr performance on one post-flight test does not necessarily correlate with performance on other post-flight tests. There are apparent relationships between individual measures and specific post-flight outcome measures; however, additional data is needed to draw conclusions. Preliminary statistical analysis indicates combining biomarkers will increase predictive power and this will be explored with future analyses. Our preliminary findings underscore the importance of a comprehensive post-flight test battery including different types of tasks with varying sensory feedback. We expect that understanding the relationships between these sensorimotor biomarkers and post-flight functional task performance will improve both our understanding of the individual variability and our strategy to optimize sensorimotor countermeasures.

S J Wood↗

NASA NTRS2024

Assessing the Relationships Between Sensorimotor Biomarkers and Post-Landing Functional Task Performance

Spaceflight drives adaptive changes in healthy individuals appropriate for sensorimotor function in a microgravity environment. These changes are maladaptive for return to earth's gravity. The inter-individual variability of sensorimotor decrements is striking, although poorly understood. The goal of this study is to identify a set of behavioral, neuroimaging and genetic measures that can be used to predict early post-flight performance on a set of sensorimotor tasks. Astronauts are recruited who previously participated in sensorimotor field tests and/or posturography soon after long-duration spaceflight. Behavioral tests include assessments of sensory dependency and adaptability. Visual dependency involves treadmill walking while viewing a moving virtual visual scene. Vestibular perceptual thresholds are measured while seated during lateral translations. Proprioception dependency is measured during one-legged stance on a horizontal air-bearing surface. Ground assessment of adaptability is performed(1) during treadmill walking with a virtual linear hallway and a moving walking surface, and (2) during multiple trials of navigating an obstacle course while wearing reversing prisms(adaptive Functional Mobility Test, aFMT). The neuroimaging tests will characterize individual differences in regional brain volumes (using Structural MRI) and white matter microstructure(using Diffusion Tensor Imaging) to serve as potential predictors of adaptive capacity. The genetic tests will utilize saliva samples to examine variations in four genes chosen because of their ability to differentiate sensorimotor adaptation ability in a normative population, including Catechol-O-methyltransferase (COMT), Dopamine Receptor D2 (DRD2), Brain-derived neurotrophic factor (BDNF) and the α2-adrenergic receptor. Twenty-seven ISS crewmembers have been tested to date, including 6 from this past year. This cohort includes 10first-time fliers, 6F, and mission durations lasting 178.6 ± 30.5 days, mean ± std. We are utilizing a combination of three post-flight functional task outcomes: tandem walk, recovery from fall and dynamic posturography. There is considerable variability among the post-flight performance outcomes for the 27participants to date. Based on a partial sample using an ordinal scale survey, 70% indicated their ability to perform functional tasks were more impacted postflight relative to inflight with 50% indicating they needed to restrict movements for a longer period postflight relative to inflight. While there is a strong association within tests obtained at different R+0 timepoints, by R+24 hr performance on one post-flight test does not necessarily correlate with performance on other post-flight tests. There are apparent relationships between individual measures and specific post-flight outcome measures, e.g., the cumulative time to complete the aFMT is significantly correlated to pre-to-post-flight changes in tandem walk (rho= 0.64, p = 0.001). Preliminary statistical analysis indicates combining biomarkers will increase predictive power and this will be explored with future analyses. Our preliminary findings underscore the importance of a comprehensive post-flight test battery including different types of tasks with varying sensory feedback. We expect that understanding the relationships between these sensorimotor biomarkers and post-flight functional task performance will improve both our understanding of the individual variability and our strategy to optimize sensorimotor countermeasures

S J Wood↗

DOE OSTI2015

Displacement Fields and Self-Energies of Circular and Polygonal Dislocation Loops in Homogeneous and Layered Anisotropic Solids

There are large classes of materials problems that involve the solutions of stress, displacement, and strain energy of dislocation loops in elastically anisotropic solids, including increasingly detailed investigations of the generation and evolution of irradiation induced defect clusters ranging in sizes from the micro- to meso-scopic length scales. Based on a two-dimensional Fourier transform and Stroh formalism that are ideal for homogeneous and layered anisotropic solids, we have developed robust and computationally efficient methods to calculate the displacement fields for circular and polygonal dislocation loops. Using the homogeneous nature of the Green tensor of order -1, we have shown that the displacement and stress fields of dislocation loops can be obtained by numerical quadrature of a line integral. In addition, it is shown that the sextuple integrals associated with the strain energy of loops can be represented by the product of a pre-factor containing elastic anisotropy effects and a universal term that is singular and equal to that for elastic isotropic case. Furthermore, we have found that the self-energy pre-factor of prismatic loops is identical to the effective modulus of normal contact, and the pre-factor of shear loops differs from the effective indentation modulus in shear by only a few percent. These results provide a convenient method for examining dislocation reaction energetic and efficient procedures for numerical computation of local displacements and stresses of dislocation loops, both of which play integral roles in quantitative defect analyses within combined experimental–theoretical investigations.

36 MATERIALS SCIENCE↗

NASA NTRS1971

Properties of hollow square gravitational wave detectors

Antennas are described with gravitational radiation properties such as polarization, directivity, and scalar and tensor characteristics. These antenna devices are mechanically resonant systems in the shape of a hollow square with fundamental frequencies in the Hz range. An antenna characteristic is that the absorption cross section for quadrupole radiation at the normal frequency of vibration is the same as that of the fundamental if the damping factors are the same. The absorption cross section of quadrupole gravitational radiation is calculated for a class-1 detector of a cylindrical cross section and for a class-2 detector shaped like a hollow square.

Douglass, D. H., Jr.↗

DOE OSTI2016

Dissipative axial inflation

We analyze in detail the background cosmological evolution of a scalar field coupled to a massless abelian gauge field through an axial term φ/ f {sub γ} F ∼ F , such as in the case of an axion. Gauge fields in this case are known to experience tachyonic growth and therefore can backreact on the background as an effective dissipation into radiation energy density ρ{sub R}, which can lead to inflation without the need of a flat potential. We analyze the system, for momenta k smaller than the cutoff f {sub γ}, including the backreaction numerically. We consider the evolution from a given static initial condition and explicitly show that, if f {sub γ} is smaller than the field excursion φ{sub 0} by about a factor of at least O (20), there is a friction effect which turns on before the field can fall down and which can then lead to a very long stage of inflation with a generic potential. In addition we find superimposed oscillations, which would get imprinted on any kind of perturbations, scalars and tensors. Such oscillations have a period of 4–5 efolds and an amplitude which is typically less than a few percent and decreases linearly with f {sub γ}. We also stress that the curvature perturbation on uniform density slices should be sensitive to slow-roll parameters related to ρ{sub R} rather than φ-dot {sup 2}/2 and we discuss the existence of friction terms acting on the perturbations, although we postpone a calculation of the power spectrum and of non-gaussianity to future work and we simply define and compute suitable slow roll parameters. Finally we stress that this scenario may be realized in the axion case, if the coupling 1/ f {sub γ} to U(1) (photons) is much larger than the coupling 1/ f {sub G} to non-abelian gauge fields (gluons), since the latter sets the range of the potential and therefore the maximal allowed φ{sub 0∼} f {sub G}.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Cholesky-based experimental design for Gaussian process and kernel-based emulation and calibration.

Gaussian processes and other kernel-based methods are used extensively to construct approximations of multivariate data sets. The accuracy of these approximations is dependent on the data used. This paper presents a computationally efficient algorithm to greedily select training samples that minimize the weighted L p error of kernel-based approximations for a given number of data. The method successively generates nested samples, with the goal of minimizing the error in high probability regions of densities specified by users. The algorithm presented is extremely simple and can be implemented using existing pivoted Cholesky factorization methods. Training samples are generated in batches which allows training data to be evaluated (labeled) in parallel. For smooth kernels, the algorithm performs comparably with the greedy integrated variance design but has significantly lower complexity. Numerical experiments demonstrate the efficacy of the approach for bounded, unbounded, multi-modal and non-tensor product densities. We also show how to use the proposed algorithm to efficiently generate surrogates for inferring unknown model parameters from data using Bayesian inference.

97 MATHEMATICS AND COMPUTING↗