# with had entries in a nice field.
return super(MatrixEuclideanJordanAlgebra, self)._charpoly_coeff(i)
else:
- n = self.natural_basis_space().nrows()
+ # If we didn't unembed first, this number would be wrong
+ # by a power-of-two factor for complex/quaternion matrices.
+ n = self.real_unembed(self.natural_basis_space().zero()).nrows()
field = self.base_ring().base_ring() # yeeeeaaaahhh
J = self.__class__(n, field, False)
(_,x,_,_) = J._charpoly_matrix_system()