Mengacu pada Subruang Krylov dan metode penyelesaian sistem persamaan linier yang memanfaatkan ruang-ruang ini.
Apa yang bisa salah ketika menggunakan metode Krylov yang dikondisikan sebelumnya dari KSP ( paket solver linier PETSc ) untuk menyelesaikan sistem linier yang jarang seperti yang diperoleh dengan mengkritisi dan membuat linierisasi persamaan diferensial parsial? Langkah apa yang bisa saya ambil...
Seperti yang saya pahami, ada dua kategori utama metode iteratif untuk menyelesaikan sistem persamaan linear: Metode Alat Tulis (Jacobi, Gauss-Seidel, SOR, Multigrid) Metode Subruang Krylov (Konjugat Gradien, GMRES, dll.) Saya mengerti bahwa sebagian besar metode stasioner bekerja dengan...
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...
Sejauh yang saya ketahui, pemecah multigrid menggunakan peranti iteratif seperti Jacobi, Gauss-Seidel, dan SOR untuk meredam kesalahan pada berbagai frekuensi. Bisakah metode ruang bagian Krylov (seperti gradien konjugat, GMRES, dll.) Digunakan sebagai gantinya? Saya tidak berpikir mereka...
Saya ingin tahu apa yang terjadi pada prekondisi polinomial. Saya tertarik pada mereka, karena mereka tampak relatif elegan dari perspektif matematika, tetapi sejauh yang saya baca dalam survei tentang metode krylov, mereka umumnya berkinerja sangat buruk sebagai prekondisi. Dalam kata-kata Saad...
Pertanyaan: Misalkan Anda memiliki dua prekondisi yang berbeda (diperhitungkan) untuk matriks pasti positif simetris : dan mana invers faktor faktor adalah mudah diaplikasikan.AAAA≈BTBA≈BTBA \approx B^TBA ≈ CTC,SEBUAH≈CTC,A \approx C^TC,B , BT, C, CTB,BT,C,CTB, B^T, C, C^T Ketika apakah...
Katakanlah saya memiliki sistem linier Ax=bAx=bA x = b , yang konvergen dengan cepat menggunakan metode Krylov yang sesuai (seperti CG atau GMRES) untuk semua bbb . Jika BBB adalah matriks dengan peringkat rendah rrr, akankah metode Krylov yang sama pada sistem (A+B)x=b(A+B)x=b(A + B) x = b juga...
Persamaan Poisson dengan semua kondisi batas Neumann memiliki ruang nol dimensi konstan tunggal. Ketika menyelesaikan melalui metode Krylov, ruang nol dapat dihapus baik dengan mengurangi rata-rata solusi setiap iterasi atau dengan menyematkan nilai satu titik. Menyematkan satu titik memiliki...
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...
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...
Sistem saya adalah masalah FE simetris dengan pengganda lagrange (misalnya aliran Stokes yang tidak dapat dikompresi): ( ABBTC)(SEBUAHBTBC)\begin{pmatrix}A & B^T \\ B & C\end{pmatrix} di mana adalah kasus khas (saya bahkan telah memastikan bahwa persamaan diberi nomor sehingga pengali Lagrange...
Dalam ilmu komputasi kita sering menemukan sistem linear besar yang harus kita pecahkan dengan beberapa cara (efisien), misalnya dengan metode langsung atau berulang. Jika kita fokus pada yang terakhir, bagaimana kita dapat menetapkan bahwa metode iteratif untuk menyelesaikan sistem linier yang...
Saya memecahkan untuk matriks A definitif positif jarang besar menggunakan metode konjugat gradien (CG). Dimungkinkan untuk menghitung penentu A menggunakan informasi yang dihasilkan selama penyelesaian?A x =
Saya tertarik dalam menghitung solusi sistem lage ODE menggunakan metode krylov seperti pada [1]. Metode semacam itu melibatkan fungsi-fungsi yang berkaitan dengan fungsi eksponensial (yang disebut fungsi ). Ini pada dasarnya terdiri dari menghitung tindakan fungsi matriks dengan membangun ruang...
Pertimbangkan situasi di mana Anda ingin menyelesaikan sistem linier menggunakan metode Krylov yang telah dikondisikan sebelumnya, tetapi menerapkan prekondisi itu sendiri melibatkan penyelesaian sistem bantu, yang dilakukan dengan metode Krylov lain yang telah dikondisikan sebelumnya. Pada satu...