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

Minor updates.

Browse Source
This commit is contained in:
rowanc1
2014-03-20 00:04:56 -07:00
parent 0a5a7aa5e3
commit 879f13df70
1 changed files with 4 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files
simpegEM/FDEM/FDEM.py
+4 -4
View File
@@ -138,13 +138,13 @@ class BaseProblemFDEM(Problem.BaseProblem):
solver = self.Solver(AT, **self.solverOpts)
for tx in self.survey.getTransmitters(freq):
P = tx.projectFieldsDeriv(self.mesh, u)
u_tx = u[tx, self.solType]
P = tx.projectFieldsDeriv(self.mesh, u)
PTv = P.T * v[tx]
fPTv = self.calcFields(PTv, freq, tx.rxList.fieldType, adjoint=True)
w = solver.solve( fPTv )
u_tx = u[tx, self.solType]
Jtv_tx = self.getADeriv(freq, u_tx, w, adjoint=True)
df_dm = self.calcFieldsDeriv(u_tx, freq, tx.rxList.fieldType, PTv, adjoint=True)
@@ -218,9 +218,9 @@ class ProblemFDEM_e(BaseProblemFDEM):
return e
elif fieldType == 'b':
if not adjoint:
b = -1./(1j*omega(freq))*self.mesh.edgeCurl*e
b = -1./(1j*omega(freq)) * ( self.mesh.edgeCurl * e )
else:
b = -1./(1j*omega(freq))*self.mesh.edgeCurl.T*e
b = -1./(1j*omega(freq)) * ( self.mesh.edgeCurl.T * e )
return b
raise NotImplementedError('fieldType "%s" is not implemented.' % fieldType)