Commit e6956b5a authored by Gordian Edenhofer's avatar Gordian Edenhofer
Browse files

correlated_fields.py: Resolve FIXMEs

Drop one and rename the other two.
parent 9d8d990b
Pipeline #92251 passed with stages
in 11 minutes and 7 seconds
...@@ -170,7 +170,7 @@ class _Normalization(Operator): ...@@ -170,7 +170,7 @@ class _Normalization(Operator):
def apply(self, x): def apply(self, x):
self._check_input(x) self._check_input(x)
spec = x.ptw("exp") spec = x.ptw("exp")
# FIXME This normalizes also the zeromode which is supposed to be left # NOTE, this normalizes also the zeromode which is supposed to be left
# untouched by this operator # untouched by this operator
return (self._specsum(spec).reciprocal()*spec).sqrt() return (self._specsum(spec).reciprocal()*spec).sqrt()
...@@ -209,7 +209,7 @@ class _AmplitudeMatern(Operator): ...@@ -209,7 +209,7 @@ class _AmplitudeMatern(Operator):
expander = VdotOperator(full(pow_spc, 1.)).adjoint expander = VdotOperator(full(pow_spc, 1.)).adjoint
k_squared = makeField(pow_spc, pow_spc.k_lengths**2) k_squared = makeField(pow_spc, pow_spc.k_lengths**2)
a = expander @ scale.log() # FIXME: look for nicer implementation a = expander @ scale.log() # TODO: look for nicer implementation
b = VdotOperator(k_squared).adjoint @ cutoff.power(-2.) b = VdotOperator(k_squared).adjoint @ cutoff.power(-2.)
c = expander.scale(-1) @ logloghalfslope c = expander.scale(-1) @ logloghalfslope
...@@ -670,7 +670,6 @@ class CorrelatedFieldMaker: ...@@ -670,7 +670,6 @@ class CorrelatedFieldMaker:
else: else:
offset = float(offset) offset = float(offset)
op = Adder(full(op.target, offset)) @ op op = Adder(full(op.target, offset)) @ op
#FIXME why does prior_info no longer works???
self.statistics_summary(prior_info) self.statistics_summary(prior_info)
return op return op
......
Supports Markdown
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