SETUP::
- sage: from mjo.basis_repr import basis_repr_of_operator
sage: from mjo.eja.eja_algebra import (JordanSpinEJA,
....: HadamardEJA,
....: RealSymmetricEJA,
sage: J = RealSymmetricEJA.random_instance()
sage: n = J.one().to_matrix().nrows()
sage: U = random_unitary_matrix(J.base_ring(), n)
- sage: f = lambda X: U.conjugate_transpose()*X*U
- sage: if n == 0:
- ....: M = matrix.zero(J.base_ring(), 0)
- ....: else:
- ....: M = basis_repr_of_operator(J.matrix_basis(), f)
+ sage: L = lambda X: J(U.conjugate_transpose()*X*U)
+ sage: columns = ( L(b).to_vector() for b in J.matrix_basis() )
+ sage: MS = MatrixSpace(J.base_ring(), J.dimension(), J.dimension())
+ sage: M = MS(columns).transpose()
sage: L = EJAOperator(J,J,M)
sage: L.is_isomorphism()
True