--- trunk/finley/test/python/test_utilOnFinley.py	2006/01/19 22:36:36	436
+++ trunk/finley/test/python/test_utilOnFinley.py	2006/01/20 00:16:58	437
@@ -2,8 +2,9 @@
 
 import unittest
 from esys.escript.test_util import Test_util_no_tagged_data as Test_util
+from esys.escript.test_util import Test_Util_SpatialFunctions
 from esys.escript import FunctionOnBoundary
-from esys.finley import Rectangle
+from esys.finley import Rectangle,Brick,JoinFaces
 import sys
 
 class Test_UtilOnFinley(Test_util):
@@ -11,9 +12,46 @@
        self.__dom =Rectangle(10,10,2)
        self.functionspace = FunctionOnBoundary(self.__dom) # due to a bug in escript python needs to hold a reference to the domain
 
+class Test_Util_SpatialFunctionsOnFinley2DOrder1(Test_Util_SpatialFunctions):
+    def setUp(self):
+        self.order=1
+        d1 = Rectangle(n0=6,n1=12,l0=0.5,order=1,useElementsOnFace=True)
+        d2 = Rectangle(n0=6,n1=12,l0=0.5,order=1,useElementsOnFace=True)
+        d2.setX(d2.getX()+[0.5,0.])
+        self.domain = JoinFaces([d1,d2])
+
+class Test_Util_SpatialFunctionsOnFinley2DOrder2(Test_Util_SpatialFunctions):
+    def setUp(self):
+        self.order=2
+        d1 = Rectangle(n0=3,n1=6,l0=0.5,order=2,useElementsOnFace=True)
+        d2 = Rectangle(n0=3,n1=6,l0=0.5,order=2,useElementsOnFace=True)
+        d2.setX(d2.getX()+[0.5,0.])
+        self.domain = JoinFaces([d1,d2])
+
+class Test_Util_SpatialFunctionsOnFinley3DOrder1(Test_Util_SpatialFunctions):
+    def setUp(self):
+        self.order=1
+        d1 = Brick(n0=6,n1=12,n2=12,l0=0.5,order=1,useElementsOnFace=True)
+        d2 = Brick(n0=6,n1=12,n2=12,l0=0.5,order=1,useElementsOnFace=True)
+        d2.setX(d2.getX()+[0.5,0.,0.])
+        self.domain = JoinFaces([d1,d2])
+
+class Test_Util_SpatialFunctionsOnFinley3DOrder2(Test_Util_SpatialFunctions):
+    def setUp(self):
+        self.order=2
+        d1 = Brick(n0=3,n1=6,n2=6,l0=0.5,order=2,useElementsOnFace=True)
+        d2 = Brick(n0=3,n1=6,n2=6,l0=0.5,order=2,useElementsOnFace=True)
+        d2.setX(d2.getX()+[0.5,0.,0.])
+        self.domain = JoinFaces([d1,d2])
+
+
 if __name__ == '__main__':
    suite = unittest.TestSuite()
    suite.addTest(unittest.makeSuite(Test_UtilOnFinley))
+   suite.addTest(unittest.makeSuite(Test_Util_SpatialFunctionsOnFinley2DOrder1))
+   suite.addTest(unittest.makeSuite(Test_Util_SpatialFunctionsOnFinley2DOrder2))
+   suite.addTest(unittest.makeSuite(Test_Util_SpatialFunctionsOnFinley3DOrder1))
+   suite.addTest(unittest.makeSuite(Test_Util_SpatialFunctionsOnFinley3DOrder2))
    s=unittest.TextTestRunner(verbosity=2).run(suite)
    if s.wasSuccessful():
      sys.exit(0)