同一发电机值的std :: normal_distribution结果?

时间:2018-11-08 01:09:53

标签: c++ random normal-distribution c++-standard-library

我注意到GNU C ++标准库的行为很奇怪。

使用具有正态分布的自定义生成器时,即使对于相同的生成器结果,该分布似乎也可能返回不同的值。

示例代码

#include <iostream>
#include <string>
#include <random>

template <int MAX>
class Generator
{
public:
  long long operator()() {return (++state) % MAX;}
  long long min() {return 0;}
  long long max() {return MAX;}
  long long state = 0;
};

int main()
{
  // 3 identical generators, 2 identical distributions
  Generator<5> gen0, gen1, gen2;
  std::normal_distribution<float> dist1(10, 1), dist2(10,1);

  // Print: iteration, generator output, distributions output
  printf("i   gen       dist1      dist2\n");    
  for (int i = 0; i < 10; i++) {
    printf("%d     %lld  %10.6f %10.6f\n", i, gen0(), dist1(gen1), dist2(gen2));
    dist2.reset();
  }

}  

这会产生(g ++ 6.2.1,与更新的相同):

i   gen       dist1      dist2
0     1   11.295982  11.295982
1     2    9.325577  10.256990
2     3   10.256990   9.672720
3     4    8.679697   9.300377
4     0    9.672720   8.957100
5     1    9.097454  11.295982
6     2    9.300377  10.256990
7     3   10.582899   9.672720
8     4    8.957100   9.300377
9     0   10.169774   8.957100

看看dist1在第0行和第5行中的结果如何不同,即使生成器产生相同的值1也是如此。对于dist2,它们与我们称为reset()的相同。

顺便说一句,原因似乎是GNU C ++库成对生成值,并在调用之间缓存它们,而不必调用生成器。这也解释了为什么dist1第4行与dist2(重置发行版)第2行相同。

我在文档中进行了搜索,但无法确定此行为是否符合标准。因此,我不知道这是否是标准库错误:)

请注意,根据我的观察,std :: normal_distribution是唯一表现出此行为的类。

1 个答案:

答案 0 :(得分:3)

符合标准。不需要生成器是无状态的。如果有的话,reset函数将毫无用处。