mirror of
https://github.com/wassname/simpeg.git
synced 2026-06-28 00:30:01 +08:00
Rowan Cockett
558e8879fe
Merge branch 'move2Simpeg' of https://github.com/simpeg/simpegmt into mt/dev Conflicts: .gitignore .travis.yml LICENSE docs/conf.py docs/index.rst requirements.txt setup.py
49 lines
1.3 KiB
Python
49 lines
1.3 KiB
Python
import unittest
|
|
from SimPEG import *
|
|
from SimPEG import MT
|
|
|
|
TOL = 1e-6
|
|
|
|
def appResPhs(freq,z):
|
|
app_res = ((1./(8e-7*np.pi**2))/freq)*np.abs(z)**2
|
|
app_phs = np.arctan2(-z.imag,z.real)*(180/np.pi)
|
|
return app_res, app_phs
|
|
|
|
def appResNorm(sigmaHalf):
|
|
nFreq = 26
|
|
|
|
m1d = Mesh.TensorMesh([[(100,5,1.5),(100.,10),(100,5,1.5)]], x0=['C'])
|
|
sigma = np.zeros(m1d.nC) + sigmaHalf
|
|
sigma[m1d.gridCC[:]>200] = 1e-8
|
|
|
|
# Calculate the analytic fields
|
|
freqs = np.logspace(4,-4,nFreq)
|
|
Z = []
|
|
for freq in freqs:
|
|
Ed, Eu, Hd, Hu = MT.Utils.getEHfields(m1d,sigma,freq,np.array([200]))
|
|
Z.append((Ed + Eu)/(Hd + Hu))
|
|
|
|
Zarr = np.concatenate(Z)
|
|
|
|
app_r, app_p = appResPhs(freqs,Zarr)
|
|
|
|
return np.linalg.norm(np.abs(app_r - np.ones(nFreq)/sigmaHalf)) / np.log10(sigmaHalf)
|
|
|
|
|
|
class TestAnalytics(unittest.TestCase):
|
|
|
|
def setUp(self):
|
|
pass
|
|
def test_appRes2en1(self):self.assertLess(appResNorm(2e-1), TOL)
|
|
def test_appRes2en2(self):self.assertLess(appResNorm(2e-2), TOL)
|
|
def test_appRes2en3(self):self.assertLess(appResNorm(2e-3), TOL)
|
|
def test_appRes2en4(self):self.assertLess(appResNorm(2e-4), TOL)
|
|
def test_appRes2en5(self):self.assertLess(appResNorm(2e-5), TOL)
|
|
def test_appRes2en6(self):self.assertLess(appResNorm(2e-6), TOL)
|
|
|
|
|
|
|
|
|
|
if __name__ == '__main__':
|
|
unittest.main()
|