boost :: random и boost :iform_real работает с двойными числами, а не с числами с плавающей точкой? - PullRequest
Новая
       34

boost :: random и boost :iform_real работает с двойными числами, а не с числами с плавающей точкой?

2 голосов
/ 24 июня 2011

Простите, если это уже обсуждалось.У меня есть функция шаблона, которая использует boost ::iform_int и boost ::iform_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и удваивается.Он отлично работает с int, отлично работает с double, но он не работает с числами с плавающей точкой.Это как если бы он либо переводил float как int, либо возникла какая-то проблема приведения.Причина, по которой я говорю это потому, что когда я делаю: 

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

, я всегда получаю int обратно.Однако, как я уже сказал, он работает с двойными.Есть ли что-то, что я делаю неправильно или отсутствует?Спасибо!

Ответы [ 3 ]

7 голосов
/ 24 июня 2011

Вы можете даже избежать написания стандартного кода с типом черт и 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(); 
};
7 голосов
/ 24 июня 2011

Аргументы 0.0,5.0 являются двойными, а не плавающими. Сделать их плавающими: 

float y = getRandom<float>(0.0f,5.0f);
5 голосов
/ 24 июня 2011

На самом деле не решение вашего вопроса, но решение:

Почему бы не использовать класс черт, чтобы получить правильный тип распределения? 

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