Search NASA⌕ Search

SEARCH · Search NASA

Results for “matrix diffusion”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

High-order matrix-free incompressible flow solvers with GPU acceleration and low-order refined preconditioners

In this work, we present a matrix-free flow solver for high-order finite element discretizations of the incompressible Navier-Stokes and Stokes equations with GPU acceleration. For high polynomial degrees, assembling the matrix for the linear systems resulting from the finite element discretization can be prohibitively expensive, both in terms of computational complexity and memory. For this reason, it is necessary to develop matrix-free operators and preconditioners, which can be used to efficiently solve these linear systems without access to the matrix entries themselves. The matrix-free operator evaluations utilize GPU-accelerated sum-factorization techniques to minimize memory movement and maximize throughput. The preconditioners developed in this work are based on a low-order refined methodology with parallel subspace corrections, as described for diffusion problems in [1]. The saddle-point Stokes system is solved using block-preconditioning techniques, which are robust in mesh size, polynomial degree, time step, and viscosity. For the incompressible Navier-Stokes equations, we make use of projection (fractional step) methods, which require Helmholtz and Poisson solves at each time step. The performance of our flow solvers is assessed on several benchmark problems in two and three spatial dimensions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Nonlocal transport in bounded two-dimensional systems: An iterative method

Here, the concept of transport mediated through the dynamics of 'jumping' particles is used to develop an iterative method for obtaining steady-state solutions to the nonlocal transport equation in two dimensions (2D). The technique is self-adjoint and capable of correctly treating spatially non-uniform, asymmetric systems. An appropriate reduced version of the iteration method is used to compare with results obtained with a self-adjoint 1D transport matrix approach. The transport 'jump' probability distribution functions (pdfs) are based on Lévy α-stable distributions. The technique can handle the entire Lévy α-parameter range from one (Lorentz distributions) to two (Gaussian distributions). Cases with α = 2 (standard diffusion) are used to establish the validity of the iterative method. The capabilities of the iterative method are demonstrated by presenting examples from systems with various source configurations, boundary shapes, boundary conditions and spatial variations in parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A two-level GPU-accelerated incomplete LU preconditioner for general sparse linear systems

This paper presents a parallel preconditioning approach based on incomplete LU (ILU) factorizations in the framework of Domain Decomposition (DD) for general sparse linear systems. We focus on distributed memory parallel architectures, specifically, those that are equipped with graphic processing units (GPUs). In addition to block-Jacobi, we present general purpose two-level ILU Schur complement-based approaches, where different strategies are presented to solve the coarse-level reduced system. These strategies are combined with modified ILU methods in the construction of the coarse-level operator, in order to effectively remove smooth errors by targeting an algebraically smooth vector. We leverage available GPU-based sparse matrix kernels to accelerate the setup and the solve phases of the proposed ILU preconditioner. We evaluate the efficiency of the proposed methods as a smoother for algebraic multigrid (AMG) and as a preconditioner for Krylov subspace methods on challenging anisotropic diffusion problems and a collection of general sparse matrices.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

CTF Void Drift Validation Study

This milestone report is a summary of work performed in support of expansion of the validation and verification (V&V) matrix for the thermal-hydraulic subchannel code, CTF. The focus of this study is on validating the void drift modeling capabilities of CTF and verifying the supporting models that impact the void drift phenomenon. CTF uses a simple turbulent-diffusion approximation to model lateral cross-flow due to turbulent mixing and void drift. The void drift component of the model is based on the Lahey and Moody model. The models are a function of two-phase mass, momentum, and energy distribution in the system; therefore, it is necessary to correctly model the ow distribution in rod bundle geometry as a first step to correctly calculating the void distribution due to void drift.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Effect of sintering time on structural, microstructural and chemical composition of Ni-doped lanthanum gallate perovskites

This work reports the effect of two different sintering times, 6 and 48 h on the structural, microstructural, and chemical features of Ni-doped La{sub 0.90}Sr{sub 0.10}GaO{sub 3.00−δ}. Independently of the sintering time, La{sub 0.90}Sr{sub 0.10}Ga{sub 1−x}Ni{sub x}O{sub 3.00−δ} (where x=0.10, and 0.20 (mol)) presents a rhombohedral symmetry with a lattice volume that decreases when NiO dopant increases. Besides the perovskite, LaSrGa{sub 3.00}O{sub 7.00} (nominal composition) is present as second phase in all cases. When the samples are doped with NiO, the peaks of this second phase are shifted with respect to the peaks of the pure phase. These shifts suggest that this second phase could admit some Ni ions in its structure. According to the XRD patterns, the amount of the latter phase is larger when sintering time is increased. Electron probe microanalysis (EPMA) indicated that the matrix of the samples sintered for 6 h is constituted by a perovskite with an experimental composition very close to the nominal one. However, when the samples are sintered for 48 h the matrix of each sample is constituted by two perovskites; both with compositional deviations with respect to their nominal one. In particular, a significant Sr depletion compensated by a La increment in the A site is observed. Those compositional deviations could be mainly due to the diffusion of the cations in the bulk and/or from the bulk to the surface of the samples. That diffusion can favour the formation, not only, of a second perovskite with a different composition in relation with the first one formed, but also, the formation of second phases. In addition, a very slight broadening of Bragg peaks of the perovskites sintered for 48 h is observed by XRD and can be related to the presence of two different perovskites in each sample according to EPMA results. By BSEM and EPMA analyses La{sub 4.00}Ga{sub 2.00}O{sub 9.00} (nominal composition) is also observed as second phase when samples are treated for 48 h. - Graphical abstract: Typical microstructures, observed by FEG-SEM, of La{sub 0.90}Sr{sub 0.10}Ga{sub 1−x}Ni{sub x}O{sub 3.00−δ} materials with x=0.00 and 0.10 (nominal compositions) and sintered at 1450 °C for 48 h after polishing and thermal etching are shown. Lighter grey grains with a bimodal grain size are observed and according to XEDS analysis correspond to two perovskites with different experimental compositions (labelled as 1 and 2 for x=0.00 and labelled as 4 and 5 for x=0.10, respectively), meanwhile the darker grey phases correspond to LaSrGa{sub 3.00}O{sub 7.00} (labelled as 3 for x=0.00) and LaSr(Ga{sub 1−y′},Ni{sub y′}){sub 3.00}O{sub 7.00} (labelled as 6 for x=0.10), respectively. - Highlights: • When NiO is introduced into La{sub 0.90}Sr{sub 0.10}GaO{sub 3.00−δ} a change in symmetry is produced. • A significant Sr depletion in the A-site of perovskites treated for 48 h is observed. • Compositional deviations could be due to a thermally-activated diffusional process. • The phase purity of these materials is deteriorated when the sintering time is extended.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2014

Finite range and upper branch effects on itinerant ferromagnetism in repulsive Fermi gases: Bethe–Goldstone ladder resummation approach

We investigate the ferromagnetic transition in repulsive Fermi gases at zero temperature with upper branch and effective range effects. Based on a general effective Lagrangian that reproduces precisely the two-body s-wave scattering phase shift, we obtain a nonperturbative expression of the energy density as a function of the polarization by using the Bethe–Goldstone ladder resummation. For hard sphere potential, the predicted critical gas parameter k{sub F}a=0.816 and the spin susceptibility agree well with the results from fixed-node diffusion Monte Carlo calculations. In general, positive and negative effective ranges have opposite effects on the critical gas parameter k{sub F}a: While a positive effective range reduces the critical gas parameter, a negative effective range increases it. For attractive potential or Feshbach resonance model, the many-body upper branch exhibits an energy maximum at k{sub F}a=α with α=1.34 from the Bethe–Goldstone ladder resummation, which is qualitatively consistent with experimental results. The many-body T-matrix has a positive-energy pole for k{sub F}a>α and it becomes impossible to distinguish the bound state and the scattering state. These positive-energy bound states become occupied and therefore the upper branch reaches an energy maximum at k{sub F}a=α. In the zero range limit, there exists a narrow window (0.86<k{sub F}a<1.56) for the ferromagnetic phase. At sufficiently large negative effective range, the ferromagnetic phase disappears. On the other hand, the appearance of positive-energy bound state resonantly enhances the two-body decay rate around k{sub F}a=α and may prevent the study of equilibrium phases and ferromagnetism of the upper branch Fermi gas. - Highlights: • Nonperturbative interaction energy is obtained within the Bethe–Goldstone ladder resummation approach. • Positive and negative effective ranges have opposite effects on the critical gas parameter. • The upper branch Fermi gas exhibits an energy maximum and reentrant ferromagnetic transition. • The ferromagnetic phase disappears for large and negative effective ranges.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Irradiation performance of U-Mo monolithic fuel

High-performance research reactors require fuel that operates at high specific power to high fission density, but at relatively low temperatures. Research reactor fuels are designed for efficient heat rejection, and are composed of assemblies of thin-plates clad in aluminum alloy. The development of low-enriched fuels to replace high-enriched fuels for these reactors requires a substantially increased uranium density in the fuel to offset the decrease in enrichment. Very few fuel phases have been identified that have the required combination of very-high uranium density and stable fuel behavior at high burnup. U-Mo alloys represent the best known tradeoff in these properties. Testing of aluminum matrix U-Mo aluminum matrix dispersion fuel revealed a pattern of breakaway swelling behavior at intermediate burnup, related to the formation of a molybdenum stabilized high aluminum intermetallic phase that forms during irradiation. In the case of monolithic fuel, this issue was addressed by eliminating, as much as possible, the interfacial area between U-Mo and aluminum. Based on scoping irradiation test data, a fuel plate system composed of solid U-10Mo fuel meat, a zirconium diffusion barrier, and Al6061 cladding was selected for development. Developmental testing of this fuel system indicates that it meets core criteria for fuel qualification, including stable and predictable swelling behavior, mechanical integrity to high burnup, and geometric stability. In addition, the fuel exhibits robust behavior during power-cooling mismatch events under irradiation at high power.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Nanostructured Oxide-Dispersion-Strengthened CoCrFeMnNi High-Entropy Alloys with High Thermal Stability

A nanostructured oxide-dispersion-strengthened (ODS) CoCrFeMnNi high-entropy alloy (HEA) is synthesized by a powder metallurgy process. The thermal stability, including the grain size and crystal structure of the HEA matrix and oxide dispersions, is carefully investigated by X-ray diffraction (XRD) and electron microscopy characterizations after annealing at 900 °C. The limited grain growth may be attributed to Zener pinning of yttria dispersions that impede the grain boundary mobility and diffusivity. The high hardness is caused by both the fine grain size and yttria dispersions, which are also retained after annealing at 900 °C. In this study, it is implied that the combination of ODS and HEA concepts may provide a new design strategy for the development of thermally stable nanostructured alloys for extreme environments.

36 MATERIALS SCIENCE↗
DOE OSTI2016

AGR-3/4 Irradiation Test Train Disassembly and Component Metrology First Look Report

The AGR-3/4 experiment was designed to study fission product transport within graphitic matrix material and nuclear-grade graphite. To this end, this experiment consisted of 12 capsules, each fueled with 4 compacts containing UCO TRISO particles as driver fuel and 20 UCO designed-to-fail (DTF) fuel particles in each compact. The DTF fuel was fabricated with a thin pyrocarbon layer which was intended to fail during irradiation and provide a source of fission products. These fission products could then migrate through the compact and into the surrounding concentric rings of graphitic matrix material and/or nuclear graphite. Through post-irradiation examination (PIE) of the rings (including physical sampling and gamma scanning) fission product concentration profiles within the rings can be determined. These data can be used to elucidate fission product transport parameters (e.g. diffusion coefficients within the test materials) which will be used to inform and refine models of fission product transport. After irradiation in the Advanced Test Reactor (ATR) had been completed in April 2014, the AGR-3/4 experiment was shipped to the Hot Fuel Examination Facility (HFEF) at the Materials and Fuels Complex (MFC) for inspection, disassembly, and metrology. The AGR-3/4 test train was received at MFC in two separate shipments between February and April 2015. Visual examinations of the test train exterior did not indicate dimensional distortion, and only two small discolored areas were observed at the bottom of Capsules 8 and 9. No corresponding discoloration was found on the inside of these capsules, however. Prior to disassembly, the two test train sections were subject to analysis via the Precision Gamma Scanner (PGS), which did not indicate that any gross fuel relocation had occurred. A series of specialized tools (including clamps, cutters, and drills) had been designed and fabricated in order to carry out test train disassembly and recovery of capsule components (graphite rings and fuel compacts). This equipment performed well for separating each capsule in the test train and extracting the capsule components. Only a few problems were encountered. In one case, the outermost ring (the sink ring) was cracked during removal of the capsule through tubes. Although the sink ring will be analyzed in order to obtain a mass balance of fission products in the experiment, these cracks do not pose a major concern because the sink ring will not be analyzed in detail to obtain the spatial distribution of fission products. In Capsules 4 and 5, the compacts could not be removed from the inner rings. Strategies for removing the compacts are being evaluated. Sampling the inner rings with the compacts in-place is also an option. Dimensional measurements were made on the compacts, inner rings, outer rings, and sink rings. The diameters of all compacts decreased by 0.5 to 2.0 %. Generally, the extent of diametric shrinkage increased linearly with increasing neutron fluence. Most compact lengths also decreased. Compact lengths decreased with increasing fluence, reaching maximum shrinkage of about 0.9 % at a fast fluence of 4.0x10 25 n/m 2 E > 0.18 MeV. Above this fluence, the extent of length shrinkage appeared to decrease with fluence, and two compacts from Capsule 7 were found to have slightly increased in length (< 0.1 %) after a fluence of 5.2x10 25 n/m 2 .

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2016

ZnO:Gd nanocrystals for fluorescent applications

Gadolinium doped ZnO crystals within the solubility limit of gadolinium in ZnO matrix were prepared by solid state reaction technique. The method is relatively less expense and enables the production in large scale. The samples were characterised by X-ray diffraction (XRD), energy dispersive X-ray spectroscopy (EDS), scanning electron microscopy (SEM), UV/Vis diffuse reflectance spectroscopy and photoluminescence techniques. Fluorescent property studies of gadolinium doped ZnO at room temperature show enhanced visible light emission due to the defects and oxygen vacancies produced via doping. This work reports the impact of gadolinium doping in the structural, optical and luminescent properties of ZnO in detail.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

A novel closed-form inversion of the convection–diffusion equation for rapid convection, diffusion, and source profile estimation

To simplify and routinize particle transport analysis in fusion devices, a novel closed form linear inversion of the 1-D convection diffusion equation to estimate diffusion and convection profiles D(r ⃗ ), v(r ⃗ ) and source distribution s(r ⃗ ), of a single species from measured data is derived and demonstrated on synthetic data. Profile estimates of D(r ⃗ ), v(r ⃗ ), s(r ⃗ ) and their uncertainties are given as a matrix expression constructed directly from the incoming density data of the transported species in space and time, as well as physics assumptions such as particle conservation and experimental geometry. The derived matrix expression can be applied to a pumped or non-pumped recycling species, or a non-recycling species that is effectively “pumped” by plasma-facing surfaces.

Hinson, Edward [ORNL] (ORCID:000000019713140X)↗
DOE OSTI2025

Robust Implicit Adaptive Low Rank Time-Stepping Methods for Matrix Differential Equations

In this work, we develop implicit rank-adaptive schemes for time-dependent matrix differential equations. The dynamic low rank approximation (DLRA) is a well-known technique to capture the dynamic low rank structure based on Dirac–Frenkel time-dependent variational principle. In recent years, it has attracted a lot of attention due to its wide applicability. Our schemes are inspired by the three-step procedure used in the rank adaptive version of the unconventional robust integrator (the so called BUG integrator) (Ceruti et al. in BIT Numer Math 62(4):1149–1174, 2022) for DLRA. First, a prediction (basis update) step is made computing the approximate column and row spaces at the next time level. Second, a Galerkin evolution step is invoked using an implicit solves for the small core matrix. Finally, a truncation is made according to a prescribed error threshold. Since the DLRA is evolving the differential equation projected on to the tangent space of the low rank manifold, the error estimate of the BUG integrator contains the tangent projection (modeling) error which cannot be easily controlled by mesh refinement. This can cause convergence issue for equations with cross terms. To address this issue, we propose a simple modification, consisting of merging the row and column spaces from the explicit step truncation method together with the BUG spaces in the prediction step. In addition, we propose an adaptive strategy where the BUG spaces are only computed if the residual for the solution obtained from the prediction space by explicit step truncation method, is too large. Here, we prove stability and estimate the local truncation error of the schemes under assumptions. We benchmark the schemes in several tests, such as anisotropic diffusion, solid body rotation and the combination of the two, to show robust convergence properties.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Influence of Alloying Elements and Effect of Stress on Anisotropic Hydrogen Diffusion in Zr-Based Alloys Predicted by Accelerated Kinetic Monte Carlo Simulations

The formation of hydrides degrades the toughness of Zr-based claddings, and thus is one of the primary concerns regarding cladding integrity during fuel operation and used fuel storage. However, the formation mechanism of hydrides, particularly the nucleation mechanism, is yet to be fully understood. To explore the nucleation process, atomic scale studies including density functional theory, molecular dynamics and kinetic Monte Carlo have been used to investigate H diffusion, segregation, and the possible nucleation paths of hydrides. Here, density functional theory calculations are used to parameterize a kinetic Monte Carlo model, which can accurately predict the effective H diffusivity in various Zr-based claddings including Zircaloy, Zirlo and Zr-2.5Nb. Assisted by density functional theory calculations, a recently developed charge-optimized-many-body potential is further improved and used to investigate the possible nucleation paths of hydrides. It’s shown that the formation of face-centered hydrides starts with coherent zeta-phase, shearing of which in the basal plane changes the hcp stacking into fcc and leads to the formation of gamma-hydride nuclei. The shearing involves a negligible barrier and negligible strain accumulation in the matrix by adopting three equivalent shear partials. The results obtained at the atomic scale are consistent with previous theory and experiments, and improves the current understanding on how and how fast hydride forms in Zr based cladding.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2017

SkyFACT: high-dimensional modeling of gamma-ray emission with adaptive templates and penalized likelihoods

We present SkyFACT (Sky Factorization with Adaptive Constrained Templates), a new approach for studying, modeling and decomposing diffuse gamma-ray emission. Like most previous analyses, the approach relies on predictions from cosmic-ray propagation codes like GALPROP and DRAGON. However, in contrast to previous approaches, we account for the fact that models are not perfect and allow for a very large number (∼> 10{sup 5}) of nuisance parameters to parameterize these imperfections. We combine methods of image reconstruction and adaptive spatio-spectral template regression in one coherent hybrid approach. To this end, we use penalized Poisson likelihood regression, with regularization functions that are motivated by the maximum entropy method. We introduce methods to efficiently handle the high dimensionality of the convex optimization problem as well as the associated semi-sparse covariance matrix, using the L-BFGS-B algorithm and Cholesky factorization. We test the method both on synthetic data as well as on gamma-ray emission from the inner Galaxy, |ℓ|<90{sup o} and | b |<20{sup o}, as observed by the Fermi Large Area Telescope. We finally define a simple reference model that removes most of the residual emission from the inner Galaxy, based on conventional diffuse emission components as well as components for the Fermi bubbles, the Fermi Galactic center excess, and extended sources along the Galactic disk. Variants of this reference model can serve as basis for future studies of diffuse emission in and outside the Galactic disk.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Collective Nanoparticle Dynamics Associated with Bridging Network Formation in Model Polymer Nanocomposites

The addition of nanoparticles (NPs) to polymers is a powerful method to improve the mechanical and other properties of macromolecular materials. Such hybrid polymer–particle systems are also rich in fundamental soft matter physics. Among several factors contributing to mechanical reinforcement, a polymer-mediated NP network is considered to be the most important in polymer nanocomposites (PNCs). Here, we present an integrated experimental–theoretical study of the collective NP dynamics in model PNCs using X-ray photon correlation spectroscopy and microscopic statistical mechanics theory. Silica NPs dispersed in unentangled or entangled poly(2-vinylpyridine) matrices over a range of NP loadings are used. Static collective structure factors of the NP subsystems at temperatures above the bulk glass transition temperature reveal the formation of a network-like microstructure via polymer-mediated bridges at high NP loadings above the percolation threshold. The NP collective relaxation times are up to 3 orders of magnitude longer than the self-diffusion limit of isolated NPs and display a rich dependence with observation wavevector and NP loading. A mode-coupling theory dynamical analysis that incorporates the static polymer-mediated bridging structure and collective motions of NPs is performed. It captures well both the observed scattering wavevector and NP loading dependences of the collective NP dynamics in the unentangled polymer matrix, with modest quantitative deviations emerging for the entangled PNC samples. Here, we identify an unusual and weak temperature dependence of collective NP dynamics, in qualitative contrast with the mechanical response. Hence, the present study has revealed key aspects of the collective motions of NPs connected by polymer bridges in contact with a viscous adsorbing polymer medium and identifies some outstanding remaining challenges for the theoretical understanding of these complex soft materials.

36 MATERIALS SCIENCE↗
DOE OSTI2016

High Temperature Fuel Cladding Chemical Interactions Between Unirradiated TRIGA Fuels and 304 Stainless Steel

High temperature fuel cladding chemical interactions (FCCI) between Training-Research, Isotopes-General-Atomics (TRIGA) fuel elements and their 304 stainless steel (304SS) cladding are of interest to develop an understanding of the fuel-behavior during transient reactor scenarios. TRIGA fuels are composed of uranium particles dispersed in a zirconium hydride (Zr H) matrix encased in a 304SS cladding. At high temperatures, the fuel and cladding have the potential experience FCCI. A large number of interactions can take place in this multi-component fuel system. Emphasis is given to the development of low melting point phases. Diffusion couples of TRIGA fuel vs. 304SS were annealed at 1000°C for 5 and 24 hours. Characterization was then carried out via scanning electron microscopy and transmission electron microscopy. The as-received and annealed microstructures of the fuel and the interactions regions between the fuel and the 304SS will be discussed.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2018

Iron Porphyrins Embedded into a Supramolecular Porous Organic Cage for Electrochemical CO 2 Reduction in Water

Abstract A porous organic cage composed of six iron tetraphenylporphyrins was used as a supramolecular catalyst for electrochemical CO 2 ‐to‐CO conversion. This strategy enhances active site exposure and substrate diffusion relative to the monomeric catalyst, resulting in CO generation with near‐quantitative Faradaic efficiency in pH 7.3 water, with activities reaching 55 250 turnovers. These results provide a starting point for the design of supramolecular catalysts that can exploit the properties of the surrounding matrix yet retain the tunability of the original molecular unit.

Smith, Peter T.↗
DOE OSTI2022

PyAlbany: A Python interface to the C++ multiphysics solver Albany

Albany is a parallel C++ finite element library for solving forward and inverse problems involving partial differential equations (PDEs). In this paper we introduce PyAlbany, a newly developed Python interface to the Albany library. PyAlbany can be used to effectively drive Albany enabling fast and easy analysis and post-processing of applications based on PDEs that are pre-implemented in Albany. PyAlbany relies on the library PyBind11 to bind Python with C++ Albany code. Here we detail the implementation of PyAlbany and showcase its capabilities through a number of examples targeting a heat-diffusion problem. In particular we consider the following: (1) the generation of samples for a Monte Carlo application, (2) a scalability study, (3) a study of parameters on the performance of a linear solver, and finally (4) a tool for performing eigenvalue decompositions of matrix-free operators for a Bayesian inference application.

97 MATHEMATICS AND COMPUTING↗