Pertanyaan yang diberi tag «linear-solver»

Mengacu pada metode untuk memecahkan sistem persamaan linear.

Apakah ada pemecah pemrograman nonlinier berkualitas tinggi untuk Python?
Saya memiliki beberapa masalah optimisasi global non-cembung yang menantang untuk dipecahkan. Saat ini saya menggunakan MATLAB's Optimization Toolbox (khusus, fmincon()dengan algoritma = 'sqp'), yang cukup efektif . Namun, sebagian besar kode saya menggunakan Python, dan saya ingin melakukan optimasi dengan Python juga. Apakah ada pemecah NLP dengan binding Python yang …

Panduan apa yang harus saya ikuti ketika memilih pemecah sistem linier yang jarang?
Sistem linear yang jarang muncul dengan meningkatnya frekuensi dalam aplikasi. Seseorang memiliki banyak rutinitas untuk dipilih untuk menyelesaikan sistem ini. Pada level tertinggi, ada batas antara metode langsung (mis. Eliminasi Gaussian yang jarang atau dekomposisi Cholesky, dengan algoritma pemesanan khusus, dan metode multirontal) dan iteratif (misalnya metode GMRES, (bi-) konjugat …

Cara memilih metode untuk menyelesaikan persamaan linear
Sepengetahuan saya, ada 4 cara untuk memecahkan sistem persamaan linear (koreksi saya jika ada lebih banyak): Jika matriks sistem adalah matriks kuadrat peringkat penuh, Anda dapat menggunakan Aturan Cramer; Hitung invers atau pseudoinverse dari matriks sistem; Gunakan metode dekomposisi matriks (eliminasi Gaussian atau Gauss-Jordan dianggap sebagai dekomposisi LU); Gunakan metode …

Memecahkan tanpa membalikkan
Saya memiliki matriks dan . jarang dan dengan sangat besar (bisa dalam urutan beberapa juta.) adalah matriks tinggi dengan agak kecil ( ) dan setiap kolom dapat hanya memiliki satu entri dengan sisanya menjadi 's, sehingga . sangat besar, jadi sangat sulit untuk membalikkan, dan saya bisa menyelesaikan sistem linear …



Apakah ada implementasi ILU multilevel berbasis sumber terbuka?
Saya sangat terkesan dengan kinerja serial prekondisi ILU berbasis bertingkat terbalik , terutama untuk Helmholtz heterogen , tetapi saya terkejut tidak dapat menemukan implementasi open source. Secara khusus, ILUPACK membuat biner tersedia secara bebas untuk akademisi, tetapi tampaknya mereka tidak merilis kode sumbernya. Apakah ini benar-benar terjadi bahwa tidak ada …


Bagaimana Multigrid yang dipercepat Krylov (menggunakan MG sebagai prasyarat) termotivasi?
Multigrid (MG) dapat digunakan untuk menyelesaikan sistem linier dengan membuat tebakan awal dan mengulangi yang berikut untuk sampai konvergensi:A x = bSEBUAHx=bAx=bx0x0x_0i = 0 , 1 ..saya=0,1 ..i=0,1.. Hitung sisarsaya= b - A xsayarsaya=b-SEBUAHxsayar_i = b-Ax_i Terapkan siklus multigrid untuk mendapatkan aproksimasi , di mana .Δ xsaya≈ esayaΔxsaya≈esaya\Delta x_i \approx …

Adalah algoritma Thomas cara tercepat untuk memecahkan sistem linier tridiagonal simetris dominan dominan diagonal
Saya bertanya-tanya apakah algoritma Thomas adalah cara tercepat (dapat dibuktikan?) Untuk menyelesaikan sistem tridiagonal simetris dominan diagonal dalam hal kompleksitas algoritmik (tidak mencari paket implementasi seperti LAPACK dll). Saya tahu bahwa kedua algoritma Thomas dan multigrid adalah kompleksitas, tetapi mungkin faktor konstan untuk multigrid kurang? Rasanya bagi saya multigrid tidak …

Presisi floating-point tunggal versus ganda
Angka floating point presisi tunggal mengambil setengah memori dan pada mesin modern (bahkan pada GPU tampaknya) operasi dapat dilakukan dengan mereka pada kecepatan hampir dua kali lipat dibandingkan dengan presisi ganda. Banyak kode FDTD yang saya temukan secara eksklusif menggunakan aritmatika dan penyimpanan presisi tunggal. Apakah ada aturan praktis tentang …

mendahului metode krylov dengan metode krylov lainnya
Dalam metode seperti gmres atau bicgstab, mungkin menarik untuk menggunakan metode krylov lain sebagai prasyarat. Bagaimanapun mereka mudah diimplementasikan dalam cara bebas-matriks dan dalam lingkungan paralel. Misalnya, satu dapat menggunakan beberapa (misalkan ~ 5) iterasi bigcstab tanpa syarat sebagai precontioner untuk gmres, atau kombinasi metode krylov lainnya. Saya tidak menemukan …

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 …


Berulang kali memecahkan dengan , berbeda
Saya menggunakan MATLAB untuk memecahkan masalah yang melibatkan penyelesaian di setiap timestep, di mana berubah seiring waktu. Saat ini, saya menyelesaikan ini menggunakan MATLAB :bAx=bAx=b\mathbf{A} \mathbf{x}=\mathbf{b}bb\mathbf{b}mldivide x = A\b Saya memiliki fleksibilitas untuk membuat sebanyak mungkin perhitungan yang diperlukan, jadi saya bertanya-tanya apakah ada metode yang lebih cepat dan / …

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.