Commit c7f8b091 authored by Philipp Arras's avatar Philipp Arras
Browse files

Fix Linearized Response

parent ea880962
Pipeline #24050 failed with stage
in 4 minutes and 7 seconds
......@@ -35,7 +35,7 @@ class NonlinearWienerFilterEnergy(Energy):
position_map = FFT.inverse_times(self.power * self.position)
self.LinearizedResponse = \
LinearizedSignalResponse(Instrument, nonlinearity, FFT, power,
position_map)
position_map, sunit)
residual = d - Instrument(sunit * nonlinearity(position_map))
self.N = N
self.S = S
......
......@@ -19,9 +19,9 @@
from ..field import exp
def LinearizedSignalResponse(Instrument, nonlinearity, FFT, power, s):
return (Instrument * nonlinearity.derivative(s) *
FFT.inverse * power)
def LinearizedSignalResponse(Instrument, nonlinearity, FFT, power, s, sunit):
return sunit * (Instrument * nonlinearity.derivative(s) * FFT.inverse * power)
def LinearizedPowerResponse(Instrument, nonlinearity, FFT, Projection, t, m):
power = exp(0.5*t)
......
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