Changed to SimPEG solver for DC problem.

2026-07-02 07:05:12 +08:00 · 2013-10-26 16:02:04 -07:00
parent 5cbc980b9c
commit 1555caa775
1 changed files with 2 additions and 7 deletions
@@ -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):