powerset function

Compute all possible combinations of a vector - ignoring the empty set