diff --git a/test/test_spaces/test_power_space.py b/test/test_spaces/test_power_space.py
index 652d37b390212ac2f2efd04a7a663168f232bfa6..4e8951ab9fb178010d0d834e7857ca42cd2f9087 100644
--- a/test/test_spaces/test_power_space.py
+++ b/test/test_spaces/test_power_space.py
@@ -109,6 +109,16 @@ class PowerSpaceInterfaceTest(unittest.TestCase):
         p = PowerSpace(r)
         assert_(isinstance(getattr(p, attribute), expected_type))
 
+class PowerSpaceConsistencyCheck(unittest.TestCase):
+    @expand(CONSTRUCTOR_CONFIGS)
+    def test_pipundexInversion(self, harmonic_partner, distribution_strategy,
+                         logarithmic, nbin, binbounds, expected):
+        #expected will not be used TODO: write expandproduct to use for this
+        p = PowerSpace(harmonic_partner=harmonic_partner,
+                           distribution_strategy=distribution_strategy,
+                           logarithmic=logarithmic, nbin=nbin,
+                           binbounds=binbounds)
+        assert_equal(p.pindex[p.pundex],np.arange(p.dim),err_msg='pundex is not right-inverse of pindex!')
 
 class PowerSpaceFunctionalityTest(unittest.TestCase):
     @expand(CONSTRUCTOR_CONFIGS)