Commit 3daaac3a authored by Martin Reinecke's avatar Martin Reinecke

temporary addition for tests by @kjako

parent bbfb9c1b
Pipeline #23774 passed with stage
in 4 minutes and 34 seconds
......@@ -69,3 +69,20 @@ class WienerFilterCurvature(EndomorphicOperator):
mock_j = self.R.adjoint_times(self.N.inverse_times(mock_data))
mock_m = self.inverse_times(mock_j)
return mock_signal - mock_m
def generate_posterior_sample2(self):
power = self.S.diagonal()
mock_signal = Field.from_random(random_type="normal",
domain=self.S.domain, dtype=noise.dtype)
mock_signal *= sqrt(power)
noise = self.N.diagonal()
mock_noise = Field.from_random(random_type="normal",
domain=self.N.domain, dtype=noise.dtype)
mock_noise *= sqrt(noise)
mock_data = self.R(mock_signal) + mock_noise
mock_j = self.R.adjoint_times(self.N.inverse_times(mock_data))
mock_m = self.inverse_times(mock_j)
return mock_signal - mock_m
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