Hi, This one looks really serious. In some ways seems to be computer dependent because my answers seem to have twice as many digits as yours.
sage: Sym = SymmetricFunctions(QQ) sage: p = Sym.p() sage: s = Sym.s() sage: s(p[2,2]) s[1, 1, 1, 1] - s[2, 1, 1] + 2*s[2, 2] - s[3, 1] + s[4] sage: g = s(p([1]*47)) sage: s(p[2,2]) s[1, 1, 1, 1] - s[2, 1, 1] - 1152912726858932132*s[2, 2] - s[3, 1] + s[4] sage: s(p[2,2]) s[1, 1, 1, 1] - s[2, 1, 1] + 3458773291570539007*s[2, 2] - s[3, 1] + s[4] sage: s(p[2,2]) s[1, 1, 1, 1] - s[2, 1, 1] - 1152912726855831925*s[2, 2] - s[3, 1] + s[4] It doesn't help to fix the problem by redefining the ring or 'reseting' sage: sage: reset() sage: Sym = SymmetricFunctions(QQ['t']) sage: s = Sym.s() sage: p = Sym.p() sage: s(p[2,2]) s[1, 1, 1, 1] - s[2, 1, 1] - 1152912726855247442*s[2, 2] - s[3, 1] + s[4] I have a suspicion that this is coming from symmetrica. If we can track the error down to it, it may require another large rewrite of the symmetric function code. -Mike On Wednesday, 29 August 2012 11:40:57 UTC-4, Anne Schilling wrote: > > On 8/29/12 8:20 AM, Franco Saliola wrote: > >>>>> Someone sent me the following sage session, which I cannot > reproduce, > >>>>> but I'm asking whether this is a known issue and whether someone can > >>>>> reproduce it: > >>>>> > >>>>> > >>>>> sage: p=SFAPower(QQ) > >>>>> sage: s=SFASchur(QQ) > >>>>> > >>>>> sage: f = p([1,1,1,1,1,1,1,1,1,1,1])/990 + p([2,2,2,2,2,1])/10 + > >>>>> p([5,5,1])/10 - p([10,1])/10 + p([9,1,1])/3 + p([3,3,3,1,1])/9 + > >>>>> 5*p([11])/11 > >>>>> > >>>>> sage: [f.scalar(s(q)) for q in Partitions(11)] > >>>>> [1, 0, 0, 0, 0, 0, 0, 1, 0, 2, 0, 2, 0, 1, 1, 1, 1, 0, 0, -1, 2, 1, > >>>>> -1, 2, 4, 3, 0, 1, 0, 2, 1, 0, 1, 1, 2, 1, 127020063634, 1, 1, 0, 0, > >>>>> 158403414193, 0, 1, -1, 1, 0, 19558976098, 0, 1, 3804836920632, > >>>>> 298744000054, 6021978496, -1, 0, 1] > >>>>> > >>>>> Expected: > >>>>> [1, 0, 0, 0, 0, 0, 0, 1, 0, 2, 0, 2, 0, 1, 1, 1, 1, 0, 0, -1, 2, 1, > >>>>> -1, 2, 4, 3, 0, 1, 0, 2, 1, 0, 1, 1, 2, 1, 2, 1, 1, 0, 0, 3, 0, 1, > -1, > >>>>> 1, 0, 2, 0, 1, 0, 0, 0, -1, 0, 1] > >> > >> I'm wondering: did this happen during a session, or is the bug > triggered > >> deterministically after the preceding sequences of command ? > > > > It sure looks likely, based on the following session (with 5.3.beta2). > > > > To begin with, let's do a change of basis in a small degree: > > > > sage: p = SymmetricFunctions(QQ).powersum() > > sage: s = SymmetricFunctions(QQ).schur() > > sage: s(p[2,2]) > > s[1, 1, 1, 1] - s[2, 1, 1] + 2*s[2, 2] - s[3, 1] + s[4] > > > > Now let's do one in a larger degree: > > > > sage: time g = s(p([1]*47)) > > Time: CPU 19.16 s, Wall: 20.91 s > > > > Now let's do that first one again: > > > > sage: s(p[2,2]) > > s[1, 1, 1, 1] - s[2, 1, 1] + 4571483302*s[2, 2] - s[3, 1] + s[4] > > > > That's not the correct answer ! And the next time you ask Sage, it > > gives different, still incorrect, answers: > > > > sage: s(p[2,2]) > > s[1, 1, 1, 1] - s[2, 1, 1] + 4614252243*s[2, 2] - s[3, 1] + s[4] > > sage: s(p[2,2]) > > s[1, 1, 1, 1] - s[2, 1, 1] + 4634718110*s[2, 2] - s[3, 1] + s[4] > > sage: s(p[2,2]) > > s[1, 1, 1, 1] - s[2, 1, 1] + 4631047636*s[2, 2] - s[3, 1] + s[4] > > I can confirm that this happens for me as well (though with different > numbers): > > sage: p = SymmetricFunctions(QQ).powersum() > sage: s = SymmetricFunctions(QQ).schur() > sage: s(p[2,2]) > s[1, 1, 1, 1] - s[2, 1, 1] + 2*s[2, 2] - s[3, 1] + s[4] > sage: time g = s(p([1]*47)) > Time: CPU 14.38 s, Wall: 14.87 s > sage: s(p[2,2]) > s[1, 1, 1, 1] - s[2, 1, 1] + 4570475192*s[2, 2] - s[3, 1] + s[4] > sage: s(p[2,2]) > s[1, 1, 1, 1] - s[2, 1, 1] + 4614292718*s[2, 2] - s[3, 1] + s[4] > > Best, > > Anne > -- You received this message because you are subscribed to the Google Groups "sage-combinat-devel" group. To view this discussion on the web visit https://groups.google.com/d/msg/sage-combinat-devel/-/hJTez9BLxXgJ. 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.
