-
Notifications
You must be signed in to change notification settings - Fork 33
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
Easier initial guess for acceleration in dynamics. #1230
Conversation
src/serac/numerics/odes.cpp
Outdated
d2u_dt2 = state_.d2u_dt2; | ||
d2u_dt2.SetSubVector(constrained_dofs, 0.0); | ||
d2U_dt2_.SetSubVectorComplement(constrained_dofs, 0.0); | ||
// d2U_dt2_.SetSubVectorComplement(constrained_dofs, 0.0); |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Looks good. Why not remove the commented line?
7b878f2
to
a1c8455
Compare
a1c8455
to
67b54f6
Compare
tests
Outdated
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
@chapman39 Is there a reason why these tests need to be in a submodule? I find it nice to always have code and tests update together in the same repo, when possible. Also, I have no idea if I even did this submodule thing correctly.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
The tests are in a submodule so that if the repo gets large with all the data that builds up, it can be cleaned out w/o having to alter the history of the serac repo. This has happened on all simulation codes I've been on
/style |
d2274fe
to
af5edb5
Compare
/style |
Sometimes the initial acceleration guess results in element inversion. This change defaults the guess to 0.0, which is more robust in a challenging problem, but could potentially be slower for some kinds of problems.