3,586 research outputs found
A Two-stage Method for Inverse Medium Scattering
We present a novel numerical method to the time-harmonic inverse medium
scattering problem of recovering the refractive index from near-field scattered
data. The approach consists of two stages, one pruning step of detecting the
scatterer support, and one resolution enhancing step with mixed regularization.
The first step is strictly direct and of sampling type, and faithfully detects
the scatterer support. The second step is an innovative application of
nonsmooth mixed regularization, and it accurately resolves the scatterer sizes
as well as intensities. The model is efficiently solved by a semi-smooth
Newton-type method. Numerical results for two- and three-dimensional examples
indicate that the approach is accurate, computationally efficient, and robust
with respect to data noise.Comment: 18 pages, 5 figure
Generalized Newton's Method based on Graphical Derivatives
This paper concerns developing a numerical method of the Newton type to solve
systems of nonlinear equations described by nonsmooth continuous functions. We
propose and justify a new generalized Newton algorithm based on graphical
derivatives, which have never been used to derive a Newton-type method for
solving nonsmooth equations. Based on advanced techniques of variational
analysis and generalized differentiation, we establish the well-posedness of
the algorithm, its local superlinear convergence, and its global convergence of
the Kantorovich type. Our convergence results hold with no semismoothness
assumption, which is illustrated by examples. The algorithm and main results
obtained in the paper are compared with well-recognized semismooth and
-differentiable versions of Newton's method for nonsmooth Lipschitzian
equations
Dynamical systems and forward-backward algorithms associated with the sum of a convex subdifferential and a monotone cocoercive operator
In a Hilbert framework, we introduce continuous and discrete dynamical
systems which aim at solving inclusions governed by structured monotone
operators , where is the subdifferential of a
convex lower semicontinuous function , and is a monotone cocoercive
operator. We first consider the extension to this setting of the regularized
Newton dynamic with two potentials. Then, we revisit some related dynamical
systems, namely the semigroup of contractions generated by , and the
continuous gradient projection dynamic. By a Lyapunov analysis, we show the
convergence properties of the orbits of these systems.
The time discretization of these dynamics gives various forward-backward
splitting methods (some new) for solving structured monotone inclusions
involving non-potential terms. The convergence of these algorithms is obtained
under classical step size limitation. Perspectives are given in the field of
numerical splitting methods for optimization, and multi-criteria decision
processes.Comment: 25 page
Nonlinear Preconditioning: How to use a Nonlinear Schwarz Method to Precondition Newton's Method
For linear problems, domain decomposition methods can be used directly as
iterative solvers, but also as preconditioners for Krylov methods. In practice,
Krylov acceleration is almost always used, since the Krylov method finds a much
better residual polynomial than the stationary iteration, and thus converges
much faster. We show in this paper that also for non-linear problems, domain
decomposition methods can either be used directly as iterative solvers, or one
can use them as preconditioners for Newton's method. For the concrete case of
the parallel Schwarz method, we show that we obtain a preconditioner we call
RASPEN (Restricted Additive Schwarz Preconditioned Exact Newton) which is
similar to ASPIN (Additive Schwarz Preconditioned Inexact Newton), but with all
components directly defined by the iterative method. This has the advantage
that RASPEN already converges when used as an iterative solver, in contrast to
ASPIN, and we thus get a substantially better preconditioner for Newton's
method. The iterative construction also allows us to naturally define a coarse
correction using the multigrid full approximation scheme, which leads to a
convergent two level non-linear iterative domain decomposition method and a two
level RASPEN non-linear preconditioner. We illustrate our findings with
numerical results on the Forchheimer equation and a non-linear diffusion
problem
- …