Introduce hjb method

[PierreMartinon]

Start with the basics, no modes.
Use a second entry point solve_global() ?
Add a warmstart from the hjb trajectory.
Maybe also in reverse, a grid definition around a local solution, to check optimally with hjb ?