calcStepvector Function

    pure function calcStepvector(ilut) result(stepVector)
        ! function to calculate stepvector of length nReps, corresponding
        ! to the ilut bit-representation, if each stepvalue is needed
        ! often within a function.
        ! there is probably a very efficient way of programming that!
        !TODO ask simon for improvements.
        integer(n_int), intent(in) :: ilut(0:GugaBits%len_orb)
        integer :: stepVector(nSpatOrbs)
        integer :: iOrb

        do iOrb = 1, nSpatOrbs
            stepVector(iOrb) = getStepvalue(ilut, iOrb)
        end do

    end function calcStepvector