Computing reduced order models via inner-outer Krylov recycling in diffuse optical tomography
From MaRDI portal
Abstract: In nonlinear imaging problems whose forward model is described by a partial differential equation (PDE), the main computational bottleneck in solving the inverse problem is the need to solve many large-scale discretized PDEs at each step of the optimization process. In the context of absorption imaging in diffuse optical tomography, one approach to addressing this bottleneck proposed recently (de Sturler, et al, 2015) reformulates the viewing of the forward problem as a differential algebraic system, and then employs model order reduction (MOR). However, the construction of the reduced model requires the solution of several full order problems (i.e. the full discretized PDE for multiple right-hand sides) to generate a candidate global basis. This step is then followed by a rank-revealing factorization of the matrix containing the candidate basis in order to compress the basis to a size suitable for constructing the reduced transfer function. The present paper addresses the costs associated with the global basis approximation in two ways. First, we use the structure of the matrix to rewrite the full order transfer function, and corresponding derivatives, such that the full order systems to be solved are symmetric (positive definite in the zero frequency case). Then we apply MOR to the new formulation of the problem. Second, we give an approach to computing the global basis approximation dynamically as the full order systems are solved. In this phase, only the incrementally new, relevant information is added to the existing global basis, and redundant information is not computed. This new approach is achieved by an inner-outer Krylov recycling approach which has potential use in other applications as well. We show the value of the new approach to approximate global basis computation on two DOT absorption image reconstruction problems.
Recommendations
- Nonlinear parametric inversion using interpolatory model reduction
- A model reduction approach for inverse problems with operator valued data
- Fast Algorithms for Hyperspectral Diffuse Optical Tomography
- Reduced-order methods in medical imaging
- A reduced basis Landweber method for nonlinear inverse problems
Cites work
- A flexible Krylov solver for shifted systems with application to oscillatory hydraulic tomography
- A model reduction approach to numerical inversion for a parabolic partial differential equation
- A reduced basis Landweber method for nonlinear inverse problems
- A regularized Gauss-Newton trust region approach to imaging in diffuse optical tomography
- A survey of projection-based model reduction methods for parametric dynamical systems
- Accelerating optimization of parametric linear systems by model order reduction
- Computational Methods for Inverse Problems
- Domain decomposition and balanced truncation model reduction for shape optimization of the Stokes system
- Extending the eigCG algorithm to nonsymmetric Lanczos for linear systems with multiple right-hand sides.
- Fast Algorithms for Hyperspectral Diffuse Optical Tomography
- Hessian-based model reduction for large-scale systems with initial-condition inputs
- scientific article; zbMATH DE number 1226426 (Why is no real title available?)
- scientific article; zbMATH DE number 2222860 (Why is no real title available?)
- Inexact solves in interpolatory model reduction
- Interpolatory projection methods for parameterized model reduction
- Large-scale topology optimization using preconditioned Krylov subspace methods with recycling
- Model Reduction for Large-Scale Systems with High-Dimensional Parametric Input Space
- Model reduction of descriptor systems by interpolatory projection methods
- Nested Krylov methods based on GCR
- Nonlinear parametric inversion using interpolatory model reduction
- On optimization techniques for solving nonlinear inverse problems
- Optical tomography in medical imaging
- Parameter and state model reduction for large-scale statistical inverse problems
- Parametric level set methods for inverse problems
- Proper orthogonal decomposition for optimality systems
- Recycling BiCGSTAB with an application to parametric model order reduction
- Recycling Krylov subspaces for efficient large-scale electrical impedance tomography
- Recycling Subspace Information for Diffuse Optical Tomography
- Reduced basis method for finite volume approximations of parametrized linear evolution equations
- Reduced order modeling based shape optimization of surface acoustic wave driven microfluidic biochips
- Solution of the time-domain inverse resistivity problem in the model reduction framework. I. One-dimensional problem with SISO data
- Subspace recycling accelerates the parametric macro-modeling of MEMS
Cited in
(9)- A multiscale method for model order reduction in PDE parameter estimation
- A model reduction approach for inverse problems with operator valued data
- Randomized approaches to accelerate MCMC algorithms for Bayesian inverse problems
- Finite element model updating for structural applications
- Compressing Large-Scale Wave Propagation Models via Phase-Preconditioned Rational Krylov Subspaces
- Reduced-order methods in medical imaging
- Nonlinear parametric inversion using interpolatory model reduction
- \( \mathcal{H}_2\)-optimal model reduction using projected nonlinear least squares
- A survey of subspace recycling iterative methods
This page was built for publication: Computing reduced order models via inner-outer Krylov recycling in diffuse optical tomography
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5738152)
