Randomized approach to nonlinear inversion combining random and optimized simultaneous sources and detectors
From MaRDI portal
(Redirected from Publication:4631407)
Stochastic approximation (62L20) Ill-posedness and regularization problems in numerical linear algebra (65F22) Numerical methods for inverse problems for initial value and initial-boundary value problems involving PDEs (65M32) Numerical methods for inverse problems for boundary value problems involving PDEs (65N21) Numerical solution of discretized equations for boundary value problems involving PDEs (65N22) Stochastic programming (90C15)
Abstract: In partial differential equations-based (PDE-based) inverse problems with many measurements, many large-scale discretized PDEs must be solved for each evaluation of the misfit or objective function. In the nonlinear case, evaluating the Jacobian requires solving an additional set of systems. This leads to a tremendous computational cost, and this is by far the dominant cost for these problems. Several authors have proposed randomization and stochastic programming techniques to drastically reduce the number of system solves by estimating the objective function using only a few appropriately chosen random linear combinations of the sources. While some have reported good solution quality at a greatly reduced cost, for our problem of interest, diffuse optical tomography, the approach often does not lead to sufficiently accurate solutions. We propose two improvements. First, to efficiently exploit Newton-type methods, we modify the stochastic estimates to include random linear combinations of detectors, drastically reducing the number of adjoint solves. Second, after solving to a modest tolerance, we compute a few simultaneous sources and detectors that maximize the Frobenius norm of the sampled Jacobian to improve the rate of convergence and obtain more accurate solutions. We complement these optimized simultaneous sources and detectors by random simultaneous sources and detectors constrained to a complementary subspace. Our approach leads to solutions of the same quality as obtained using all sources and detectors but at a greatly reduced computational cost, as the number of large-scale linear systems to be solved is significantly reduced.
Recommendations
- Stochastic algorithms for inverse problems involving PDEs and many measurements
- A general framework for nonlinear multigrid inversion
- On unifying randomized methods for inverse problems
- Data completion and stochastic algorithms for PDE inversion problems with many measurements
- Randomize-then-optimize: a method for sampling from posterior distributions in nonlinear inverse problems
Cites work
- A Multilinear Singular Value Decomposition
- A regularized Gauss-Newton trust region approach to imaging in diffuse optical tomography
- A Stochastic Estimator of the Trace of the Influence Matrix for Laplacian Smoothing Splines
- Adaptive sampling strategies for stochastic optimization
- An effective method for parameter estimation with PDE constraints with multiple right-hand sides
- Computational Methods for Inverse Problems
- Finding structure with randomness: probabilistic algorithms for constructing approximate matrix decompositions
- scientific article; zbMATH DE number 194139 (Why is no real title available?)
- scientific article; zbMATH DE number 1460605 (Why is no real title available?)
- scientific article; zbMATH DE number 852532 (Why is no real title available?)
- Lectures on Stochastic Programming
- Nonlinear parametric inversion using interpolatory model reduction
- On the Best Rank-1 and Rank-(R1 ,R2 ,. . .,RN) Approximation of Higher-Order Tensors
- Optical tomography in medical imaging
- Parametric level set methods for inverse problems
- Principal component analysis of three-mode data by means of alternating least squares algorithms
- Recycling Subspace Information for Diffuse Optical Tomography
- Sample size selection in optimization methods for machine learning
- Stochastic algorithms for inverse problems involving PDEs and many measurements
- Sub-sampled Newton methods
- Tensor Decompositions and Applications
Cited in
(7)- Randomized approaches to accelerate MCMC algorithms for Bayesian inverse problems
- Data completion and stochastic algorithms for PDE inversion problems with many measurements
- Stochastic algorithms for inverse problems involving PDEs and many measurements
- Simultaneous-shot inversion for PDE-constrained optimization problems with missing data
- Randomized Discrete Empirical Interpolation Method for Nonlinear Model Reduction
- A-optimal encoding weights for nonlinear inverse problems, with application to the Helmholtz inverse problem
- Parametric level-sets enhanced to improve reconstruction (PaLEnTIR)
This page was built for publication: Randomized approach to nonlinear inversion combining random and optimized simultaneous sources and detectors
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q4631407)
