This is the mail archive of the
gsl-discuss@sourceware.org
mailing list for the GSL project.
Re: unsymmetric eigenvalue problem
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