boost :: random和boost:uniform_real可以使用双倍而不是浮点数吗?

时间:2011-06-24 15:28:21

标签: c++ boost floating-point uniform boost-random

如果已经讨论过,请原谅我。我有一个模板函数,它根据模板参数使用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转换为int,或者存在一些转换问题。原因我说这是因为当我这样做时:

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

我总是得到一个int。 但是,就像我说的那样,它适用于双打。 有什么我做错了或错过了吗? 谢谢!

3 个答案:

答案 0 :(得分:7)

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

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

答案 1 :(得分:7)

您甚至可以避免使用类型特征和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(); 
};

答案 2 :(得分:5)

本身并没有真正解决你的问题,而是一个解决方案:

为什么不使用traits类来获得正确的分发类型?

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(); 
};