Untuk pertanyaan tentang desain dan implementasi prekondisi untuk menyelesaikan sistem linier.
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...
Untuk solusi sistem linear besar menggunakan metode iteratif, sering kali menarik untuk memperkenalkan prakondisi, misalnya menyelesaikan , di mana sini digunakan untuk prakondisi kiri sistem. . Biasanya, kita harus memiliki M ^ {- 1} \ approx A ^ {- 1} dan memberikan dasar untuk (lebih banyak)...
Ini sebagian besar ditujukan untuk PDE elips lebih dari domain cembung, sehingga saya bisa mendapatkan gambaran yang baik dari dua
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...
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...
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...
Saya ingin menyelesaikan masalah non-linear dengan kendala kesetaraan non-linear dan saya menggunakan Lagrangian yang diperbesar dengan istilah regularisasi penalti yang, sebagaimana diketahui, merusak jumlah kondisi sistem linierisasi saya (pada setiap iterasi Newton yang saya maksud) . Semakin...
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...
Metode Jacobian-Free Newton-Krylov (JFNK), dan metode Krylov secara umum, dapat sangat berguna karena mereka tidak memerlukan penyimpanan eksplisit atau konstruksi matriks, hanya hasil dari produk vektor-matriks. Jika Anda benar-benar membentuk sistem sparse, ada banyak prasyarat di luar sana untuk...
Saya memiliki sistem linier dengan matriks yang nilai eigennya terdistribusi secara seragam pada satuan lingkaran seperti ini: Apakah mungkin untuk menyelesaikan sistem semacam ini secara efektif dengan metode iteratif, mungkin dengan beberapa
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...