# Diff of /trunk/downunder/test/python/run_forward.py

revision 5554 by gross, Mon Mar 23 06:54:53 2015 UTC revision 5555 by caltinay, Tue Mar 24 06:31:43 2015 UTC
# Line 599  class TestMT2DModelTMMode(unittest.TestC Line 599  class TestMT2DModelTMMode(unittest.TestC
599          self.assertRaises(ValueError, MT2DModelTMMode, domain, omega, [(6.7,5)], Z_XY, eta=[1.,1.], w0=[2.,3.], Hx_bottom=complex(4.5,6) )          self.assertRaises(ValueError, MT2DModelTMMode, domain, omega, [(6.7,5)], Z_XY, eta=[1.,1.], w0=[2.,3.], Hx_bottom=complex(4.5,6) )
600
601      def test_PDE(self):      def test_PDE(self):
602          omega=1.          omega=10.
603          mu0=0.123          mu0=0.123
604          RHO=0.15          RHO=0.15
605          k=cmath.sqrt(1j*omega*mu0/RHO)  # Hx=exp(k*z)          k=cmath.sqrt(1j*omega*mu0/RHO)  # Hx=exp(k*z)
# Line 682  class TestMT2DModelTMMode(unittest.TestC Line 682  class TestMT2DModelTMMode(unittest.TestC
682          acw=MT2DModelTMMode(domain, omega, x, Z_XY, eta, mu=mu0, tol=1e-9,  directSolver=True)          acw=MT2DModelTMMode(domain, omega, x, Z_XY, eta, mu=mu0, tol=1e-9,  directSolver=True)
683
684          # this is the base line:          # this is the base line:
685          RHO0=1. # was 100.15          RHO0=2. # was 100.15
686          args0=acw.getArguments(RHO0)          args0=acw.getArguments(RHO0)
687          d0=acw.getDefect(RHO0, *args0)          d0=acw.getDefect(RHO0, *args0)