temporary fix for probing

ChangeLog

@@ -5,6 +5,7 @@ Interface change in EndomorphicOperator.draw_sample()
 
 This method now requires a `dtype` argument to be passed.
 As a consequence, `dtype` moves to the first place of the argument list.
+(This of course applies to all derived classes as well.)
 
 MPI parallelisation over samples in MetricGaussianKL
 ====================================================

nifty6/probing.py

@@ -97,10 +97,12 @@ def probe_with_posterior_samples(op, post_op, nprobes):
             raise ValueError
     sc = StatCalculator()
     for i in range(nprobes):
+        # FIXME which dtype should we use here?
+        import numpy as np
         if post_op is None:
-            sc.add(op.draw_sample(from_inverse=True))
+            sc.add(op.draw_sample(dtype=np.float64, from_inverse=True))
         else:
-            sc.add(post_op(op.draw_sample(from_inverse=True)))
+            sc.add(post_op(op.draw_sample(dtype=np.float64, from_inverse=True)))
 
     if nprobes == 1:
         return sc.mean, None