c-std :: piecewise_linear_distribution在VS2012下不工作？

我决定知道c 11<随机>更好,所以我写了这样一个代码：

std::mt19937 gen(10);
std::piecewise_linear_distribution<> d(Range.begin(),Range.end(),RangeValues.begin());

std::map<int,unsigned int> hist;
for (int n = 0; ++n != iterations;)
    ++hist[std::round(d(gen))];

for (auto p : hist)
   std::cout << std::setw(2) << p.first << ": "
             << std::string(p.second/(iterations/200),'*') << '\n';

由于某些原因std :: random_device似乎不适用于Coliru,所以我输入了一个常量样本种子.我认为,这是UB,因此IIRC它基于硬件,这就是为什么它不支持Coliru(纠正我,如果我错了).在Ideone它的工作,但是.

然后我把它移植到VS2012上运行,唯一的区别就是自己执行的std :: round：

return number < 0.0 ? std::ceil(number - 0.5) : std::floor(number + 0.5);

它在Coliru上完美正确,但是当我编译并在VS2012,the output is just wrong上运行它.

任何想法如何纠正这一点,更重要的是,为什么会发生这种情况？
我做某事迟缓了,还是VS2012在这里不聪明？

解决方法

看来这是一个Visual Studio的问题.我已经尝试了下面的程序(从OP改编),由GCC 4.7.2,Clang 3.2和Intel 13.1.0生成的输出是非常合理的,而Visual Studio 2012年11月生成的CTP完全不同.

概率密度是分段线性的,由数组x和p以下列方式定义.构建了i = 0,…,N(其中N = x.size()-1)的点(x [i],p [i])的分段线性函数.然后将该函数归一化(通过将其除以其积分)以获得概率密度.

#include <iostream>
#include <iomanip>
#include <string>
#include <random>
#include <array>

int main() {

    std::mt19937 gen(10);

    std::array<double,3> x = {{0,20,40}};
    std::array<double,3> p = {{0,1,0}};
    std::piecewise_linear_distribution<> dist(x.begin(),x.end(),p.begin());

    std::array<int,40> hist = {{0}};

    for (size_t i = 0; i < 200000; ++i)
        ++hist[static_cast<size_t>(dist(gen))];

    for (size_t n = 0; n < hist.size(); ++n)
        std::cout << std::setfill('0') << std::setw(2) << n << ' ' << 
          std::string(hist[n] / 200,'*') << std::endl;

    std::cout << "\nValues in interval [20,21[ : " << hist[20] << std::endl;
}

在我们的例子中,多边形函数连接(0,0),(20,1)和(40,0).因此,其形状是基座40和高度1的等腰三角形,其面积为20.因此,概率密度f连接(0,1 / 20)和(40,0).这意味着在间隔[20,21 [我们可以期待f(20)*(21 – 20)= 1/20 * 1 = 1/20的结果.总共我们画出20万个价值,然后我们可以预期在[20,21]中有大约10,000点.

GCC,Clang和Intel在[20,21]中报告了9734点,并显示了与等腰三角形非常相似的图案：

00 *
01 ***
02 *****
03 ********
04 ***********
05 **************
06 ***************
07 ******************
08 ********************
09 ************************
10 **************************
11 ****************************
12 *******************************
13 *********************************
14 ***********************************
15 ***************************************
16 *****************************************
17 ******************************************
18 **********************************************
19 ************************************************
20 ************************************************
21 *********************************************
22 *******************************************
23 *****************************************
24 **************************************
25 ************************************
26 **********************************
27 ******************************
28 ****************************
29 **************************
30 ***********************
31 ********************
32 ******************
33 ****************
34 *************
35 ***********
36 *********
37 ******
38 ***
39 *

Values in interval [20,21[ : 9734

不幸的是,Visual Studio 2012年11月CTP给出了这一点：

00 ********************************************** [truncated]
01 **********************************************
02 ***********************************
03 *****************************
04 **************************
05 ***********************
06 *********************
07 ********************
08 *******************
09 ******************
10 *****************
11 ****************
12 ***************
13 **************
14 **************
15 **************
16 *************
17 *************
18 *************
19 ************
20 ************
21 *************
22 *************
23 *************
24 *************
25 **************
26 ***************
27 ***************
28 ****************
29 *****************
30 ******************
31 *******************
32 *******************
33 *********************
34 ***********************
35 **************************
36 *****************************
37 ***********************************
38 **********************************************
39 ********************************************** [truncated]

Values in interval [20,21[ : 2496

笔记：

>我截断了Visual Studio输出以更好地显示.对于[20,21 [20万*(0.5 *(f(20)f(21)))*(21 – 20)= 100,000 *(1/20 1/20 – 1) / 400)= 10,000-250 = 9750.

原文地址：https://www.jb51.cc/c/114456.html