wassname/simpeg
1
0
Fork 0
mirror of https://github.com/wassname/simpeg.git synced 2026-06-27 18:25:42 +08:00
Code Issues Packages Projects Releases Wiki Activity

Remove data class for Jv (possible memory leakage)

Browse Source
This commit is contained in:
seogi_macbook
2016-07-03 23:22:04 -07:00
parent 3f3b2e9232
commit 2d992f32d4
1 changed files with 8 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files
SimPEG/EM/Static/IP/ProblemIP.py
+8 -4
View File
@@ -45,7 +45,8 @@ class BaseIPProblem(BaseEMProblem):
self.curModel = m
Jv = self.dataPair(self.survey) #same size as the data
# Jv = self.dataPair(self.survey) #same size as the data
Jv = []
A = self.getA()
@@ -58,13 +59,16 @@ class BaseIPProblem(BaseEMProblem):
for rx in src.rxList:
df_dmFun = getattr(f, '_%sDeriv'%rx.projField, None)
df_dm_v = df_dmFun(src, du_dm_v, v, adjoint=False)
Jv[src, rx] = rx.evalDeriv(src, self.mesh, f, df_dm_v)
# Jv[src, rx] = rx.evalDeriv(src, self.mesh, f, df_dm_v)
Jv.append(rx.evalDeriv(src, self.mesh, f, df_dm_v))
# Conductivity (d u / d log sigma)
if self._formulation is 'EB':
return -Utils.mkvc(Jv)
# return -Utils.mkvc(Jv)
return -np.hstack(Jv)
# Conductivity (d u / d log rho)
if self._formulation is 'HJ':
return Utils.mkvc(Jv)
# return Utils.mkvc(Jv)
return np.hstack(Jv)
def Jtvec(self, m, v, f=None):
if f is None: