boost::random 和 boost:uniform_real 适用于双精度而不适用于浮点数？

Question

如果已经讨论过这个问题，请原谅我。我有一个模板函数，它根据模板参数使用 boost::uniform_int 和 boost::uniform_real ，并且应该返回相同的类型：

template <typename N> N getRandom(int min, int max)
{
  timeval t;
  gettimeofday(&t,NULL);
  boost::mt19937 seed((int)t.tv_sec);
  boost::uniform_int<> dist(min, max);
  boost::variate_generator<boost::mt19937&, boost::uniform_int<> > random(seed, dist);
  return random(); 
}
//! partial specialization for real numbers
template <typename N> N getRandom(N min, N max)
{
  timeval t;
  gettimeofday(&t,NULL);
  boost::mt19937 seed( (int)t.tv_sec );
  boost::uniform_real<> dist(min,max);
  boost::variate_generator<boost::mt19937&, boost::uniform_real<> > random(seed,dist);
  return random(); 
}

现在我已经用 int、float 和 double 测试了该函数。 它适用于 int，适用于 double，但不适用于 float。就好像它要么将 float 转换为 int，要么存在一些转换问题。我这么说的原因是因为当我这样做时：

float y = getRandom<float>(0.0,5.0);

我总是得到一个 int 返回。 然而，就像我说的，它适用于双打。 我做错了什么或遗漏了什么吗？ 谢谢 ！

Answer 1

参数 0.0,5.0 是双精度数，而不是浮点数。让它们漂浮：

float y = getRandom<float>(0.0f,5.0f);

Answer 2

您甚至可以避免编写具有类型特征和 MPL 的样板代码：

template <typename N>
N getRandom(N min, N max)
{
  typedef typename boost::mpl::if_<
    boost::is_floating_point<N>, // if we have a floating point type
    boost::uniform_real<>,       // use this, or
    boost::uniform_int<>         // else use this one
  >::type distro_type;

  timeval t;
  gettimeofday(&t,NULL);
  boost::mt19937 seed( (int)t.tv_sec );
  distro_type dist(min,max);
  boost::variate_generator<boost::mt19937&, distro_type > random(seed,dist);
  return random(); 
};

Answer 3

本身并没有真正解决您的问题，而是解决方案：

为什么不使用特征类来获得正确的分布类型？

template<class T>
struct distribution
{ // general case, assuming T is of integral type
  typedef boost::uniform_int<> type;
};

template<>
struct distribution<float>
{ // float case
  typedef boost::uniform_real<> type;
};

template<>
struct distribution<double>
{ // double case
  typedef boost::uniform_real<> type;
};

有了这一组，您就可以拥有一个通用功能：

template <typename N> N getRandom(N min, N max)
{
  typedef typename distribution<N>::type distro_type;

  timeval t;
  gettimeofday(&t,NULL);
  boost::mt19937 seed( (int)t.tv_sec );
  distro_type dist(min,max);
  boost::variate_generator<boost::mt19937&, distro_type > random(seed,dist);
  return random(); 
};