Hierarchical off-diagonal low-rank approximation of Hessians in inverse problems, with application to ice sheet model initialization
From MaRDI portal
Publication:6166000
Abstract: Obtaining lightweight and accurate approximations of Hessian applies in inverse problems governed by partial differential equations (PDEs) is an essential task to make both deterministic and Bayesian statistical large-scale inverse problems computationally tractable. The computational complexity of dense linear algebraic routines such as that needed for sampling from Gaussian proposal distributions and Newton solves by direct linear methods, can be reduced to log-linear complexity by utilizing hierarchical off-diagonal low-rank (HODLR) matrix approximations. In this work, we show that a class of Hessians that arise from inverse problems governed by PDEs are well approximated by the HODLR matrix format. In particular, we study inverse problems governed by PDEs that model the instantaneous viscous flow of ice sheets. In these problems, we seek a spatially distributed basal sliding parameter field such that the flow predicted by the ice sheet model is consistent with ice sheet surface velocity observations. We demonstrate the use of HODLR approximation by efficiently generating Hessian approximations that allow fast generation of samples from a Gaussianized posterior proposal distribution. Computational studies are performed which illustrate ice sheet problem regimes for which the Gauss-Newton data-misfit Hessian is more efficiently approximated by the HODLR matrix format than the low-rank (LR) format. We then demonstrate that HODLR approximations can be favorable, when compared to global low-rank approximations, for large-scale problems by studying the data-misfit Hessian associated to inverse problems governed by the Stokes flow model on the Humboldt glacier and Greenland ice sheets.
Cites work
- scientific article; zbMATH DE number 1350773 (Why is no real title available?)
- scientific article; zbMATH DE number 467276 (Why is no real title available?)
- scientific article; zbMATH DE number 3206496 (Why is no real title available?)
- scientific article; zbMATH DE number 5060482 (Why is no real title available?)
- A computational framework for infinite-dimensional Bayesian inverse problems. I: The linearized case, with application to global seismic inversion
- A computational framework for infinite-dimensional Bayesian inverse problems. II: stochastic Newton MCMC with application to ice sheet flow inverse problems
- A fast randomized algorithm for computing a hierarchically semiseparable representation of a matrix
- A fast randomized eigensolver with structured LDL factorization update
- A sparse matrix arithmetic based on \({\mathfrak H}\)-matrices. I: Introduction to \({\mathfrak H}\)-matrices
- Algorithms for Kullback-Leibler approximation of probability measures in infinite dimensions
- An \(\mathcal O(N\log N)\) fast direct solver for partial hierarchically semi-separable matrices. With application to radial basis function interpolation
- An explicit link between Gaussian fields and Gaussian Markov random fields: the stochastic partial differential equation approach
- Compressing Rank-Structured Matrices via Randomized Sampling
- Data-sparse approximation by adaptive \({\mathcal H}^2\)-matrices
- Fast algorithms for Bayesian uncertainty quantification in large-scale linear inverse problems based on low-rank partial Hessian approximations
- Fast construction of hierarchical matrix representation from matrix-vector multiplication
- Finding structure with randomness: probabilistic algorithms for constructing approximate matrix decompositions
- Hierarchical Matrix Approximations of Hessians Arising in Inverse Problems Governed by PDEs
- Inverse Problem Theory and Methods for Model Parameter Estimation
- Inverse problems: a Bayesian perspective
- Likelihood approximation with hierarchical matrices for large spatial datasets
- Mitigating the influence of the boundary on PDE-based covariance operators
- Monte Carlo sampling methods using Markov chains and their applications
- Nonstationary inverse problems and state estimation
- On a generalization of the preconditioned Crank-Nicolson metropolis algorithm
- Optimal low-rank approximations of Bayesian linear inverse problems
- PyAlbany: a Python interface to the C++ multiphysics solver Albany
- Randomized GPU Algorithms for the Construction of Hierarchical Matrices from Matrix-Vector Operations
- Robust and accurate stopping criteria for adaptive randomized sampling in matrix-free hierarchically semiseparable construction
- Scalable Gaussian Process Computations Using Hierarchical Matrices
- Scalable Matrix-Free Adaptive Product-Convolution Approximation for Locally Translation-Invariant Operators
- 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
Cited in
(3)
This page was built for publication: Hierarchical off-diagonal low-rank approximation of Hessians in inverse problems, with application to ice sheet model initialization
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6166000)
