We currently initialize joint positions as zeros by default:
https://github.com/ami-iit/jaxsim/blob/9d55da69f3ed8a1db28dfe2fb63371b7482d79d8/src/jaxsim/api/data.py#L125-L133
However, some joints have limits that don't include zero (for example, stickbot shoulder).
We may want to initialize them differently, for example in the mean of limits.