Computing Reduced Order Models via Inner-Outer Krylov Recycling in Diffuse Optical Tomography
Meghan O'Connell,Misha E. Kilmer,Eric de Sturler,Serkan Gugercin
DOI: https://doi.org/10.1137/16m1062880
IF: 2.968
2017-01-01
SIAM Journal on Scientific Computing
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, de Sturler et al., [SIAM J. Sci. Compute., 37 (2015), pp. B495--B517] addres this bottleneck by using parametrized reduced order models, since the forward problem, typically posed in the frequency domain, is essentially the transfer function of a high-dimensional differential algebraic system. Although this approach drastically reduces the number of large linear systems to be solved, the construction of a candidate global basis for parametrized reduced models still requires the solution of many full order problems, the discretized PDE for multiple right-hand sides, and multiple parameter vectors. This step is followed by a rank-revealing factorization to compress the candidate global basis, in order to produce a much smaller number of vectors suitable for efficient construction of reduced transfer functions. These two steps are still expensive. First, a reduced model must be computed for every reconstruction of a distinct anomaly (say, for every patient), in contrast to applications where a reduced model is computed once and used forever. Moreover, the computation of the reduced model must be efficient on commonly used hardware. 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, and positive definite in the zero frequency case. We apply model order reduction to the new formulation of the problem. Second, we give an approach to computing the global basis approximation directly as the full order systems are solved. This new approach (1) avoids the expensive rank-revealing factorization and (2) substantially reduces both the number of large linear systems to be solved and the number of iterations for each solve, leading to huge savings. Only the incrementally new, relevant information is solved for and added to the existing global basis, avoiding the computation of redundant information. This new approach is achieved by an inner-outer Krylov recycling approach, which has wider potential use in other applications and in contexts beyond reduced order modeling. We show the value of the new approach to approximate global basis computation on two diffuse optical tomography absorption image reconstruction problems.
mathematics, applied