Menu

#3 EigenVect: complex eigenvalues issue

1.0
open
None
2026-01-01
2026-01-01
No

Hello Glorfin,

first of all, thank you for this very useful library, which is a real gem in the Free Pascal world.
I think there is an issue in the EigenVect procedure:

unit ueigvec;
...
procedure EigenVect(A      : TMatrix;
                    Lb, Ub : Integer;
                    Lambda : TCompVector;
                    V      : TMatrix);
...

If the eigenvalues (Lambda) are complex numbers, the eigenvectors (V) should
be complex vectors, not real vectors as shown in the definition.

Best regards
Melchiorre

Discussion


Log in to post a comment.