wassname/simpeg
1
0
Fork 0
mirror of https://github.com/wassname/simpeg.git synced 2026-07-07 05:36:02 +08:00
Code Issues Packages Projects Releases Wiki Activity

make sure factors of JtVec get cleaned

Browse Source
This commit is contained in:
Lindsey Heagy
2016-06-28 11:25:17 -07:00
parent 43d0b0369f
commit e179b5e71d
1 changed files with 2 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files
SimPEG/EM/TDEM/TDEM.py
+2 -2
View File
@@ -175,7 +175,7 @@ class BaseTDEMProblem(Problem.BaseTimeProblem, BaseEMProblem):
for src in self.survey.srcList:
PT_v = Fields_Derivs(self.mesh, self.survey) # initialize storage for PT_v (don't need to preserve over sources)
# initialize size
df_duT_v[src, '%sDeriv'%self._fieldType, :] = np.zeros_like(f[src, self._fieldType, :])
df_duT_v[src, '%sDeriv'%self._fieldType, :] = sp.zeros_like(f[src, self._fieldType, :])
for rx in src.rxList:
PT_v[src,'%sDeriv'%rx.projField,:] = rx.evalDeriv(src, self.mesh, self.timeMesh, Utils.mkvc(v[src,rx]), adjoint=True) # this is +=
@@ -272,7 +272,7 @@ class BaseTDEMProblem(Problem.BaseTimeProblem, BaseEMProblem):
# JTv = JTv - Utils.mkvc(self.getAdiagDeriv(0, f[src, ftype, tInd], v, adjoint = True))
# # JTv = JTv + self.getInitialFieldsDeriv(Utils.mkvc(df_duT_v[src,'%sDeriv'%self._fieldType,0] - Asubdiag.T * Utils.mkvc(ATinv_df_duT_v[isrc,:])), adjoint=True)
AdiagTinv.clean()
return Utils.mkvc(JTv).astype(float)