From 1555caa7758aa46d5771a3221391c171fbfd645e Mon Sep 17 00:00:00 2001
From: Rowan Cockett <rowan@visiblegeology.com>
Date: Sat, 26 Oct 2013 16:02:04 -0700
Subject: [PATCH] Changed to SimPEG solver for DC problem.

---
 SimPEG/forward/DCProblem.py | 9 ++-------
 1 file changed, 2 insertions(+), 7 deletions(-)

diff --git a/SimPEG/forward/DCProblem.py b/SimPEG/forward/DCProblem.py
index b8e8248d..c2a615ed 100644
--- a/SimPEG/forward/DCProblem.py
+++ b/SimPEG/forward/DCProblem.py
@@ -40,13 +40,8 @@ class DCProblem(ModelTransforms.LogModel, Problem):
 
     def field(self, m):
         A = self.createMatrix(m)
-        solve = linalg.factorized(A)
-
-        nRHSs = self.RHS.shape[1]  # Number of RHSs
-        phi = np.zeros((self.mesh.nC, nRHSs)) + np.nan
-        for ii in range(nRHSs):
-            phi[:,ii] = solve(self.RHS[:,ii])
-
+        solve = Solver(A)
+        phi = solve.solve(self.RHS)
         return phi
 
     def J(self, m, v, u=None):