There do seem to be slight variations between results from the Windows and
Linux versions of LAPACK. They are mathematically insignificant, but have
an effect if you are erroneously relying on order of results. Here's my
revised code for calculating a principal eigenvector using LAPACK.
require '~addons/lapack/lapack.ijs'
require '~addons/lapack/dgeev.ijs'
la=:3 : 0
'l c r'=.dgeev_jlapack_ y.
({. \: c){"1 r
)
und=: >(1,1,5);(1,1,3);((%5),(%3),1) NB. "Understandability
(%+/)power und
0.48064 0.405388 0.113972
%+/)la und
0.48064 0.405388 0.113972
Best wishes,
John
----------------------------------------------------------------------
For information about J forums see http://www.jsoftware.com/forums.htm
