在类中随机使用 C++0x TR1，以降低开销

Question

我正在使用 VC 2010，并试图通过将随机定义放置在每个类实例的构造函数中，然后根据需要从那里进行调用，来降低某些函数的开销和重复代码。我现在所拥有的，简化的是：

#include <random>
#include <Windows.h>
mt19937 eng(GetTickCount());

class Cycles {
    int line;
    normal_distribution<> rand_norm;
    variate_generator<mt19937&,normal_distribution<>> r_norm;
public:
    Cycles() 
    : rand_norm(0.85,0.05),
      r_norm(eng,rand_norm) { 
        line=0; 
    }
}

不幸的是，这不起作用，我最终得到了这个错误：

\vc\include\random(513): 错误 C2248: 'std::tr1::_Ewrap<_Engine,_Tgt_type>::operator =' : 无法访问类 'std::tr1::_Ewrap< 中声明的私有成员;_Engine,_Tgt_type>'

\vc\include\random(446) ：参见 'std::tr1::_Ewrap<_Engine,_Tgt_type>::operator =' 的声明

此诊断发生在编译器生成的函数“std::tr1::variate_generator<_Engine,_Distrib>”中&std::tr1::variate_generator<_Engine,_Distrib>::operator =(const std::tr1::variate_generator<_Engine,_Distrib> &)'

我知道这些应该在构造函数打开之前初始化，否则它会因为缺少默认构造函数而出错，但我不明白为什么会失败。我的C++ fu已经很生疏了。

我见过的每个示例都显示分配器和生成器在调用它的函数中全局或本地初始化，这对我来说似乎很愚蠢，因为我有几个将使用在紧密循环中调用的 r_norm 的成员函数。它严重未通过气味测试。没有人知道我错过了什么吗？

Answer 1

c++0x 中使用 bind() 而不是 variate_generator：

#include <functional> // bind
#include <iostream>
#include <vector>
#include <random>

namespace {
  std::random_device generate_seed;
  std::mt19937 eng(generate_seed());    
  class Cycles {
    int line;
    std::normal_distribution<> rand_norm;
    std::function<double()> r_norm;
  public:
    Cycles(): line(0), rand_norm(0.85, 0.05), r_norm(std::bind(rand_norm, eng))
    {
      // print distribution
      std::vector<size_t> freq(200);
      for (int i = 0; i < 900; ++i) ++freq.at(std::round(r_norm()*100));

      size_t line = 0, max_line = freq.size()-1;
      for ( ; not freq[line] and line <= max_line;  ++line);
      for ( ; not freq[max_line] and max_line >= line; --max_line);
      for ( ; line <= max_line; ++line) {
        std::cout << line << '\t';
        for (size_t j = 0; j < freq[line]; ++j) std::cout << '*';
        std::cout << std::endl;
      }
    }
  };
}    
int main() {
  Cycles c;
}

该图的灵感来自 C++0xFAQ。

输出

72 **
73 ***
74 ******
75 ********
76 *****************
77 ********************
78 *******************************
79 ************************************************
80 *************************************************
81 **********************************************************
82 ************************************************************
83 ***********************************************************
84 ************************************************************************
85 ********************************************************************************
86 *********************************************************************
87 ************************************************************
88 ****************************************************
89 *************************************
90 **********************************
91 **************************************
92 *************************
93 ******************
94 ********************
95 ************
96 ****
97 **********
98 ***
99 **
100 **
101 *

Answer 2

手头没有 MSVC，但是使用 gcc，如果您替换

variate_generator<mt19937&,normal_distribution<>> r_norm;

为

variate_generator<mt19937, normal_distribution<> > r_norm;

（注意缺少 &），

则可以编译。顺便说一句， std::random_device 可能是比 GetTickCount() 更好的随机源，尽管我不知道它是如何在 MSVC 中实现的。

Answer 3

mt19937 eng(GetTickCount());

不能使用非静态表达式初始化全局变量。