temporary renaming to make things less confusing

nifty4/library/nonlinear_power_energy.py

@@ -70,7 +70,7 @@ class NonlinearPowerEnergy(Energy):
             if samples is None or samples == 0:
                 xi_sample_list = [xi]
             else:
-                xi_sample_list = [D.draw_sample() + xi
+                xi_sample_list = [D.draw_inverse_sample() + xi
                                   for _ in range(samples)]
         self.xi_sample_list = xi_sample_list
         self.inverter = inverter

nifty4/library/wiener_filter_curvature.py

@@ -60,7 +60,7 @@ class WienerFilterCurvature(EndomorphicOperator):
     def apply(self, x, mode):
         return self._op.apply(x, mode)
 
-    def draw_sample(self, dtype=np.float64):
+    def draw_inverse_sample(self, dtype=np.float64):
         n = self.N.draw_sample(dtype)
         s = self.S.draw_sample(dtype)
 

nifty4/probing/utils.py

@@ -51,7 +51,7 @@ class StatCalculator(object):
 def probe_with_posterior_samples(op, post_op, nprobes):
     sc = StatCalculator()
     for i in range(nprobes):
-        sample = post_op(op.draw_sample())
+        sample = post_op(op.draw_inverse_sample())
         sc.add(sample)
 
     if nprobes == 1:

test/test_energies/test_noise.py

@@ -84,7 +84,7 @@ class Noise_Energy_Tests(unittest.TestCase):
             S=S,
             inverter=inverter).curvature
 
-        res_sample_list = [d - R(f(ht(C.draw_sample() + xi)))
+        res_sample_list = [d - R(f(ht(C.draw_inverse_sample() + xi)))
                            for _ in range(10)]
 
         energy0 = ift.library.NoiseEnergy(eta0, alpha, q, res_sample_list)