Files in this item



application/pdfRecycling Subsp ... use Optical Tomography.pdf (1MB)
(no description provided)PDF


Title:Recycling Subspace Information for Diffuse Optical Tomography
Author(s):Kilmer, Misha; de Sturler, Eric
Subject(s):Scientific Computing
Abstract:We discuss the efficient solution of a large sequence of slowly varying linear systems arising in computations for diffuse optical tomographic imaging. In particular, we analyze a number of strategies for recycling Krylov subspace information for the most efficient solution. We reconstruct three-dimensional absorption and scattering information by matching computed solutions from a parameterized model to measured data. For this nonlinear least squares problem we use the Gauss-Newton method with a line search. This algorithm requires the solution of a large sequence of linear systems. Each choice of parameters in the nonlinear least squares algorithm results in a different matrix describing the optical properties of the medium. These matrices change slowly from one step to the next, but may change significantly over many steps. For each matrix we must solve a set of linear systems involving both multiple shifts and multiple right-hand sides. We discuss strategies that minimize the overall solution time. In particular, we show how we can tune the linear solver for both the nonlinear optimization algorithm and the underlying application. Although we focus on a particular application and optimization algorithm, we feel that our approach is applicable generally to problems where many linear systems must be solved. We describe extensions to the GCRO algorithm to deal efficiently with symmetric problems and to combine subspace recycling with solving for multiple shifts using a single Krylov subspace. We provide results for two sets of numerical experiments to demonstrate the effectiveness of the resulting method.
Issue Date:2004-06
Genre:Technical Report
Rights Information:You are granted permission for the non-commercial reproduction, distribution, display, and performance of this technical report in any format, BUT this permission is only for a period of 45 (forty-five) days from the most recent time that you verified that this technical report is still available from the University of Illinois at Urbana-Champaign Computer Science Department under terms that include this permission. All other rights are reserved by the author(s).
Date Available in IDEALS:2009-04-14

This item appears in the following Collection(s)

Item Statistics