passagemath-modules 10.5.32__cp310-cp310-macosx_14_0_arm64.whl → 10.6.20__cp310-cp310-macosx_14_0_arm64.whl

sage/modules/ore_module_element.py

@@ -0,0 +1,178 @@
+# sage_setup: distribution = sagemath-modules
+r"""
+Elements in Ore modules
+
+AUTHOR:
+
+- Xavier Caruso (2024-10)
+"""
+
+# ***************************************************************************
+#    Copyright (C) 2024 Xavier Caruso <xavier.caruso@normalesup.org>
+#
+#    This program is free software: you can redistribute it and/or modify
+#    it under the terms of the GNU General Public License as published by
+#    the Free Software Foundation, either version 2 of the License, or
+#    (at your option) any later version.
+#                  https://www.gnu.org/licenses/
+# ***************************************************************************
+
+from sage.misc.repr import repr_lincomb
+from sage.modules.free_module_element import FreeModuleElement_generic_dense
+
+
+class OreModuleElement(FreeModuleElement_generic_dense):
+    r"""
+    A generic element of a Ore module.
+    """
+    def _repr_(self):
+        r"""
+        Return a string representation of this element.
+
+        EXAMPLES::
+
+            sage: K.<t> = Frac(QQ['t'])
+            sage: S.<X> = OrePolynomialRing(K, K.derivation())
+            sage: M.<v,w> = S.quotient_module((t+1)*X^2 + 1)
+
+            sage: v - w  # indirect doctest
+            v - w
+            sage: w - v  # indirect doctest
+            -v + w
+            sage: X^5*v  # indirect doctest
+            ((-4*t+2)/(t^4+4*t^3+6*t^2+4*t+1))*v + ((t-5)/(t^3+3*t^2+3*t+1))*w
+        """
+        parent = self.parent()
+        names = parent._names
+        if parent._names is None:
+            return self.parent()._repr_element(self)
+        else:
+            return repr_lincomb([(names[i], self[i]) for i in range(len(names))])
+
+    def _latex_(self):
+        r"""
+        Return a LaTeX representation of this element.
+
+        EXAMPLES::
+
+            sage: K.<t> = Frac(QQ['t'])
+            sage: S.<X> = OrePolynomialRing(K, K.derivation())
+            sage: M.<v,w> = S.quotient_module((t+1)*X^2 + 1)
+
+            sage: latex(v - w)
+            v - w
+            sage: latex(w - v)
+            -v + w
+            sage: latex(X^5*v)
+            \left(\frac{-4 t + 2}{t^{4} + 4 t^{3} + 6 t^{2} + 4 t + 1}\right) v + \left(\frac{t - 5}{t^{3} + 3 t^{2} + 3 t + 1}\right) w
+        """
+        parent = self.parent()
+        if parent._names is None:
+            return self.parent()._latex_element(self)
+        else:
+            names = parent._latex_names
+            return repr_lincomb([(names[i], self[i]) for i in range(len(names))], is_latex=True)
+
+    def is_mutable(self) -> bool:
+        r"""
+        Always return ``False`` since elements in Ore modules
+        are all immutable.
+
+        EXAMPLES::
+
+            sage: K.<t> = Frac(QQ['t'])
+            sage: S.<X> = OrePolynomialRing(K, K.derivation())
+            sage: M = S.quotient_module(X^2 + t)
+
+            sage: v, w = M.basis()
+            sage: v
+            (1, 0)
+            sage: v.is_mutable()
+            False
+            sage: v[1] = 1
+            Traceback (most recent call last):
+            ...
+            ValueError: vectors in Ore modules are immutable
+        """
+        return False
+
+    def __setitem__(self, i, v):
+        r"""
+        Always raise an error since elements in Ore modules are
+        all immutable.
+
+        TESTS::
+
+            sage: K.<t> = Frac(QQ['t'])
+            sage: S.<X> = OrePolynomialRing(K, K.derivation())
+            sage: M.<v,w> = S.quotient_module(X^2 + t)
+            sage: w[1] = 0
+            Traceback (most recent call last):
+            ...
+            ValueError: vectors in Ore modules are immutable
+        """
+        raise ValueError("vectors in Ore modules are immutable")
+
+    def __hash__(self):
+        r"""
+        Return a hash of this element.
+
+        TESTS::
+
+            sage: K.<t> = Frac(QQ['t'])
+            sage: S.<X> = OrePolynomialRing(K, K.derivation())
+            sage: M.<v,w> = S.quotient_module(X^2 + t)
+            sage: hash(v)  # random
+            -5164621852614943976
+            sage: hash(w)  # random
+            -1950498447580522560
+        """
+        return hash(tuple(self))
+
+    def vector(self):
+        r"""
+        Return the coordinates vector of this element.
+
+        EXAMPLES::
+
+            sage: K.<t> = Frac(QQ['t'])
+            sage: S.<X> = OrePolynomialRing(K, K.derivation())
+            sage: M.<v,w> = S.quotient_module(X^2 + t)
+            sage: v.vector()
+            (1, 0)
+
+        We underline that this vector is not an element of the
+        Ore module; it lives in `K^2`. Compare::
+
+            sage: v.parent()
+            Ore module <v, w> over Fraction Field of Univariate Polynomial Ring in t over Rational Field twisted by d/dt
+            sage: v.vector().parent()
+            Vector space of dimension 2 over Fraction Field of Univariate Polynomial Ring in t over Rational Field
+        """
+        V = self.parent().module()
+        return V(self.list())
+
+    def image(self):
+        r"""
+        Return the image of this element by the pseudomorphism
+        defining the action of the Ore variable on this Ore module.
+
+        EXAMPLES::
+
+            sage: K.<t> = Frac(QQ['t'])
+            sage: S.<X> = OrePolynomialRing(K, K.derivation())
+            sage: M.<v,w> = S.quotient_module(X^2 + t)
+            sage: v.image()
+            w
+            sage: w.image()
+            -t*v
+
+        TESTS:
+
+        We check that this corresponds to the action of `X`::
+
+            sage: x = M.random_element()
+            sage: x.image() == X*x
+            True
+        """
+        return self.parent()._pseudohom(self)

sage/modules/ore_module_homspace.py

@@ -0,0 +1,147 @@
+# sage_setup: distribution = sagemath-modules
+# sage.doctest: needs sage.rings.finite_rings
+r"""
+Space of morphisms between Ore modules
+
+AUTHOR:
+
+- Xavier Caruso (2024-10)
+"""
+
+# ***************************************************************************
+#    Copyright (C) 2024 Xavier Caruso <xavier.caruso@normalesup.org>
+#
+#    This program is free software: you can redistribute it and/or modify
+#    it under the terms of the GNU General Public License as published by
+#    the Free Software Foundation, either version 2 of the License, or
+#    (at your option) any later version.
+#                  https://www.gnu.org/licenses/
+# ***************************************************************************
+
+from sage.structure.unique_representation import UniqueRepresentation
+from sage.categories.homset import HomsetWithBase
+from sage.matrix.matrix_space import MatrixSpace
+
+from sage.modules.ore_module import OreModule
+from sage.modules.ore_module_morphism import OreModuleMorphism
+
+
+class OreModule_homspace(UniqueRepresentation, HomsetWithBase):
+    r"""
+    Class for hom spaces between Ore modules.
+    """
+    Element = OreModuleMorphism
+
+    def __init__(self, domain, codomain, category=None):
+        r"""
+        Initialize this homspace.
+
+        INPUT:
+
+        - ``domain`` -- a Ore module
+
+        - ``codomain`` -- a Ore module
+
+        - ``category`` (default: ``None``) -- the category in which
+          the morphisms are
+
+        TESTS::
+
+            sage: K.<z> = GF(7^2)
+            sage: S.<X> = OrePolynomialRing(K, K.frobenius_endomorphism())
+            sage: M = S.quotient_module(X + z)
+            sage: N = S.quotient_module(X^2 + z)
+            sage: Hom(M, N)  # indirect doctest
+            Set of Morphisms
+            from Ore module of rank 1 over Finite Field in z of size 7^2 twisted by z |--> z^7
+            to Ore module of rank 2 over Finite Field in z of size 7^2 twisted by z |--> z^7
+            in Category of enumerated finite dimensional Ore modules with basis over Finite Field in z of size 7^2 twisted by z |--> z^7
+            sage: End(M)     # indirect doctest
+            Set of Morphisms
+            from Ore module of rank 1 over Finite Field in z of size 7^2 twisted by z |--> z^7
+            to Ore module of rank 1 over Finite Field in z of size 7^2 twisted by z |--> z^7
+            in Category of enumerated finite dimensional Ore modules with basis over Finite Field in z of size 7^2 twisted by z |--> z^7
+
+        ::
+
+            sage: V = M.module()
+            sage: Hom(M, V)
+            Traceback (most recent call last):
+            ...
+            ValueError: codomain must be a Ore module
+        """
+        if not isinstance(domain, OreModule):
+            raise ValueError("domain must be a Ore module")
+        if not isinstance(codomain, OreModule):
+            raise ValueError("codomain must be a Ore module")
+        if domain.ore_ring(action=False) is not codomain.ore_ring(action=False):
+            raise ValueError("domain and codomain must be defined over the same ring with same twisting maps")
+        super().__init__(domain, codomain, category)
+        base = domain.base_ring()
+        self._matrix_space = MatrixSpace(base, domain.dimension(), codomain.dimension())
+
+    def _element_constructor_(self, *args, **kwds):
+        r"""
+        Return a morphism in this homspace.
+
+        TESTS::
+
+            sage: K.<z> = GF(7^2)
+            sage: S.<X> = OrePolynomialRing(K, K.frobenius_endomorphism())
+            sage: M = S.quotient_module(X + z)
+            sage: H = End(M)
+            sage: H(2)
+            Ore module endomorphism of Ore module of rank 1 over Finite Field in z of size 7^2 twisted by z |--> z^7
+        """
+        return self.element_class(self, *args, **kwds)
+
+    def matrix_space(self):
+        r"""
+        Return the matrix space used to represent the
+        morphisms in this homspace.
+
+        EXAMPLES::
+
+            sage: K.<z> = GF(7^2)
+            sage: S.<X> = OrePolynomialRing(K, K.frobenius_endomorphism())
+            sage: M = S.quotient_module(X^3 + z*X + 1)
+            sage: End(M).matrix_space()
+            Full MatrixSpace of 3 by 3 dense matrices over Finite Field in z of size 7^2
+
+        ::
+
+            sage: N = S.quotient_module(X^2 + z)
+            sage: Hom(M, N).matrix_space()
+            Full MatrixSpace of 3 by 2 dense matrices over Finite Field in z of size 7^2
+        """
+        return self._matrix_space
+
+    def identity(self):
+        r"""
+        Return the identity morphism in this homspace.
+
+        EXAMPLES::
+
+            sage: K.<z> = GF(7^2)
+            sage: S.<X> = OrePolynomialRing(K, K.frobenius_endomorphism())
+            sage: M = S.quotient_module(X^3 + z*X + 1)
+            sage: End(M).identity()
+            Ore module endomorphism of Ore module of rank 3 over Finite Field in z of size 7^2 twisted by z |--> z^7
+        """
+        one = self.base_ring().one()
+        return self(one, check=False)
+
+    def zero(self):
+        r"""
+        Return the zero morphism in this homspace.
+
+        EXAMPLES::
+
+            sage: K.<z> = GF(7^2)
+            sage: S.<X> = OrePolynomialRing(K, K.frobenius_endomorphism())
+            sage: M = S.quotient_module(X^3 + z*X + 1)
+            sage: End(M).zero()
+            Ore module endomorphism of Ore module of rank 3 over Finite Field in z of size 7^2 twisted by z |--> z^7
+        """
+        zero = self.base_ring().zero()
+        return self(zero, check=False)