X-Git-Url: http://gitweb.michael.orlitzky.com/?a=blobdiff_plain;f=mjo%2Feja%2Feuclidean_jordan_algebra.py;h=1426d5e16be4b6acc7c68a4494e2c6f1c4d61819;hb=4156fccce1265f500fe432b0f5567e43fbbc23d6;hp=747e198f8922e3d6c195c4c7552829526463d405;hpb=ec37d0e9ac40f4345da8e52a1b3dd35599e161fb;p=sage.d.git

diff --git a/mjo/eja/euclidean_jordan_algebra.py b/mjo/eja/euclidean_jordan_algebra.py
index 747e198..1426d5e 100644
--- a/mjo/eja/euclidean_jordan_algebra.py
+++ b/mjo/eja/euclidean_jordan_algebra.py
@@ -124,6 +124,18 @@ class FiniteDimensionalEuclideanJordanAlgebra(FiniteDimensionalAlgebra):
                 sage: (x*x)*(x*x*x) == x^5
                 True
 
+            We also know that powers operator-commute (Koecher, Chapter
+            III, Corollary 1)::
+
+                sage: set_random_seed()
+                sage: x = random_eja().random_element()
+                sage: m = ZZ.random_element(0,10)
+                sage: n = ZZ.random_element(0,10)
+                sage: Lxm = (x^m).matrix()
+                sage: Lxn = (x^n).matrix()
+                sage: Lxm*Lxn == Lxn*Lxm
+                True
+
             """
             A = self.parent()
             if n == 0:
@@ -188,7 +200,7 @@ class FiniteDimensionalEuclideanJordanAlgebra(FiniteDimensionalAlgebra):
             Example 11.11::
 
                 sage: set_random_seed()
-                sage: n = ZZ.random_element(1,10).abs()
+                sage: n = ZZ.random_element(1,10)
                 sage: J = JordanSpinSimpleEJA(n)
                 sage: x = J.random_element()
                 sage: while x.is_zero():
@@ -355,7 +367,7 @@ class FiniteDimensionalEuclideanJordanAlgebra(FiniteDimensionalAlgebra):
             aren't multiples of the identity are regular::
 
                 sage: set_random_seed()
-                sage: n = ZZ.random_element(1,10).abs()
+                sage: n = ZZ.random_element(1,10)
                 sage: J = JordanSpinSimpleEJA(n)
                 sage: x = J.random_element()
                 sage: x == x.coefficient(0)*J.one() or x.degree() == 2
@@ -454,7 +466,7 @@ class FiniteDimensionalEuclideanJordanAlgebra(FiniteDimensionalAlgebra):
             identity::
 
                 sage: set_random_seed()
-                sage: n = ZZ.random_element(2,10).abs()
+                sage: n = ZZ.random_element(2,10)
                 sage: J = JordanSpinSimpleEJA(n)
                 sage: y = J.random_element()
                 sage: while y == y.coefficient(0)*J.one():
@@ -499,7 +511,7 @@ class FiniteDimensionalEuclideanJordanAlgebra(FiniteDimensionalAlgebra):
             Alizadeh's Example 11.12::
 
                 sage: set_random_seed()
-                sage: n = ZZ.random_element(1,10).abs()
+                sage: n = ZZ.random_element(1,10)
                 sage: J = JordanSpinSimpleEJA(n)
                 sage: x = J.random_element()
                 sage: x_vec = x.vector()
@@ -549,7 +561,7 @@ class FiniteDimensionalEuclideanJordanAlgebra(FiniteDimensionalAlgebra):
 
             Property 6:
 
-                sage: k = ZZ.random_element(1,10).abs()
+                sage: k = ZZ.random_element(1,10)
                 sage: actual = (x^k).quadratic_representation()
                 sage: expected = (x.quadratic_representation())^k
                 sage: actual == expected
@@ -790,7 +802,7 @@ def random_eja():
         Euclidean Jordan algebra of degree...
 
     """
-    n = ZZ.random_element(1,5).abs()
+    n = ZZ.random_element(1,5)
     constructor = choice([eja_rn,
                           JordanSpinSimpleEJA,
                           RealSymmetricSimpleEJA,
@@ -825,7 +837,7 @@ def _complex_hermitian_basis(n, field=QQ):
     TESTS::
 
         sage: set_random_seed()
-        sage: n = ZZ.random_element(1,5).abs()
+        sage: n = ZZ.random_element(1,5)
         sage: all( M.is_symmetric() for M in _complex_hermitian_basis(n) )
         True
 
@@ -1002,7 +1014,7 @@ def RealSymmetricSimpleEJA(n, field=QQ):
     The degree of this algebra is `(n^2 + n) / 2`::
 
         sage: set_random_seed()
-        sage: n = ZZ.random_element(1,5).abs()
+        sage: n = ZZ.random_element(1,5)
         sage: J = RealSymmetricSimpleEJA(n)
         sage: J.degree() == (n^2 + n)/2
         True
@@ -1026,7 +1038,7 @@ def ComplexHermitianSimpleEJA(n, field=QQ):
     The degree of this algebra is `n^2`::
 
         sage: set_random_seed()
-        sage: n = ZZ.random_element(1,5).abs()
+        sage: n = ZZ.random_element(1,5)
         sage: J = ComplexHermitianSimpleEJA(n)
         sage: J.degree() == n^2
         True