This should just require implementing some abs() function on AlgebraicVector, but needs to be checked with details [here](https://www.boost.org/doc/libs/1_82_0/libs/numeric/odeint/doc/html/boost_numeric_odeint/odeint_in_detail/state_types__algebras_and_operations.html)
