- constexpr double turn_speed = 10.0;
- // TODO: this is crap
- state.dir = glm::normalize(state.dir + f.turn * turn_speed * dt);
+ if (length2(state.vel) > 0.000001) {
+ glm::dvec3 nvel(normalize(state.vel));
+ double ang = angle(nvel, state.dir);
+ double turn_rate = PI * dt;
+ if (ang < turn_rate) {
+ state.dir = normalize(state.vel);
+ } else if (std::abs(ang - PI) < 0.001) {
+ state.dir = rotate(state.dir, turn_rate, world::Planet::SurfaceNormal(situation.Surface()));
+ } else {
+ state.dir = rotate(state.dir, turn_rate, normalize(cross(state.dir, nvel)));
+ }
+ }
+