Mittwoch, 16. Mai 2012

Schmetterlingseffekt

Differentialgleichungen numerisch zu lösen macht mehr Spaß, als man erwarten würde, wenn man es hört. Und sobald man den ersten Runge-Kutta-Algorithmus in einer kommerziellen Interpretersprache geschrieben hat, bemerkt man, dass dieses Skript doch recht lange braucht.
Für dieses Problem gibt es zwei Lösungen: Entweder wird man zum Guru und wendet irgendeine okkulte Code-Magie an, um das Programm schneller laufen zu lassen, oder man schreibt das Programm in einer schönen Sprache neu. In C zum Beispiel. (Quellen hier)

Ich habe mich für den einfachen Weg entschieden und wenig überraschend eine Tempoverbesserung von Faktor ~140 festgestellt. Jedenfalls für diesen Lorenz-Attraktor. Geplottet habe ich die Werte dann mit Python und matplotlib

Warum ich den Titel "Schmetterlingseffekt" gewählt habe? Naja, das Bild hier sieht ein wenig nach einem Schmetterling aus :)
Davon abgesehen, wurde der Schmetterlingseffekt nach diesen Differentialgleichungen benannt.