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

Fixing dtypes in the MT1Danalytic

Browse Source
This commit is contained in:
GudniRos
2016-03-09 14:24:15 -08:00
parent 9a739d8380
commit fc444a345f
1 changed files with 2 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files
SimPEG/MT/Utils/MT1Danalytic.py
+2 -2
View File
@@ -38,9 +38,9 @@ def getEHfields(m1d,sigma,freq,zd,scaleUD=True):
# Build the propagation matrix
# Convert fields to down/up going components in layer below current layer
Pj1 = np.array([[1,1],[yp1,-yp1]])
Pj1 = np.array([[1,1],[yp1[0],-yp1[0]]],dtype=complex)
# Convert fields to down/up going components in current layer
Pjinv = 1./2*np.array([[1,zp],[1,-zp]])
Pjinv = 1./2*np.array([[1,zp[0]],[1,-zp[0]]],dtype=complex)
# Propagate down and up components through the current layer
elamh = np.array([[np.exp(-1j*k[lnr+1]*h),0],[0,np.exp(1j*k[lnr+1]*h)]])