{
gmx::RVec r{ 0, 0, 0 };
pbc_dx(pbc, fr.x[A], fr.x[B], r.as_vec());
-// return r.norm(); // TODO * gmx::c_nm2A; if not PDB, i guess????
+ return r.norm(); // TODO * gmx::c_nm2A; if not PDB, i guess????
return r.norm() * gmx::c_nm2A; // TODO * gmx::c_nm2A; if not PDB, i guess????
}
{
gmx::RVec r{ 0, 0, 0 };
pbc_dx(pbc, A, fr.x[B], r.as_vec());
-// return r.norm(); // TODO * gmx::c_nm2A; if not PDB, i guess????
+ return r.norm(); // TODO * gmx::c_nm2A; if not PDB, i guess????
return r.norm() * gmx::c_nm2A; // TODO * gmx::c_nm2A; if not PDB, i guess????
}