[sage-support] Re: Jordan Normal Form Problem

Hi Jason, thanks, you hit the nail on the head - I had not downloaded the latest version..it works fine now, and I've learned a valuable lesson! Ciaran On Jan 19, 11:28 pm, Jason Grout wrote: > c mullan wrote: > > Hi all, > > > by general theory I know that an in

[sage-support] Jordan Normal Form Problem

Hi all, by general theory I know that an invertible transformation matrix P exists such that A = ~P*J*P where J is the Jordan Normal Form of a square matrix A. When I try to calculate P, some strange things happen.. M=MatrixSpace(GF(2),7) A=M.random_element() f=A.charpoly() d = lcm([p.degree() f

[sage-support] Re: Jordan Normal Form Block Sizes

Excellent, yes, this solution will work for me! Thanks for your swift reply! On Jan 15, 11:12 pm, Jason Grout wrote: > c mullan wrote: > > Hi, > > > Suppose I compute the Jordan Normal Form of a matrix, > > > A.jordan_form() > > > Then in the output I can se

[sage-support] Jordan Normal Form Block Sizes

Hi, Suppose I compute the Jordan Normal Form of a matrix, A.jordan_form() Then in the output I can see that the block sizes are indicated (by subdivide='True'), but I can't extract this information. I would like a list of block sizes, (e.g. [2,2,1,1,1] for a 7x7 matrix). I cannot figure out ho