> > 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. ________________________________________________________________________ _______________________________________________ Axiom-mail mailing list [email protected] http://lists.nongnu.org/mailman/listinfo/axiom-mail
