reals with an orthonormal basis::
sage: set_random_seed()
- sage: x = random_eja(AA).random_element()
+ sage: x = random_eja().random_element()
sage: A = x.subalgebra_generated_by(orthonormalize_basis=True)
sage: x = A.random_element()
sage: A.one()*x == x and x*A.one() == x
the algebraic reals with an orthonormal basis::
sage: set_random_seed()
- sage: x = random_eja(AA).random_element()
+ sage: x = random_eja().random_element()
sage: A = x.subalgebra_generated_by(orthonormalize_basis=True)
sage: actual = A.one().operator().matrix()
sage: expected = matrix.identity(A.base_ring(), A.dimension())