[Fortran] 纯文本查看 复制代码
subroutine integralk(numk, numw, w, gamma, sigma)
use ABA
implicit none
integer, intent(in) :: numw, numk
real(8),intent(in) :: w(numw)
real(8), intent(in) :: gamma
complex(8),intent(out) ::sigma(2, 2, numw)
integer :: N = 6, i1, i3, i4
real(8) :: k(2), en(6), ds
complex(8) :: t_sigma(2, 2, numw) , sigmak(2, 2, numw)
complex(8) :: H(6,6), dagger_H(6,6), v(2, 6, 6), r(2, 6, 6)
ds = 8d0*(PI**2)/((a0**2)*sqrt(3.0)*(numk**2))
t_sigma = 0.d0
do i3 = 1,numk
do i4 = 1,numk
k(1) = real(i3)/numk, k(2) = real(i4)/numk