diff --git a/tests/test_qpms_p.py b/tests/test_qpms_p.py
new file mode 100644
index 0000000..2b6b354
--- /dev/null
+++ b/tests/test_qpms_p.py
@@ -0,0 +1,29 @@
+import unittest
+import qpms
+import numpy as np
+
+# Some constants go here.
+maxx = 20 # The maximum argument of the Bessel's functions, i.e. maximum wave number times the distance
+lMax = 25 # To which order we decompose the waves
+
+
+class PlaneWaveDecompositionTests(unittest.TestCase):
+
+    def testRandomInc(self):
+        pass
+
+    def testCornerCases(self):
+        pass
+
+
+class SphericalWaveTranslationTests(unittest.TestCase):
+
+    def sometest(self):
+        pass
+
+def main():
+    unittest.main()
+
+if __name__ == '__main__':
+    main()
+