wassname/simpeg
1
0
Fork 0
mirror of https://github.com/wassname/simpeg.git synced 2026-08-11 11:26:01 +08:00
Code Issues Packages Projects Releases Wiki Activity

Adjoint for all 4 formulations with Zero Identity

Browse Source
This commit is contained in:
Lindsey Heagy
2015-11-08 14:09:18 -08:00
parent 0af1f59891
commit 639a6e7a06
2 changed files with 18 additions and 18 deletions
Download Patch File Download Diff File Expand all files Collapse all files
SimPEG/EM/FDEM/FDEM.py
+16 -16
View File
@@ -130,30 +130,30 @@ class BaseFDEMProblem(BaseEMProblem):
df_duTFun = getattr(f, '_%sDeriv_u'%rx.projField, None)
df_duT = df_duTFun(src, PTv, adjoint=True)
if df_duT is not None:
dA_duIT = ATinv * df_duT
else:
dA_duIT = ATinv * PTv
# if df_duT is not None:
dA_duIT = ATinv * df_duT
# else:
# dA_duIT = ATinv * PTv
dA_dmT = self.getADeriv_m(freq, u_src, dA_duIT, adjoint=True)
dRHS_dmT = self.getRHSDeriv_m(src, dA_duIT, adjoint=True)
dRHS_dmT = self.getRHSDeriv_m(freq,src, dA_duIT, adjoint=True)
if dRHS_dmT is None:
du_dmT = - dA_dmT
else:
du_dmT = -dA_dmT + dRHS_dmT
# if dRHS_dmT is None:
# du_dmT = - dA_dmT
# else:
du_dmT = -dA_dmT + dRHS_dmT
df_dmFun = getattr(f, '_%sDeriv_m'%rx.projField, None)
dfT_dm = df_dmFun(src, PTv, adjoint=True)
if dfT_dm is not None:
du_dmT += dfT_dm
# if dfT_dm is not None:
du_dmT += dfT_dm
real_or_imag = rx.projComp
if real_or_imag == 'real':
Jtv += du_dmT.real
elif real_or_imag == 'imag':
Jtv += - du_dmT.real
if real_or_imag is 'real':
Jtv += np.array(du_dmT,dtype=complex).real
elif real_or_imag is 'imag':
Jtv += - np.array(du_dmT,dtype=complex).real
else:
raise Exception('Must be real or imag')
@@ -399,7 +399,7 @@ class ProblemFDEM_b(BaseFDEMProblem):
if not adjoint:
SrcDeriv = S_mDeriv + C * (self.MeSigmaI * S_eDeriv)
elif adjoint:
SrcDeriv = S_mDeriv + Self.MeSigmaI.T * ( C.T * S_eDeriv)
SrcDeriv = S_mDeriv + self.MeSigmaI.T * ( C.T * S_eDeriv)
# elif S_mDeriv is not None:
# SrcDeriv = S_mDeriv
# elif S_eDeriv is not None:
SimPEG/EM/FDEM/SurveyFDEM.py
+2 -2
View File
@@ -118,10 +118,10 @@ class SrcFDEM(Survey.BaseSrc):
return None
def S_m(self, prob):
return Zero()
return None
def S_e(self, prob):
return Zero()
return None
def S_mDeriv(self, prob, v, adjoint = False):
return Zero()