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. Implement spectral projector decomposition for EJA operators
+ using jordan_form() or eigenmatrix_right(). I suppose we can
+ ignore the problem of base rings for now and just let it crash
+ if we're not using AA as our base field.
-7. If we factor out a "matrix algebra" class, then it would make sense
- to replace the custom embedding/unembedding functions with static
- _real_embedding() and _real_unembedding() methods.
+7. Do we really need to orthonormalize the basis in a subalgebra?
+ So long as we can decompose the operator (which is invariant
+ under changes of basis), who cares?