格林函数数值计算
假設(shè)能譜已經(jīng)算出來(lái),那么算格林函數(shù)是非常容易的
可以直接用向量化編程實(shí)現(xiàn),這樣的計(jì)算效率最高。
nk = 1024 nw = 1024 delta = 0.01 ks = np.linspace(-pi,pi,nk) omega = np.linspace(-2,2,nw) Ks, Omega = np.meshgrid(ks,omega) Ek = -2*np.cos(Ks) Gk = np.zeros((nw,nk),dtype='double') Gk = np.power((Omega+1j*delta-Ek),-1) Ak = -Gk.imag/pi plt.pcolormesh(Ks, Omega, Ak) plt.colorbar()這里積分只對(duì)k進(jìn)行,我們令\omga = 0,這個(gè)公式就可以寫(xiě)成兩個(gè)格林函數(shù)的卷積
如何用卷積的方式解決這個(gè)計(jì)算。
chi0 = signal.fftconvolve(Gk,Gk,'same')/nkx?
總結(jié)
- 上一篇: 子网掩码,CIDR前缀法表示掩码
- 下一篇: 多重引导