I have a problem where a part of my state vector is a quaternion [q0,q1,q2,q3]. After each integration step I would like to trigger a normalization (sucha that q0^2 +q1^2 +q2^2 +q3^2 = 1). Is this possible with ode-solvers? Or is there some other ODE library in Rust which supports this?
I have a problem where a part of my state vector is a quaternion [q0,q1,q2,q3]. After each integration step I would like to trigger a normalization (sucha that q0^2 +q1^2 +q2^2 +q3^2 = 1). Is this possible with ode-solvers? Or is there some other ODE library in Rust which supports this?