Hello Florent!
On Wed, Aug 29, 2012 at 11:00 AM, Florent Hivert <florent.hiv...@lri.fr> wrote:
> Hi Franco,
>
>> >> 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]
Franco
--
--
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.
