Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Complex tridiagonal skewhermitian eigensolver #56

Closed
smataigne opened this issue Aug 16, 2022 · 2 comments
Closed

Complex tridiagonal skewhermitian eigensolver #56

smataigne opened this issue Aug 16, 2022 · 2 comments

Comments

@smataigne
Copy link
Collaborator

Find a method to provide an eigensolver for SkewHermtridiagonal complex matrices

@smataigne
Copy link
Collaborator Author

smataigne commented Aug 16, 2022

A SkewHermTridiagonal matrix can be transformed into a complex SymTridiagonal by unitary transformation but LAPACK provides no solver for complex symmetric tridiagonal matrices. Maybe implementing my own QR-algorithm would be worth since the matrix is already tridiagonal. It could certainly outperform hermitian eigensolvers.

@smataigne
Copy link
Collaborator Author

smataigne commented Aug 17, 2022

SkewHermTridiagonal can be multiplied by $i$ to be turned into a Hermitian matrix. It can then be turned into real SymTridiagonal using "scalar " givens rotations. Then efficient solvers can be used. The orthogonal transformation will be diagonal, so very cheap to store and compute.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant