Thank you for the clear report.
The documentation is correct, but there was a bug in do_dgeev in that it forgot
to transpose vl and vr matrix because J is row major.
A fix has been pushed to jal. Please download and try again.
A script to reproduce the issue
load '~addons/math/lapack2/example/dgeev.ijs'
echo t0=: 0 0 0 1 1, 0 0 1 0 1, 0 1 0 0 1, 1 0 0 0 1,:1 0 0 1 0
echo 0 2 { do_dgeev t0
echo <.@(+&0.5)&.(100&*) >2 { do_dgeev t0
> On 5 Feb 2022, at 6:48 PM, R.E. Boss <[email protected]> wrote:
>
> Thanks for the assistance offered, but I think the chance is epsilon I will
> suggest a change in the future, and I saw you as an author (in the github
> already), plus the suggested change is minimal, so I think it's most
> efficient if you would be so kind to submit this pull request.
> Thnx in advance.
>
>
> R.E. Boss
>
>
> -----Original Message-----
> From: Programming <[email protected]> On Behalf Of
> Raul Miller
> Sent: vrijdag 4 februari 2022 15:43
> To: Programming forum <[email protected]>
> Subject: Re: [Jprogramming] j903/addons/math/lapack2/example/dgeev.ijs
>
> I would recommend a small pull request against
> https://github.com/jsoftware/math_lapack2
>
> (If you are uncomfortable with the mechanics of pull requests, let me know
> and I can walk you through it. Possibly in the chat forum, possibly on a wiki
> page, possibly in private email, etc..)
>
> Thanks,
>
> --
> Raul
>
> On Fri, Feb 4, 2022 at 4:41 AM R.E. Boss <[email protected]> wrote:
>>
>> This script says:
>>
>> NB. [out] VR If the j-th eigenvalue is real, then v(j) =
>> VR(:,j),
>> NB. the j-th column of VR.
>>
>> but it appears to be "the j-th ROW of VR"
>>
>>
>> R.E. Boss
>>
>> ----------------------------------------------------------------------
>> For information about J forums see http://www.jsoftware.com/forums.htm
> ----------------------------------------------------------------------
> For information about J forums see http://www.jsoftware.com/forums.htm
> ----------------------------------------------------------------------
> For information about J forums see http://www.jsoftware.com/forums.htm
----------------------------------------------------------------------
For information about J forums see http://www.jsoftware.com/forums.htm
