Hot-keys on this page
r m x p toggle line displays
j k next/prev highlighted chunk
0 (zero) top of page
1 (one) first highlighted chunk
""" Symmetric functions defined by orthogonality and triangularity.
One characterization of Schur functions is that they are upper triangularly related to the monomial symmetric functions and orthogonal with respect to the Hall scalar product. We can use the class SymmetricFunctionAlgebra_orthotriang to obtain the Schur functions from this definition.
::
sage: from sage.combinat.sf.sfa import zee sage: from sage.combinat.sf.orthotriang import SymmetricFunctionAlgebra_orthotriang sage: Sym = SymmetricFunctions(QQ) sage: m = Sym.m() sage: s = SymmetricFunctionAlgebra_orthotriang(Sym, m, zee, 's', 'Schur functions') sage: s([2,1])^2 s[2, 2, 1, 1] + s[2, 2, 2] + s[3, 1, 1, 1] + 2*s[3, 2, 1] + s[3, 3] + s[4, 1, 1] + s[4, 2]
::
sage: s2 = SymmetricFunctions(QQ).s() sage: s2([2,1])^2 s[2, 2, 1, 1] + s[2, 2, 2] + s[3, 1, 1, 1] + 2*s[3, 2, 1] + s[3, 3] + s[4, 1, 1] + s[4, 2] """ from __future__ import absolute_import #***************************************************************************** # Copyright (C) 2007 Mike Hansen <mhansen@gmail.com> # 2012 Mike Zabrocki <mike.zabrocki@gmail.com> # # Distributed under the terms of the GNU General Public License (GPL) # # This code is distributed in the hope that it will be useful, # but WITHOUT ANY WARRANTY; without even the implied warranty of # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU # General Public License for more details. # # The full text of the GPL is available at: # # http://www.gnu.org/licenses/ #*****************************************************************************
from . import sfa from sage.categories.morphism import SetMorphism from sage.categories.homset import Hom
class SymmetricFunctionAlgebra_orthotriang(sfa.SymmetricFunctionAlgebra_generic):
class Element(sfa.SymmetricFunctionAlgebra_generic.Element): pass
def __init__(self, Sym, base, scalar, prefix, basis_name, leading_coeff=None): r""" Initialization of the symmetric function algebra defined via orthotriangular rules.
INPUT:
- ``self`` -- a basis determined by an orthotriangular definition - ``Sym`` -- ring of symmetric functions - ``base`` -- an instance of a basis of the ring of symmetric functions (e.g. the Schur functions) - ``scalar`` -- a function ``zee`` on partitions. The function ``zee`` determines the scalar product on the power sum basis with normalization `\langle p_{\mu}, p_{\mu} \rangle = \mathrm{zee}(\mu)`. - ``prefix`` -- the prefix used to display the basis - ``basis_name`` -- the name used for the basis
.. NOTE::
The base ring is required to be a `\QQ`-algebra for this method to be useable, since the scalar product is defined by its values on the power sum basis.
EXAMPLES::
sage: from sage.combinat.sf.sfa import zee sage: from sage.combinat.sf.orthotriang import SymmetricFunctionAlgebra_orthotriang sage: Sym = SymmetricFunctions(QQ) sage: m = Sym.m() sage: s = SymmetricFunctionAlgebra_orthotriang(Sym, m, zee, 's', 'Schur'); s Symmetric Functions over Rational Field in the Schur basis
TESTS::
sage: TestSuite(s).run(elements = [s[1,1]+2*s[2], s[1]+3*s[1,1]]) sage: TestSuite(s).run(skip = ["_test_associativity", "_test_prod"]) # long time (7s on sage.math, 2011)
Note: ``s.an_element()`` is of degree 4; so we skip ``_test_associativity`` and ``_test_prod`` which involve (currently?) expensive calculations up to degree 12. """
def _base_to_self(self, x): """ Coerce a symmetric function in base ``x`` into ``self``.
INPUT:
- ``self`` -- a basis determined by an orthotriangular definition - ``x`` -- an element of the basis `base` to which ``self`` is triangularly related
OUTPUT:
- an element of ``self`` equivalent to ``x``
EXAMPLES::
sage: from sage.combinat.sf.sfa import zee sage: from sage.combinat.sf.orthotriang import SymmetricFunctionAlgebra_orthotriang sage: Sym = SymmetricFunctions(QQ) sage: m = Sym.m() sage: s = SymmetricFunctionAlgebra_orthotriang(Sym, m, zee, 's', 'Schur functions') sage: s._base_to_self(m([2,1])) -2*s[1, 1, 1] + s[2, 1] """
def _self_to_base(self, x): """ Coerce a symmetric function in ``self`` into base ``x``.
INPUT:
- ``self`` -- a basis determined by an orthotriangular definition - ``x`` -- an element of ``self`` as a basis of the ring of symmetric functions.
OUTPUT:
- the element ``x`` expressed in the basis to which ``self`` is triangularly related
EXAMPLES::
sage: from sage.combinat.sf.sfa import zee sage: from sage.combinat.sf.orthotriang import SymmetricFunctionAlgebra_orthotriang sage: Sym = SymmetricFunctions(QQ) sage: m = Sym.m() sage: s = SymmetricFunctionAlgebra_orthotriang(Sym, m, zee, 's', 'Schur functions') sage: s._self_to_base(s([2,1])) 2*m[1, 1, 1] + m[2, 1] """
def _base_cache(self, n): """ Computes the change of basis between ``self`` and base for the homogeneous component of size ``n``
INPUT:
- ``self`` -- a basis determined by an orthotriangular definition - ``n`` -- a nonnegative integer
EXAMPLES::
sage: from sage.combinat.sf.sfa import zee sage: from sage.combinat.sf.orthotriang import SymmetricFunctionAlgebra_orthotriang sage: Sym = SymmetricFunctions(QQ) sage: m = Sym.m() sage: s = SymmetricFunctionAlgebra_orthotriang(Sym, m, zee, 's', 'Schur functions') sage: s._base_cache(2) sage: l = lambda c: [ (i[0],[j for j in sorted(i[1].items())]) for i in sorted(c.items())] sage: l(s._base_to_self_cache[2]) [([1, 1], [([1, 1], 1)]), ([2], [([1, 1], -1), ([2], 1)])] sage: l(s._self_to_base_cache[2]) [([1, 1], [([1, 1], 1)]), ([2], [([1, 1], 1), ([2], 1)])] """ return else:
leading_coeff=self._leading_coeff, upper_triangular=True) self._base_to_self_cache, to_other_function = self._to_base)
def _to_base(self, part): """ Returns a function which takes in a partition `\mu` and returns the coefficient of a partition in the expansion of ``self`` `(part)` in base.
INPUT:
- ``self`` -- a basis determined by an orthotriangular definition - ``part`` -- a partition
.. note::
We assume that self._gram_schmidt has been called before self._to_base is called.
OUTPUT:
- a function which accepts a partition ``mu`` and returns the coefficients in the expansion of ``self(part)`` in the triangularly related basis.
EXAMPLES::
sage: from sage.combinat.sf.sfa import zee sage: from sage.combinat.sf.orthotriang import SymmetricFunctionAlgebra_orthotriang sage: Sym = SymmetricFunctions(QQ) sage: m = Sym.m() sage: s = SymmetricFunctionAlgebra_orthotriang(Sym, m, zee, 's', 'Schur functions') sage: m(s([2,1])) 2*m[1, 1, 1] + m[2, 1] sage: f = s._to_base(Partition([2,1])) sage: [f(p) for p in Partitions(3)] [0, 1, 2] """
def _multiply(self, left, right): """ Returns ``left`` * ``right`` by converting both to the base and then converting back to ``self``.
INPUT:
- ``self`` -- a basis determined by an orthotriangular definition - ``left``, ``right`` -- elements in ``self``
OUTPUT:
- the expansion of the product of ``left`` and ``right`` in the basis ``self``.
EXAMPLES::
sage: from sage.combinat.sf.sfa import zee sage: from sage.combinat.sf.orthotriang import SymmetricFunctionAlgebra_orthotriang sage: Sym = SymmetricFunctions(QQ) sage: m = Sym.m() sage: s = SymmetricFunctionAlgebra_orthotriang(Sym, m, zee, 's', 'Schur functions') sage: s([1])*s([2,1]) #indirect doctest s[2, 1, 1] + s[2, 2] + s[3, 1] """
# Backward compatibility for unpickling from sage.structure.sage_object import register_unpickle_override register_unpickle_override('sage.combinat.sf.orthotriang', 'SymmetricFunctionAlgebraElement_orthotriang', SymmetricFunctionAlgebra_orthotriang.Element) |