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

Same as #3, but targeting the develop branch instead of master.