From adf51d7a568de0ae03aef5f774502fe2e993484c Mon Sep 17 00:00:00 2001
From: rowanc1 <rowan@visiblegeology.com>
Date: Fri, 16 May 2014 15:39:09 -0700
Subject: [PATCH] Remove digs function and replace with sdiag which is backward
 compatible

---
 SimPEG/Mesh/CylMesh.py | 4 +++-
 1 file changed, 3 insertions(+), 1 deletion(-)

diff --git a/SimPEG/Mesh/CylMesh.py b/SimPEG/Mesh/CylMesh.py
index 574ae454..3214a9ed 100644
--- a/SimPEG/Mesh/CylMesh.py
+++ b/SimPEG/Mesh/CylMesh.py
@@ -241,8 +241,10 @@ class CylMesh(BaseTensorMesh, InnerProducts):
             Dr =  sp.kron(sp.identity(self.nNz), dr)
             Dz = -sp.kron(dz, sp.identity(self.nCx))
 
+            A = self.area
+            E = self.edge
             #Edge curl operator
-            self._edgeCurl = sp.diags(1/self.area,0)*sp.vstack((Dz, Dr))*sp.diags(self.edge,0)
+            self._edgeCurl = sdiag(1/A)*sp.vstack((Dz, Dr))*sdiag(E)
         return self._edgeCurl
 
     # @property