Let e[j] be the j-th standard unit basis vector.
CanonicalSimplex(n) gives the simplex with columns being
vertices of the canonical simplex in n-dimensions:
the n-dim. simplex with vertices (0,0,...,0) and e[1],...,e[n]. A vector
(u[1],...,u[n]) is in the canonical simplex if 0