From aab304b5646b8c0a5e6acf1c08d7ddab41de71fd Mon Sep 17 00:00:00 2001
From: Philipp Arras <parras@mpa-garching.mpg.de>
Date: Mon, 16 Aug 2021 12:22:13 +0200
Subject: [PATCH] 12/n
---
teaser_critical_filter.py | 6 ++----
1 file changed, 2 insertions(+), 4 deletions(-)
diff --git a/teaser_critical_filter.py b/teaser_critical_filter.py
index 1c13b77..0b57738 100644
--- a/teaser_critical_filter.py
+++ b/teaser_critical_filter.py
@@ -96,10 +96,8 @@ plot_WF('teaser_unknown_power', ground_truth, data, m=mean,
samples=posterior_samples)
# Plot the reconstruction of the power spectrum
-mysterious_spectrum = lambda k: 5/((7**2 - k**2)**2 + 3**2*k**2)
-# FIXME There is a simpler way, isn't it?
-ground_truth_spectrum = ift.makeField(power_space,
- mysterious_spectrum(power_space.k_lengths))
+ground_truth_spectrum = ift.PS_field(power_space,
+ lambda k: 5/((7**2 - k**2)**2 + 3**2*k**2))
posterior_power_samples = [pspec.force(KL.position+samp) for samp in KL.samples]
power_mean = reduce(add, posterior_power_samples)
power_plot('teaser_power_reconstruction', ground_truth_spectrum, power_mean,
--
GitLab