diff --git a/SimPEG/Solver.py b/SimPEG/Solver.py
index 3e9c2d02..60ca5558 100644
--- a/SimPEG/Solver.py
+++ b/SimPEG/Solver.py
@@ -156,7 +156,7 @@ class Solver(object):
         if len(b.shape) == 1 or b.shape[1] == 1:
             # Just one RHS
             if factorize:
-                return self.dsolve(b)
+                return self.dsolve(b.flatten())
             else:
                 return linalg.dsolve.spsolve(self.A, b)