From 2ce694e9d5f020fb26cafd0fbc99a091176c5d82 Mon Sep 17 00:00:00 2001
From: Rowan Cockett <rowan@3ptscience.com>
Date: Mon, 9 Nov 2015 10:36:59 -0800
Subject: [PATCH] test edgeCurl

---
 tests/mesh/test_pointerMesh.py | 15 +++++++++++++++
 1 file changed, 15 insertions(+)

diff --git a/tests/mesh/test_pointerMesh.py b/tests/mesh/test_pointerMesh.py
index 4147a0a9..7f8d7e3f 100644
--- a/tests/mesh/test_pointerMesh.py
+++ b/tests/mesh/test_pointerMesh.py
@@ -103,6 +103,21 @@ class TestOperatorsOcTree(unittest.TestCase):
         assert (Mr.faceDiv - M.permuteCC*M.faceDiv*M.permuteF.T).nnz == 0
 
 
+    def test_edgeCurl(self):
+
+        hx, hy, hz = np.r_[1.,2,3,4], np.r_[5.,6,7,8], np.r_[9.,10,11,12]
+        M = Tree([hx, hy, hz], levels=2)
+        M.refine(lambda xc:2)
+        # M.plotGrid(showIt=True)
+        Mr = Mesh.TensorMesh([hx, hy, hz])
+
+        # plt.subplot(211).spy(Mr.faceDiv)
+        # plt.subplot(212).spy(M.permuteCC.T*M.faceDiv*M.permuteF)
+        # plt.show()
+
+        assert (Mr.edgeCurl - M.permuteF*M.edgeCurl*M.permuteE.T).nnz == 0
+
+
 
 if __name__ == '__main__':
     unittest.main()