constr,
do_log,
do_ene);
- upd->finish_update(*ir, mdatoms, state, wcycle, constr != nullptr);
+ upd->finish_update(
+ *ir, mdatoms->havePartiallyFrozenAtoms, mdatoms->homenr, state, wcycle, constr != nullptr);
if (ir->eI == IntegrationAlgorithm::VVAK)
{
* to numerical errors, or are they important
* physically? I'm thinking they are just errors, but not completely sure.
* For now, will call without actually constraining, constr=NULL*/
- upd->finish_update(*ir, mdatoms, state, wcycle, false);
+ upd->finish_update(*ir, mdatoms->havePartiallyFrozenAtoms, mdatoms->homenr, state, wcycle, false);
}
/* this factor or 2 correction is necessary
because half of the constraint force is removed