Generates MOLS by cyclic permuation of a basic Latin square L0 constructed
from a vector of ordered elements X of a prime-power finite field of size
p**q and a unit vector 1 of length p**q (see Chapter 1 of Raghavarao 1971).
\(L0 = crossprod(X,t(1)) + crosssprod(1,t(X))\)
The primitive polynomials for the MOLS generated by this package were
extracted from the Table of Primitive Polynomials given in the Supplement to Hansen and Mullen (1992).
The output is a single data frame for a \(p**q x p**q\) square classified by rows and columns
with a separate column for the allocation of treatments to each individual square.