Mehrgittermethoden

Seminar

Robert Speck & Dieter Moser, Sommersemester 2016

</center>


Mehrgitter für nicht-lineare Funktionen

Implementieren Sie sowohl Newton-Mehrgitter und FAS-Mehrgitter in Python. Sie können dafür die Struktur von pyMG benutzen, müssen dies aber nicht tun. Dokumentieren Sie den Code angemessen und schreiben Sie geeignete Tests für nose. Teilen Sie die Arbeit geeignet auf: Ziel dieser erweiterten Übungsaufgabe ist es auch, gemeinschaftliche Software-Entwicklung zu lernen.

Vergleichen Sie beide Ansätze für folgendes Modell-Problem:

$$-\Delta u(x,y) + \gamma u(x,y) e^{u(x,y)} = f(x,y),\quad (x,y)\in[0,1]^2,$$$$u(x,y) = 0,\quad (x,y)\in\partial[0,1]^2$$

Folgen Sie bei den numerischen Tests von Kapitel 5 aus "MULTIGRID METHODS FOR NONLINEAR PROBLEMS: AN OVERVIEW", siehe hier. Wählen Sie insbes. den Quellterm $f$ gemäß Gleichung (5.4).

Erstellen Sie einen Release und reichen Sie den Link zum Release ein. Bereiten Sie außerdem eine ca. 30-minütige Vorstellung Ihres Programms vor, bei dem Sie uns Ihre Implementierung und die Ergebnisse der numerischen Tests erläutern.