]> gitweb.michael.orlitzky.com - sage.d.git/commitdiff
mjo/eja/eja_operator.py: simplify an is_isomorphism() test
authorMichael Orlitzky <michael@orlitzky.com>
Thu, 28 Nov 2024 16:43:11 +0000 (11:43 -0500)
committerMichael Orlitzky <michael@orlitzky.com>
Thu, 28 Nov 2024 16:43:11 +0000 (11:43 -0500)
mjo/eja/eja_operator.py

index be22c4a88b438b10454b2a2992f9315b7ee3affb..c2a8d9c782ac5b5e10d7ba6f7c8a1089631cddb1 100644 (file)
@@ -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