Devon McCormick wrote:
> This is odd. Even in a new session, I get
>
> NB. [Session started: 2006 3 1 18 21 37.817]
> load '~addons/lapack/lapack.ijs'
> load '~addons/lapack/dgeev.ijs'
> ]rod=. >(1,(%3),%9);(3,1,%3);9 3 1 NB. "Replication of detail", p. 7
> 1 0.33333333 0.11111111
> 3 1 0.33333333
> 9 3 1
> (%+/){."1 >2{dgeev_jlapack_ rod
> _0.2 0.3 0.9
> (9!:12 '');9!:14 ''
> +-+---------------------+
> |6|j504/2005-03-30/13:35|
> +-+---------------------+
>
> Is this a J version or machine difference?
>
> Devon
>
Looks like a machine problem. I can replicate your results with both J504
and J601 on Windows XP, and get the correct result with either version of
J on Linux32. I don't know enough about how the LAPACK library interface
works to speculate on the cause.
Best wishes,
John
> On 3/1/06, John Randall <[EMAIL PROTECTED]> wrote:
>>
>> Devon McCormick wrote:
>>
>> > ]rod=. >(1,(%3),%9);(3,1,%3);9 3 1 NB. "Replication of detail", p.
>> 7
>> > 1 0.33333333 0.11111111
>> > 3 1 0.33333333
>> > 9 3 1
>> > (%+/){."1 >2{dgeev_jlapack_ rod
>> > _0.2 0.3 0.9
>> > NB. This differs... (vs. 0.077 0.231 0.692)
>>
>> Devon:
>>
>> I get
>>
>> (%+/){."1 >2{dgeev_jlapack_ rod
>> 0.0769231 0.230769 0.692308
>>
>> which looks OK.
>>
>> Best wishes,
>>
>> John
>>
> ...
> --
> Devon McCormick
> ^me^ at acm.
> org is my
> preferred e-mail
> ----------------------------------------------------------------------
> For information about J forums see http://www.jsoftware.com/forums.htm
>
----------------------------------------------------------------------
For information about J forums see http://www.jsoftware.com/forums.htm
