From fc444a345f9d4c62071f59e8d9a9cd25aba93ba4 Mon Sep 17 00:00:00 2001 From: GudniRos Date: Wed, 9 Mar 2016 14:24:15 -0800 Subject: [PATCH] Fixing dtypes in the MT1Danalytic --- SimPEG/MT/Utils/MT1Danalytic.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/SimPEG/MT/Utils/MT1Danalytic.py b/SimPEG/MT/Utils/MT1Danalytic.py index a28777cb..6083aa9e 100644 --- a/SimPEG/MT/Utils/MT1Danalytic.py +++ b/SimPEG/MT/Utils/MT1Danalytic.py @@ -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)]])