unsymmetric eigenvalue problem

Brian Gough bjg@network-theory.co.uk
Tue May 23 18:35:00 GMT 2006


You can find my test program here
http://www.briangough.ukfsn.org/unsymm/
It's not totally finished, example output below.

-- 
Brian Gough

Network Theory Ltd,
Publishing Free Software Manuals --- http://www.network-theory.co.uk/

$ ./a.out -n 3 -i -l -1 -u 1
testing N=3 incrementally on element range [-1,1]
==== CASE 17 ==========================

 A = [ 
         0,         1,         0;
        -1,        -1,        -1;
        -1,        -1,        -1;
]

=== eval - lapack dgees ===
-1.813985266369727769e-32  0.000000000000000000e+00
-9.999999999999996669e-01  1.950115636758225308e-08
-9.999999999999996669e-01 -1.950115636758225308e-08

=== eval - alken ===
-6.631056343480628004e-33  0.000000000000000000e+00
-9.999999999999997780e-01  1.633303231307825499e-08
-9.999999999999997780e-01 -1.633303231307825499e-08

max abserr = 3.16812e-09  max relerr = 0.162458



More information about the Gsl-discuss mailing list