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)