On 1/30/07, David Joyner <[EMAIL PROTECTED]> wrote:
>
> Hello all:
>
> In trying to compute the eigenvalues and eigenvectors of the
> discrete Fourier transformation for Z/5Z for a class I'm teaching,
> I came across the following:
>
> sage: MS = MatrixSpace(CyclotomicField(5),5,5)
> sage: z = CyclotomicField(5).gen()
> sage: F5 = 
> MS([[1,1,1,1,1],[1,z,z^2,z^3,z^4],[1,z^2,z^4,z^6,z^8],[1,z^3,z^6,z^9,z^(12)],[1,z^4,z^8,z^(12),z^(16)]])
> sage: F5.eigenspaces()
> ---------------------------------------------------------------------------
> <class 'gen.PariError'>                   Traceback (most recent call last)
>
> /mnt/hd200/sagefiles/sage-1.5/<ipython console> in <module>()
>
> /mnt/hd200/sagefiles/sage-1.5/matrix2.pyx in matrix2.Matrix.eigenspaces()
>
> /mnt/hd200/sagefiles/sage-1.5/matrix2.pyx in matrix2.Matrix.kernel()
>
> /mnt/hd200/sagefiles/sage-1.5/sage_object.pyx in 
> sage_object.SageObject._pari_()
>
> /mnt/hd200/sagefiles/sage-1.5/gen.pyx in gen._pari_trap()
>
> <class 'gen.PariError'>: forbidden (45)

This seems to have been corrected as of 2.0, so I would suggest you upgrade :) :

sage: F5.eigenspaces()
[
(2*zeta5^3 + 2*zeta5^2 + 1, [
(1, 1/2*zeta5^3 + 1/2*zeta5^2, 1/2*zeta5^3 + 1/2*zeta5^2, 1/2*zeta5^3
+ 1/2*zeta5^2, 1/2*zeta5^3 + 1/2*zeta5^2)
]),
(-2*zeta5^3 - 2*zeta5^2 - 1, [
(1, 0, -zeta5^3 - zeta5^2 - 1, -zeta5^3 - zeta5^2 - 1, 0),
(0, 1, -1, -1, 1)
]),
(a2, [

])
]

--~--~---------~--~----~------------~-------~--~----~
To post to this group, send email to sage-support@googlegroups.com
To unsubscribe from this group, send email to [EMAIL PROTECTED]
For more options, visit this group at 
http://groups.google.com/group/sage-support
URLs: http://sage.math.washington.edu/sage/ and http://sage.scipy.org/sage/
-~----------~----~----~----~------~----~------~--~---

Reply via email to