Commit 81ec4f36 authored by Nicolas Docquier's avatar Nicolas Docquier
Browse files

Update state vector with qdd[ind_qu] rather than udd

parent 4e730968
......@@ -391,16 +391,8 @@ void mbs_fct_dirdyn(double tsim, double y[], double dydt[], MbsData *s, MbsAux *
// Update state vector
for(i=1;i<=s->nqu;i++)
{
dydt[i-1] = s->qd[s->qu[i]];
if (opts->accelred)
{
dydt[i+s->nqu-1] = s->qdd[i];
}
else
{
dydt[i+s->nqu-1] = s->udd[i];
}
dydt[i-1] = s->qd [s->qu[i]];
dydt[i+s->nqu-1] = s->qdd[s->qu[i]];
}
for(i=1;i<=s->Nux;i++)
{
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment