Replace std::numeric_limits<T>::epsilon calls with boost::math::tools::epsilon<T>

Same as #3, but targeted at the develop branch.