critical_power_curvature.py

from ...operators.endomorphic_operator import EndomorphicOperator
from ...operators.invertible_operator_mixin import InvertibleOperatorMixin
from ...operators.diagonal_operator import DiagonalOperator


class CriticalPowerCurvature(InvertibleOperatorMixin, EndomorphicOperator):
    """The curvature of the CriticalPowerEnergy.

    This operator implements the second derivative of the
    CriticalPowerEnergy used in some minimization algorithms or
    for error estimates of the powerspectrum.


    Parameters
    ----------
    theta: Field,
        The map and inverse gamma prior contribution to the curvature.
    T : SmoothnessOperator,
        The smoothness prior contribution to the curvature.
    """

    # ---Overwritten properties and methods---

    def __init__(self, theta, T, inverter=None, preconditioner=None, **kwargs):

        self.theta = DiagonalOperator(theta.domain, diagonal=theta)
        self.T = T
        if preconditioner is None:
            preconditioner = self.theta.inverse_times
        self._domain = self.theta.domain
        super(CriticalPowerCurvature, self).__init__(
                                                 inverter=inverter,
                                                 preconditioner=preconditioner,
                                                 **kwargs)

    def _times(self, x, spaces):
        return self.T(x) + self.theta(x)

    # ---Mandatory properties and methods---

    @property
    def domain(self):
        return self._domain

    @property
    def self_adjoint(self):
        return True

    @property
    def unitary(self):
        return False