Hi there,
Sage is producing a wrong answer when asked to compute a Frobenius
normal form of a matrix, with the transformation matrix:
sage: sage: A=matrix(ZZ,8,[[6,0,-2,4,0,0,0,-2],[14,-1,0,6,0,-1,-1,1],\
[2,2,0,1,0,0,1,0],[-12,0,5,-8,0,0,0,4],[0,4,0,0,0,0,4,0],\
[0,0,0,0,1,0,0,0],[-14,2,0,-6,0,2,2,-1],[-4,0,2,-4,0,0,0,4]])
sage: F,K=A.frobenius(2)
sage: ~K*F*K==a
False
And the answer should be True.
As Sage directly wraps PARI's frobenius form routine, I guess the bug is
there.
Before I post a ticket, has anyone heard about/met this bug before? Is
someone already working on it?
Cheers,
Clément
--
To post to this group, send an email to sage-devel@googlegroups.com
To unsubscribe from this group, send an email to
sage-devel+unsubscr...@googlegroups.com
For more options, visit this group at http://groups.google.com/group/sage-devel
URL: http://www.sagemath.org
