gamma = (vol / (2* pi)**2* sqrt(pi / alpha) * self.kd.nbzkpts)forG2_Ginpd.G2_qG:ifG2_G[0] <1e-7: G2_G = G2_G[1:] gamma -= np.dot(np.exp(-alpha * G2_G), G2_G**-1)returngamma / self.qstride_c.prod() 開發者ID:ryancoleman,項目名稱:lotsofcoresbook2code...