eja: begin refactoring to allow noncanonical inner products.

[sage.d.git] / mjo / eja / eja_subalgebra.py

diff --git a/mjo/eja/eja_subalgebra.py b/mjo/eja/eja_subalgebra.py
index 110b049573bbc8d8aeaaab768e9cebceadbf3c7e..6a9d10f65b7164627394c5ddb32bd682c323c5b2 100644 (file)
--- a/mjo/eja/eja_subalgebra.py
+++ b/mjo/eja/eja_subalgebra.py
@@ -177,6 +177,10 @@ class FiniteDimensionalEuclideanJordanSubalgebra(FiniteDimensionalEuclideanJorda
 
         n = len(basis)
         mult_table = [[W.zero() for i in range(n)] for j in range(n)]
+        ip_table = [ [ self._superalgebra.inner_product(basis[i],basis[j])
+                       for i in range(n) ]
+                     for j in range(n) ]
+
         for i in range(n):
             for j in range(n):
                 product = basis[i]*basis[j]
@@ -187,6 +191,7 @@ class FiniteDimensionalEuclideanJordanSubalgebra(FiniteDimensionalEuclideanJorda
                 product_vector = V.from_vector(product.to_vector())
                 mult_table[i][j] = W.coordinate_vector(product_vector)
 
+        self._inner_product_matrix = matrix(field, ip_table)
         natural_basis = tuple( b.natural_representation() for b in basis )