diff --git a/nifty/field.py b/nifty/field.py
index d48e5a7ba16a33b0f03a92ba30dd305c2dcf7a37..3fca09e0437cab467432106d362b8f627ee464d9 100644
--- a/nifty/field.py
+++ b/nifty/field.py
@@ -605,19 +605,26 @@ class Field(Loggable, Versionable, object):
         if flipped_Q:
             h *= np.sqrt(2)
             a *= np.sqrt(2)
-            fixed_points = [[fp] if fp is None else fp for fp in fixed_points]
-            for product_point in itertools.product(*fixed_points):
-                slice_object = np.array((slice(None), )*len(val.shape),
-                                        dtype=np.object)
-                for i, sp in enumerate(spaces):
-                    point_component = product_point[i]
-                    if point_component is None:
-                        point_component = slice(None)
-                    slice_object[list(domain_axes[sp])] = point_component
-
-                slice_object = tuple(slice_object)
-                h[slice_object] /= np.sqrt(2)
-                a[slice_object] /= np.sqrt(2)
+
+            # in principle one must not correct the variance for the fixed
+            # points of the hermitianization. However, for a complex field
+            # the input field looses half of its power at its fixed points
+            # in the `hermitian` part. Hence, here a factor of sqrt(2) is
+            # also necessary!
+
+#            fixed_points = [[fp] if fp is None else fp for fp in fixed_points]
+#            for product_point in itertools.product(*fixed_points):
+#                slice_object = np.array((slice(None), )*len(val.shape),
+#                                        dtype=np.object)
+#                for i, sp in enumerate(spaces):
+#                    point_component = product_point[i]
+#                    if point_component is None:
+#                        point_component = slice(None)
+#                    slice_object[list(domain_axes[sp])] = point_component
+#
+#                slice_object = tuple(slice_object)
+#                h[slice_object] /= np.sqrt(2)
+#                a[slice_object] /= np.sqrt(2)
 
         return (h, a)
 
diff --git a/nifty/spaces/lm_space/lm_space.py b/nifty/spaces/lm_space/lm_space.py
index 5cadf785bdb0bd89947d704242fb0fea21db56ce..0cf01177e084d0cbbd8e5703be5f6372cde7626a 100644
--- a/nifty/spaces/lm_space/lm_space.py
+++ b/nifty/spaces/lm_space/lm_space.py
@@ -97,8 +97,8 @@ class LMSpace(Space):
         anti_hermitian_part[:] = x.imag * 1j
         return (hermitian_part, anti_hermitian_part)
 
-    def hermitian_fixed_points(self):
-        return None
+#    def hermitian_fixed_points(self):
+#        return None
 
     # ---Mandatory properties and methods---
 
diff --git a/nifty/spaces/rg_space/rg_space.py b/nifty/spaces/rg_space/rg_space.py
index 887e52de8431296ed9d8de35e8f7dd5259f5afde..e06307d0511c85d5d95ecbb74fcd15cf57857bd4 100644
--- a/nifty/spaces/rg_space/rg_space.py
+++ b/nifty/spaces/rg_space/rg_space.py
@@ -120,17 +120,17 @@ class RGSpace(Space):
 
         return (hermitian_part, anti_hermitian_part)
 
-    def hermitian_fixed_points(self):
-        shape = self.shape
-        mid_index = np.array(shape)//2
-        ndlist = [2 if (shape[i] % 2 == 0) else 1 for i in xrange(len(shape))]
-        ndlist = tuple(ndlist)
-        odd_axes_list = np.array([1 if (shape[i] % 2 == 1) else 0
-                                  for i in xrange(len(shape))])
-        fixed_points = []
-        for i in np.ndindex(ndlist):
-            fixed_points += [tuple((i+odd_axes_list) * mid_index)]
-        return fixed_points
+#    def hermitian_fixed_points(self):
+#        shape = self.shape
+#        mid_index = np.array(shape)//2
+#        ndlist = [2 if (shape[i] % 2 == 0) else 1 for i in xrange(len(shape))]
+#        ndlist = tuple(ndlist)
+#        odd_axes_list = np.array([1 if (shape[i] % 2 == 1) else 0
+#                                  for i in xrange(len(shape))])
+#        fixed_points = []
+#        for i in np.ndindex(ndlist):
+#            fixed_points += [tuple((i+odd_axes_list) * mid_index)]
+#        return fixed_points
 
     def _hermitianize_correct_variance(self, hermitian_part,
                                        anti_hermitian_part, axes):