https://github.com/mseri/ocaml-cviode

Contact variational integrators - native ocaml implementation
https://github.com/mseri/ocaml-cviode

geometric-integrator hacktoberfest solvers

Last synced: 11 months ago
JSON representation

Contact variational integrators - native ocaml implementation

• Host: GitHub
• URL: https://github.com/mseri/ocaml-cviode
• Owner: mseri
• License: mit
• Created: 2019-02-01T09:34:23.000Z (over 7 years ago)
• Default Branch: master
• Last Pushed: 2020-10-04T09:35:35.000Z (over 5 years ago)
• Last Synced: 2025-01-22T00:36:01.476Z (over 1 year ago)
• Topics: geometric-integrator, hacktoberfest, solvers
• Language: OCaml
• Homepage:
• Size: 254 KB
• Stars: 3
• Watchers: 3
• Forks: 0
• Open Issues: 0
• Metadata Files:
  • Readme: README.md
  • Changelog: CHANGES.md
  • License: LICENSE

Awesome Lists containing this project

README

          
# Contact variational integrators - native ocaml implementation

This is a collection of geometric solvers for initial value problems derived from contact Lagrangians.

The provided solvers concern Lagrangians of the form

```latex

$$                                         

L(x, \dot{x}, z, t) = \frac12|\dot{x}|^2 + g_1(x)g_2(z) + h(z) + f(t)x

$$

```

For further information refer to [_Vermeeren, Bravetti, Seri: Contact Variational Integrators (2019)_](https://arxiv.org/abs/1902.00436).

Running the integrators on the damped oscillator with and without forcing (see [example/damped.ml](example/damped.ml)) produces the following output:

![no forcing, critical damping](img/damped.png)

![no forcing, critical damping](img/forced.png)

## TODO

- [ ] Add mli with documentation

- [ ] Add implementation with support for $g_2(z)$ as per description

- [ ] Figure out how to test