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

smaller test mesh

Browse Source
This commit is contained in:
Lindsey Heagy
2015-12-08 21:48:01 -08:00
parent adda7a43dc
commit b620c4286b
1 changed files with 7 additions and 7 deletions
Download Patch File Download Diff File Expand all files Collapse all files
SimPEG/EM/Utils/testingUtils.py
+7 -7
View File
@@ -6,17 +6,17 @@ from scipy.constants import mu_0
def getFDEMProblem(fdemType, comp, SrcList, freq, verbose=False):
cs = 10.
ncx, ncy, ncz = 2, 2, 2
ncx, ncy, ncz = 0, 0, 0
npad = 8
hx = [(cs,npad,-1.5), (cs,ncx), (cs,npad,1.5)]
hy = [(cs,npad,-1.5), (cs,ncy), (cs,npad,1.5)]
hz = [(cs,npad,-1.5), (cs,ncz), (cs,npad,1.5)]
hx = [(cs,npad,-1.3), (cs,ncx), (cs,npad,1.3)]
hy = [(cs,npad,-1.3), (cs,ncy), (cs,npad,1.3)]
hz = [(cs,npad,-1.3), (cs,ncz), (cs,npad,1.3)]
mesh = Mesh.TensorMesh([hx,hy,hz],['C','C','C'])
mapping = Maps.ExpMap(mesh)
x = np.array([np.linspace(-5.5*cs,-2.5*cs,4),np.linspace(3.5*cs,2.5*cs,4)]) #don't sample right by the source
XYZ = Utils.ndgrid(x,x,np.linspace(-2.25*cs,2.25*cs,5))
x = np.array([np.linspace(-5.*cs,-2.*cs,3),np.linspace(5.*cs,2.*cs,3)]) #don't sample right by the source
XYZ = Utils.ndgrid(x,x,np.linspace(-2.*cs,2.*cs,5))
Rx0 = EM.FDEM.Rx(XYZ, comp)
Src = []
@@ -70,6 +70,6 @@ def getFDEMProblem(fdemType, comp, SrcList, freq, verbose=False):
from pymatsolver import MumpsSolver
prb.Solver = MumpsSolver
except ImportError, e:
pass
prb.Solver = SolverLU
return prb