A force for determining the CG direction was not communicated
correctly with domain decomposition to recent refactoring change
5ed943bc.
Fixes #3112
Change-Id: Ic33062f4aad0ccd7c1296c7351bc71d735ebba24
rvec *fmg;
snew(fmg, natoms);
- gmx::ArrayRef<const int> indicesMin = s_b->s.cg_gl;
+ gmx::ArrayRef<const int> indicesMin = s_min->s.cg_gl;
int i = 0;
for (int a : indicesMin)
{
copy_rvec(fm[i], fmg[a]);
+ i++;
}
gmx_sum(top_global->natoms*3, fmg[0], cr);