Module regpy.solvers.nonlinear.landweber

Classes

class Landweber (setting, data, init, stepsize=None, op_norm_method='lanczos')
Expand source code 
class Landweber(RegSolver):
    r"""The Landweber method. Solves the potentially non-linear, ill-posed equation
    \[
        F(x) = g^\delta,
    \]
    where \(T)\ is a Frechet-differentiable operator, by gradient descent for the residual
    \[
        \Vert F(x) - g^\delta\Vert^2,
    \]
    where \(\Vert\cdot\Vert)\ is the Hilbert space norm in the codomain, and gradients are computed with
    respect to the Hilbert space structure on the domain.

    The number of iterations is effectively the regularization parameter and needs to be picked
    carefully.

    Parameters
    ----------
    setting : regpy.solvers.RegularizationSetting
        The setting of the forward problem.
    data : array-like
        The measured data/right hand side.
    init : array-like
        The initial guess.
    stepsize : float, optional
        The step length; must be chosen not too large. If omitted, it is guessed from the norm of
        the derivative at the initial guess.
    """

    def __init__(self, setting, data, init, stepsize=None, op_norm_method = "lanczos"):
        super().__init__(setting)
        self.rhs = data
        """The right hand side gets initialized with the measured data."""
        self.x = init
        self.y, deriv = self.op.linearize(self.x)
        self.deriv = deriv
        """The derivative at the current iterate."""

        self.stepsize = stepsize or 0.9 / setting.op_norm(op=self.deriv, method = op_norm_method)**2
        """The stepsize."""

    def _next(self):
        self._residual = self.y - self.rhs
        self._gy_residual = self.h_codomain.gram(self._residual)
        self._update = self.deriv.adjoint(self._gy_residual)
        self.x -= self.stepsize * self.h_domain.gram_inv(self._update)
        self.y, self.deriv = self.op.linearize(self.x)

        if self.log.isEnabledFor(logging.INFO):
            norm_residual = sqrt((self.op.codomain.vdot(self._residual, self._gy_residual)).real)
            self.log.info('|residual| = {}'.format(norm_residual))

The Landweber method. Solves the potentially non-linear, ill-posed equation F(x) = g^\delta, where (T)\ is a Frechet-differentiable operator, by gradient descent for the residual \Vert F(x) - g^\delta\Vert^2, where (\Vert\cdot\Vert)\ is the Hilbert space norm in the codomain, and gradients are computed with respect to the Hilbert space structure on the domain.

The number of iterations is effectively the regularization parameter and needs to be picked carefully.

Parameters

setting : RegularizationSetting
The setting of the forward problem.
data : array-like
The measured data/right hand side.
init : array-like
The initial guess.
stepsize : float, optional
The step length; must be chosen not too large. If omitted, it is guessed from the norm of the derivative at the initial guess.

Ancestors

Instance variables

var rhs

The right hand side gets initialized with the measured data.

var deriv

The derivative at the current iterate.

var stepsize

The stepsize.

Inherited members