X-Git-Url: http://gitweb.michael.orlitzky.com/?a=blobdiff_plain;f=mjo%2Feja%2Feja_algebra.py;h=4947fe82bb981401a09a039306516b4b0631113a;hb=f911ca107b6ad2a92547a255311ef7f16978feac;hp=d0e7b074ce41a00f66d6f6dcd487653f3a8b1674;hpb=11e681d6320f0b7ddbb834931845b6f4a745da93;p=sage.d.git

diff --git a/mjo/eja/eja_algebra.py b/mjo/eja/eja_algebra.py
index d0e7b07..4947fe8 100644
--- a/mjo/eja/eja_algebra.py
+++ b/mjo/eja/eja_algebra.py
@@ -60,6 +60,9 @@ class FiniteDimensionalEuclideanJordanAlgebra(CombinatorialFreeModule):
         self._rank = rank
         self._natural_basis = natural_basis
 
+        # TODO: HACK for the charpoly.. needs redesign badly.
+        self._basis_normalizers = None
+
         if category is None:
             category = MagmaticAlgebras(field).FiniteDimensional()
             category = category.WithBasis().Unital()
@@ -224,6 +227,7 @@ class FiniteDimensionalEuclideanJordanAlgebra(CombinatorialFreeModule):
         return V.span_of_basis(b)
 
 
+
     @cached_method
     def _charpoly_coeff(self, i):
         """
@@ -234,6 +238,19 @@ class FiniteDimensionalEuclideanJordanAlgebra(CombinatorialFreeModule):
         store the trace/determinant (a_{r-1} and a_{0} respectively)
         separate from the entire characteristic polynomial.
         """
+        if self._basis_normalizers is not None:
+             # Must be a matrix class?
+             # WARNING/TODO: this whole mess is mis-designed.
+             n = self.natural_basis_space().nrows()
+             field = self.base_ring().base_ring() # yeeeeaaaahhh
+             J = self.__class__(n, field, False)
+             (_,x,_,_) = J._charpoly_matrix_system()
+             p = J._charpoly_coeff(i)
+             # p might be missing some vars, have to substitute "optionally"
+             pairs = zip(x.base_ring().gens(), self._basis_normalizers)
+             substitutions = { v: v*c for (v,c) in pairs }
+             return p.subs(substitutions)
+
         (A_of_x, x, xr, detA) = self._charpoly_matrix_system()
         R = A_of_x.base_ring()
         if i >= self.rank():
@@ -924,7 +941,7 @@ def _complex_hermitian_basis(n, field):
 
 
 
-def _quaternion_hermitian_basis(n, field, normalize):
+def _quaternion_hermitian_basis(n, field):
     """
     Returns a basis for the space of quaternion Hermitian n-by-n matrices.
 
@@ -942,7 +959,7 @@ def _quaternion_hermitian_basis(n, field, normalize):
 
         sage: set_random_seed()
         sage: n = ZZ.random_element(1,5)
-        sage: B = _quaternion_hermitian_basis(n, QQ, False)
+        sage: B = _quaternion_hermitian_basis(n, QQ)
         sage: all( M.is_symmetric() for M in B )
         True
 
@@ -1339,10 +1356,10 @@ class RealSymmetricEJA(FiniteDimensionalEuclideanJordanAlgebra):
             if p.is_irreducible():
                 field = NumberField(p, 'sqrt2', embedding=RLF(2).sqrt())
             S = [ s.change_ring(field) for s in S ]
-            self._basis_denormalizers = tuple(
-                self.__class__.natural_inner_product(s,s).sqrt()
+            self._basis_normalizers = tuple(
+                ~(self.__class__.natural_inner_product(s,s).sqrt())
                 for s in S )
-            S = tuple( s/c for (s,c) in zip(S,self._basis_denormalizers) )
+            S = tuple( s*c for (s,c) in zip(S,self._basis_normalizers) )
 
         Qs = _multiplication_table_from_matrix_basis(S)
 
@@ -1439,10 +1456,10 @@ class ComplexHermitianEJA(FiniteDimensionalEuclideanJordanAlgebra):
             if p.is_irreducible():
                 field = NumberField(p, 'sqrt2', embedding=RLF(2).sqrt())
             S = [ s.change_ring(field) for s in S ]
-            self._basis_denormalizers = tuple(
-                self.__class__.natural_inner_product(s,s).sqrt()
+            self._basis_normalizers = tuple(
+                ~(self.__class__.natural_inner_product(s,s).sqrt())
                 for s in S )
-            S = tuple( s/c for (s,c) in zip(S,self._basis_denormalizers) )
+            S = tuple( s*c for (s,c) in zip(S,self._basis_normalizers) )
 
         Qs = _multiplication_table_from_matrix_basis(S)