---------------------------------------------------------------------------
KeyboardInterrupt Traceback (most recent call last)
<ipython-input-41-e70b86582598> in <module>()
----> 1 mopt = inv.run(m0)
/home/seogi/Documents/simpeg/SimPEG/Utils/CounterUtils.pyc in wrapper(self, *args, **kwargs)
90 counter = getattr(self,'counter',None)
91 if type(counter) is Counter: counter.countTic(self.__class__.__name__+'.'+f.__name__)
---> 92 out = f(self,*args,**kwargs)
93 if type(counter) is Counter: counter.countToc(self.__class__.__name__+'.'+f.__name__)
94 return out
/home/seogi/Documents/simpeg/SimPEG/Inversion.pyc in run(self, m0)
60 self.invProb.startup(m0)
61 self.directiveList.call('initialize')
---> 62 self.m = self.opt.minimize(self.invProb.evalFunction, self.invProb.curModel)
63 self.directiveList.call('finish')
64
/home/seogi/Documents/simpeg/SimPEG/Utils/CounterUtils.pyc in wrapper(self, *args, **kwargs)
90 counter = getattr(self,'counter',None)
91 if type(counter) is Counter: counter.countTic(self.__class__.__name__+'.'+f.__name__)
---> 92 out = f(self,*args,**kwargs)
93 if type(counter) is Counter: counter.countToc(self.__class__.__name__+'.'+f.__name__)
94 return out
/home/seogi/Documents/simpeg/SimPEG/Optimization.pyc in minimize(self, evalFunction, x0)
179 self.printIter()
180 if self.stoppingCriteria(): break
--> 181 self.searchDirection = self.findSearchDirection()
182 del self.H #: Doing this saves memory, as it is not needed in the rest of the computations.
183 p = self.scaleSearchDirection(self.searchDirection)
/home/seogi/Documents/simpeg/SimPEG/Utils/CounterUtils.pyc in wrapper(self, *args, **kwargs)
90 counter = getattr(self,'counter',None)
91 if type(counter) is Counter: counter.countTic(self.__class__.__name__+'.'+f.__name__)
---> 92 out = f(self,*args,**kwargs)
93 if type(counter) is Counter: counter.countToc(self.__class__.__name__+'.'+f.__name__)
94 return out
/home/seogi/Documents/simpeg/SimPEG/Optimization.pyc in findSearchDirection(self)
974
975 # Form product Hessian*pc.
--> 976 Hp = self.H*pc
977 Hp = (1-Active)*Hp
978
/usr/local/lib/python2.7/dist-packages/scipy/sparse/linalg/interface.pyc in __mul__(self, x)
201
202 if x.ndim == 1 or x.ndim == 2 and x.shape[1] == 1:
--> 203 return self.matvec(x)
204 elif x.ndim == 2:
205 return self.matmat(x)
/usr/local/lib/python2.7/dist-packages/scipy/sparse/linalg/interface.pyc in matvec(self, x)
132 raise ValueError('dimension mismatch')
133
--> 134 y = self._matvec(x)
135
136 if isinstance(x, np.matrix):
/home/seogi/Documents/simpeg/SimPEG/InvProblem.pyc in H_fun(v)
132 if return_H:
133 def H_fun(v):
--> 134 phi_d2Deriv = self.dmisfit.eval2Deriv(m, v, u=u)
135 phi_m2Deriv = self.reg.eval2Deriv(m, v=v)
136
/home/seogi/Documents/simpeg/SimPEG/Utils/CounterUtils.pyc in wrapper(self, *args, **kwargs)
90 counter = getattr(self,'counter',None)
91 if type(counter) is Counter: counter.countTic(self.__class__.__name__+'.'+f.__name__)
---> 92 out = f(self,*args,**kwargs)
93 if type(counter) is Counter: counter.countToc(self.__class__.__name__+'.'+f.__name__)
94 return out
/home/seogi/Documents/simpeg/SimPEG/DataMisfit.pyc in eval2Deriv(self, m, v, u)
135 prob = self.prob
136 if u is None: u = prob.fields(m)
--> 137 return prob.Jtvec_approx(m, self.Wd * (self.Wd * prob.Jvec_approx(m, v, u=u)), u=u)
/home/seogi/Documents/simpeg/SimPEG/Utils/CounterUtils.pyc in wrapper(self, *args, **kwargs)
90 counter = getattr(self,'counter',None)
91 if type(counter) is Counter: counter.countTic(self.__class__.__name__+'.'+f.__name__)
---> 92 out = f(self,*args,**kwargs)
93 if type(counter) is Counter: counter.countToc(self.__class__.__name__+'.'+f.__name__)
94 return out
/home/seogi/Documents/simpeg/SimPEG/Problem.pyc in Jtvec_approx(self, m, v, u)
399 :return: JTv
400 """
--> 401 return self.Jtvec(m, v, u)
402
403 def fields(self, m):
/home/seogi/Documents/simpeg/SimPEG/Utils/CounterUtils.pyc in wrapper(self, *args, **kwargs)
90 counter = getattr(self,'counter',None)
91 if type(counter) is Counter: counter.countTic(self.__class__.__name__+'.'+f.__name__)
---> 92 out = f(self,*args,**kwargs)
93 if type(counter) is Counter: counter.countToc(self.__class__.__name__+'.'+f.__name__)
94 return out
/home/seogi/Documents/simpegpf/simpegPF/Magnetics.pyc in Jtvec(self, m, v, u)
328 Atemp = Utils.sdiag(self.MfMu0*B0)*(dMfMuI * (dmudm))
329 Btemp = Utils.sdiag(Div.T*u)*(dMfMuI* (dmudm))
--> 330 Jtv = Atemp.T*(P.T*v) - Btemp.T*(P.T*v) - Ctv
331
332 return Utils.mkvc(Jtv)
/usr/local/lib/python2.7/dist-packages/scipy/sparse/base.pyc in __mul__(self, other)
294 # Fast path for the most common case
295 if other.shape == (N,):
--> 296 return self._mul_vector(other)
297 elif other.shape == (N, 1):
298 return self._mul_vector(other.ravel()).reshape(M, 1)
/usr/local/lib/python2.7/dist-packages/scipy/sparse/compressed.pyc in _mul_vector(self, other)
445 # csr_matvec or csc_matvec
446 fn = getattr(_sparsetools,self.format + '_matvec')
--> 447 fn(M, N, self.indptr, self.indices, self.data, other, result)
448
449 return result
KeyboardInterrupt: