serious bug in eigenvalues and eigenvectors
- Subject: serious bug in eigenvalues and eigenvectors
- From: go_furuya at infoseek
- Date: Tue, 11 May 2004 14:05:20 +0900
hi all
5.9.0.1cvs version cannnot calc eigenvalues for
matrix([3,1,0,0],[-4,-1,0,0],[7,1,2,1],[-17,-6,-1,0]),
but 5.9.0 can calc this. [[1],[4]] is OK.
Though share/matrix/eigen.mac was not changed.
I detected this bug,when I checked my diag,dem.
ver 5.9.0 can do . for example batch(2diag.dem")
(C7) A4 : MATRIX([3, 1, 0, 0], [- 4, - 1, 0, 0], [7, 1, 2, 1],[- 17, - 6, - 1, 0])
(C33) jordan(A4)
(D33) [[1, 2, 2]]
(C34) dispJordan(jordan(A4))
[ 1 1 0 0 ]
[ ]
[ 0 1 0 0 ]
(D34) [ ]
[ 0 0 1 1 ]
[ ]
[ 0 0 0 1 ]
C35) mm : ModeMatrix(A4, jordan(A4))
[ 2 1 1 0 ]
[ ]
[ - 4 0 - 2 1 ]
(D35) [ ]
[ 7 0 1 0 ]
[ ]
[ - 17 0 - 6 0 ]
<- 1>
(C36) mm . (A4 . mm)
[ 1 1 0 0 ]
[ ]
[ 0 1 0 0 ]
(D36) [ ]
[ 0 0 1 1 ]
[ ]
[ 0 0 0 1 ]
this is correct.
but in 5.9.0.1 cvs
(C6) eigenvalues(a4);
PART fell off end.
#0: EIGENVALUES(mat=MATRIX([3,1,0,0],[-4,-1,0,0],[7,1,2,1],[-17,-6,-1,0]))(eigen.mac line 92)
-- an error. Quitting. To debug this try DEBUGMODE(TRUE);)
G.F
------------------------------------------------------------------------
$BO"5YCf!!$"$N;v7o!"$=$N8e!!$I$&$J$C$?!)!!(Bby infoseek
http://ap.infoseek.co.jp/news6.html