Hierarchical Matrix Approximations of Hessians Arising in Inverse Problems Governed by PDEs
From MaRDI portal
Publication:5132022
DOI10.1137/19M1270367zbMath1453.65289arXiv2003.10173OpenAlexW3093797476WikidataQ114074247 ScholiaQ114074247MaRDI QIDQ5132022
Georg Stadler, Stefano Zampini, Tan Bui-Thanh, Wajih Boukaram, Ilona Ambartsumyan, Omar Ghattas, George M. Turkiyyah, David E. Keyes
Publication date: 9 November 2020
Published in: SIAM Journal on Scientific Computing (Search for Journal in Brave)
Full work available at URL: https://arxiv.org/abs/2003.10173
inverse problemsNewton methodsmatrix compressionhierarchical matricesPDE-constrained optimizationGPUHessianslog-linear complexitylow rank updatesNewton-Schulz
Iterative numerical methods for linear systems (65F10) Parallel numerical computation (65Y05) Discrete approximations in optimal control (49M25) Numerical methods for inverse problems for initial value and initial-boundary value problems involving PDEs (65M32) Inverse problems (including inverse scattering) in optics and electromagnetic theory (78A46) Inverse problems in optimal control (49N45) Numerical linear algebra (65F99) Numerical algorithms for specific classes of architectures (65Y10) PDEs in connection with control and optimization (35Q93)
Related Items
An Offline-Online Decomposition Method for Efficient Linear Bayesian Goal-Oriented Optimal Experimental Design: Application to Optimal Sensor Placement ⋮ Learning physics-based models from data: perspectives from inverse problems and model reduction ⋮ H2Opus: a distributed-memory multi-GPU software package for non-local operators ⋮ A Fast and Scalable Computational Framework for Large-Scale High-Dimensional Bayesian Optimal Experimental Design ⋮ Sparse Recovery of Elliptic Solvers from Matrix-Vector Products ⋮ Hierarchical off-diagonal low-rank approximation of Hessians in inverse problems, with application to ice sheet model initialization ⋮ Enabling Hyper-Differential Sensitivity Analysis for Ill-Posed Inverse Problems ⋮ Scalable Physics-Based Maximum Likelihood Estimation Using Hierarchical Matrices ⋮ KSPHPDDM and PCHPDDM: extending PETSc with advanced Krylov methods and robust multilevel overlapping Schwarz preconditioners
Uses Software
Cites Work
- Finding structure with randomness: Probabilistic algorithms for constructing approximate matrix decompositions
- Tensor Decompositions and Applications
- Matrix probing: a randomized preconditioner for the wave-equation Hessian
- On Bayesian A- and D-optimal experimental designs in infinite dimensions
- Fast construction of hierarchical matrix representation from matrix-vector multiplication
- Automated solution of differential equations by the finite element method. The FEniCS book
- Scalable and efficient algorithms for the propagation of uncertainty from data through inference to prediction for large-scale problems, with application to flow of the antarctic ice sheet
- Data-sparse approximation of non-local operators by \(\mathcal H^2\)-matrices
- Approximate iterations for structured matrices
- Approximation of solution operators of elliptic partial differential equations by \({\mathcal H}\)- and \({\mathcal H}^2\)-matrices
- A perfectly matched layer for the absorption of electromagnetic waves
- Existence of \(\mathcal H\)-matrix approximants to the inverse FE-matrix of elliptic operators with \(L^\infty\)-coefficients
- Construction and arithmetics of \(\mathcal H\)-matrices
- Geometric MCMC for infinite-dimensional inverse problems
- An efficient computation of generalized inverse of a matrix
- MFEM: a modular finite element methods library
- Dimension-independent likelihood-informed MCMC
- Analysis of the Hessian for inverse scattering problems. III: Inverse medium scattering of electromagnetic waves in three dimensions.
- An \(\mathcal O(N\log N)\) fast direct solver for partial hierarchically semi-separable matrices. With application to radial basis function interpolation
- Low-Rank Approximation of Elliptic Boundary Value Problems with High-Contrast Coefficients
- PCBDDC: A Class of Robust Dual-Primal Methods in PETSc
- Analysis of the Hessian for inverse scattering problems: I. Inverse shape scattering of acoustic waves
- A Stochastic Newton MCMC Method for Large-Scale Statistical Inverse Problems with Application to Seismic Inversion
- A Computational Framework for Infinite-Dimensional Bayesian Inverse Problems, Part II: Stochastic Newton MCMC with Application to Ice Sheet Flow Inverse Problems
- Solving large-scale PDE-constrained Bayesian inverse problems with Riemann manifold Hamiltonian Monte Carlo
- Hierarchical Matrices: Algorithms and Analysis
- Fast algorithms for hierarchically semiseparable matrices
- Fast Algorithms for Bayesian Uncertainty Quantification in Large-Scale Linear Inverse Problems Based on Low-Rank Partial Hessian Approximations
- Computational Optimization of Systems Governed by Partial Differential Equations
- A Fast Randomized Algorithm for Computing a Hierarchically Semiseparable Representation of a Matrix
- An Improved Newton Iteration for the Generalized Inverse of a Matrix, with Applications
- A Newton-CG method for large-scale three-dimensional elastic full-waveform seismic inversion
- Discontinuous Galerkin Methods for Solving Elliptic and Parabolic Equations
- Optimization with PDE Constraints
- Discontinuous Galerkin Methods for Advection-Diffusion-Reaction Problems
- Hessian-based model reduction for large-scale systems with initial-condition inputs
- A Nonlinear Primal-Dual Method for Total Variation-Based Image Restoration
- Unified Analysis of Discontinuous Galerkin Methods for Elliptic Problems
- Directional ‐matrix compression for high‐frequency problems
- A comparative study of structural similarity and regularization for joint inverse problems governed by PDEs
- Mesh Independence for Nonlinear Least Squares Problems with Norm Constraints
- Accuracy and Stability of Numerical Algorithms
- Choosing the Forcing Terms in an Inexact Newton Method
- Balancing Domain Decomposition by Constraints Algorithms for Curl-Conforming Spaces of Arbitrary Order
- Randomized GPU Algorithms for the Construction of Hierarchical Matrices from Matrix-Vector Operations
- Scalable Matrix-Free Adaptive Product-Convolution Approximation for Locally Translation-Invariant Operators
- Hierarchical Matrix Operations on GPUs
- Adaptive BDDC Deluxe Methods for H(curl)
- A Data Scalable Augmented Lagrangian KKT Preconditioner for Large-Scale Inverse Problems
- A Computational Framework for Infinite-Dimensional Bayesian Inverse Problems Part I: The Linearized Case, with Application to Global Seismic Inversion
- On the Numerical Properties of an Iterative Method for Computing the Moore–Penrose Generalized Inverse
- Compressing Rank-Structured Matrices via Randomized Sampling
- Unnamed Item
