Re: Problem with function geev

>>>>> "Thibault" == Thibault Langlois <[email protected]> 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



-------------------------------------------------------
SF.Net email is Sponsored by the Better Software Conference & EXPO September
19-22, 2005 * San Francisco, CA * Development Lifecycle Practices
Agile & Plan-Driven Development * Managing Projects & Teams * Testing & QA
Security * Process Improvement & Measurement * http://www.sqe.com/bsce5sf