cartesian joins

Julia:

sample_vv = [[1, 2], [2, 3, 4], [5, 1], [7, 9]]

product(vv)=foldl((v1, v2)->vec(hcat.(v1, v2')), vv)

product(sample_vv)