Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Resolution of Arnold's equation with minimal testing #2351

Merged
merged 28 commits into from
Oct 6, 2019

Conversation

pleroy
Copy link
Member

@pleroy pleroy commented Oct 6, 2019

No description provided.

physics/euler_solver.hpp Outdated Show resolved Hide resolved
physics/euler_solver_body.hpp Outdated Show resolved Hide resolved
Bivector<Variation<AngularMomentum>, PrincipalAxesFrame> const ṁ =
Commutator(m, ω) / Radian;

// If ṁ is constant in the principal axes frame, pick 𝒫ₜ = identity.
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

If m is constant, or If ṁ is 0.
But how is this correct? If 𝒫ₜ is the identity in coordinates, 𝒫ₜ(m_normalized) is e₃ (in ℬₜ) only if m_normalized is e₃ (in PrincipalAxesFrame); but m can be constant for at least 6 different values of m_normalized (±eᵢ in PrincipalAxesFrame), and for many values if some moments of inertia are equal.

Copy link
Member Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Done, it gets a bit ugly.

if (m.x < AngularMomentum()) {
B₁₃_ = -B₁₃_;
σB₁₃_ = -B₁₃_;
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Add a TODO to investigate the behaviour with respect to nothing's sign bit.

@eggrobin eggrobin added the LGTM label Oct 6, 2019
@pleroy pleroy merged commit 129e850 into mockingbirdnest:master Oct 6, 2019
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Projects
None yet
Development

Successfully merging this pull request may close these issues.

None yet

2 participants