This is the mail archive of the gsl-discuss@sourceware.org mailing list for the GSL project.


Index Nav: [Date Index] [Subject Index] [Author Index] [Thread Index]
Message Nav: [Date Prev] [Date Next] [Thread Prev] [Thread Next]
Other format: [Raw text]

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


Index Nav: [Date Index] [Subject Index] [Author Index] [Thread Index]
Message Nav: [Date Prev] [Date Next] [Thread Prev] [Thread Next]