Commit 5f91d098 authored by Reimar Heinrich Leike's avatar Reimar Heinrich Leike

now exporting the linear part of the amplitude operator

parent eac959bd
Pipeline #46361 failed with stages
in 3 minutes and 46 seconds
...@@ -74,7 +74,7 @@ from .minimization.metric_gaussian_kl import MetricGaussianKL ...@@ -74,7 +74,7 @@ from .minimization.metric_gaussian_kl import MetricGaussianKL
from .sugar import * from .sugar import *
from .plot import Plot from .plot import Plot
from .library.smooth_linear_amplitude import SLAmplitude, CepstrumOperator from .library.smooth_linear_amplitude import SLAmplitude, LinearSLAmplitude, CepstrumOperator
from .library.inverse_gamma_operator import InverseGammaOperator from .library.inverse_gamma_operator import InverseGammaOperator
from .library.los_response import LOSResponse from .library.los_response import LOSResponse
from .library.dynamic_operator import (dynamic_operator, from .library.dynamic_operator import (dynamic_operator,
......
...@@ -111,7 +111,6 @@ def CepstrumOperator(target, a, k0): ...@@ -111,7 +111,6 @@ def CepstrumOperator(target, a, k0):
return sym @ qht @ makeOp(cepstrum.sqrt()) return sym @ qht @ makeOp(cepstrum.sqrt())
def SLAmplitude(*, target, n_pix, a, k0, sm, sv, im, iv, keys=['tau', 'phi']): def SLAmplitude(*, target, n_pix, a, k0, sm, sv, im, iv, keys=['tau', 'phi']):
'''Operator for parametrizing smooth amplitudes (square roots of power '''Operator for parametrizing smooth amplitudes (square roots of power
spectra). spectra).
...@@ -169,6 +168,16 @@ def SLAmplitude(*, target, n_pix, a, k0, sm, sv, im, iv, keys=['tau', 'phi']): ...@@ -169,6 +168,16 @@ def SLAmplitude(*, target, n_pix, a, k0, sm, sv, im, iv, keys=['tau', 'phi']):
which returns on its target a power spectrum which consists out of a which returns on its target a power spectrum which consists out of a
smooth and a linear part. smooth and a linear part.
''' '''
return LinearSLAmplitude(target, n_pix, a, k0, sm, sv, im, iv, keys).exp()
def LinearSLAmplitude(*, target, n_pix, a, k0, sm, sv, im, iv, keys=['tau', 'phi']):
'''LinearOperator for parametrizing smooth log-amplitudes (square roots of power
spectra).
Logarithm of SLAmplitude
See documentation of SLAmplitude for more details
'''
if not (isinstance(n_pix, int) and isinstance(target, PowerSpace)): if not (isinstance(n_pix, int) and isinstance(target, PowerSpace)):
raise TypeError raise TypeError
...@@ -196,4 +205,4 @@ def SLAmplitude(*, target, n_pix, a, k0, sm, sv, im, iv, keys=['tau', 'phi']): ...@@ -196,4 +205,4 @@ def SLAmplitude(*, target, n_pix, a, k0, sm, sv, im, iv, keys=['tau', 'phi']):
loglog_ampl = 0.5*(smooth + linear) loglog_ampl = 0.5*(smooth + linear)
# Go from loglog-space to linear-linear-space # Go from loglog-space to linear-linear-space
return et @ loglog_ampl.exp() return et @ loglog_ampl
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment