fmmax = nfr - 1;
}
snew(mcount, fmmax);
- nbin = static_cast<int>(rmax*invbin + 0.5);
+ nbin = gmx::roundToInt(rmax*invbin);
if (sbin == 0)
{
mat_nx = fmmax + 1;
}
else
{
- mbin = static_cast<int>(std::sqrt(fbin*dt)*invsbin + 0.5);
+ mbin = gmx::roundToInt(std::sqrt(fbin*dt)*invsbin);
}
for (i = 0; i < isize; i++)
{
d2 = distance2(sx[f][i], sx[ff][i]);
if (mbin < mat_nx && d2 < rmax2)
{
- bin = static_cast<int>(std::sqrt(d2)*invbin + 0.5);
+ bin = gmx::roundToInt(std::sqrt(d2)*invbin);
if (bin < nbin)
{
mat[mbin][bin] += 1;
for (i = 0; i < isize; i++)
{
d2 = distance2(sx[f][i], sx[ff][i]);
- bin = static_cast<int>(std::sqrt(d2)*invbin + 0.5);
+ bin = gmx::roundToInt(std::sqrt(d2)*invbin);
if (bin >= nalloc)
{
nallocn = 10*(bin/10) + 11;