#18453: Infinite affine crystals should use extended weight lattice
-------------------------------------------------+-------------------------
Reporter: bump | Owner:
Type: defect | Status: closed
Priority: major | Milestone: sage-6.8
Component: combinatorics | Resolution: fixed
Keywords: crystals, days65 | Merged in:
Authors: Ben Salisbury, Anne | Reviewers: Dan Bump
Schilling, Travis Scrimshaw | Work issues:
Report Upstream: N/A | Commit:
Branch: | Stopgaps:
ee398efc22755947ea98f0a8250333d21352f2db |
Dependencies: |
-------------------------------------------------+-------------------------
Changes (by nthiery):
* commit: ee398efc22755947ea98f0a8250333d21352f2db =>
Comment:
Replying to [comment:19 aschilling]:
> Replying to [comment:17 nthiery]:
> > Replying to [comment:8 aschilling]:
> > > I had to fix the Weyl dimension formula (which I guess was written
by Nicolas and only worked in the ambient space).
> >
> > By Dan, if I recall correctly.
>
> {{{
> 7d693534 (Nicolas M. Thiery 2012-03-19 21:38:26 +0100 859) def
weyl_dimension(self, highest_weight):
> 7d693534 (Nicolas M. Thiery 2012-03-19 21:38:26 +0100 860)
"""
> 7d693534 (Nicolas M. Thiery 2012-03-19 21:38:26 +0100 861)
EXAMPLES::
> 7d693534 (Nicolas M. Thiery 2012-03-19 21:38:26 +0100 862)
> 7d693534 (Nicolas M. Thiery 2012-03-19 21:38:26 +0100 863)
sage: RootSystem(['A',3]).ambient_lattice().weyl_dimension([2,1,0,0])
> 7d693534 (Nicolas M. Thiery 2012-03-19 21:38:26 +0100 864)
20
> 7d693534 (Nicolas M. Thiery 2012-03-19 21:38:26 +0100 865)
> 7d693534 (Nicolas M. Thiery 2012-03-19 21:38:26 +0100 866)
sage:
type(RootSystem(['A',3]).ambient_lattice().weyl_dimension([2,1,0,0]))
> 7d693534 (Nicolas M. Thiery 2012-03-19 21:38:26 +0100 867)
<type 'sage.rings.integer.Integer'>
> 7d693534 (Nicolas M. Thiery 2012-03-19 21:38:26 +0100 868)
"""
> 7d693534 (Nicolas M. Thiery 2012-03-19 21:38:26 +0100 869)
highest_weight = self(highest_weight)
> f562ca2b (Travis Scrimshaw 2014-10-01 17:16:31 -0700 870)
if not highest_weight.is_dominant():
> f562ca2b (Travis Scrimshaw 2014-10-01 17:16:31 -0700 871)
raise ValueError("the highest weight must be dominant")
> 7d693534 (Nicolas M. Thiery 2012-03-19 21:38:26 +0100 872)
rho = self.rho()
> 7d693534 (Nicolas M. Thiery 2012-03-19 21:38:26 +0100 873)
n = prod([(rho+highest_weight).dot_product(x) for x in
self.positive_roots()])
> 7d693534 (Nicolas M. Thiery 2012-03-19 21:38:26 +0100 874)
d = prod([ rho.dot_product(x) for x in self.positive_roots()])
> 7d693534 (Nicolas M. Thiery 2012-03-19 21:38:26 +0100 875)
from sage.rings.integer import Integer
> 7d693534 (Nicolas M. Thiery 2012-03-19 21:38:26 +0100 876)
return Integer(n/d)
> }}}
This got me curious about my memory failing. Digging back in history ended
up in a mercurial patch from early 2008 from which one can't extract the
author. As far as I remember, I refactored this code at different
occasions (in the above to move it from
weight_lattice_realization to weight_lattice_realizations when
creating the categories), but the original code was by Dan in the
`WeylDim` function. Anyway just ranting this has no relevance.
--
Ticket URL: <http://trac.sagemath.org/ticket/18453#comment:61>
Sage <http://www.sagemath.org>
Sage: Creating a Viable Open Source Alternative to Magma, Maple, Mathematica,
and MATLAB
--
You received this message because you are subscribed to the Google Groups
"sage-trac" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
To post to this group, send email to [email protected].
Visit this group at http://groups.google.com/group/sage-trac.
For more options, visit https://groups.google.com/d/optout.