"""
-Euclidean Jordan Algebras. These are formally-real Jordan Algebras;
-specifically those where u^2 + v^2 = 0 implies that u = v = 0. They
-are used in optimization, and have some additional nice methods beyond
-what can be supported in a general Jordan Algebra.
-
+Representations and constructions for Euclidean Jordan algebras.
+
+A Euclidean Jordan algebra is a Jordan algebra that has some
+additional properties:
+
+ 1. It is finite-dimensional.
+ 2. Its scalar field is the real numbers.
+ 3a. An inner product is defined on it, and...
+ 3b. That inner product is compatible with the Jordan product
+ in the sense that `<x*y,z> = <y,x*z>` for all elements
+ `x,y,z` in the algebra.
+
+Every Euclidean Jordan algebra is formally-real: for any two elements
+`x` and `y` in the algebra, `x^{2} + y^{2} = 0` implies that `x = y =
+0`. Conversely, every finite-dimensional formally-real Jordan algebra
+can be made into a Euclidean Jordan algebra with an appropriate choice
+of inner-product.
+
+Formally-real Jordan algebras were originally studied as a framework
+for quantum mechanics. Today, Euclidean Jordan algebras are crucial in
+symmetric cone optimization, since every symmetric cone arises as the
+cone of squares in some Euclidean Jordan algebra.
+
+It is known that every Euclidean Jordan algebra decomposes into an
+orthogonal direct sum (essentially, a Cartesian product) of simple
+algebras, and that moreover, up to Jordan-algebra isomorphism, there
+are only five families of simple algebras. We provide constructions
+for these simple algebras:
+
+ * :class:`BilinearFormEJA`
+ * :class:`RealSymmetricEJA`
+ * :class:`ComplexHermitianEJA`
+ * :class:`QuaternionHermitianEJA`
+
+Missing from this list is the algebra of three-by-three octononion
+Hermitian matrices, as there is (as of yet) no implementation of the
+octonions in SageMath. In addition to these, we provide two other
+example constructions,
+
+ * :class:`HadamardEJA`
+ * :class:`TrivialEJA`
+
+The Jordan spin algebra is a bilinear form algebra where the bilinear
+form is the identity. The Hadamard EJA is simply a Cartesian product
+of one-dimensional spin algebras. And last but not least, the trivial
+EJA is exactly what you think. Cartesian products of these are also
+supported using the usual ``cartesian_product()`` function; as a
+result, we support (up to isomorphism) all Euclidean Jordan algebras
+that don't involve octonions.
SETUP::
sage: random_eja()
Euclidean Jordan algebra of dimension...
-
"""
from itertools import repeat
from sage.rings.all import (ZZ, QQ, AA, QQbar, RR, RLF, CLF,
PolynomialRing,
QuadraticField)
-from mjo.eja.eja_element import (CartesianProductEJAElement,
- FiniteDimensionalEJAElement)
+from mjo.eja.eja_element import FiniteDimensionalEJAElement
from mjo.eja.eja_operator import FiniteDimensionalEJAOperator
-from mjo.eja.eja_utils import _mat2vec
+from mjo.eja.eja_utils import _all2list, _mat2vec
class FiniteDimensionalEJA(CombinatorialFreeModule):
r"""
INPUT:
- - basis -- a tuple of basis elements in "matrix form," which
- must be the same form as the arguments to ``jordan_product``
- and ``inner_product``. In reality, "matrix form" can be either
- vectors, matrices, or a Cartesian product (ordered tuple)
- of vectors or matrices. All of these would ideally be vector
- spaces in sage with no special-casing needed; but in reality
- we turn vectors into column-matrices and Cartesian products
- `(a,b)` into column matrices `(a,b)^{T}` after converting
- `a` and `b` themselves.
-
- - jordan_product -- function of two elements (in matrix form)
- that returns their jordan product in this algebra; this will
- be applied to ``basis`` to compute a multiplication table for
- the algebra.
-
- - inner_product -- function of two elements (in matrix form) that
- returns their inner product. This will be applied to ``basis`` to
- compute an inner-product table (basically a matrix) for this algebra.
+ - ``basis`` -- a tuple; a tuple of basis elements in "matrix
+ form," which must be the same form as the arguments to
+ ``jordan_product`` and ``inner_product``. In reality, "matrix
+ form" can be either vectors, matrices, or a Cartesian product
+ (ordered tuple) of vectors or matrices. All of these would
+ ideally be vector spaces in sage with no special-casing
+ needed; but in reality we turn vectors into column-matrices
+ and Cartesian products `(a,b)` into column matrices
+ `(a,b)^{T}` after converting `a` and `b` themselves.
+
+ - ``jordan_product`` -- a function; afunction of two ``basis``
+ elements (in matrix form) that returns their jordan product,
+ also in matrix form; this will be applied to ``basis`` to
+ compute a multiplication table for the algebra.
+
+ - ``inner_product`` -- a function; a function of two ``basis``
+ elements (in matrix form) that returns their inner
+ product. This will be applied to ``basis`` to compute an
+ inner-product table (basically a matrix) for this algebra.
+
+ - ``field`` -- a subfield of the reals (default: ``AA``); the scalar
+ field for the algebra.
+
+ - ``orthonormalize`` -- boolean (default: ``True``); whether or
+ not to orthonormalize the basis. Doing so is expensive and
+ generally rules out using the rationals as your ``field``, but
+ is required for spectral decompositions.
+
+ SETUP::
+
+ sage: from mjo.eja.eja_algebra import random_eja
+
+ TESTS:
+
+ We should compute that an element subalgebra is associative even
+ if we circumvent the element method::
+
+ sage: set_random_seed()
+ sage: J = random_eja(field=QQ,orthonormalize=False)
+ sage: x = J.random_element()
+ sage: A = x.subalgebra_generated_by(orthonormalize=False)
+ sage: basis = tuple(b.superalgebra_element() for b in A.basis())
+ sage: J.subalgebra(basis, orthonormalize=False).is_associative()
+ True
"""
Element = FiniteDimensionalEJAElement
inner_product,
field=AA,
orthonormalize=True,
- associative=False,
+ associative=None,
cartesian_product=False,
check_field=True,
check_axioms=True,
prefix='e'):
+ # Keep track of whether or not the matrix basis consists of
+ # tuples, since we need special cases for them damned near
+ # everywhere. This is INDEPENDENT of whether or not the
+ # algebra is a cartesian product, since a subalgebra of a
+ # cartesian product will have a basis of tuples, but will not
+ # in general itself be a cartesian product algebra.
+ self._matrix_basis_is_cartesian = False
+ n = len(basis)
+ if n > 0:
+ if hasattr(basis[0], 'cartesian_factors'):
+ self._matrix_basis_is_cartesian = True
+
if check_field:
if not field.is_subring(RR):
# Note: this does return true for the real algebraic
# If the basis given to us wasn't over the field that it's
# supposed to be over, fix that. Or, you know, crash.
- basis = tuple( b.change_ring(field) for b in basis )
+ if not cartesian_product:
+ # The field for a cartesian product algebra comes from one
+ # of its factors and is the same for all factors, so
+ # there's no need to "reapply" it on product algebras.
+ if self._matrix_basis_is_cartesian:
+ # OK since if n == 0, the basis does not consist of tuples.
+ P = basis[0].parent()
+ basis = tuple( P(tuple(b_i.change_ring(field) for b_i in b))
+ for b in basis )
+ else:
+ basis = tuple( b.change_ring(field) for b in basis )
+
if check_axioms:
# Check commutativity of the Jordan and inner-products.
category = MagmaticAlgebras(field).FiniteDimensional()
- category = category.WithBasis().Unital()
+ category = category.WithBasis().Unital().Commutative()
+
+ if associative is None:
+ # We should figure it out. As with check_axioms, we have to do
+ # this without the help of the _jordan_product_is_associative()
+ # method because we need to know the category before we
+ # initialize the algebra.
+ associative = all( jordan_product(jordan_product(bi,bj),bk)
+ ==
+ jordan_product(bi,jordan_product(bj,bk))
+ for bi in basis
+ for bj in basis
+ for bk in basis)
+
if associative:
# Element subalgebras can take advantage of this.
category = category.Associative()
# Call the superclass constructor so that we can use its from_vector()
# method to build our multiplication table.
- n = len(basis)
- super().__init__(field,
- range(n),
- prefix=prefix,
- category=category,
- bracket=False)
+ CombinatorialFreeModule.__init__(self,
+ field,
+ range(n),
+ prefix=prefix,
+ category=category,
+ bracket=False)
# Now comes all of the hard work. We'll be constructing an
# ambient vector space V that our (vectorized) basis lives in,
degree = 0
if n > 0:
- # Works on both column and square matrices...
- degree = len(basis[0].list())
+ degree = len(_all2list(basis[0]))
# Build an ambient space that fits our matrix basis when
# written out as "long vectors."
# Save a copy of the un-orthonormalized basis for later.
# Convert it to ambient V (vector) coordinates while we're
# at it, because we'd have to do it later anyway.
- deortho_vector_basis = tuple( V(b.list()) for b in basis )
+ deortho_vector_basis = tuple( V(_all2list(b)) for b in basis )
from mjo.eja.eja_utils import gram_schmidt
basis = tuple(gram_schmidt(basis, inner_product))
# Now create the vector space for the algebra, which will have
# its own set of non-ambient coordinates (in terms of the
# supplied basis).
- vector_basis = tuple( V(b.list()) for b in basis )
+ vector_basis = tuple( V(_all2list(b)) for b in basis )
W = V.span_of_basis( vector_basis, check=check_axioms)
if orthonormalize:
# The jordan product returns a matrixy answer, so we
# have to convert it to the algebra coordinates.
elt = jordan_product(q_i, q_j)
- elt = W.coordinate_vector(V(elt.list()))
+ elt = W.coordinate_vector(V(_all2list(elt)))
self._multiplication_table[i][j] = self.from_vector(elt)
if not orthonormalize:
def product_on_basis(self, i, j):
+ r"""
+ Returns the Jordan product of the `i` and `j`th basis elements.
+
+ This completely defines the Jordan product on the algebra, and
+ is used direclty by our superclass machinery to implement
+ :meth:`product`.
+
+ SETUP::
+
+ sage: from mjo.eja.eja_algebra import random_eja
+
+ TESTS::
+
+ sage: set_random_seed()
+ sage: J = random_eja()
+ sage: n = J.dimension()
+ sage: ei = J.zero()
+ sage: ej = J.zero()
+ sage: ei_ej = J.zero()*J.zero()
+ sage: if n > 0:
+ ....: i = ZZ.random_element(n)
+ ....: j = ZZ.random_element(n)
+ ....: ei = J.gens()[i]
+ ....: ej = J.gens()[j]
+ ....: ei_ej = J.product_on_basis(i,j)
+ sage: ei*ej == ei_ej
+ True
+
+ """
# We only stored the lower-triangular portion of the
# multiplication table.
if j <= i:
sage: y = J.random_element()
sage: (n == 1) or (x.inner_product(y) == (x*y).trace()/2)
True
+
"""
B = self._inner_product_matrix
return (B*x.to_vector()).inner_product(y.to_vector())
+ def is_associative(self):
+ r"""
+ Return whether or not this algebra's Jordan product is associative.
+
+ SETUP::
+
+ sage: from mjo.eja.eja_algebra import ComplexHermitianEJA
+
+ EXAMPLES::
+
+ sage: J = ComplexHermitianEJA(3, field=QQ, orthonormalize=False)
+ sage: J.is_associative()
+ False
+ sage: x = sum(J.gens())
+ sage: A = x.subalgebra_generated_by(orthonormalize=False)
+ sage: A.is_associative()
+ True
+
+ """
+ return "Associative" in self.category().axioms()
+
def _is_commutative(self):
r"""
Whether or not this algebra's multiplication table is commutative.
this algebra was constructed with ``check_axioms=False`` and
passed an invalid multiplication table.
"""
- return all( self.product_on_basis(i,j) == self.product_on_basis(i,j)
- for i in range(self.dimension())
- for j in range(self.dimension()) )
+ return all( x*y == y*x for x in self.gens() for y in self.gens() )
def _is_jordanian(self):
r"""
for i in range(self.dimension())
for j in range(self.dimension()) )
+ def _jordan_product_is_associative(self):
+ r"""
+ Return whether or not this algebra's Jordan product is
+ associative; that is, whether or not `x*(y*z) = (x*y)*z`
+ for all `x,y,x`.
+
+ This method should agree with :meth:`is_associative` unless
+ you lied about the value of the ``associative`` parameter
+ when you constructed the algebra.
+
+ SETUP::
+
+ sage: from mjo.eja.eja_algebra import (random_eja,
+ ....: RealSymmetricEJA,
+ ....: ComplexHermitianEJA,
+ ....: QuaternionHermitianEJA)
+
+ EXAMPLES::
+
+ sage: J = RealSymmetricEJA(4, orthonormalize=False)
+ sage: J._jordan_product_is_associative()
+ False
+ sage: x = sum(J.gens())
+ sage: A = x.subalgebra_generated_by()
+ sage: A._jordan_product_is_associative()
+ True
+
+ ::
+
+ sage: J = ComplexHermitianEJA(2,field=QQ,orthonormalize=False)
+ sage: J._jordan_product_is_associative()
+ False
+ sage: x = sum(J.gens())
+ sage: A = x.subalgebra_generated_by(orthonormalize=False)
+ sage: A._jordan_product_is_associative()
+ True
+
+ ::
+
+ sage: J = QuaternionHermitianEJA(2)
+ sage: J._jordan_product_is_associative()
+ False
+ sage: x = sum(J.gens())
+ sage: A = x.subalgebra_generated_by()
+ sage: A._jordan_product_is_associative()
+ True
+
+ TESTS:
+
+ The values we've presupplied to the constructors agree with
+ the computation::
+
+ sage: set_random_seed()
+ sage: J = random_eja()
+ sage: J.is_associative() == J._jordan_product_is_associative()
+ True
+
+ """
+ R = self.base_ring()
+
+ # Used to check whether or not something is zero.
+ epsilon = R.zero()
+ if not R.is_exact():
+ # I don't know of any examples that make this magnitude
+ # necessary because I don't know how to make an
+ # associative algebra when the element subalgebra
+ # construction is unreliable (as it is over RDF; we can't
+ # find the degree of an element because we can't compute
+ # the rank of a matrix). But even multiplication of floats
+ # is non-associative, so *some* epsilon is needed... let's
+ # just take the one from _inner_product_is_associative?
+ epsilon = 1e-15
+
+ for i in range(self.dimension()):
+ for j in range(self.dimension()):
+ for k in range(self.dimension()):
+ x = self.gens()[i]
+ y = self.gens()[j]
+ z = self.gens()[k]
+ diff = (x*y)*z - x*(y*z)
+
+ if diff.norm() > epsilon:
+ return False
+
+ return True
+
def _inner_product_is_associative(self):
r"""
Return whether or not this algebra's inner product `B` is
this algebra was constructed with ``check_axioms=False`` and
passed an invalid Jordan or inner-product.
"""
+ R = self.base_ring()
- # Used to check whether or not something is zero in an inexact
- # ring. This number is sufficient to allow the construction of
- # QuaternionHermitianEJA(2, field=RDF) with check_axioms=True.
- epsilon = 1e-16
+ # Used to check whether or not something is zero.
+ epsilon = R.zero()
+ if not R.is_exact():
+ # This choice is sufficient to allow the construction of
+ # QuaternionHermitianEJA(2, field=RDF) with check_axioms=True.
+ epsilon = 1e-15
for i in range(self.dimension()):
for j in range(self.dimension()):
z = self.gens()[k]
diff = (x*y).inner_product(z) - x.inner_product(y*z)
- if self.base_ring().is_exact():
- if diff != 0:
- return False
- else:
- if diff.abs() > epsilon:
- return False
+ if diff.abs() > epsilon:
+ return False
return True
SETUP::
- sage: from mjo.eja.eja_algebra import (JordanSpinEJA,
+ sage: from mjo.eja.eja_algebra import (random_eja,
+ ....: JordanSpinEJA,
....: HadamardEJA,
....: RealSymmetricEJA)
...
ValueError: not an element of this algebra
+ Tuples work as well, provided that the matrix basis for the
+ algebra consists of them::
+
+ sage: J1 = HadamardEJA(3)
+ sage: J2 = RealSymmetricEJA(2)
+ sage: J = cartesian_product([J1,J2])
+ sage: J( (J1.matrix_basis()[1], J2.matrix_basis()[2]) )
+ e(0, 1) + e(1, 2)
+
TESTS:
- Ensure that we can convert any element of the two non-matrix
- simple algebras (whose matrix representations are columns)
- back and forth faithfully::
+ Ensure that we can convert any element back and forth
+ faithfully between its matrix and algebra representations::
sage: set_random_seed()
- sage: J = HadamardEJA.random_instance()
- sage: x = J.random_element()
- sage: J(x.to_vector().column()) == x
- True
- sage: J = JordanSpinEJA.random_instance()
+ sage: J = random_eja()
sage: x = J.random_element()
- sage: J(x.to_vector().column()) == x
+ sage: J(x.to_matrix()) == x
True
+ We cannot coerce elements between algebras just because their
+ matrix representations are compatible::
+
+ sage: J1 = HadamardEJA(3)
+ sage: J2 = JordanSpinEJA(3)
+ sage: J2(J1.one())
+ Traceback (most recent call last):
+ ...
+ ValueError: not an element of this algebra
+ sage: J1(J2.zero())
+ Traceback (most recent call last):
+ ...
+ ValueError: not an element of this algebra
"""
msg = "not an element of this algebra"
- if elt == 0:
- # The superclass implementation of random_element()
- # needs to be able to coerce "0" into the algebra.
- return self.zero()
- elif elt in self.base_ring():
+ if elt in self.base_ring():
# Ensure that no base ring -> algebra coercion is performed
# by this method. There's some stupidity in sage that would
# otherwise propagate to this method; for example, sage thinks
raise ValueError(msg)
try:
+ # Try to convert a vector into a column-matrix...
elt = elt.column()
except (AttributeError, TypeError):
- # Try to convert a vector into a column-matrix
+ # and ignore failure, because we weren't really expecting
+ # a vector as an argument anyway.
pass
if elt not in self.matrix_space():
# closure whereas the base ring of the 3-by-3 identity matrix
# could be QQ instead of QQbar.
#
+ # And, we also have to handle Cartesian product bases (when
+ # the matrix basis consists of tuples) here. The "good news"
+ # is that we're already converting everything to long vectors,
+ # and that strategy works for tuples as well.
+ #
# We pass check=False because the matrix basis is "guaranteed"
# to be linearly independent... right? Ha ha.
- V = VectorSpace(self.base_ring(), elt.nrows()*elt.ncols())
- W = V.span_of_basis( (_mat2vec(s) for s in self.matrix_basis()),
+ elt = _all2list(elt)
+ V = VectorSpace(self.base_ring(), len(elt))
+ W = V.span_of_basis( (V(_all2list(s)) for s in self.matrix_basis()),
check=False)
try:
- coords = W.coordinate_vector(_mat2vec(elt))
+ coords = W.coordinate_vector(V(elt))
except ArithmeticError: # vector is not in free module
raise ValueError(msg)
# And to each subsequent row, prepend an entry that belongs to
# the left-side "header column."
- M += [ [self.gens()[i]] + [ self.product_on_basis(i,j)
+ M += [ [self.gens()[i]] + [ self.gens()[i]*self.gens()[j]
for j in range(n) ]
for i in range(n) ]
we think of them as matrices (including column vectors of the
appropriate size).
- Generally this will be an `n`-by-`1` column-vector space,
+ "By default" this will be an `n`-by-`1` column-matrix space,
except when the algebra is trivial. There it's `n`-by-`n`
(where `n` is zero), to ensure that two elements of the matrix
- space (empty matrices) can be multiplied.
+ space (empty matrices) can be multiplied. For algebras of
+ matrices, this returns the space in which their
+ real embeddings live.
+
+ SETUP::
+
+ sage: from mjo.eja.eja_algebra import (ComplexHermitianEJA,
+ ....: JordanSpinEJA,
+ ....: QuaternionHermitianEJA,
+ ....: TrivialEJA)
+
+ EXAMPLES:
+
+ By default, the matrix representation is just a column-matrix
+ equivalent to the vector representation::
+
+ sage: J = JordanSpinEJA(3)
+ sage: J.matrix_space()
+ Full MatrixSpace of 3 by 1 dense matrices over Algebraic
+ Real Field
+
+ The matrix representation in the trivial algebra is
+ zero-by-zero instead of the usual `n`-by-one::
+
+ sage: J = TrivialEJA()
+ sage: J.matrix_space()
+ Full MatrixSpace of 0 by 0 dense matrices over Algebraic
+ Real Field
+
+ The matrix space for complex/quaternion Hermitian matrix EJA
+ is the space in which their real-embeddings live, not the
+ original complex/quaternion matrix space::
+
+ sage: J = ComplexHermitianEJA(2,field=QQ,orthonormalize=False)
+ sage: J.matrix_space()
+ Full MatrixSpace of 4 by 4 dense matrices over Rational Field
+ sage: J = QuaternionHermitianEJA(1,field=QQ,orthonormalize=False)
+ sage: J.matrix_space()
+ Full MatrixSpace of 4 by 4 dense matrices over Rational Field
- Matrix algebras override this with something more useful.
"""
if self.is_trivial():
return MatrixSpace(self.base_ring(), 0)
if not c.is_idempotent():
raise ValueError("element is not idempotent: %s" % c)
- from mjo.eja.eja_subalgebra import FiniteDimensionalEJASubalgebra
-
# Default these to what they should be if they turn out to be
# trivial, because eigenspaces_left() won't return eigenvalues
# corresponding to trivial spaces (e.g. it returns only the
# eigenspace corresponding to lambda=1 if you take the
# decomposition relative to the identity element).
- trivial = FiniteDimensionalEJASubalgebra(self, ())
+ trivial = self.subalgebra(())
J0 = trivial # eigenvalue zero
J5 = VectorSpace(self.base_ring(), 0) # eigenvalue one-half
J1 = trivial # eigenvalue one
J5 = eigspace
else:
gens = tuple( self.from_vector(b) for b in eigspace.basis() )
- subalg = FiniteDimensionalEJASubalgebra(self,
- gens,
- check_axioms=False)
+ subalg = self.subalgebra(gens, check_axioms=False)
if eigval == 0:
J0 = subalg
elif eigval == 1:
return len(self._charpoly_coefficients())
+ def subalgebra(self, basis, **kwargs):
+ r"""
+ Create a subalgebra of this algebra from the given basis.
+ """
+ from mjo.eja.eja_subalgebra import FiniteDimensionalEJASubalgebra
+ return FiniteDimensionalEJASubalgebra(self, basis, **kwargs)
+
+
def vector_space(self):
"""
Return the vector space that underlies this algebra.
return self.zero().to_vector().parent().ambient_vector_space()
- Element = FiniteDimensionalEJAElement
class RationalBasisEJA(FiniteDimensionalEJA):
r"""
if not all( all(b_i in QQ for b_i in b.list()) for b in basis ):
raise TypeError("basis not rational")
+ super().__init__(basis,
+ jordan_product,
+ inner_product,
+ field=field,
+ check_field=check_field,
+ **kwargs)
+
self._rational_algebra = None
if field is not QQ:
# There's no point in constructing the extra algebra if this
jordan_product,
inner_product,
field=QQ,
+ associative=self.is_associative(),
orthonormalize=False,
check_field=False,
check_axioms=False)
- super().__init__(basis,
- jordan_product,
- inner_product,
- field=field,
- check_field=check_field,
- **kwargs)
-
@cached_method
def _charpoly_coefficients(self):
r"""
In theory, our "field" can be any subfield of the reals::
- sage: RealSymmetricEJA(2, field=RDF)
+ sage: RealSymmetricEJA(2, field=RDF, check_axioms=True)
Euclidean Jordan algebra of dimension 3 over Real Double Field
- sage: RealSymmetricEJA(2, field=RR)
+ sage: RealSymmetricEJA(2, field=RR, check_axioms=True)
Euclidean Jordan algebra of dimension 3 over Real Field with
53 bits of precision
# if the user passes check_axioms=True.
if "check_axioms" not in kwargs: kwargs["check_axioms"] = False
- super(RealSymmetricEJA, self).__init__(self._denormalized_basis(n),
- self.jordan_product,
- self.trace_inner_product,
- **kwargs)
+ associative = False
+ if n <= 1:
+ associative = True
+
+ super().__init__(self._denormalized_basis(n),
+ self.jordan_product,
+ self.trace_inner_product,
+ associative=associative,
+ **kwargs)
# TODO: this could be factored out somehow, but is left here
# because the MatrixEJA is not presently a subclass of the
True
"""
- super(ComplexMatrixEJA,cls).real_embed(M)
+ super().real_embed(M)
n = M.nrows()
# We don't need any adjoined elements...
True
"""
- super(ComplexMatrixEJA,cls).real_unembed(M)
+ super().real_unembed(M)
n = ZZ(M.nrows())
d = cls.dimension_over_reals()
F = cls.complex_extension(M.base_ring())
In theory, our "field" can be any subfield of the reals::
- sage: ComplexHermitianEJA(2, field=RDF)
+ sage: ComplexHermitianEJA(2, field=RDF, check_axioms=True)
Euclidean Jordan algebra of dimension 4 over Real Double Field
- sage: ComplexHermitianEJA(2, field=RR)
+ sage: ComplexHermitianEJA(2, field=RR, check_axioms=True)
Euclidean Jordan algebra of dimension 4 over Real Field with
53 bits of precision
# if the user passes check_axioms=True.
if "check_axioms" not in kwargs: kwargs["check_axioms"] = False
- super(ComplexHermitianEJA, self).__init__(self._denormalized_basis(n),
- self.jordan_product,
- self.trace_inner_product,
- **kwargs)
+ associative = False
+ if n <= 1:
+ associative = True
+
+ super().__init__(self._denormalized_basis(n),
+ self.jordan_product,
+ self.trace_inner_product,
+ associative=associative,
+ **kwargs)
# TODO: this could be factored out somehow, but is left here
# because the MatrixEJA is not presently a subclass of the
# FDEJA class that defines rank() and one().
True
"""
- super(QuaternionMatrixEJA,cls).real_embed(M)
+ super().real_embed(M)
quaternions = M.base_ring()
n = M.nrows()
True
"""
- super(QuaternionMatrixEJA,cls).real_unembed(M)
+ super().real_unembed(M)
n = ZZ(M.nrows())
d = cls.dimension_over_reals()
In theory, our "field" can be any subfield of the reals::
- sage: QuaternionHermitianEJA(2, field=RDF)
+ sage: QuaternionHermitianEJA(2, field=RDF, check_axioms=True)
Euclidean Jordan algebra of dimension 6 over Real Double Field
- sage: QuaternionHermitianEJA(2, field=RR)
+ sage: QuaternionHermitianEJA(2, field=RR, check_axioms=True)
Euclidean Jordan algebra of dimension 6 over Real Field with
53 bits of precision
# if the user passes check_axioms=True.
if "check_axioms" not in kwargs: kwargs["check_axioms"] = False
- super(QuaternionHermitianEJA, self).__init__(self._denormalized_basis(n),
- self.jordan_product,
- self.trace_inner_product,
- **kwargs)
+ associative = False
+ if n <= 1:
+ associative = True
+
+ super().__init__(self._denormalized_basis(n),
+ self.jordan_product,
+ self.trace_inner_product,
+ associative=associative,
+ **kwargs)
+
# TODO: this could be factored out somehow, but is left here
# because the MatrixEJA is not presently a subclass of the
# FDEJA class that defines rank() and one().
if "check_axioms" not in kwargs: kwargs["check_axioms"] = False
column_basis = tuple( b.column() for b in FreeModule(ZZ, n).basis() )
- super().__init__(column_basis, jordan_product, inner_product, **kwargs)
+ super().__init__(column_basis,
+ jordan_product,
+ inner_product,
+ associative=True,
+ **kwargs)
self.rank.set_cache(n)
if n == 0:
n = B.nrows()
column_basis = tuple( b.column() for b in FreeModule(ZZ, n).basis() )
- super(BilinearFormEJA, self).__init__(column_basis,
- jordan_product,
- inner_product,
- **kwargs)
+
+ # TODO: I haven't actually checked this, but it seems legit.
+ associative = False
+ if n <= 2:
+ associative = True
+
+ super().__init__(column_basis,
+ jordan_product,
+ inner_product,
+ associative=associative,
+ **kwargs)
# The rank of this algebra is two, unless we're in a
# one-dimensional ambient space (because the rank is bounded
# But also don't pass check_field=False here, because the user
# can pass in a field!
- super(JordanSpinEJA, self).__init__(B, **kwargs)
+ super().__init__(B, **kwargs)
@staticmethod
def _max_random_instance_size():
if "orthonormalize" not in kwargs: kwargs["orthonormalize"] = False
if "check_axioms" not in kwargs: kwargs["check_axioms"] = False
- super(TrivialEJA, self).__init__(basis,
- jordan_product,
- inner_product,
- **kwargs)
+ super().__init__(basis,
+ jordan_product,
+ inner_product,
+ associative=True,
+ **kwargs)
+
# The rank is zero using my definition, namely the dimension of the
# largest subalgebra generated by any element.
self.rank.set_cache(0)
sage: J.rank() == J1.rank() + J2.rank()
True
+ The product algebra will be associative if and only if all of its
+ components are associative::
+
+ sage: J1 = HadamardEJA(2)
+ sage: J1.is_associative()
+ True
+ sage: J2 = HadamardEJA(3)
+ sage: J2.is_associative()
+ True
+ sage: J3 = RealSymmetricEJA(3)
+ sage: J3.is_associative()
+ False
+ sage: CP1 = cartesian_product([J1,J2])
+ sage: CP1.is_associative()
+ True
+ sage: CP2 = cartesian_product([J1,J3])
+ sage: CP2.is_associative()
+ False
+
TESTS:
All factors must share the same base field::
...
ValueError: all factors must share the same base field
- The "cached" Jordan and inner products are the componentwise
- ones::
-
- sage: set_random_seed()
- sage: J1 = random_eja()
- sage: J2 = random_eja()
- sage: J = cartesian_product([J1,J2])
- sage: x,y = J.random_elements(2)
- sage: x*y == J.cartesian_jordan_product(x,y)
- True
- sage: x.inner_product(y) == J.cartesian_inner_product(x,y)
- True
-
The cached unit element is the same one that would be computed::
sage: set_random_seed() # long time
True
"""
- def __init__(self, modules, **kwargs):
+ Element = FiniteDimensionalEJAElement
+
+
+ def __init__(self, algebras, **kwargs):
CombinatorialFreeModule_CartesianProduct.__init__(self,
- modules,
+ algebras,
**kwargs)
- field = modules[0].base_ring()
- if not all( J.base_ring() == field for J in modules ):
+ field = algebras[0].base_ring()
+ if not all( J.base_ring() == field for J in algebras ):
raise ValueError("all factors must share the same base field")
- basis = tuple( b.to_vector().column() for b in self.basis() )
+ associative = all( m.is_associative() for m in algebras )
- # Define jordan/inner products that operate on the basis.
- def jordan_product(x_mat,y_mat):
- x = self.from_vector(_mat2vec(x_mat))
- y = self.from_vector(_mat2vec(y_mat))
- return self.cartesian_jordan_product(x,y).to_vector().column()
+ # The definition of matrix_space() and self.basis() relies
+ # only on the stuff in the CFM_CartesianProduct class, which
+ # we've already initialized.
+ Js = self.cartesian_factors()
+ m = len(Js)
+ MS = self.matrix_space()
+ basis = tuple(
+ MS(tuple( self.cartesian_projection(i)(b).to_matrix()
+ for i in range(m) ))
+ for b in self.basis()
+ )
- def inner_product(x_mat, y_mat):
- x = self.from_vector(_mat2vec(x_mat))
- y = self.from_vector(_mat2vec(y_mat))
- return self.cartesian_inner_product(x,y)
+ # Define jordan/inner products that operate on that matrix_basis.
+ def jordan_product(x,y):
+ return MS(tuple(
+ (Js[i](x[i])*Js[i](y[i])).to_matrix() for i in range(m)
+ ))
+
+ def inner_product(x, y):
+ return sum(
+ Js[i](x[i]).inner_product(Js[i](y[i])) for i in range(m)
+ )
- # Use whatever category the superclass came up with. Usually
- # some join of the EJA and Cartesian product
- # categories. There's no need to check the field since it
- # already came from an EJA.
+ # There's no need to check the field since it already came
+ # from an EJA. Likewise the axioms are guaranteed to be
+ # satisfied, unless the guy writing this class sucks.
#
# If you want the basis to be orthonormalized, orthonormalize
# the factors.
inner_product,
field=field,
orthonormalize=False,
+ associative=associative,
cartesian_product=True,
check_field=False,
check_axioms=False)
- ones = tuple(J.one() for J in modules)
+ ones = tuple(J.one() for J in algebras)
self.one.set_cache(self._cartesian_product_of_elements(ones))
- self.rank.set_cache(sum(J.rank() for J in modules))
-
- # Now that everything else is ready, we clobber our computed
- # matrix basis with the "correct" one consisting of ordered
- # tuples. Since we didn't orthonormalize our basis, we can
- # create these from the basis that was handed to us; that is,
- # we don't need to use the one that the earlier __init__()
- # method came up with.
- m = len(self.cartesian_factors())
- MS = self.matrix_space()
- self._matrix_basis = tuple(
- MS(tuple( self.cartesian_projection(i)(b).to_matrix()
- for i in range(m) ))
- for b in self.basis()
- )
+ self.rank.set_cache(sum(J.rank() for J in algebras))
+
+ def _monomial_to_generator(self, mon):
+ r"""
+ Convert a monomial index into a generator index.
+
+ SETUP::
+
+ sage: from mjo.eja.eja_algebra import random_eja()
+
+ TESTS::
+
+ sage: J1 = random_eja(field=QQ, orthonormalize=False)
+ sage: J2 = random_eja(field=QQ, orthonormalize=False)
+ sage: J = cartesian_product([J1,J2])
+ sage: all( J.monomial(m)
+ ....: ==
+ ....: J.gens()[J._monomial_to_generator(m)]
+ ....: for m in J.basis().keys() )
+
+ """
+ # The superclass method indexes into a matrix, so we have to
+ # turn the tuples i and j into integers. This is easy enough
+ # given that the first coordinate of i and j corresponds to
+ # the factor, and the second coordinate corresponds to the
+ # index of the generator within that factor.
+ try:
+ factor = mon[0]
+ except TypeError: # 'int' object is not subscriptable
+ return mon
+ idx_in_factor = self._monomial_to_generator(mon[1])
+
+ offset = sum( f.dimension()
+ for f in self.cartesian_factors()[:factor] )
+ return offset + idx_in_factor
+
+ def product_on_basis(self, i, j):
+ r"""
+ Return the product of the monomials indexed by ``i`` and ``j``.
+
+ This overrides the superclass method because here, both ``i``
+ and ``j`` will be ordered pairs.
+
+ SETUP::
+
+ sage: from mjo.eja.eja_algebra import (HadamardEJA,
+ ....: JordanSpinEJA,
+ ....: QuaternionHermitianEJA,
+ ....: RealSymmetricEJA,)
+
+ EXAMPLES::
+
+ sage: J1 = JordanSpinEJA(2, field=QQ)
+ sage: J2 = RealSymmetricEJA(2, field=QQ, orthonormalize=False)
+ sage: J3 = HadamardEJA(1, field=QQ)
+ sage: K1 = cartesian_product([J1,J2])
+ sage: K2 = cartesian_product([K1,J3])
+ sage: list(K2.basis())
+ [e(0, (0, 0)), e(0, (0, 1)), e(0, (1, 0)), e(0, (1, 1)),
+ e(0, (1, 2)), e(1, 0)]
+ sage: sage: g = K2.gens()
+ sage: (g[0] + 2*g[3]) * (g[1] - 4*g[2])
+ e(0, (0, 1)) - 4*e(0, (1, 1))
+
+ TESTS::
+
+ sage: J1 = RealSymmetricEJA(1,field=QQ)
+ sage: J2 = QuaternionHermitianEJA(1,field=QQ)
+ sage: J = cartesian_product([J1,J2])
+ sage: x = sum(J.gens())
+ sage: x == J.one()
+ True
+ sage: x*x == x
+ True
+
+ """
+ l = self._monomial_to_generator(i)
+ m = self._monomial_to_generator(j)
+ return FiniteDimensionalEJA.product_on_basis(self, l, m)
def matrix_space(self):
r"""
return FiniteDimensionalEJAOperator(Ji,self,Ei.matrix())
- def cartesian_jordan_product(self, x, y):
- r"""
- The componentwise Jordan product.
-
- We project ``x`` and ``y`` onto our factors, and add up the
- Jordan products from the subalgebras. This may still be useful
- after (if) the default Jordan product in the Cartesian product
- algebra is overridden.
-
- SETUP::
-
- sage: from mjo.eja.eja_algebra import (HadamardEJA,
- ....: JordanSpinEJA)
-
- EXAMPLE::
-
- sage: J1 = HadamardEJA(3)
- sage: J2 = JordanSpinEJA(3)
- sage: J = cartesian_product([J1,J2])
- sage: x1 = J1.from_vector(vector(QQ,(1,2,1)))
- sage: y1 = J1.from_vector(vector(QQ,(1,0,2)))
- sage: x2 = J2.from_vector(vector(QQ,(1,2,3)))
- sage: y2 = J2.from_vector(vector(QQ,(1,1,1)))
- sage: z1 = J.from_vector(vector(QQ,(1,2,1,1,2,3)))
- sage: z2 = J.from_vector(vector(QQ,(1,0,2,1,1,1)))
- sage: (x1*y1).to_vector()
- (1, 0, 2)
- sage: (x2*y2).to_vector()
- (6, 3, 4)
- sage: J.cartesian_jordan_product(z1,z2).to_vector()
- (1, 0, 2, 6, 3, 4)
-
- """
- m = len(self.cartesian_factors())
- projections = ( self.cartesian_projection(i) for i in range(m) )
- products = ( P(x)*P(y) for P in projections )
- return self._cartesian_product_of_elements(tuple(products))
-
- def cartesian_inner_product(self, x, y):
- r"""
- The standard componentwise Cartesian inner-product.
-
- We project ``x`` and ``y`` onto our factors, and add up the
- inner-products from the subalgebras. This may still be useful
- after (if) the default inner product in the Cartesian product
- algebra is overridden.
-
- SETUP::
-
- sage: from mjo.eja.eja_algebra import (HadamardEJA,
- ....: QuaternionHermitianEJA)
- EXAMPLE::
+FiniteDimensionalEJA.CartesianProduct = CartesianProductEJA
- sage: J1 = HadamardEJA(3,field=QQ)
- sage: J2 = QuaternionHermitianEJA(2,field=QQ,orthonormalize=False)
- sage: J = cartesian_product([J1,J2])
- sage: x1 = J1.one()
- sage: x2 = x1
- sage: y1 = J2.one()
- sage: y2 = y1
- sage: x1.inner_product(x2)
- 3
- sage: y1.inner_product(y2)
- 2
- sage: z1 = J._cartesian_product_of_elements((x1,y1))
- sage: z2 = J._cartesian_product_of_elements((x2,y2))
- sage: J.cartesian_inner_product(z1,z2)
- 5
+class RationalBasisCartesianProductEJA(CartesianProductEJA,
+ RationalBasisEJA):
+ r"""
+ A separate class for products of algebras for which we know a
+ rational basis.
- """
- m = len(self.cartesian_factors())
- projections = ( self.cartesian_projection(i) for i in range(m) )
- return sum( P(x).inner_product(P(y)) for P in projections )
+ SETUP::
+ sage: from mjo.eja.eja_algebra import (JordanSpinEJA,
+ ....: RealSymmetricEJA)
- def _element_constructor_(self, elt):
- r"""
- Construct an element of this algebra from an ordered tuple.
+ EXAMPLES:
- We just apply the element constructor from each of our factors
- to the corresponding component of the tuple, and package up
- the result.
+ This gives us fast characteristic polynomial computations in
+ product algebras, too::
- SETUP::
- sage: from mjo.eja.eja_algebra import (HadamardEJA,
- ....: RealSymmetricEJA)
+ sage: J1 = JordanSpinEJA(2)
+ sage: J2 = RealSymmetricEJA(3)
+ sage: J = cartesian_product([J1,J2])
+ sage: J.characteristic_polynomial_of().degree()
+ 5
+ sage: J.rank()
+ 5
- EXAMPLES::
+ """
+ def __init__(self, algebras, **kwargs):
+ CartesianProductEJA.__init__(self, algebras, **kwargs)
- sage: J1 = HadamardEJA(3)
- sage: J2 = RealSymmetricEJA(2)
- sage: J = cartesian_product([J1,J2])
- sage: J( (J1.matrix_basis()[1], J2.matrix_basis()[2]) )
- e(0, 1) + e(1, 2)
- """
- m = len(self.cartesian_factors())
- try:
- z = tuple( self.cartesian_factors()[i](elt[i]) for i in range(m) )
- return self._cartesian_product_of_elements(z)
- except:
- raise ValueError("not an element of this algebra")
+ self._rational_algebra = None
+ if self.vector_space().base_field() is not QQ:
+ self._rational_algebra = cartesian_product([
+ r._rational_algebra for r in algebras
+ ])
- Element = CartesianProductEJAElement
+RationalBasisEJA.CartesianProduct = RationalBasisCartesianProductEJA
-FiniteDimensionalEJA.CartesianProduct = CartesianProductEJA
random_eja = ConcreteEJA.random_instance
+
+# def random_eja(*args, **kwargs):
+# J1 = ConcreteEJA.random_instance(*args, **kwargs)
+
+# # This might make Cartesian products appear roughly as often as
+# # any other ConcreteEJA.
+# if ZZ.random_element(len(ConcreteEJA.__subclasses__()) + 1) == 0:
+# # Use random_eja() again so we can get more than two factors.
+# J2 = random_eja(*args, **kwargs)
+# J = cartesian_product([J1,J2])
+# return J
+# else:
+# return J1
projects / sage.d.git / blobdiff