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
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.
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
On Fri, 1 Jul 2022, Luis Mochan wrote:
> # read all data, skip first row
> my ($f, $M11, $A11, $M21, $A21, $M12, $A12, $M22, $A22)=rcols "rem.s2p",
> {LINES=>"1:-1:1"};
> my $rows=$f->nelem;
> my $M=pdl($M11, $M12, $M21, $M22); # Indices: filerow, position
> my $R=pdl($A11, $A12, $A21,