calc_pgen_mol_guga_single_guga Function

    function calc_pgen_mol_guga_single_guga(ilutI, nI, ilutJ, nJ, excitInfo) result(pgen)
        ! this function calculates the guga branching part of the
        ! generation probability!
        integer(n_int), intent(in) :: ilutI(0:niftot), ilutJ(0:niftot)
        integer, intent(in) :: nI(nel), nJ(nel)
        type(ExcitationInformation_t), intent(in) :: excitInfo
        real(dp) :: pgen
        character(*), parameter :: this_routine = "calc_pgen_mol_guga_single_guga"

        call stop_all(this_routine, "TODO")

        ! todo
        unused_var(ilutI)
        unused_var(nI)
        unused_var(ilutJ)
        unused_var(nJ)
        unused_var(excitInfo)

        pgen = 0.0_dp

    end function calc_pgen_mol_guga_single_guga