diff --git a/nifty/operators/invertible_operator_mixin/invertible_operator_mixin.py b/nifty/operators/invertible_operator_mixin/invertible_operator_mixin.py
index af952177e2e46948e0e68ef0c293fca91a2d3e00..483840cc60be9df46e2f3fba114da783a128b1a8 100644
--- a/nifty/operators/invertible_operator_mixin/invertible_operator_mixin.py
+++ b/nifty/operators/invertible_operator_mixin/invertible_operator_mixin.py
@@ -71,7 +71,8 @@ class InvertibleOperatorMixin(object):
 
     def _times(self, x, spaces, x0=None):
         if x0 is None:
-            x0 = Field(self.target, val=0., dtype=x.dtype)
+            x0 = Field(self.target, val=0., dtype=x.dtype,
+                       distribution_strategy=x.distribution_strategy)
 
         (result, convergence) = self.__inverter(A=self.inverse_times,
                                                 b=x,
@@ -80,7 +81,8 @@ class InvertibleOperatorMixin(object):
 
     def _adjoint_times(self, x, spaces, x0=None):
         if x0 is None:
-            x0 = Field(self.domain, val=0., dtype=x.dtype)
+            x0 = Field(self.domain, val=0., dtype=x.dtype,
+                       distribution_strategy=x.distribution_strategy)
 
         (result, convergence) = self.__inverter(A=self.adjoint_inverse_times,
                                                 b=x,
@@ -89,7 +91,8 @@ class InvertibleOperatorMixin(object):
 
     def _inverse_times(self, x, spaces, x0=None):
         if x0 is None:
-            x0 = Field(self.domain, val=0., dtype=x.dtype)
+            x0 = Field(self.domain, val=0., dtype=x.dtype,
+                       distribution_strategy=x.distribution_strategy)
 
         (result, convergence) = self.__inverter(A=self.times,
                                                 b=x,
@@ -98,7 +101,8 @@ class InvertibleOperatorMixin(object):
 
     def _adjoint_inverse_times(self, x, spaces, x0=None):
         if x0 is None:
-            x0 = Field(self.target, val=0., dtype=x.dtype)
+            x0 = Field(self.target, val=0., dtype=x.dtype,
+                       distribution_strategy=x.distribution_strategy)
 
         (result, convergence) = self.__inverter(A=self.adjoint_times,
                                                 b=x,