get_umat_el_hub Function

public pure function get_umat_el_hub(i, j, k, l) result(hel)

Arguments

Type IntentOptional Attributes Name
integer, intent(in) :: i
integer, intent(in) :: j
integer, intent(in) :: k
integer, intent(in) :: l

Return Value real(kind=dp)


Source Code

    pure function get_umat_el_hub(i, j, k, l) result(hel)
        integer, intent(in) :: i, j, k, l
        HElement_t(dp) :: hel
#ifdef DEBUG_
        character(*), parameter :: this_routine = "get_umat_el_hub"
#endif

        if (i == j .and. i == k .and. i == l) then
            hel = h_cast(uhub)
        else
            hel = h_cast(0.0_dp)
        end if

        ASSERT(i > 0)
        ASSERT(i <= nbasis / 2)
        ASSERT(j > 0)
        ASSERT(j <= nbasis / 2)
        ASSERT(k > 0)
        ASSERT(k <= nbasis / 2)
        ASSERT(l > 0)
        ASSERT(l <= nbasis / 2)

    end function get_umat_el_hub