Bagaimana tepatnya algoritma multigrid * penuh * berjalan?

Jadi saya mengerti (atau setidaknya saya percaya saya tahu) bagaimana siklus-V berjalan. Saya telah menulis dalam Matlab versi 1-D, versi rekursif dari siklus-V. Namun, ketika saya menjalankan kode saya untuk FMG, solusi saya tidak konvergen. Saya percaya masalah saya terletak pada pemahaman saya tentang bagian FMG yang sebenarnya. Yang saya tahu saat ini adalah:

1. Hanya sebelum FMG interpolasi, saya sudah santai solusi saya $u$
2. Interpolasi baik kesalahan dan (?)$u$
3. Jalankan v-cycle 2-grid, meneruskan kesalahan ke dalam v-cycle (?)
4. Santai kesalahan (di grid 2 kasar)
5. Interpolasi dan kesalahannya$u$
6. Perbarui dengan menambahkan kesalahan padanya.$u$
7. Jalankan siklus-v, lalu ulangi dari langkah 4.

Saya tidak yakin tentang pesanan, tetapi saya juga bisa salah tentang apa yang sebenarnya saya interpolasi dan lulus ke dalam v-cycle saya. Jika saya kehilangan sesuatu dari algoritme, beri tahu saya.

Di mana Anda menemukan interpolasi "kesalahan"? (Dan bagaimana Anda mengukur kesalahan?)

$u$$u^h \gets \mathbb{I}_H^h u^H$$\mathbb{I}_h^H = I_h^H$

$r^h = A^h u^h - b^h$

$\tilde u^H \gets \hat I_h^H \tilde u^h$$A$

$b^H$$\tau_h^H$$u^{h*}$$A^H \hat I_h^H u^{h*} = b^H + \tau_h^{H*}$$u^h \gets \tilde u^h + I_H^h (u^H - \hat I_h^H \tilde u^h)$