Dear Dan and Nicolas,
Thank you for your comments! I posted a revised patch for #10632 on the
sage-combinat
queue and the trac server which incorporates the suggestions below as well as
Nicolas' suggestions regarding the latex options/printing for crystals.
The new patch now depends on #10723 which is available on the sage-combinat
server.
Best wishes,
Anne
On 2/2/11 4:59 PM, Daniel Bump wrote:
I'm looking at #10632.
Comment 1.
After the patch, I observe the following. Line 757 of tensor_product.py
contains:
shapes = tuple( tuple(shape) for shape in shapes )
I think here the making of a tuple is important because a tuple is immutable,
a list is not.
To illustrate, suppose take this line out. Then
CrystalOfTableaux("B3",shape=[1,0,0]) produces a
crash. Something about a list not being an unhashable type.
Now the line 764, which is new code, contains
shapes = [[int(i) for i in sh] for sh in shapes].
Then shapes is potentially further modified in Type D.
It seems to me that this should also be turned into
tuples. This would have to be done immediately before
B is defined. Then when B is defined, one could call
CrystalOfTableauxGeneric instead of CrystalOfTableaux.
Comment 2.
Currently we have:
sage: CrystalOfTableaux("D4",shape=[1,0,0,0])
The crystal of tableaux of type ['D', 4] and shape(s) [[1, 0, 0, 0]]
sage: CrystalOfTableaux("D4",shape=[1/2,1/2,1/2,-1/2])
The tensor product of the crystals (The minus crystal of spins for type ['D',
4], The crystal of tableaux of type ['D', 4] and shape(s) [[0, 0, 0, 0]])
Wouldn't it be better to get:
sage: CrystalOfTableaux("D4",shape=[1/2,1/2,1/2,-1/2])
The crystal of tableaux of type ['D', 4] and shape(s) [[1/2, 1/2, 1/2, -1/2]]
Dan
--
You received this message because you are subscribed to the Google Groups
"sage-combinat-devel" group.
To post to this group, send email to sage-combinat-devel@googlegroups.com.
To unsubscribe from this group, send email to
sage-combinat-devel+unsubscr...@googlegroups.com.
For more options, visit this group at
http://groups.google.com/group/sage-combinat-devel?hl=en.
