symAllowed Function

public function symAllowed(a, b) result(allowed)

Arguments

Type IntentOptional Attributes Name
integer, intent(in) :: a
integer, intent(in) :: b

Return Value logical


Source Code

    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