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

cleaned up Jvec

Browse Source
This commit is contained in:
Lindsey Heagy
2016-03-08 19:56:14 -08:00
parent ceff861413
commit 1b401feb54
1 changed files with 5 additions and 5 deletions
Download Patch File Download Diff File Expand all files Collapse all files
SimPEG/EM/TDEM/TDEM.py
+5 -5
View File
@@ -93,22 +93,22 @@ class BaseTDEMProblem(Problem.BaseTimeProblem, BaseEMProblem):
Asubdiag = self.getAsubdiag(tInd)
for i, src in enumerate(self.survey.srcList):
# compute next du_dm_v for next timestep
un_src = u[src,ftype,tInd+1]
# rhs_v = self.getJRHS(tInd, src, un_src, v)
dA_dm_v = self.getAdiagDeriv(tInd, un_src, v)
dRHS_dm_v = self.getRHSDeriv(tInd, src, v)
# dAsubdiag_dm_v = 0
JRHS = - dA_dm_v - Asubdiag * dun_dm_v[:,i] + dRHS_dm_v
JRHS = dRHS_dm_v - dA_dm_v # - dAsubdiag_dm_v (which is zero)
for rx in src.rxList:
df_dmFun = getattr(u, '_%sDeriv'%rx.projField, None)
df_dm_v[src, '%sDeriv'%rx.projField , tInd] = df_dmFun(tInd, src, dun_dm_v[:,i], v)
# over-write with this time-steps (if not on last timestep, no need to do the extra solve)
# step in time and overwrite
if tInd != len(self.timeSteps):
dun_dm_v[:,i] = Ainv * JRHS
dun_dm_v[:,i] = Ainv * (JRHS - Asubdiag * dun_dm_v[:,i])
for src in self.survey.srcList:
for rx in src.rxList: