Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
integer(kind=n_int), | intent(in) | :: | ilut(0:GugaBits%len_orb) | |||
integer, | intent(in) | :: | s |
pure function getSpatialOccupation(iLut, s) result(nOcc)
integer(n_int), intent(in) :: ilut(0:GugaBits%len_orb)
integer, intent(in) :: s
real(dp) :: nOcc
if (isZero(ilut, s)) then
nOcc = 0.0_dp
else if (isThree(ilut, s)) then
nOcc = 2.0_dp
else
nOcc = 1.0_dp
end if
end function getSpatialOccupation