[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Axiom-mail] The Cayley-Hamilton theorem and finite fields - a small
From: |
Themos Tsikas |
Subject: |
Re: [Axiom-mail] The Cayley-Hamilton theorem and finite fields - a small problem |
Date: |
Mon, 16 Jul 2007 11:07:20 +0000 |
User-agent: |
KMail/1.9.1 |
>
> What's going on, and why?
>
Don't know! I suspect that the interpreter is constructing some coercion that
is doomed to failure. You can do the problem if you are more careful with
your types
)cl c
D:=FFP(PF 2,x^4+x+1)
M:=matrix([[random()$D for i in 1..4] for j in 1..4])
p:=characteristicPolynomial(M,y)
SM:=SquareMatrix(4,D)
sp:=map(z+->(squareMatrix$SM)diagonalMatrix([z,z,z,z]),p)
sm:=squareMatrix(M)$SM
eval(sp,y=sm)
)cl c
D:=INT
M:=matrix([[random()$D for i in 1..4] for j in 1..4])
p:=characteristicPolynomial(M,y)
SM:=SquareMatrix(4,D)
sp:=map(z+->(squareMatrix$SM)diagonalMatrix([z,z,z,z]),p)
sm:=squareMatrix(M)$SM
eval(sp,y=sm)
)cl c
D:=PF 7
M:=matrix([[random()$D for i in 1..4] for j in 1..4])
p:=characteristicPolynomial(M,y)
SM:=SquareMatrix(4,D)
sp:=map(z+->(squareMatrix$SM)diagonalMatrix([z,z,z,z]),p)
sm:=squareMatrix(M)$SM
eval(sp,y=sm)
)cl c
D:=FF(2,4)
M:=matrix([[random()$D for i in 1..4] for j in 1..4])
p:=characteristicPolynomial(M,y)
SM:=SquareMatrix(4,D)
sp:=map(z+->(squareMatrix$SM)diagonalMatrix([z,z,z,z]),p)
sm:=squareMatrix(M)$SM
eval(sp,y=sm)
________________________________________________________________________
The Numerical Algorithms Group Ltd is a company registered in England
and Wales with company number 1249803. The registered office is:
Wilkinson House, Jordan Hill Road, Oxford OX2 8DR, United Kingdom.
This e-mail has been scanned for all viruses by Star. The service is
powered by MessageLabs.
________________________________________________________________________