Commit eaaed9cc authored by Lukas Platz's avatar Lukas Platz

cosmetics

parent 2943d105
......@@ -192,17 +192,17 @@ def LinearSLAmplitude(*, target, n_pix, a, k0, sm, sv, im, iv,
et = ExpTransform(target, n_pix)
dom = et.domain[0]
# Smooth component
dct = {'a': a, 'k0': k0}
smooth = CepstrumOperator(dom, **dct).ducktape(keys[0])
# Linear component
sl = SlopeOperator(dom)
mean = np.array([sm, im + sm*dom.t_0[0]])
sig = np.array([sv, iv])
mean = Field.from_global_data(sl.domain, mean)
sig = Field.from_global_data(sl.domain, sig)
linear = sl @ Adder(mean) @ makeOp(sig).ducktape(keys[1])
sigma = np.array([sv, iv])
sigma = Field.from_global_data(sl.domain, sigma)
linear = sl @ Adder(mean) @ makeOp(sigma).ducktape(keys[1])
# Smooth deviations from linear component
dct = {'a': a, 'k0': k0}
smooth = CepstrumOperator(dom, **dct).ducktape(keys[0])
# Combine linear and smooth component
loglog_ampl = 0.5*(smooth + linear)
......
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