from sage.rings.all import (ZZ, QQ, AA, QQbar, RR, RLF, CLF,
PolynomialRing,
QuadraticField)
-from mjo.eja.eja_element import FiniteDimensionalEJAElement
+from mjo.eja.eja_element import (CartesianProductEJAElement,
+ FiniteDimensionalEJAElement)
from mjo.eja.eja_operator import FiniteDimensionalEJAOperator
-from mjo.eja.eja_utils import _all2list, _mat2vec
+from mjo.eja.eja_utils import _all2list
def EuclideanJordanAlgebras(field):
r"""
#
# Of course, matrices aren't vectors in sage, so we have to
# appeal to the "long vectors" isometry.
- oper_vecs = [ _mat2vec(g.operator().matrix()) for g in self.gens() ]
+
+ V = VectorSpace(self.base_ring(), self.dimension()**2)
+ oper_vecs = [ V(g.operator().matrix().list()) for g in self.gens() ]
# Now we use basic linear algebra to find the coefficients,
# of the matrices-as-vectors-linear-combination, which should
# We used the isometry on the left-hand side already, but we
# still need to do it for the right-hand side. Recall that we
# wanted something that summed to the identity matrix.
- b = _mat2vec( matrix.identity(self.base_ring(), self.dimension()) )
+ b = V( matrix.identity(self.base_ring(), self.dimension()).list() )
# Now if there's an identity element in the algebra, this
# should work. We solve on the left to avoid having to
check_field=False,
check_axioms=False)
+ def rational_algebra(self):
+ # Using None as a flag here (rather than just assigning "self"
+ # to self._rational_algebra by default) feels a little bit
+ # more sane to me in a garbage-collected environment.
+ if self._rational_algebra is None:
+ return self
+ else:
+ return self._rational_algebra
+
@cached_method
def _charpoly_coefficients(self):
r"""
Algebraic Real Field
"""
- if self._rational_algebra is None:
- # There's no need to construct *another* algebra over the
- # rationals if this one is already over the
- # rationals. Likewise, if we never orthonormalized our
- # basis, we might as well just use the given one.
+ if self.rational_algebra() is self:
+ # Bypass the hijinks if they won't benefit us.
return super()._charpoly_coefficients()
# Do the computation over the rationals. The answer will be
# the same, because all we've done is a change of basis.
# Then, change back from QQ to our real base ring
a = ( a_i.change_ring(self.base_ring())
- for a_i in self._rational_algebra._charpoly_coefficients() )
+ for a_i in self.rational_algebra()._charpoly_coefficients() )
# Otherwise, convert the coordinate variables back to the
# deorthonormalized ones.
# if the user passes check_axioms=True.
if "check_axioms" not in kwargs: kwargs["check_axioms"] = False
-
super().__init__(self._denormalized_basis(matrix_space),
self.jordan_product,
self.trace_inner_product,
return cls(n, **kwargs)
def __init__(self, n, field=AA, **kwargs):
- # We know this is a valid EJA, but will double-check
- # if the user passes check_axioms=True.
- if "check_axioms" not in kwargs: kwargs["check_axioms"] = False
-
A = MatrixSpace(field, n)
super().__init__(A, **kwargs)
from mjo.eja.eja_cache import real_symmetric_eja_coeffs
a = real_symmetric_eja_coeffs(self)
if a is not None:
- if self._rational_algebra is None:
- self._charpoly_coefficients.set_cache(a)
- else:
- self._rational_algebra._charpoly_coefficients.set_cache(a)
+ self.rational_algebra()._charpoly_coefficients.set_cache(a)
"""
def __init__(self, n, field=AA, **kwargs):
- # We know this is a valid EJA, but will double-check
- # if the user passes check_axioms=True.
- if "check_axioms" not in kwargs: kwargs["check_axioms"] = False
-
from mjo.hurwitz import ComplexMatrixAlgebra
A = ComplexMatrixAlgebra(n, scalars=field)
super().__init__(A, **kwargs)
from mjo.eja.eja_cache import complex_hermitian_eja_coeffs
a = complex_hermitian_eja_coeffs(self)
if a is not None:
- if self._rational_algebra is None:
- self._charpoly_coefficients.set_cache(a)
- else:
- self._rational_algebra._charpoly_coefficients.set_cache(a)
+ self.rational_algebra()._charpoly_coefficients.set_cache(a)
@staticmethod
def _max_random_instance_size(max_dimension):
"""
def __init__(self, n, field=AA, **kwargs):
- # We know this is a valid EJA, but will double-check
- # if the user passes check_axioms=True.
- if "check_axioms" not in kwargs: kwargs["check_axioms"] = False
-
from mjo.hurwitz import QuaternionMatrixAlgebra
A = QuaternionMatrixAlgebra(n, scalars=field)
super().__init__(A, **kwargs)
from mjo.eja.eja_cache import quaternion_hermitian_eja_coeffs
a = quaternion_hermitian_eja_coeffs(self)
if a is not None:
- if self._rational_algebra is None:
- self._charpoly_coefficients.set_cache(a)
- else:
- self._rational_algebra._charpoly_coefficients.set_cache(a)
+ self.rational_algebra()._charpoly_coefficients.set_cache(a)
from mjo.eja.eja_cache import octonion_hermitian_eja_coeffs
a = octonion_hermitian_eja_coeffs(self)
if a is not None:
- if self._rational_algebra is None:
- self._charpoly_coefficients.set_cache(a)
- else:
- self._rational_algebra._charpoly_coefficients.set_cache(a)
+ self.rational_algebra()._charpoly_coefficients.set_cache(a)
class AlbertEJA(OctonionHermitianEJA):
sage: from mjo.eja.eja_algebra import (random_eja,
....: CartesianProductEJA,
+ ....: ComplexHermitianEJA,
....: HadamardEJA,
....: JordanSpinEJA,
....: RealSymmetricEJA)
| b2 || 0 | 0 | b2 |
+----++----+----+----+
+ The "matrix space" of a Cartesian product always consists of
+ ordered pairs (or triples, or...) whose components are the
+ matrix spaces of its factors::
+
+ sage: J1 = HadamardEJA(2)
+ sage: J2 = ComplexHermitianEJA(2)
+ sage: J = cartesian_product([J1,J2])
+ sage: J.matrix_space()
+ The Cartesian product of (Full MatrixSpace of 2 by 1 dense
+ matrices over Algebraic Real Field, Module of 2 by 2 matrices
+ with entries in Algebraic Field over the scalar ring Algebraic
+ Real Field)
+ sage: J.one().to_matrix()[0]
+ [1]
+ [1]
+ sage: J.one().to_matrix()[1]
+ +---+---+
+ | 1 | 0 |
+ +---+---+
+ | 0 | 1 |
+ +---+---+
+
TESTS:
All factors must share the same base field::
sage: expected = J.one() # long time
sage: actual == expected # long time
True
-
"""
+ Element = CartesianProductEJAElement
def __init__(self, factors, **kwargs):
m = len(factors)
if m == 0:
if associative: category = category.Associative()
category = category.join([category, category.CartesianProducts()])
- # Compute my matrix space. This category isn't perfect, but
- # is good enough for what we need to do.
+ # Compute my matrix space. We don't simply use the
+ # ``cartesian_product()`` functor here because it acts
+ # differently on SageMath MatrixSpaces and our custom
+ # MatrixAlgebras, which are CombinatorialFreeModules. We
+ # always want the result to be represented (and indexed) as an
+ # ordered tuple. This category isn't perfect, but is good
+ # enough for what we need to do.
MS_cat = MagmaticAlgebras(field).FiniteDimensional().WithBasis()
MS_cat = MS_cat.Unital().CartesianProducts()
MS_factors = tuple( J.matrix_space() for J in factors )
self._inner_product_matrix = matrix.block_diagonal(
[J._inner_product_matrix for J in factors]
)
+ self._inner_product_matrix._cache = {'hermitian': True}
+ self._inner_product_matrix.set_immutable()
# Building the multiplication table is a bit more tricky
# because we have to embed the entries of the factors'
ones = tuple(J.one().to_matrix() for J in factors)
self.one.set_cache(self(ones))
+ def _sets_keys(self):
+ r"""
+
+ SETUP::
+
+ sage: from mjo.eja.eja_algebra import (ComplexHermitianEJA,
+ ....: RealSymmetricEJA)
+
+ TESTS:
+
+ The superclass uses ``_sets_keys()`` to implement its
+ ``cartesian_factors()`` method::
+
+ sage: J1 = RealSymmetricEJA(2,
+ ....: field=QQ,
+ ....: orthonormalize=False,
+ ....: prefix="a")
+ sage: J2 = ComplexHermitianEJA(2,field=QQ,orthonormalize=False)
+ sage: J = cartesian_product([J1,J2])
+ sage: x = sum(i*J.gens()[i] for i in range(len(J.gens())))
+ sage: x.cartesian_factors()
+ (a1 + 2*a2, 3*b0 + 4*b1 + 5*b2 + 6*b3)
+
+ """
+ # Copy/pasted from CombinatorialFreeModule_CartesianProduct,
+ # but returning a tuple instead of a list.
+ return tuple(range(len(self.cartesian_factors())))
+
def cartesian_factors(self):
# Copy/pasted from CombinatorialFreeModule_CartesianProduct.
return self._sets
return cartesian_product.symbol.join("%s" % factor
for factor in self._sets)
- def matrix_space(self):
- r"""
- Return the space that our matrix basis lives in as a Cartesian
- product.
-
- We don't simply use the ``cartesian_product()`` functor here
- because it acts differently on SageMath MatrixSpaces and our
- custom MatrixAlgebras, which are CombinatorialFreeModules. We
- always want the result to be represented (and indexed) as
- an ordered tuple.
-
- SETUP::
-
- sage: from mjo.eja.eja_algebra import (ComplexHermitianEJA,
- ....: HadamardEJA,
- ....: OctonionHermitianEJA,
- ....: RealSymmetricEJA)
-
- EXAMPLES::
-
- sage: J1 = HadamardEJA(1)
- sage: J2 = RealSymmetricEJA(2)
- sage: J = cartesian_product([J1,J2])
- sage: J.matrix_space()
- The Cartesian product of (Full MatrixSpace of 1 by 1 dense
- matrices over Algebraic Real Field, Full MatrixSpace of 2
- by 2 dense matrices over Algebraic Real Field)
-
- ::
-
- sage: J1 = ComplexHermitianEJA(1)
- sage: J2 = ComplexHermitianEJA(1)
- sage: J = cartesian_product([J1,J2])
- sage: J.one().to_matrix()[0]
- +---+
- | 1 |
- +---+
- sage: J.one().to_matrix()[1]
- +---+
- | 1 |
- +---+
-
- ::
-
- sage: J1 = OctonionHermitianEJA(1)
- sage: J2 = OctonionHermitianEJA(1)
- sage: J = cartesian_product([J1,J2])
- sage: J.one().to_matrix()[0]
- +----+
- | e0 |
- +----+
- sage: J.one().to_matrix()[1]
- +----+
- | e0 |
- +----+
-
- """
- return super().matrix_space()
-
@cached_method
def cartesian_projection(self, i):
SETUP::
- sage: from mjo.eja.eja_algebra import (HadamardEJA,
+ sage: from mjo.eja.eja_algebra import (FiniteDimensionalEJA,
+ ....: HadamardEJA,
....: JordanSpinEJA,
- ....: OctonionHermitianEJA,
....: RealSymmetricEJA)
EXAMPLES:
The ``cartesian_product()`` function only uses the first factor to
decide where the result will live; thus we have to be careful to
- check that all factors do indeed have a `_rational_algebra` member
- before we try to access it::
-
- sage: J1 = OctonionHermitianEJA(1) # no rational basis
- sage: J2 = HadamardEJA(2)
- sage: cartesian_product([J1,J2])
- Euclidean Jordan algebra of dimension 1 over Algebraic Real Field
- (+) Euclidean Jordan algebra of dimension 2 over Algebraic Real Field
- sage: cartesian_product([J2,J1])
- Euclidean Jordan algebra of dimension 2 over Algebraic Real Field
- (+) Euclidean Jordan algebra of dimension 1 over Algebraic Real Field
+ check that all factors do indeed have a ``rational_algebra()`` method
+ before we construct an algebra that claims to have a rational basis::
+
+ sage: J1 = HadamardEJA(2)
+ sage: jp = lambda X,Y: X*Y
+ sage: ip = lambda X,Y: X[0,0]*Y[0,0]
+ sage: b1 = matrix(QQ, [[1]])
+ sage: J2 = FiniteDimensionalEJA((b1,), jp, ip)
+ sage: cartesian_product([J2,J1]) # factor one not RationalBasisEJA
+ Euclidean Jordan algebra of dimension 1 over Algebraic Real
+ Field (+) Euclidean Jordan algebra of dimension 2 over Algebraic
+ Real Field
+ sage: cartesian_product([J1,J2]) # factor one is RationalBasisEJA
+ Traceback (most recent call last):
+ ...
+ ValueError: factor not a RationalBasisEJA
"""
def __init__(self, algebras, **kwargs):
+ if not all( hasattr(r, "rational_algebra") for r in algebras ):
+ raise ValueError("factor not a RationalBasisEJA")
+
CartesianProductEJA.__init__(self, algebras, **kwargs)
- self._rational_algebra = None
- if self.vector_space().base_field() is not QQ:
- if all( hasattr(r, "_rational_algebra") for r in algebras ):
- self._rational_algebra = cartesian_product([
- r._rational_algebra for r in algebras
- ])
+ @cached_method
+ def rational_algebra(self):
+ if self.base_ring() is QQ:
+ return self
+
+ return cartesian_product([
+ r.rational_algebra() for r in self.cartesian_factors()
+ ])
RationalBasisEJA.CartesianProduct = RationalBasisCartesianProductEJA