...
 
Commits (3)
......@@ -61,7 +61,7 @@ int main(int argc, char const *argv[])
/* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
/* COORDINATE PARTITIONING *
/* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * */
mbs_data->process = 1;
mbs_part = mbs_new_part(mbs_data);
mbs_part->options->rowperm=1;
mbs_part->options->verbose = 1;
......@@ -73,7 +73,7 @@ int main(int argc, char const *argv[])
/* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
/* EQUILIBRIUM *
/* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * */
mbs_data->process = 2;
mbs_equil = mbs_new_equil(mbs_data);
// dirdyn options (see documentations for additional options)
mbs_equil->options->method = 1;
......
......@@ -82,7 +82,7 @@ int invdynared(MbsAux *mbs_aux, MbsData *s)
}
}
copy_dvec_0(mbs_aux->Rred+1, mbs_aux->b+1, s->nqa);
copy_dvec_0(mbs_aux->Rred+1, mbs_aux->b+1, s->nqu);
// Switch case between over-iso-under actuated system
if (s->nqa == s->nqu) { // solution of the original equations: Qact(ind_a) = A\b;
......