Hi Bruce, If you are only interested in the classical weight, then you can do
sage: C = KirillovReshetikhinCrystal(['A',1,1],1,2) sage: B = C.classical_decomposition() sage: T = TensorProductOfCrystals(*[B]*3) sage: L = CombinatorialFreeModule(ZZ,B.weight_lattice_realization()) sage: hw = [a for a in T if a.is_highest_weight()] sage: sum( L.term( a.weight(), 1 ) for a in hw ) B[(3, 3)] + 3*B[(4, 2)] + 2*B[(5, 1)] + B[(6, 0)] Best wishes, Anne On 11/23/13 2:13 AM, Bruce wrote: > Sorry about being difficult. > > Here is one command: > > lie.p_tensor(3,[2],'A1') > > Here is an alternative: > > C = KirillovReshetikhinCrystal(['A',1,1],1,2) > L = CombinatorialFreeModule(ZZ,C.weight_lattice_realization()) > T = TensorProductOfCrystals(*[C]*3) > hw = [ a for a in T if a.e(1) == None ] > sum( L.term( a.weight(), 1 ) for a in hw ) > > I would like to convince sage (in this simplified example) that these are > "the same". > > Thank you for your patience. > > On Saturday, November 23, 2013 7:28:26 AM UTC, Nicolas M. Thiery wrote: > > On Fri, Nov 22, 2013 at 08:09:25AM -0800, Bruce wrote: > > Here is the result of the first calculation (using the KR crystal) > and its > > parent: > > Please, not the result but the command (or a simplified version) > producing the result! Otherwise one has to reconstruct the command to > play with the objects :-) > > Cheers, > Nicolas -- You received this message because you are subscribed to the Google Groups "sage-combinat-devel" group. To unsubscribe from this group and stop receiving emails from it, send an email to sage-combinat-devel+unsubscr...@googlegroups.com. To post to this group, send email to sage-combinat-devel@googlegroups.com. Visit this group at http://groups.google.com/group/sage-combinat-devel. For more options, visit https://groups.google.com/groups/opt_out.
