+0. Add tests for orthogonality in the Peirce decomposition.
+
1. Add CartesianProductEJA.
2. Check the axioms in the constructor when check != False?
5. Factor out the unit-norm basis (and operator symmetry) tests once
all of the algebras pass.
-6. Refactor the current ungodly fast charpoly hack (relies on the
- theory to ensure that the charpolys are equal.)
+6. Can we make the minimal and characteristic polynomial tests work
+ for trivial algebras, too? Then we wouldn't need the "nontrivial"
+ argument to random_eja(). This would also allow us to potentially
+ pass a "size" of zero to the other constructors in random_eja() --
+ a code path that is basically untested right now.