LargeScale Optimization Techniques for the Regularization of IllPosed
Problems
Marielba Rojas
Department of Mathematics, Wake Forest University
Thursday, May 1, 2003 15:0016:00, Room 145, NIST North (820) Gaithersburg Thursday, May 1, 2003 13:0014:00, Room 4550 Boulder
Abstract:
We will discuss the use of optimization techniques, in particular,
trustregion techniques,
in the solution of largescale inverse problems. We will describe the
matrixfree method
LSTRS for the largescale trustregion subproblem (TRS) of minimizing a
quadratic functional
subject to a norm constraint. The method is based on a reformulation of
the TRS as a parameterized
eigenvalue problem. The strategy consists of an iterative procedure that
drives the parameter
toward its optimal value. The solution to the TRS is then recovered from
the solution of the
eigenvalue problem corresponding to the optimal parameter. A largescale
eigenvalue
problem must be solved at each iteration. This is accomplished by means
of the Implicitly
Restarted Lanczos Method.
We will describe the method, discuss the issues associated with
illposed problems, and
present numerical results on largescale inverse problems. The results
were obtained with
a MATLAB implementation of LSTRS. A MATLAB 6 version of LSTRS will be
publicly
available shortly. The problems discussed include largescale seismic
inversion problems
with field data.
Contact: A. J. KearsleyNote: Visitors from outside NIST must contact
Robin Bickel; (301) 9753668;
at least 24 hours in advance.
