From: Michael Orlitzky Date: Thu, 28 Nov 2024 16:43:11 +0000 (-0500) Subject: mjo/eja/eja_operator.py: simplify an is_isomorphism() test X-Git-Url: https://gitweb.michael.orlitzky.com/?a=commitdiff_plain;h=b6dd8cba5f819b5710177407575e61444cdd2006;p=sage.d.git mjo/eja/eja_operator.py: simplify an is_isomorphism() test --- diff --git a/mjo/eja/eja_operator.py b/mjo/eja/eja_operator.py index be22c4a..c2a8d9c 100644 --- a/mjo/eja/eja_operator.py +++ b/mjo/eja/eja_operator.py @@ -723,7 +723,6 @@ class EJAOperator(Map): SETUP:: - sage: from mjo.basis_repr import basis_repr_of_operator sage: from mjo.eja.eja_algebra import (JordanSpinEJA, ....: HadamardEJA, ....: RealSymmetricEJA, @@ -776,11 +775,10 @@ class EJAOperator(Map): 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