Saya memiliki masalah di mana saya harus menemukan semua positif (seperti dalam nilai eigen positif) pasang eigen dari matriks nonsimetrik kecil (biasanya lebih kecil dari 60x60). Saya bisa berhenti menghitung ketika nilai eigen lebih kecil dari ambang tertentu. Saya tahu bahwa nilai eigennya nyata. Adakah saran tentang algoritma yang dapat saya gunakan untuk mencoba memeras kinerja terbaik? Saya harus melakukan beberapa ribu dekomposisi ini, jadi kecepatan itu penting.
Terima kasih sebelumnya.
EDIT: Saya perlu melakukan ini pada GPU dalam memori bersama. Matriks juga tidak harus berukuran sama. Saya tidak mengetahui ada perpustakaan yang melakukan ini saat ini. Saran algoritma yang cocok untuk masalah akan dihargai.
performance
eigensystem
gpu
Kantoku
sumber
sumber
Jawaban:
Tanpa melakukan banyak pencarian, saya sarankan Anda untuk melihat perpustakaan MAGMA . Kode tersedia secara bebas dengan dukungan terus menerus. NVIDIA mengakui MAGMA sebagai "Terobosan dalam Solver untuk Masalah Nilai Eigen".
Ada juga perpustakaan CULA , yang umumnya merupakan produk komersial, meskipun baru-baru ini telah dibuat gratis untuk penggunaan akademis (lihat detailnya di sini ).
sumber
Gunakan fungsi di LAPACK, tidak mungkin Anda bisa mengalahkan mereka dalam implementasi Anda sendiri.
sumber