EncodeBitDet Subroutine

    pure subroutine EncodeBitDet(nI, iLut)
        integer, intent(in) :: nI(:)
        integer(kind=n_int), intent(out) :: iLut(0:NIfTot)
        integer :: i, pos

        iLut(:) = 0_n_int

        !Decode determinant
        do i = 1, size(nI)
            pos = (nI(i) - 1) / bits_n_int
            iLut(pos) = ibset(iLut(pos), mod(nI(i) - 1, bits_n_int))
        end do

    end subroutine EncodeBitDet