change i^-n to i^n in cylindrical radial function so that k-vector points the incoming direction

Author: narahahn

examples/modal_beamforming_open_circular_array.py

@@ -26,7 +26,7 @@
 D = micarray.modal.radial.circ_diagonal_mode_mat(Dn)
 
 # compute microphone signals for an incident broad-band plane wave
-p = np.exp(-1j * k[:, np.newaxis]*r * np.cos(pol - pw_angle))
+p = np.exp(1j * k[:, np.newaxis]*r * np.cos(pol - pw_angle))
 # compute plane wave decomposition
 A_pwd = np.matmul(np.matmul(np.conj(Psi_q.T), D), Psi_p)
 q_pwd = np.squeeze(np.matmul(A_pwd, np.expand_dims(p, 2)))

micarray/modal/radial.py

@@ -260,7 +260,7 @@ def circular_pw(N, k, r, setup):
 
     bn = circ_radial_weights(N, kr, setup)
     for i, x in enumerate(kr):
-        bn[i, :] = bn[i, :] * (1j)**(-n)
+        bn[i, :] = bn[i, :] * (1j)**(n)
     return bn