From 5c8fba424251b066d4893a827ea15089e9f0e068 Mon Sep 17 00:00:00 2001
From: Lindsey Heagy <lheagy@eos.ubc.ca>
Date: Sat, 27 Feb 2016 15:44:07 -0800
Subject: [PATCH] changing btwn e,j and h,b depends on which formulation we are
 going between and is now in the fields object

---
 SimPEG/EM/Utils/EMUtils.py  | 33 ---------------------------------
 SimPEG/EM/Utils/__init__.py |  3 ---
 2 files changed, 36 deletions(-)

diff --git a/SimPEG/EM/Utils/EMUtils.py b/SimPEG/EM/Utils/EMUtils.py
index 660ef117..e7dbf441 100644
--- a/SimPEG/EM/Utils/EMUtils.py
+++ b/SimPEG/EM/Utils/EMUtils.py
@@ -13,37 +13,4 @@ def k(freq, sigma, mu=mu_0, eps=epsilon_0):
     beta = w * np.sqrt( mu*eps/2 * ( np.sqrt(1. + (sigma / (eps*w))**2 ) - 1) ) 
     return alp - 1j*beta
 
-# # Constitutive relations
-# def e_from_j(prob,j):
-#     eqLocs = prob._eqLocs
-#     if eqLocs is 'FE':
-#         MSigmaI = prob.MeSigmaI
-#     elif eqLocs is 'EF':
-#         MSigmaI = prob.MfRho
-#     return MSigmaI*j
-
-# def j_from_e(prob,e):
-#     eqLocs = prob._eqLocs
-#     if eqLocs is 'FE':
-#         MSigma = prob.MeSigma
-#     elif eqLocs is 'EF':
-#         MSigma = prob.MfRhoI
-#     return MSigma*e
-
-# def b_from_h(prob,h):
-#     eqLocs = prob._eqLocs
-#     if eqLocs is 'FE':
-#         MMu = prob.MfMuiI
-#     elif eqLocs is 'EF':
-#         MMu = prob.MeMu
-#     return MMu*h
-
-# def h_from_b(prob,b):
-#     eqLocs = prob._eqLocs
-#     if eqLocs is 'FE':
-#         MMuI = prob.MfMui
-#     elif eqLocs is 'EF':
-#         MMuI = prob.MeMuI
-#     return MMuI*b 
-
 
diff --git a/SimPEG/EM/Utils/__init__.py b/SimPEG/EM/Utils/__init__.py
index a60d9701..ef779042 100644
--- a/SimPEG/EM/Utils/__init__.py
+++ b/SimPEG/EM/Utils/__init__.py
@@ -1,5 +1,2 @@
-# import Sources
-# import Ana
-# import Solver
 from EMUtils import omega, k
 from AnalyticUtils import MagneticDipoleFields, MagneticDipoleVectorPotential, MagneticLoopVectorPotential
\ No newline at end of file