Fixed bug in the solver. Rotation must be around the cm

iagoleal · iagoleal · commit ea3417736cda · 2020-03-13T21:03:25.000-03:00
diff --git a/src/model.jl b/src/model.jl
@@ -64,7 +64,7 @@ function solve!(m::Model; reltol=1e-8, abstol=1e-8, solver=ODE.Tsit5())
     # Evolution of angular momentum
     m.momentum = t -> solution(t)[5:end]
     # Store solution on model
-    m.inertialframe = t -> map(x -> rotate(x, m.rotation(t)), m.bodyframe(t))
+    m.inertialframe = t -> map(x -> rotate(x, m.rotation(t), center_of_mass(m.bodyframe(t))), m.bodyframe(t))
     return m.inertialframe, m.rotation, m.momentum
 end
 
