Beispiel - SAM - ETH Zürich

–y0: Anfangszustandy 0 ,
–h0: Schrittweiteh 0 für den ersten Zeitschritt
Numerische
Mathemtik
–reltol,abstol: Relative and absolute Toleranzen, siehe oben,
–hmin: minimale Zeitschrittweite, Verfahren bricht ab, fallsh k < h min , was für das Erkennen
von Blow-ups und Kollaps wichtig ist.
• Zeile 3: Überprüfe, ob der Endzeitpunkt erreicht ist oder das Verfahren steckengeblieben ist
(h k < h min ).
• Zeile 4, 5: Propagiere den aktuellen Zustand mit Hilfe beider Einschrittverfahren.
• Zeile 6: Berechne die Norm des geschätzten Fehlers, siehe (2.6.3).
• Zeile 8: Vergleich, um zu entscheiden, ob der aktuelle Schritt akzeptiert oder verworfen werden
sollte.
• Zeile 9, 10: Schritt akzeptiert: Aktualisiere den Zustand und schlage 1.1 mal die aktuelle
Schrittweite für den nächsten Schritt vor.
R. Hiptmair
rev 35327,
25. April
2011
• Zeile 11 Schritt verworfen: Versuche es nochmal mit der halben Zeitschrittweite.
• Rückgabewerte
–t: Zeitgittert 0 < t 1 < t 2 < ... < t N < T , wobeit N < T auf vorzeitigen Abbruch hinweist
(Kollaps, Blow-up),
–y: Folge von Zuständen(y k ) N k=0 .
2.6
p. 297