349
edits
Line 78: | Line 78: | ||
Details on the methods to be implemented can be found in [http://dx.doi.org/10.1145/502800.502801 this paper] on bvp4c and [http://www.jnaiam.net/new/uploads/files/014dde86eef73328e7ab674d1a32aa9c.pdf this paper] on bvp5c. Further details are available in [http://books.google.it/books/about/Nonlinear_two_point_boundary_value_probl.html?id=s_pQAAAAMAAJ&redir_esc=y this book]. | Details on the methods to be implemented can be found in [http://dx.doi.org/10.1145/502800.502801 this paper] on bvp4c and [http://www.jnaiam.net/new/uploads/files/014dde86eef73328e7ab674d1a32aa9c.pdf this paper] on bvp5c. Further details are available in [http://books.google.it/books/about/Nonlinear_two_point_boundary_value_probl.html?id=s_pQAAAAMAAJ&redir_esc=y this book]. | ||
=== | === Geometric integrators for Hamiltonian Systems === | ||
The | [http://openlibrary.org/books/OL9056139M/Geometric_Numerical_Integration | Geometric (AKA Symplectic) integrators] are useful for | ||
multi-dimensional classical mechanics problems and for molecular dynamics simulations. | |||
The odepkg package has a number of solvers for ODE, DAE and DDE problems but none of them is currently | |||
specifically suited for second order problems in general and Hamiltonian systems in particular. | |||
Therefore a new package for geometric integrators would be a useful contribution. | |||
This could be created as new package or added as a set of new functions for odepkg. | |||
The function interface should be consistent throughout the package and should be modeled to follow | |||
that of other functions in odepkg (or that of DASPK and LSODE) but will need specific extensions to accommodate for specific options that only make sense for this specific class of solvers. | |||
An initial list of methods to be implemented includes (but is not limited to) | |||
* Symplectic Euler methods | |||
* Störmer-Verlet method | |||
* Velocity Verlet method | |||
* Spectral Variational Integrator methods, see [http://arxiv.org/abs/1211.4534 here] or [http://www.math.ucsd.edu/~mleok/pdf/HaLe2012_SVI.pdf here] | |||
For this latter there is an existing code which is already working but needs to be improved, posted on the patch tracker (TODO) | |||
Furthermore, methods to implement solutions of problems with rigid constraints should be implemented, e.g. | |||
* SHAKE | |||
* RATTLE | |||
== GUI == | == GUI == |
edits