Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
integer, | intent(in) | :: | a | |||
integer, | intent(in) | :: | b |
function symAllowed(a, b) result(allowed)
! Check if a transition from a to be is symmetry-allowed
! Input: a,b - orbitals to check
! Output: allowed - true if and only if a,b have the same symmetries
implicit none
integer, intent(in) :: a, b
logical :: allowed
allowed = same_spin(a, b) .and. (G1(a)%Sym%s == G1(b)%Sym%s)
end function symAllowed