Saya harus menyelesaikan masalah nilai eigen umum mana dan keduanya tridiagonal, adalah simetris positif pasti dan nyata, tetapi hanya simetris yang kompleks (tidak pasti atau Hermitian). Selain itu, saya perlu komposisi eigend lengkap. Saat ini saya hanya memanggil eigensolver umum dari Lapack, tetapi saya bertanya-tanya apakah ada metode yang lebih baik untuk masalah khusus yang sangat terstruktur ini. Secara khusus, memiliki kode yang tersedia secara bebas (C ++) akan menjadi yang terbaik.ZGGEV
linear-algebra
eigensystem
Victor Liu
sumber
sumber
Jawaban:
Metode Pole EXpansion dan Selected Inversion ( PEXSI ) mungkin jawabannya. Saya belum pernah menggunakan metode ini, tetapi ia menawarkan inversi rutin untuk matriks simetris yang kompleks. Ini tidak spesifik untuk matriks tridiagonal, tetapi memanfaatkan sparsity.
sumber