wassname/simpeg
1
0
Fork 0
mirror of https://github.com/wassname/simpeg.git synced 2026-07-07 07:44:07 +08:00
Code Issues Packages Projects Releases Wiki Activity

More analytic tests.

Browse Source
This commit is contained in:
rowanc1
2014-04-17 11:54:51 -07:00
parent 7b8cbabbb3
commit 44ec85cfaa
1 changed files with 23 additions and 12 deletions
Download Patch File Download Diff File Expand all files Collapse all files
simpegEM/Tests/test_TDEM_forward_Analytic.py
+23 -12
View File
@@ -56,24 +56,35 @@ def halfSpaceProblemAnaDiff(meshType, sig_half=1e-2, rxOffset=50., bounds=[1e-5,
if showIt == True:
plt.loglog(prb.times[bz_calc>0], bz_calc[bz_calc>0], 'r', prb.times[bz_calc<0], -bz_calc[bz_calc<0], 'r--')
plt.loglog(prb.times, abs(bz_ana), 'b*')
plt.title('sig_half = %e'%sig_half)
plt.show()
return log10diff
class TDEM_bTests(unittest.TestCase):
def test_analitic_p2_CYL(self):
self.assertTrue(halfSpaceProblemAnaDiff('CYL',sig_half=1e+2) < 0.01)
def test_analitic_p1_CYL(self):
self.assertTrue(halfSpaceProblemAnaDiff('CYL',sig_half=1e+1) < 0.01)
def test_analitic_p0_CYL(self):
self.assertTrue(halfSpaceProblemAnaDiff('CYL',sig_half=1e+0) < 0.01)
def test_analitic_m1_CYL(self):
self.assertTrue(halfSpaceProblemAnaDiff('CYL',sig_half=1e-1) < 0.01)
def test_analitic_m2_CYL(self):
self.assertTrue(halfSpaceProblemAnaDiff('CYL',sig_half=1e-2) < 0.01)
def test_analitic_m3_CYL(self):
self.assertTrue(halfSpaceProblemAnaDiff('CYL',sig_half=1e-3) < 0.02)
def test_analitic_p2_CYL_50m(self):
self.assertTrue(halfSpaceProblemAnaDiff('CYL', rxOffset=50., sig_half=1e+2) < 0.01)
def test_analitic_p1_CYL_50m(self):
self.assertTrue(halfSpaceProblemAnaDiff('CYL', rxOffset=50., sig_half=1e+1) < 0.01)
def test_analitic_p0_CYL_50m(self):
self.assertTrue(halfSpaceProblemAnaDiff('CYL', rxOffset=50., sig_half=1e+0) < 0.01)
def test_analitic_m1_CYL_50m(self):
self.assertTrue(halfSpaceProblemAnaDiff('CYL', rxOffset=50., sig_half=1e-1) < 0.01)
def test_analitic_m2_CYL_50m(self):
self.assertTrue(halfSpaceProblemAnaDiff('CYL', rxOffset=50., sig_half=1e-2) < 0.01)
def test_analitic_m3_CYL_50m(self):
self.assertTrue(halfSpaceProblemAnaDiff('CYL', rxOffset=50., sig_half=1e-3) < 0.02)
def test_analitic_p0_CYL_1m(self):
self.assertTrue(halfSpaceProblemAnaDiff('CYL', rxOffset=1.0, sig_half=1e+0) < 0.01)
def test_analitic_m1_CYL_1m(self):
self.assertTrue(halfSpaceProblemAnaDiff('CYL', rxOffset=1.0, sig_half=1e-1) < 0.01)
def test_analitic_m2_CYL_1m(self):
self.assertTrue(halfSpaceProblemAnaDiff('CYL', rxOffset=1.0, sig_half=1e-2) < 0.01)
def test_analitic_m3_CYL_1m(self):
self.assertTrue(halfSpaceProblemAnaDiff('CYL', rxOffset=1.0, sig_half=1e-3) < 0.02)