critical filter seems now to work for space volumes other than one as well

demos/critical_filtering.py

@@ -59,7 +59,7 @@ if __name__ == "__main__":
     distribution_strategy = 'not'
 
     # Set up position space
-    dist = 1/128. *10
+    dist = 1/128. *0.1
     s_space = ift.RGSpace([128, 128], distances=[dist,dist])
     # s_space = ift.HPSpace(32)
 
@@ -93,7 +93,7 @@ if __name__ == "__main__":
     # Add a harmonic transformation to the instrument
     R = AdjointFFTResponse(fft, Instrument)
 
-    noise = .1
+    noise = 100.
     ndiag = ift.Field(s_space, noise).weight(1)
     N = ift.DiagonalOperator(s_space, ndiag)
     n = ift.Field.from_random(domain=s_space,

nifty/sugar.py

@@ -78,7 +78,7 @@ def create_power_operator(domain, power_spectrum, dtype=None,
         f = f.real
 
     # f **= 2
-    diag = Field(domain,f)#.weight()
+    diag = Field(domain,f).weight()
     return DiagonalOperator(domain, diag)
 
 
@@ -109,7 +109,7 @@ def generate_posterior_sample(mean, covariance):
     R = covariance.R
     N = covariance.N
 
-    power = sqrt(S.diagonal())
+    power = sqrt(S.diagonal().weight(-1))
     mock_signal = Field.from_random(random_type="normal", domain=S.domain,
                                    dtype=power.dtype)
     mock_signal = power * mock_signal