- bOK = constrain_lincs(bLog || bEner, ir, step, lincsd, md, cr, ms, x, xprime, min_proj, box,
- pbc_null, lambda, dvdlambda, invdt, v.unpaddedArrayRef(),
- vir != nullptr, vir_r_m_dr, econq, nrnb, maxwarn, &warncount_lincs);
+ bOK = constrain_lincs(bLog || bEner, ir, step, lincsd, md.invmass, cr, ms, x, xprime,
+ min_proj, box, pbc_null, md.nMassPerturbed != 0, lambda, dvdlambda,
+ invdt, v.unpaddedArrayRef(), vir != nullptr, vir_r_m_dr, econq, nrnb,
+ maxwarn, &warncount_lincs);