diff --git a/src/library/correlated_fields_simple.py b/src/library/correlated_fields_simple.py index ea1614888c28813967be1c703fa10fc5907199f3..e1a0560d7cca5e57bef868fa69fa0db8ff43c7fd 100644 --- a/src/library/correlated_fields_simple.py +++ b/src/library/correlated_fields_simple.py @@ -87,13 +87,13 @@ class SimpleCorrelatedField(Operator): xi = ducktape(dom, None, prefix + 'spectrum') smooth = xi*sig_flex*(Adder(shift) @ sig_asp).ptw("sqrt") smooth = _SlopeRemover(pspace, 0) @ twolog @ smooth - op = _Normalization(pspace, 0) @ (slope + smooth) + a = _Normalization(pspace, 0) @ (slope + smooth) maskzm = np.ones(pspace.shape) maskzm[0] = 0 maskzm = makeOp(makeField(pspace, maskzm)) insert = ValueInserter(pspace, (0,)) - a = (maskzm @ ((ps_expander @ fluct)*op)) + insert(zm) + a = (maskzm @ ((ps_expander @ fluct)*a)) + insert(zm) self._a = a.scale(target.total_volume) ht = HarmonicTransformOperator(harmonic_partner, target)