diff --git a/simpegEM/Base.py b/simpegEM/Base.py
index 82a508d1..00b628ca 100644
--- a/simpegEM/Base.py
+++ b/simpegEM/Base.py
@@ -124,5 +124,5 @@ class BaseEMProblem(Problem.BaseProblem):
 
     def fields(self, m):
         self.curModel = m
-        F = self.forward(m, self.getRHS)
+        F = self.forward(m)
         return F
diff --git a/simpegEM/FDEM/FDEM.py b/simpegEM/FDEM/FDEM.py
index d350c18d..0eb5962f 100644
--- a/simpegEM/FDEM/FDEM.py
+++ b/simpegEM/FDEM/FDEM.py
@@ -19,13 +19,13 @@ class BaseFDEMProblem(BaseEMProblem):
     surveyPair = SurveyFDEM
     fieldsPair = FieldsFDEM
 
-    def forward(self, m, RHS):
+    def forward(self, m):
 
         F = self.fieldsPair(self.mesh, self.survey)
 
         for freq in self.survey.freqs:
             A = self.getA(freq)
-            rhs = RHS(freq)
+            rhs = self.getRHS(freq)
             Ainv = self.Solver(A, **self.solverOpts)
             sol = Ainv * rhs
             Srcs = self.survey.getSrcByFreq(freq)
diff --git a/simpegEM/FDEM/FieldsFDEM.py b/simpegEM/FDEM/FieldsFDEM.py
index 3120f7a7..1917a379 100644
--- a/simpegEM/FDEM/FieldsFDEM.py
+++ b/simpegEM/FDEM/FieldsFDEM.py
@@ -22,7 +22,9 @@ class FieldsFDEM_e(FieldsFDEM):
         self._edgeCurl = self.survey.prob.mesh.edgeCurl
 
     def _b_sec(self, e, src):
-        return - 1./(1j*omega(src.freq)) * (self._edgeCurl * e)
+        C = self._edgeCurl
+        b_sec =  - 1./(1j*omega(src.freq))*(C * e)
+        return b_sec
 
     def _b_secDeriv(self, e, src, v, adjoint=False): 
         return None