Re: [Yanl-develop] Contribution

[Karsten A. <kar...@gm...>]

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
>