Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

bug in border case in highest weight vectors of tensor product of crystals #30493

Closed
mantepse opened this issue Sep 2, 2020 · 7 comments
Closed

Comments

@mantepse
Copy link
Contributor

mantepse commented Sep 2, 2020

sage: CW = CartanType("G", 2)
sage: C = crystals.Letters(CW)
sage: C.highest_weight_vectors()
(1,)
sage: T = crystals.TensorProduct(C)
sage: T.highest_weight_vectors()
---------------------------------------------------------------------------
IndexError                                Traceback (most recent call last)
<ipython-input-33-feba30ccc2e6> in <module>()
----> 1 T.highest_weight_vectors()

/home/martin/sage-develop/local/lib/python3.7/site-packages/sage/misc/cachefunc.pyx in sage.misc.cachefunc.CachedMethodCallerNoArgs.__call__ (build/cythonized/sage/misc/cachefunc.c:12717)()
   2308         if self.cache is None:
   2309             f = self.f
-> 2310             self.cache = f(self._instance)
   2311         return self.cache
   2312 

/home/martin/sage-develop/local/lib/python3.7/site-packages/sage/categories/highest_weight_crystals.py in highest_weight_vectors(self)
    713                      [[[-2]], [[1, 1], [2, 2]]])
    714                 """
--> 715                 return tuple(self.highest_weight_vectors_iterator())
    716 
    717             def highest_weight_vectors_iterator(self):

/home/martin/sage-develop/local/lib/python3.7/site-packages/sage/categories/highest_weight_crystals.py in highest_weight_vectors_iterator(self)
    786                     T_phi = [None]*(m-1) + [{i: b.phi(i) for i in I}]
    787                     while T_pos < m:
--> 788                         if T_cur[T_pos] == T_len[T_pos]:
    789                             T_cur[T_pos] = 0
    790                             T_pos += 1

IndexError: list index out of range
sage: T = crystals.TensorProduct(C, C)
sage: T.highest_weight_vectors()
([1, 1], [2, 1], [0, 1], [-1, 1])

CC: @tscrim

Component: combinatorics

Author: Travis Scrimshaw

Branch/Commit: 86b067d

Reviewer: Martin Rubey

Issue created by migration from https://trac.sagemath.org/ticket/30493

@mantepse mantepse added this to the sage-9.2 milestone Sep 2, 2020
@tscrim
Copy link
Collaborator

tscrim commented Sep 3, 2020

comment:1

Simple fix to handle the corner case.


New commits:

86b067dFixing corner case of highest weight vectors iterator.

@tscrim
Copy link
Collaborator

tscrim commented Sep 3, 2020

Commit: 86b067d

@tscrim
Copy link
Collaborator

tscrim commented Sep 3, 2020

@tscrim
Copy link
Collaborator

tscrim commented Sep 3, 2020

Author: Travis Scrimshaw

@mantepse
Copy link
Contributor Author

mantepse commented Sep 3, 2020

comment:2

Thank you!

@mantepse
Copy link
Contributor Author

mantepse commented Sep 3, 2020

Reviewer: Martin Rubey

@vbraun
Copy link
Member

vbraun commented Sep 23, 2020

Changed branch from public/crystals/tensor_product_hw_elts-30493 to 86b067d

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

No branches or pull requests

3 participants