Re: [Matlisp-users] Problem with function geev

[Raymond T. <ray...@er...>]

>>>>> "Thibault" == Thibault Langlois <tl...@di...> writes:

    Thibault> #<The MATLISP-USER package, 1/9 internal, 0/2 external>
    Thibault> * (setq m (make-real-matrix '((1 1 2) (2 2 4) (3 3 6))))

    Thibault> Warning:  Declaring M special.

    Thibault> #<REAL-MATRIX  3 x 3
    Thibault>     1.0000      1.0000      2.0000     
    Thibault>     2.0000      2.0000      4.0000     
    Thibault>     3.0000      3.0000      6.0000     >
    Thibault> * (geev m)

Just tried this (finally) using cvs matlisp and cmucl 2005-07.
(geev m) returns

#<REAL-MATRIX  3 x 1
    0.0000     
    9.0000     
   -5.24537E-16 >

which is pretty close to what matlab returns.

I'll have to look at the code some more.  Perhaps it's doing something
wrong.  (There used to be a memory corruption error in some eigenvalue
routine, but that's been fixed.)

Ray