QPMS
/
qpms
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

Fix variable names in scatter_plane_wave_rectarray 

Former-commit-id: 1766eab4d6fd6bd98ecba6a188a25c89a58e5561
This commit is contained in:
Marek Nečada 2016-07-18 11:46:48 +03:00
parent 550901a18a
commit ab6c9fc267
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
qpms/qpms_p.py
View File

@ -1073,7 +1073,7 @@ def scatter_plane_wave_rectarray(omega, epsilon_b, xN, yN, xd, yd, TMatrices, k_
if ((1 == k_dirs.ndim) and (1 == E_0s.ndim)): if ((1 == k_dirs.ndim) and (1 == E_0s.ndim)):
k_cart = k_dirs * k_out # wave vector of the incident plane wave k_cart = k_dirs * k_out # wave vector of the incident plane wave
pq_0 = np.zeros((N,2,nelem), dtype=np.complex_) pq_0 = np.zeros((N,2,nelem), dtype=np.complex_)
p_y0, q_y0 = plane_pq_y(nmax, k_cart, E_0s) p_y0, q_y0 = plane_pq_y(lMax, k_cart, E_0s)
pq_0[:,0,:] = np.exp(1j*np.sum(k_cart[ň,:]*cart_lattice,axis=-1))[:, ň] * p_y0[ň, :] pq_0[:,0,:] = np.exp(1j*np.sum(k_cart[ň,:]*cart_lattice,axis=-1))[:, ň] * p_y0[ň, :]
pq_0[:,1,:] = np.exp(1j*np.sum(k_cart[ň,:]*cart_lattice,axis=-1))[:, ň] * q_y0[ň, :] pq_0[:,1,:] = np.exp(1j*np.sum(k_cart[ň,:]*cart_lattice,axis=-1))[:, ň] * q_y0[ň, :]
if (return_pq_0): if (return_pq_0):
@ -1105,7 +1105,7 @@ def scatter_plane_wave_rectarray(omega, epsilon_b, xN, yN, xd, yd, TMatrices, k_
for ki in range(K): for ki in range(K):
k_cart = k_dirs[ki] * k_out k_cart = k_dirs[ki] * k_out
pq_0 = np.zeros((N,2,nelem), dtype=np.complex_) pq_0 = np.zeros((N,2,nelem), dtype=np.complex_)
p_y0, q_y0 = plane_pq_y(nmax, k_cart, E_0s[ki]) p_y0, q_y0 = plane_pq_y(lMax, k_cart, E_0s[ki])
pq_0[:,0,:] = np.exp(1j*np.sum(k_cart[ň,:]*cart_lattice,axis=-1))[:, ň] * p_y0[ň, :] pq_0[:,0,:] = np.exp(1j*np.sum(k_cart[ň,:]*cart_lattice,axis=-1))[:, ň] * p_y0[ň, :]
pq_0[:,1,:] = np.exp(1j*np.sum(k_cart[ň,:]*cart_lattice,axis=-1))[:, ň] * q_y0[ň, :] pq_0[:,1,:] = np.exp(1j*np.sum(k_cart[ň,:]*cart_lattice,axis=-1))[:, ň] * q_y0[ň, :]
if (return_pq_0): if (return_pq_0):