A = OctonionMatrixAlgebra(n, scalars=field)
super().__init__(A, **kwargs)
+ if n == 3:
+ from mjo.eja.eja_cache import albert_eja_coeffs
+ a = albert_eja_coeffs(self.coordinate_polynomial_ring())
+ if self._rational_algebra is None:
+ self._charpoly_coefficients.set_cache(a)
+ else:
+ self._rational_algebra._charpoly_coefficients.set_cache(a)
+
class AlbertEJA(OctonionHermitianEJA):
r"""