Re: [Pdl-devel] meigen sort criteria for eigenvalue and eigenvectors.

2022-07-05 Thread Ed .
Thank you Luis! Guillermo, if you believe it would help to document this a bit better (i.e. more explicitly regarding input and output vectors/matrices) in PDL::LinearAlgebra, please open an issue there and spell out your thoughts. Best regards, Ed From: Guillermo P. Ortiz

Re: [Pdl-devel] meigen sort criteria for eigenvalue and eigenvectors.

2022-07-05 Thread Guillermo P. Ortiz
Yes, Dear Luis, your are fast and right ! just confirm that using meigen() to obtain eigenvectors (right) if you perform matrix multiplication (x) the vector are in the rows but if you want to see they you need the transpose Thanks ! El mar, 5 jul 2022 a las 16:10, Luis Mochan () escribió: > > De

Re: [Pdl-devel] meigen sort criteria for eigenvalue and eigenvectors.

2022-07-05 Thread Luis Mochan
Dear Guillermo, I guess the eigenvectors are the rows, not the columns, and they are normalized! Regards, Luis On Tue, Jul 05, 2022 at 03:05:16PM -0300, Guillermo P. Ortiz wrote: > Dear PDL developers, > > Trying to use meigen() for eigenvalue problem with matrix 4x4 without > threading for now,

Re: [Pdl-devel] meigen sort criteria for eigenvalue and eigenvectors.

2022-07-05 Thread Ed .
Hi Guillermo, https://metacpan.org/pod/PDL::LinearAlgebra#meigen says it uses https://metacpan.org/pod/PDL::LinearAlgebra::Real#geev. That obviously uses the relevant LAPACK function. The double-precision version, dgeev, is documented at https://www.netlib.org/lapack/explore-html/d9/d8e/group__

[Pdl-devel] meigen sort criteria for eigenvalue and eigenvectors.

2022-07-05 Thread Guillermo P. Ortiz
Dear PDL developers, Trying to use meigen() for eigenvalue problem with matrix 4x4 without threading for now, I realize that eigenvalues seems did not have any ordered output, is not it? In the first example below the second eigenvalue is larger than the first === Matrix: [ [