Re: [Yanl-develop] Contribution
Brought to you by:
karstenahnert
|
From: Karsten A. <kar...@gm...> - 2009-10-30 14:23:25
|
Hi Constantine, Thanks for your reply. Maybe, first of all it would be nice, if you could contribute some examples to odeint (similar to lorenz.cpp in libs/numeric/odeint/examples): 1. Numerical integration of the harmonic oscillator [1] and comparison with the analytical solution. The ODE is d^2x/dt^2 = - x - gamma dx / dt, and the analyitcal solution has the form x(t) = A cos( omega t + phi ) E^( - kappa t ) with appropriate constants A, omega, phi, kappa. It would be great here, if it could be shown that the precision of Euler is much smaller than that of RK4. 2. Speed test for the same ODE, with different containers. For example, it would be nice to see, which Container in lorenz.cpp is the fastest: vector<double>, list<double>, array<double,3>? If you need any help or have any questions, just ask. Maybe, you have also some suggestion on the interface of the odeint classes. Thank you a lot, Karsten [1] http://en.wikipedia.org/wiki/Harmonic_oscillator , the universal oscillator and the damped oscillator section. Constantine wrote: > Message body follows: > > > Karsten, > > Thanks for the reply. I have checked out the code from > https://boost.org/svn/boost/sandbox/odeint. As far as ODEs > and TMP, I > am comfortable with both. I would be happy to jump into > either or both > on this project. > > Constantine > > -- > This message has been sent to you, a registered SourceForge.net user, > by another site user, through the SourceForge.net site. This message > has been delivered to your SourceForge.net mail alias. You may reply > to this message using the "Reply" feature of your email client, or > using the messaging facility of SourceForge.net at: > https://sourceforge.net/sendmessage.php?touser=2670242 > |