Numerische Bibliothek in C++
-
Hallo headmyshoulder!
Ich habe erst gestern das entsprechende Kapitel in den Num. Recipes gelesen und bin heute auf deine Idee von dem C++ Odeint Framework gestoßen. Ich hoffe das Projekt läuft noch!
Ich selber habe Mathematik studiert. Mit der ODE-Solvern und Schrittweitensteurungen beschäftige mich erst seit Kurzem. Auch in C++ bin ich noch relativ unerfahren. Aber ich werde muss mich sowieso auf beiden Gebieten in nächster Zeit einarbeiten. Über jeden Austausch und gegebenenfalls eine Zusammenarbeit würde ich mich freuen!
Viele Grüße
FelixKlein
-
Hi FelixKlein,
das odeint Projekt läuft noch und ist auch schon recht fortgeschritten, du hast quasi eine menge aufzuholen. Zu finden ist das Ganze in der boost sandbox
https://www.boost.org/svn/boost/sandbox/odeint
Akutell könnten wir die Dormand-Prince Methoden gebrauchen. Wenn Du Bock kannst Du dir ja ansehen wie unsere lib funktioneirt und dann die Methoden implementieren.
Viele Grüße,
headmyshoulder
-
Also DOPRI interessiert mich persönlich weniger.
Habt ihr euch eigentlich schon Gedanken über Extrapolationsverfahren gemacht?In jedem Fall muss ich erst einmal noch einiges lesen. Erst wenn ich die theoretischen Grundlagen durchblicke will ich mir eure Implementierungen im Detail ansehen und erst wenn ich die verstehe mache ich mich selber ans Werk.
Soweit ich das bisher abschätzen kann ist euer Projekt sehr strukturiert und sauber programmiert. Ob ich mit dem Code tatsächlich etwas anfangen kann weiß ich noch nicht. Toll jedenfalls, dass ihr das als Open Source ins Netz stellt.
Ich würde gerne weiter in Kontakt bleiben. Gerne auch persönlicher als nur über das Forum (E-Mail o.Ä.)
Gruß
FelixKlein
-
FelixKlein schrieb:
Also DOPRI interessiert mich persönlich weniger.
Ok, kein Problem. Es gibt noch ne Menge anderer Sachen...
FelixKlein schrieb:
Habt ihr euch eigentlich schon Gedanken über Extrapolationsverfahren gemacht?
Ja, haben wir, und das steht auch auf unserer Wunschliste. Allerdings gibt es nicht für jede Methode ein Extrapolationsverfahren, glaube ich zumindest. Deshalb haben wir das erstmal nicht mit aufgenommen. Implizite Methoden wäre prinzipiell auch sehr interessant. Allerdings braucht man dann lineare Algebra (LU-Zerlegung,...) was nicht trivial ist. Und symplektische Verfahren mit Fehlerabschätzung haben wir auch noch nicht implementiert.
FelixKlein schrieb:
In jedem Fall muss ich erst einmal noch einiges lesen. Erst wenn ich die theoretischen Grundlagen durchblicke will ich mir eure Implementierungen im Detail ansehen und erst wenn ich die verstehe mache ich mich selber ans Werk.
Welche Bücher hast Du dazu? Die von Hairer kann ich empfehlen.
FelixKlein schrieb:
Soweit ich das bisher abschätzen kann ist euer Projekt sehr strukturiert und sauber programmiert. Ob ich mit dem Code tatsächlich etwas anfangen kann weiß ich noch nicht. Toll jedenfalls, dass ihr das als Open Source ins Netz stellt.
Wir sind im Moment eigentlich zu zweit. Und das Ziel ist es, das Ganze in die boost bibliotheken zu submitten. Zu den technischen Details haben wir einen Artikel geschrieben:
http://www.codeproject.com/KB/recipes/odeint.aspx, obwohl der schon leicht veraltet ist. Ich werde den aber demnächst aktualisieren.FelixKlein schrieb:
Ich würde gerne weiter in Kontakt bleiben. Gerne auch persönlicher als nur über das Forum (E-Mail o.Ä.)
Jo, ich schick Dir mal eine pn. Wir haben auch eine Mailingliste zu diesem Thema bei sf: https://lists.sourceforge.net/lists/listinfo/yanl-develop. Wenn Du willst können wir auch darüber kommunizieren.
Viele Grüße,
headmyshoulder
-
headmyshoulder schrieb:
Welche Bücher hast Du dazu? Die von Hairer kann ich empfehlen.
Momentan arbeite ich mich durch "Numerische Mathematik (2)" von Bulirsch/Stoer. Ich habe aber auch schon andere Lektüre gesichtet. Unter anderem "Solving ODEs I/II" von Hairer/Norsett/Wanner. Meintest du das oder andere Bücher von ihm?
headmyshoulder schrieb:
Zu den technischen Details haben wir einen Artikel geschrieben:
http://www.codeproject.com/KB/recipes/odeint.aspx, obwohl der schon leicht veraltet ist. Ich werde den aber demnächst aktualisieren.Den Artikel hatte ich schon entdeckt. Ist interessant. Ich freue mich auf die Aktualisierung.
headmyshoulder schrieb:
Wir haben auch eine Mailingliste zu diesem Thema bei sf: https://lists.sourceforge.net/lists/listinfo/yanl-develop. Wenn Du willst können wir auch darüber kommunizieren.
Was genau ist der Zweck der Mailingliste? Ich möchte weder vollgespammt werden noch, dass meine E-Mails öffentlich gemacht werden.
Ich melde mich heute Abend wieder per E-Mail.
Viele Grüße
FelixKlein
-
FelixKlein schrieb:
Momentan arbeite ich mich durch "Numerische Mathematik (2)" von Bulirsch/Stoer. Ich habe aber auch schon andere Lektüre gesichtet. Unter anderem "Solving ODEs I/II" von Hairer/Norsett/Wanner. Meintest du das oder andere Bücher von ihm?
Jepp, die Solving ODEs I/II" mein ich. Für symplektische Verfahren haben wir auch mit "B. Leimkuhler and S. Reich, Simulating Hamiltonian Dynamics." gearbeitet.
FelixKlein schrieb:
Was genau ist der Zweck der Mailingliste? Ich möchte weder vollgespammt werden noch, dass meine E-Mails öffentlich gemacht werden.
Soviel passiert im Moment nicht auf der Liste, also keine Angst vor Spam. Man kann da halt diskutieren, alles was odeint betrifft. Ist wichtig, falls mehrere Leute programmieren, daß alle halbwegs auf dem Laufenden sind. Einige Sachen sollte man bei solch einem Projekt öffentlich ankündigen.
Wenns um Details geht, kann man natürlich auch direkt via Email kommunizieren.
Bis dahin,
headmyshoulder
-
Nur mal als kleines Update: Was hier als Bibliothek zum Lösen von Differentialgleichungen (odeint) angefangen hat ist jetzt akzeptiert als Boost-Bibliothek. Hat ziemlich lange gedauert, und es war eine Menge Arbeit, aber es hat sich gelohnt. Und natürlich ist der Umfang enorm gewachsen. Ein Hauptvorteil ist die Unabhängigkeit der Algorithmen von arithmetischen Operationen und Containern, so dass man mit relativ einfachem Aufwand auch Differentialgleichungen auf Grafikkarten lösen kann.
-
ihr freaks
ich gratuliere, schoen mal ein erfolgreiches projekt zu sehen!
-
-
Ich habe die Entwicklung aus der ferne verfolgt. Habt ja auch einige Vortraege gehalten. Sieht gut aus ... wenn ich nur mehr Zeit haette.
-
Wo finde ich die Bibliothek?
-
Kellerautomat schrieb:
Wo finde ich die Bibliothek?
Die Bibliothek liegt unter odeint.com. Es wird wahrscheinlich bis Release 1.53 dauern bis die direkt in der Boost ist.