12/n

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,