diff --git a/SimPEG/__init__.py b/SimPEG/__init__.py
index e8946f88..a17eaa74 100644
--- a/SimPEG/__init__.py
+++ b/SimPEG/__init__.py
@@ -1,4 +1,6 @@
 import mesh
 import utils
 import inverse
+import forward
+import regularization
 from Solver import Solver
diff --git a/SimPEG/forward/DCProblem/DCProblem.py b/SimPEG/forward/DCProblem.py
similarity index 82%
rename from SimPEG/forward/DCProblem/DCProblem.py
rename to SimPEG/forward/DCProblem.py
index 45b8a2c1..5aad45bd 100644
--- a/SimPEG/forward/DCProblem/DCProblem.py
+++ b/SimPEG/forward/DCProblem.py
@@ -3,8 +3,8 @@ from SimPEG.forward import Problem, SyntheticProblem
 from SimPEG.tests import checkDerivative
 from SimPEG.utils import ModelBuilder, sdiag
 import numpy as np
+import scipy as sp
 import scipy.sparse.linalg as linalg
-import DCutils
 
 class DCProblem(Problem):
     """
@@ -134,7 +134,7 @@ if __name__ == '__main__':
     elecend = 0.5+spacelec*(nelec-1)
     elecLocR = np.linspace(elecini, elecend, nelec)
     rxmidLoc = (elecLocR[0:nelec-1]+elecLocR[1:nelec])*0.5
-    q, Q, rxmidloc = DCutils.genTxRxmat(nelec, spacelec, surfloc, elecini, mesh)
+    q, Q, rxmidloc = genTxRxmat(nelec, spacelec, surfloc, elecini, mesh)
     P = Q.T
 
     # Create some data
@@ -179,3 +179,32 @@ if __name__ == '__main__':
     print v.dot(problem.Jt(mSynth, w, u=u))
 
 
+
+
+
+def genTxRxmat(nelec, spacelec, surfloc, elecini, mesh):
+    """ Generate projection matrix (Q) and """
+    elecend = 0.5+spacelec*(nelec-1)
+    elecLocR = np.linspace(elecini, elecend, nelec)
+    elecLocT = elecLocR+1
+    nrx = nelec-1
+    ntx = nelec-1
+    q = np.zeros((mesh.nC, ntx))
+    Q = np.zeros((mesh.nC, nrx))
+
+    for i in range(nrx):
+
+        rxind1 = np.argwhere((mesh.gridCC[:,0]==surfloc) & (mesh.gridCC[:,1]==elecLocR[i]))
+        rxind2 = np.argwhere((mesh.gridCC[:,0]==surfloc) & (mesh.gridCC[:,1]==elecLocR[i+1]))
+
+        txind1 = np.argwhere((mesh.gridCC[:,0]==surfloc) & (mesh.gridCC[:,1]==elecLocT[i]))
+        txind2 = np.argwhere((mesh.gridCC[:,0]==surfloc) & (mesh.gridCC[:,1]==elecLocT[i+1]))
+
+        q[txind1,i] = 1
+        q[txind2,i] = -1
+        Q[rxind1,i] = 1
+        Q[rxind2,i] = -1
+
+    Q = sp.csr_matrix(Q)
+    rxmidLoc = (elecLocR[0:nelec-1]+elecLocR[1:nelec])*0.5
+    return q, Q, rxmidLoc
diff --git a/SimPEG/forward/DCProblem/DCutils.py b/SimPEG/forward/DCProblem/DCutils.py
deleted file mode 100644
index f3445096..00000000
--- a/SimPEG/forward/DCProblem/DCutils.py
+++ /dev/null
@@ -1,29 +0,0 @@
-import numpy as np
-import scipy.sparse as sp
-
-def genTxRxmat(nelec, spacelec, surfloc, elecini, mesh):
-    """ Generate projection matrix (Q) and """
-    elecend = 0.5+spacelec*(nelec-1)
-    elecLocR = np.linspace(elecini, elecend, nelec)
-    elecLocT = elecLocR+1
-    nrx = nelec-1
-    ntx = nelec-1
-    q = np.zeros((mesh.nC, ntx))
-    Q = np.zeros((mesh.nC, nrx))
-
-    for i in range(nrx):
-
-        rxind1 = np.argwhere((mesh.gridCC[:,0]==surfloc) & (mesh.gridCC[:,1]==elecLocR[i]))
-        rxind2 = np.argwhere((mesh.gridCC[:,0]==surfloc) & (mesh.gridCC[:,1]==elecLocR[i+1]))
-
-        txind1 = np.argwhere((mesh.gridCC[:,0]==surfloc) & (mesh.gridCC[:,1]==elecLocT[i]))
-        txind2 = np.argwhere((mesh.gridCC[:,0]==surfloc) & (mesh.gridCC[:,1]==elecLocT[i+1]))
-
-        q[txind1,i] = 1
-        q[txind2,i] = -1
-        Q[rxind1,i] = 1
-        Q[rxind2,i] = -1
-
-    Q = sp.csr_matrix(Q)
-    rxmidLoc = (elecLocR[0:nelec-1]+elecLocR[1:nelec])*0.5
-    return q, Q, rxmidLoc
diff --git a/SimPEG/forward/DCProblem/__init__.py b/SimPEG/forward/DCProblem/__init__.py
deleted file mode 100644
index a868cf80..00000000
--- a/SimPEG/forward/DCProblem/__init__.py
+++ /dev/null
@@ -1,2 +0,0 @@
-from DCProblem import *
-from DCutils import *
diff --git a/SimPEG/forward/__init__.py b/SimPEG/forward/__init__.py
index fe849d41..5e326a59 100644
--- a/SimPEG/forward/__init__.py
+++ b/SimPEG/forward/__init__.py
@@ -1,2 +1,2 @@
 from Problem import *
-import DCProblem
+from DCProblem import DCProblem