Merged in boundaryConditions (pull request #25)

Cell Gradient and Boundary Conditions

SimPEG/utils/__init__.py

@@ -4,7 +4,7 @@ import lomutils
 import interputils
 import ModelBuilder
 from matutils import getSubArray, mkvc, ndgrid, ind2sub, sub2ind
-from sputils import spzeros, kron3, speye, sdiag, ddx, av
+from sputils import spzeros, kron3, speye, sdiag, ddx, av, avExtrap
 from lomutils import volTetra, faceInfo, inv2X2BlockDiagonal, inv3X3BlockDiagonal, indexCube, exampleLomGird
 from interputils import interpmat
 from ipythonUtils import easyAnimate as animate

SimPEG/utils/sputils.py

@@ -29,5 +29,10 @@ def ddx(n):
 
 
 def av(n):
-    """Define 1D averaging operator from cell-centres to nodes."""
+    """Define 1D averaging operator from nodes to cell-centers."""
     return sp.spdiags((0.5*np.ones((n+1, 1))*[1, 1]).T, [0, 1], n, n+1, format="csr")
+
+def avExtrap(n):
+    """Define 1D averaging operator from cell-centers to nodes."""
+    Av = sp.spdiags((0.5*np.ones((n, 1))*[1, 1]).T, [-1, 0], n+1, n, format="csr") + sp.csr_matrix(([0.5,0.5],([0,n],[0,n-1])),shape=(n+1,n))
+    return Av