Skip to content

eig to return eigenvalues in order of descending absolute value #324

Closed
@mhauru

Description

Eigenvalue decomposition of a positive-definite Hermitian matrix is a special case of SVD. However, svd returns the singular values in descending order, where as eig returns them in ascending order. The order of singular values is fixed popular convention, so I think for consistency eig should adhere to the same order.

Further more, I would prefer the ordering for eig to be descending absolute value. This would also make it consistent with the default behavior of eigs and with the fact that singular values of Hermitian matrices are the absolute values of the eigenvalues.

I don't know how this would affect perfomance, i.e. can LAPACK do this for us naturally or would we need to sort (and permute the eigenvectors) afterwards.

Activity

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Metadata

Assignees

No one assigned

    Labels

    docsThis change adds or pertains to documentation

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions