def epsfun(v):
om = v
omsi = -ky
omdi = 2 * omsi * LnbyR
za = -om / omdi
zb = -kpar / omdi * np.sqrt(2)
b = ky**2
i10 = gp.Inm(za, zb, b, 1, 0)
i12 = gp.Inm(za, zb, b, 1, 2)
i30 = gp.Inm(za, zb, b, 3, 0)
eps = 1 + 1 / tau + (i10 * (om - omsi * (1 - 1.5 * etai)) - omsi * etai *
(i12 + i30)) / omdi
return eps
def epsfun(v):
om=v[0]+1j*v[1]
omsi=-ky
omdi=2*omsi*LnbyR
za=-om/omdi
zb=-np.sqrt(2)*kpar/omdi
b=ky**2
i10=gp.Inm(za,zb,b,1,0)
i12=gp.Inm(za,zb,b,1,2)
i30=gp.Inm(za,zb,b,3,0)
eps=1+1/tau+(i10*(om-omsi*(1-1.5*etai))-omsi*etai*(i12+i30))/omdi
res=[float(np.real(eps)),float(np.imag(eps))]
return res
wdts[0] = 2.0
wdts[5] = 2.0
wdts[10] = 4.0
wdts[15] = 2.0
wdts[20] = 2.0
vcb = np.arange(-12, 14, 2)
ii = 0
xx, yy = np.meshgrid(np.arange(-6, 6, 0.1), np.arange(-6, 6, 0.1))
za = xx + 1j * yy
zb = 0.0
b = 0.09
for ns in nlist:
ii = ii + 1
[n, m] = ns
print('computing I' + str(n) + str(m) + ' ...\n')
inm = gp.Inm(za, zb, b, n, m)
plt.subplot(3, 3, ii)
plt.pcolormesh(np.real(za),
np.imag(za),
np.imag(inm),
shading='gouraud',
rasterized=True)
plt.clim(-10, 10)
cs = plt.contour(np.real(za),
np.imag(za),
np.imag(inm),
cnts,
colors='k',
linewidths=wdts * 0.5)
aa = r'$Im[I_{' + str(n) + str(m) + r'}(\zeta_\alpha,\zeta_\beta,b)]$'
plt.text(-1, 3, aa, fontsize=14)