genidxthc.f90

Go to the documentation of this file.

! Copyright (C) 2024 J. K. Dewhurst and S. Sharma.
! This file is distributed under the terms of the GNU General Public License.
! See the file COPYING for license details.

subroutine genidxthc
use modmain
use modtdhfc
implicit none
! local variables
integer ik,jk,nst,nmax,ist
! determine the maximum number of states within energy window over all k-points
nmax=0
do ik=1,nkpt
  nst=0
  do ist=1,nstsv
    if (abs(evalsv(ist,ik)-efermi) < ecutthc) nst=nst+1
  end do
  nmax=max(nmax,nst)
end do
if (nmax == 0) then
  write(*,*)
  write(*,'("Error(genidxthc): no states within energy window ecutthc")')
  write(*,*)
  stop
end if
! allocate global arrays
if (allocated(nthck)) deallocate(nthck)
allocate(nthck(nkpt))
if (allocated(idxthc)) deallocate(idxthc)
allocate(idxthc(nmax,nkpt))
if (allocated(istthc)) deallocate(istthc)
allocate(istthc(nmax,nkptnr))
! determine the number of and index to used states
do ik=1,nkpt
  nst=0
  do ist=1,nstsv
    if (abs(evalsv(ist,ik)-efermi) < ecutthc) then
      nst=nst+1
      idxthc(nst,ik)=ist
    end if
  end do
  nthck(ik)=nst
end do
! calculate the index to and total number of TDHFC states
nthc=0
do ik=1,nkptnr
  jk=ivkik(ivk(1,ik),ivk(2,ik),ivk(3,ik))
  do ist=1,nthck(jk)
    nthc=nthc+1
    istthc(ist,ik)=nthc
  end do
end do
end subroutine