From 71a5c8b2c1743cf8a892c3d3d59404c46bc4c947 Mon Sep 17 00:00:00 2001 From: Lindsey Heagy Date: Tue, 24 Nov 2015 11:01:51 -0800 Subject: [PATCH] fix sizes in Fields object (make an array not a vector) --- SimPEG/EM/FDEM/FieldsFDEM.py | 5 ++--- 1 file changed, 2 insertions(+), 3 deletions(-) diff --git a/SimPEG/EM/FDEM/FieldsFDEM.py b/SimPEG/EM/FDEM/FieldsFDEM.py index f90dbfef..8f6fafe9 100644 --- a/SimPEG/EM/FDEM/FieldsFDEM.py +++ b/SimPEG/EM/FDEM/FieldsFDEM.py @@ -159,7 +159,7 @@ class Fields_b(Fields): Me = Me.T w = self._edgeCurl.T * (self._MfMui * bSolution) - w += -Utils.mkvc(Me * S_e,2) + w = w - Utils.mkvc(Me * S_e,2) if not adjoint: de_dm = self._MeSigmaIDeriv(w) * v @@ -169,8 +169,7 @@ class Fields_b(Fields): _, S_eDeriv = src.evalDeriv(self.prob, adjoint) Se_Deriv = S_eDeriv(v) - if Se_Deriv is not None: - de_dm = de_dm-self._MeSigmaI * Se_Deriv + de_dm = de_dm - self._MeSigmaI * Se_Deriv return de_dm