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

model versus sigma bug

Browse Source
This commit is contained in:
rowanc1
2014-03-17 10:49:06 -07:00
parent 95e5a09523
commit 3fd840201d
1 changed files with 4 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files
simpegEM/FDEM/FDEM.py
+4 -2
View File
@@ -97,6 +97,7 @@ class ProblemFDEM_e(Problem.BaseProblem):
u = self.fields(m)
Jv = self.dataPair(self.survey)
sig = self.model.transform(m)
dsig_dm = self.model.transformDeriv(m)
for i, freq in enumerate(self.survey.freqs):
@@ -105,7 +106,7 @@ class ProblemFDEM_e(Problem.BaseProblem):
solver = Solver(A, options=self.solveOpts)
for txi, tx in enumerate(self.survey.getTransmitters(freq)):
dMe_dsig = self.mesh.getEdgeInnerProductDeriv(m, v=e[:,txi])
dMe_dsig = self.mesh.getEdgeInnerProductDeriv(sig, v=e[:,txi])
P = tx.projectFieldsDeriv(self.mesh, u)
b = 1j*omega(freq) * ( dMe_dsig * ( dsig_dm * v ) )
@@ -125,6 +126,7 @@ class ProblemFDEM_e(Problem.BaseProblem):
Jtv = np.zeros(self.model.nP, dtype=complex)
sig = self.model.transform(m)
dsig_dm = self.model.transformDeriv(m)
for i, freq in enumerate(self.survey.freqs):
@@ -133,7 +135,7 @@ class ProblemFDEM_e(Problem.BaseProblem):
solver = Solver(AT, options=self.solveOpts)
for txi, tx in enumerate(self.survey.getTransmitters(freq)):
dMe_dsig = self.mesh.getEdgeInnerProductDeriv(m, v=e[:,txi])
dMe_dsig = self.mesh.getEdgeInnerProductDeriv(sig, v=e[:,txi])
P = tx.projectFieldsDeriv(self.mesh, u)
w = solver.solve(P.T * v[tx])