X-Git-Url: http://gitweb.michael.orlitzky.com/?a=blobdiff_plain;f=mjo%2Feja%2Feja_algebra.py;h=06f6f531ac46305da59ec2b4290ba808602cea50;hb=0fd07263cc543e345f3cd7668938f8a0de70641f;hp=a9bb6c7de4fc829f881978de63d1e78597a51d73;hpb=208bc5d64bd206684e59f757d932768552e9f7ba;p=sage.d.git diff --git a/mjo/eja/eja_algebra.py b/mjo/eja/eja_algebra.py index a9bb6c7..06f6f53 100644 --- a/mjo/eja/eja_algebra.py +++ b/mjo/eja/eja_algebra.py @@ -12,7 +12,7 @@ from sage.matrix.constructor import matrix from sage.misc.cachefunc import cached_method from sage.misc.prandom import choice from sage.misc.table import table -from sage.modules.free_module import VectorSpace +from sage.modules.free_module import FreeModule, VectorSpace from sage.rings.integer_ring import ZZ from sage.rings.number_field.number_field import QuadraticField from sage.rings.polynomial.polynomial_ring_constructor import PolynomialRing @@ -207,8 +207,13 @@ class FiniteDimensionalEuclideanJordanAlgebra(CombinatorialFreeModule): determinant). """ z = self._a_regular_element() - V = self.vector_space() - V1 = V.span_of_basis( (z**k).to_vector() for k in range(self.rank()) ) + # Don't use the parent vector space directly here in case this + # happens to be a subalgebra. In that case, we would be e.g. + # two-dimensional but span_of_basis() would expect three + # coordinates. + V = VectorSpace(self.base_ring(), self.vector_space().dimension()) + basis = [ (z**k).to_vector() for k in range(self.rank()) ] + V1 = V.span_of_basis( basis ) b = (V1.basis() + V1.complement().basis()) return V.span_of_basis(b) @@ -259,16 +264,40 @@ class FiniteDimensionalEuclideanJordanAlgebra(CombinatorialFreeModule): r = self.rank() n = self.dimension() - # Construct a new algebra over a multivariate polynomial ring... + # Turn my vector space into a module so that "vectors" can + # have multivatiate polynomial entries. names = tuple('X' + str(i) for i in range(1,n+1)) R = PolynomialRing(self.base_ring(), names) - # Hack around the fact that our multiplication table is in terms of - # algebra elements but the constructor wants it in terms of vectors. - vmt = [ tuple(map(lambda x: x.to_vector(), ls)) - for ls in self._multiplication_table ] - J = FiniteDimensionalEuclideanJordanAlgebra(R, tuple(vmt), r) - idmat = matrix.identity(J.base_ring(), n) + # Using change_ring() on the parent's vector space doesn't work + # here because, in a subalgebra, that vector space has a basis + # and change_ring() tries to bring the basis along with it. And + # that doesn't work unless the new ring is a PID, which it usually + # won't be. + V = FreeModule(R,n) + + # Now let x = (X1,X2,...,Xn) be the vector whose entries are + # indeterminates... + x = V(names) + + # And figure out the "left multiplication by x" matrix in + # that setting. + lmbx_cols = [] + monomial_matrices = [ self.monomial(i).operator().matrix() + for i in range(n) ] # don't recompute these! + for k in range(n): + ek = self.monomial(k).to_vector() + lmbx_cols.append( + sum( x[i]*(monomial_matrices[i]*ek) + for i in range(n) ) ) + Lx = matrix.column(R, lmbx_cols) + + # Now we can compute powers of x "symbolically" + x_powers = [self.one().to_vector(), x] + for d in range(2, r+1): + x_powers.append( Lx*(x_powers[-1]) ) + + idmat = matrix.identity(R, n) W = self._charpoly_basis_space() W = W.change_ring(R.fraction_field()) @@ -288,18 +317,10 @@ class FiniteDimensionalEuclideanJordanAlgebra(CombinatorialFreeModule): # We want the middle equivalent thing in our matrix, but use # the first equivalent thing instead so that we can pass in # standard coordinates. - x = J.from_vector(W(R.gens())) - - # Handle the zeroth power separately, because computing - # the unit element in J is mathematically suspect. - x0 = W.coordinate_vector(self.one().to_vector()) - l1 = [ x0.column() ] - l1 += [ W.coordinate_vector((x**k).to_vector()).column() - for k in range(1,r) ] - l2 = [idmat.column(k-1).column() for k in range(r+1, n+1)] - A_of_x = matrix.block(R, 1, n, (l1 + l2)) - xr = W.coordinate_vector((x**r).to_vector()) - return (A_of_x, x, xr, A_of_x.det()) + x_powers = [ W.coordinate_vector(xp) for xp in x_powers ] + l2 = [idmat.column(k-1) for k in range(r+1, n+1)] + A_of_x = matrix.column(R, n, (x_powers[:r] + l2)) + return (A_of_x, x, x_powers[r], A_of_x.det()) @cached_method @@ -394,6 +415,29 @@ class FiniteDimensionalEuclideanJordanAlgebra(CombinatorialFreeModule): return x.trace_inner_product(y) + def is_trivial(self): + """ + Return whether or not this algebra is trivial. + + A trivial algebra contains only the zero element. + + SETUP:: + + sage: from mjo.eja.eja_algebra import ComplexHermitianEJA + + EXAMPLES:: + + sage: J = ComplexHermitianEJA(3) + sage: J.is_trivial() + False + sage: A = J.zero().subalgebra_generated_by() + sage: A.is_trivial() + True + + """ + return self.dimension() == 0 + + def multiplication_table(self): """ Return a visual representation of this algebra's multiplication @@ -536,6 +580,15 @@ class FiniteDimensionalEuclideanJordanAlgebra(CombinatorialFreeModule): return self.linear_combination(zip(self.gens(), coeffs)) + def random_element(self): + # Temporary workaround for https://trac.sagemath.org/ticket/28327 + if self.is_trivial(): + return self.zero() + else: + s = super(FiniteDimensionalEuclideanJordanAlgebra, self) + return s.random_element() + + def rank(self): """ Return the rank of this EJA. @@ -649,14 +702,21 @@ class RealCartesianProductEJA(FiniteDimensionalEuclideanJordanAlgebra): sage: e2*e2 e2 + TESTS: + + We can change the generator prefix:: + + sage: RealCartesianProductEJA(3, prefix='r').gens() + (r0, r1, r2) + """ - def __init__(self, n, field=QQ): + def __init__(self, n, field=QQ, **kwargs): V = VectorSpace(field, n) mult_table = [ [ V.gen(i)*(i == j) for j in range(n) ] for i in range(n) ] fdeja = super(RealCartesianProductEJA, self) - return fdeja.__init__(field, mult_table, rank=n) + return fdeja.__init__(field, mult_table, rank=n, **kwargs) def inner_product(self, x, y): return _usual_ip(x,y) @@ -1135,8 +1195,13 @@ class RealSymmetricEJA(FiniteDimensionalEuclideanJordanAlgebra): sage: J(expected) == x*y True + We can change the generator prefix:: + + sage: RealSymmetricEJA(3, prefix='q').gens() + (q0, q1, q2, q3, q4, q5) + """ - def __init__(self, n, field=QQ): + def __init__(self, n, field=QQ, **kwargs): S = _real_symmetric_basis(n, field=field) Qs = _multiplication_table_from_matrix_basis(S) @@ -1144,7 +1209,8 @@ class RealSymmetricEJA(FiniteDimensionalEuclideanJordanAlgebra): return fdeja.__init__(field, Qs, rank=n, - natural_basis=S) + natural_basis=S, + **kwargs) def inner_product(self, x, y): return _matrix_ip(x,y) @@ -1187,8 +1253,13 @@ class ComplexHermitianEJA(FiniteDimensionalEuclideanJordanAlgebra): sage: J(expected) == x*y True + We can change the generator prefix:: + + sage: ComplexHermitianEJA(2, prefix='z').gens() + (z0, z1, z2, z3) + """ - def __init__(self, n, field=QQ): + def __init__(self, n, field=QQ, **kwargs): S = _complex_hermitian_basis(n) Qs = _multiplication_table_from_matrix_basis(S) @@ -1196,7 +1267,8 @@ class ComplexHermitianEJA(FiniteDimensionalEuclideanJordanAlgebra): return fdeja.__init__(field, Qs, rank=n, - natural_basis=S) + natural_basis=S, + **kwargs) def inner_product(self, x, y): @@ -1247,8 +1319,13 @@ class QuaternionHermitianEJA(FiniteDimensionalEuclideanJordanAlgebra): sage: J(expected) == x*y True + We can change the generator prefix:: + + sage: QuaternionHermitianEJA(2, prefix='a').gens() + (a0, a1, a2, a3, a4, a5) + """ - def __init__(self, n, field=QQ): + def __init__(self, n, field=QQ, **kwargs): S = _quaternion_hermitian_basis(n) Qs = _multiplication_table_from_matrix_basis(S) @@ -1256,7 +1333,8 @@ class QuaternionHermitianEJA(FiniteDimensionalEuclideanJordanAlgebra): return fdeja.__init__(field, Qs, rank=n, - natural_basis=S) + natural_basis=S, + **kwargs) def inner_product(self, x, y): # Since a+bi+cj+dk on the diagonal is represented as @@ -1303,8 +1381,13 @@ class JordanSpinEJA(FiniteDimensionalEuclideanJordanAlgebra): sage: e2*e3 0 + We can change the generator prefix:: + + sage: JordanSpinEJA(2, prefix='B').gens() + (B0, B1) + """ - def __init__(self, n, field=QQ): + def __init__(self, n, field=QQ, **kwargs): V = VectorSpace(field, n) mult_table = [[V.zero() for j in range(n)] for i in range(n)] for i in range(n): @@ -1325,7 +1408,7 @@ class JordanSpinEJA(FiniteDimensionalEuclideanJordanAlgebra): # one-dimensional ambient space (because the rank is bounded by # the ambient dimension). fdeja = super(JordanSpinEJA, self) - return fdeja.__init__(field, mult_table, rank=min(n,2)) + return fdeja.__init__(field, mult_table, rank=min(n,2), **kwargs) def inner_product(self, x, y): return _usual_ip(x,y)