Playing with 2.5D DC inversion

TODOs:

     Parallelize ky ... (currently pretty slow)

SimPEG/EM/Static/DC/ProblemDC_2D.py

@@ -50,6 +50,7 @@ class BaseDCProblem_2D(BaseEMProblem):
         dky = np.r_[dky[0], dky]
         y = 0.
 
+        #TODO: this loop is pretty slow .. (Parellize)
         for iky in range(self.nky):
             ky = self.kys[iky]
             A = self.getA(ky)
@@ -94,6 +95,7 @@ class BaseDCProblem_2D(BaseEMProblem):
             for rx in src.rxList:
                 Jtv_temp1 = np.zeros(m.size)
                 Jtv_temp0 = np.zeros(m.size)
+                #TODO: this loop is pretty slow .. (Parellize)
                 for iky in range(self.nky):
                     u_src = f[src, self._solutionType, iky]
                     ky = self.kys[iky]

SimPEG/EM/Static/Utils/StaticUtils.py

@@ -213,7 +213,6 @@ def gen_DCIPsurvey(endl, mesh, stype, a, b, n):
     # Pole-dipole: Moving pole on one end -> [A a MN1 a MN2 ... MNn a B]
     SrcList = []
 
-
     if stype != 'gradient':
 
         for ii in range(0, int(nstn)-1):
@@ -262,6 +261,7 @@ def gen_DCIPsurvey(endl, mesh, stype, a, b, n):
                     srcClass = DC.Src.Pole([rxClass], M[ii,:])
             SrcList.append(srcClass)
 
+    # TODO: test gradient array!
     elif stype == 'gradient':
 
         # Gradient survey only requires Tx at end of line and creates a square