A randomized algorithm for computing the affine dimension of a join of toric varieties using Terracini's Lemma.
Each input matrix defines a parameterization of the variety. For each variety, a vector of parameter values is chosen at random from a large finite field. The dimension of the sum of the tangent spaces at those points is computed.
This algorithm is much much faster than computing the join variety.




All input matrices must have the same number of columns.
The object toricJoinDim is a method function.