Pertanyaan yang diberi tag «lapack»

LAPACK (Linear Algebra PACKage) adalah pustaka subrutin yang umum digunakan untuk tugas-tugas aljabar linear numerik, termasuk solusi dari himpunan persamaan linier, kuadrat terkecil linier, masalah nilai eigen, dan dekomposisi nilai singular. Rutinitas LAPACK dapat digunakan dengan fortran, C dan kerabat dan berbagai bahasa lainnya.


3
Mengapa Penskalaan Multiplikasi Matriks-Vektor saya tidak?
Maaf untuk posting lama tapi saya ingin memasukkan semua yang saya pikir relevan pada langkah pertama. Apa yang saya inginkan Saya menerapkan versi paralel dari Metode Subruang Krylov untuk Matriks Padat.Terutama GMRES, QMR dan CG. Saya menyadari (setelah membuat profil) bahwa rutinitas DGEMV saya menyedihkan. Jadi saya memutuskan untuk berkonsentrasi …

4
Cepat menentukan apakah matriks padat peringkat rendah atau tidak
Dalam proyek perangkat lunak yang saya kerjakan, perhitungan tertentu jauh lebih mudah untuk matriks peringkat rendah yang padat. Beberapa contoh masalah melibatkan matriks peringkat rendah yang padat, tetapi mereka diberikan kepada saya secara penuh, bukan sebagai faktor, jadi saya harus memeriksa peringkat dan faktor matriks jika saya ingin mengambil keuntungan …


3
Memahami bagaimana Numpy melakukan SVD
Saya telah menggunakan metode yang berbeda untuk menghitung peringkat matriks dan solusi sistem persamaan matriks. Saya menemukan fungsi linalg.svd. Membandingkan ini dengan upaya saya sendiri untuk memecahkan sistem dengan Gaussian Elimination, tampaknya lebih cepat dan lebih tepat. Saya mencoba memahami bagaimana ini mungkin. Sejauh yang saya tahu, fungsi linalg.svd menggunakan …

2
selesaikan
Saya porting kode yang ada dari MATLAB ke C ++ dan memiliki sistem linier untuk menyelesaikan (daripada bentuk yang lebih khas A x = b )xA=bxA=bxA=bAx=bAx=bAx=b Matriks padat, dan berbentuk umum, tetapi tidak lebih besar dari 1000x1000. Jadi dalam MATLAB, solusinya ditemukan oleh fungsi atau notasi forward-slashAAAmrdivide(b,A)x = b/A; Bagaimana …

2
Apa fungsi LAPACK yang sesuai di belakang Matlab [Q, R, E] = qr (A)?
Saat ini saya mencoba untuk murah menghitung perkiraan yang baik peringkat untuk matriks . Oleh karena itu saya menghitung penggunaan QR decompostion columnt pivotingSEBUAHAA [Q,R,E]=qr(A) di Matlab. Saya memperkirakan peringkat menggunakanSEBUAHAA tol = size(A,n)*eps*norm(A,'fro'); r = sum(abs(diag(R))>tol) Ini berfungsi dengan baik dan plot pada semua entri diagonal R terlihat seperti: …
12 matlab  lapack 


4
Matriks eksponensial dari matriks miring-Hermitian dengan fortran 95 dan LAPACK
Saya baru saja dimasukkan ke fortran 95 untuk beberapa simulasi mekanika kuantum. Jujur, saya telah dimanjakan oleh Oktaf jadi saya telah mengambil matriks eksponensial begitu saja. Dengan kemiringan (kecil, ) - Matriks Jerman dengan ukuran , apa cara paling efisien menggunakan LAPACK untuk menyelesaikan masalah ini? Saya tidak menggunakan pembungkus …
11 fortran  lapack 



3
Matriks eksponensial dari matriks asimetris nyata dengan Fortran 95 dan LAPACK
Baru-baru ini saya mengajukan pertanyaan di sepanjang baris yang sama untuk matriks miring-Hermitian. Terinspirasi oleh keberhasilan pertanyaan itu, dan setelah membenturkan kepalaku ke dinding selama beberapa jam, aku melihat matriks eksponensial dari matriks asimetris yang nyata. Rute untuk menemukan nilai eigen dan vektor eigen tampaknya agak berbelit-belit, dan saya khawatir …

3
Bagaimana cara mulai menggunakan LAPACK di c ++?
Saya baru mengenal ilmu komputasi dan saya sudah mempelajari metode dasar untuk integrasi, interpolasi, metode seperti RK4, Numerov dll pada c ++ tetapi baru-baru ini profesor saya meminta saya untuk belajar bagaimana menggunakan LAPACK untuk menyelesaikan masalah yang berkaitan dengan matriks. Seperti misalnya menemukan nilai eigen dari matriks kompleks. Saya …
10 lapack 

3
Apa alasan bahwa LAPACK menggunakan
QR rutin LAPACK menyimpan Q sebagai reflektor rumah tangga. Ini skala vektor refleksi dengan , jadi elemen pertama hasilnya menjadi , sehingga tidak harus disimpan. Dan itu menyimpan vektor terpisah , yang berisi faktor skala yang diperlukan. Jadi matriks reflektor adalah seperti ini:vvv1/v11/v11/v_1111ττ\tauH= Saya- τv vT,H=I−τvvT,H=I-\tau v v^T, di mana …

1
Algoritma Penyeimbang Matriks
Saya telah menulis toolbox sistem kontrol dari awal dan murni di Python3 (plug shameless:) harold. Dari penelitian saya sebelumnya, saya selalu mengeluh tentang pemecah Riccati care.mkarena alasan teknis / tidak relevan. Karenanya, saya telah menulis serangkaian rutinitas saya sendiri. Satu hal yang saya tidak dapat menemukan jalan keluar adalah untuk …

Dengan menggunakan situs kami, Anda mengakui telah membaca dan memahami Kebijakan Cookie dan Kebijakan Privasi kami.
Licensed under cc by-sa 3.0 with attribution required.