Hi Bruce, You can also do
sage: C = KirillovReshetikhinCrystal(['A',1,1],1,2) sage: B = C.classical_decomposition() sage: L = CombinatorialFreeModule(ZZ,B.weight_lattice_realization()) sage: T = TensorProductOfCrystals(*[C]*3) sage: hw = [a for a in T if a.is_highest_weight(index_set=[1])] sage: sum( L.term( sum(a.lift().weight() for a in b), 1 ) for b in hw ) B[(3, 3)] + 3*B[(4, 2)] + 2*B[(5, 1)] + B[(6, 0)] Then you can still compute the energy function in T. Best wishes, Anne On 11/24/13 8:42 AM, Bruce wrote: > Although I do not have an answer to the question I originally asked, I can > now do the calculation I wanted to do using Anne's suggestion of using > C.classical_decomposition(). Many thanks to Nicolas > and Anne. -- 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.