Jak dokładnie działa * pełny * algorytm wielosiatkowy?

Więc rozumiem (a przynajmniej tak mi się wydaje), jak przebiega cykl V. Napisałem w Matlabie 1-D, rekurencyjną wersję cyklu V. Jednak kiedy uruchomiłem kod dla FMG, moje rozwiązanie nie było zbieżne. Wierzę, że mój problem polega na zrozumieniu faktycznej części FMG. To, co obecnie wiem, to:

1. Tuż przed FMG interpolacji Mam rozluźnił moje rozwiązanie $u$
2. Interpolate zarówno błąd i (?)$u$
3. Wykonaj cykl V z 2 siatkami, przekazując błąd do cyklu V (?)
4. Zrelaksuj błąd (na 2. najgrubszej siatce)
5. Interpoluj i błąd$u$
6. Zaktualizuj , dodając do niego błąd.$u$
7. Uruchom cykl V, a następnie powtórz od kroku 4.

Nie jestem pewien co do kolejności, ale mogę się również mylić, co dokładnie interpoluję i przechodzę do mojego cyklu V. Jeśli czegoś brakuje w algorytmie, daj mi znać.

Skąd pomysł na interpolację „błędu”? (A jak mierzysz błąd?)

$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)$