diff --git a/test/test_field.py b/test/test_field.py
index 590f8ba54ccb1e0bcd8923a7587629e1bc49f074..c089e1242687eaeb179c4caf7aeb0dfe573a90b9 100644
--- a/test/test_field.py
+++ b/test/test_field.py
@@ -73,6 +73,25 @@ class Test_Interface(unittest.TestCase):
         assert_almost_equal(h1.get_full_data(),h3.get_full_data())
         assert_almost_equal(a1.get_full_data(),a3.get_full_data())
 
+    @expand(product([False,True],[False,True]))
+    def test_hermitian_decomposition1(self, complexdata, preserve):
+        s0=(1,)
+        s1=(56,25)
+        r0 = RGSpace(s0, harmonic=True)
+        r1 = RGSpace(s1, harmonic=True)
+        ra = RGSpace(s0+s1, harmonic=True)
+        v = np.random.random(s1)
+        if (complexdata):
+            v = v + 1j*np.random.random(s1)
+        f1=Field(r1,val=v,copy=True)
+        f2=Field((r0,r1),val=v,copy=True)
+        h1,a1 = Field._hermitian_decomposition((r1,),f1.val,(0,),((0,1,),),preserve)
+        h2,a2 = Field._hermitian_decomposition((r0,r1),f2.val,(0,1),((0,),(1,2)),preserve)
+        h2=h2[0,:,:]
+        a2=a2[0,:,:]
+        assert_almost_equal(h1.get_full_data(),h2.get_full_data())
+        assert_almost_equal(a1.get_full_data(),a2.get_full_data())
+
 #class Test_Initialization(unittest.TestCase):
 #
 #    @parameterized.expand(