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