Rc,
*krf,
*crf,
- ONE_4PI_EPS0 / eps_r);
+ gmx::c_one4PiEps0 / eps_r);
// Make sure we don't lose resolution in pow() by casting real arg to double
real rmin = gmx::invcbrt(static_cast<double>(*krf * 2.0));
fprintf(fplog, "The electrostatics potential has its minimum at r = %g\n", rmin);