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

Avoid using cross product for "v x (a v)" expecting the result is zero #229

Open
jslee02 opened this issue Aug 5, 2014 · 2 comments
Open
Assignees

Comments

@jslee02
Copy link
Member

jslee02 commented Aug 5, 2014

The result of v x (a v) is not zero.

double a = 10.0;
Eigen::Vector3d v = Eigen::Vector3d::Random();

Eigen::Vector3d res1 = v.cross(v);  // res1 is zero.
Eigen::Vector3d res2 = v.cross(a * v);  // res2 is not zero!

This case occurs in math::dad(V, GV).

@jslee02 jslee02 added this to the Release DART (14.12) milestone Aug 5, 2014
@jslee02 jslee02 changed the title Avoid using cross product for "v x (a v)" where v is 3-dim vector and a is scalar Avoid using cross product for "v x (a v)" where "v" is a 3-dim vector and "a" is a scalar Aug 5, 2014
@scpeters
Copy link
Collaborator

scpeters commented Aug 5, 2014

Is this related to the rigid body rotation errors?

@jslee02
Copy link
Member Author

jslee02 commented Aug 5, 2014

This is related to the angular momentum error but not to other errors such as linear velocity error.

@jslee02 jslee02 changed the title Avoid using cross product for "v x (a v)" where "v" is a 3-dim vector and "a" is a scalar Avoid using cross product for "v x (a v)" expecting the result is zero Dec 8, 2014
@jslee02 jslee02 modified the milestones: Release DART (15.06), Release DART 4.3 Jan 4, 2015
@jslee02 jslee02 removed this from the Release DART 5.0 milestone Jun 11, 2015
@jslee02 jslee02 added this to the DART 6.7.0 milestone Sep 5, 2018
@jslee02 jslee02 self-assigned this Sep 5, 2018
@jslee02 jslee02 removed this from the DART 6.7.0 milestone Jan 10, 2019
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants