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 ekstrim, Anda bisa menjalankan penyelesaian batin untuk konvergensi dalam setiap langkah pemecahan luar.
Di sisi lain, Anda tidak bisa menyelesaikan batin sama sekali, tetapi menggantinya dengan prekondisi batin.
Di suatu tempat di tengah, Anda dapat memotong loop Krylov bagian dalam setelah beberapa iterasi tetap, atau setelah toleransi tertentu tercapai.
Secara empiris, saya telah menemukan situasi di mana ekstrim pertama lebih baik, dan situasi berbeda di mana ekstrim kedua lebih baik (dalam hal total biaya). Namun, saya tidak dapat menemukan alasan yang jelas mengapa situasi tertentu mendukung satu strategi daripada yang lain.
Apakah ada petunjuk atau teori tentang kapan strategi yang berbeda ini lebih disukai?
sumber
Jawaban:
Pertanyaan ini telah terbuka untuk waktu yang lama, tetapi saya pikir masih layak untuk dijawab.
Ini berbeda dengan banyak metode lain yang digunakan untuk prasyarat: misalnya, satu langkah SSOR adalah operasi linear pada vektor yang Anda terapkan, seperti juga semua metode lain yang menerapkan satu langkah dari iterasi titik tetap.
Masalah mendasar sekarang adalah bahwa sebagian besar metode ruang Krylov memang mengharuskan pendahulunya adalah operator linier. Mereka tidak akan konvergen jika preconditioner tidak linier, menjelaskan pengamatan Anda. Di sisi lain, ada variasi beberapa metode ruang Krylov - biasanya diawali dengan kata "Fleksibel", seperti F-GMRES di "Fleksibel GMRES" - yang mengatasi hal ini dan yang dapat menangani prekondisi yang tidak linier operator. Varian fleksibel dari metode asli ini masih akan menyatu, dan seringkali merupakan metode yang kuat ketika digabungkan dengan prekondisi yang baik (tetapi nonlinier).
sumber