From b6dd8cba5f819b5710177407575e61444cdd2006 Mon Sep 17 00:00:00 2001 From: Michael Orlitzky Date: Thu, 28 Nov 2024 11:43:11 -0500 Subject: [PATCH] mjo/eja/eja_operator.py: simplify an is_isomorphism() test --- mjo/eja/eja_operator.py | 10 ++++------ 1 file changed, 4 insertions(+), 6 deletions(-) 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 -- 2.49.0