c++ - 指数ジェネレーターは時々「奇妙な」結果を与える

Question

私は C++ でシミュレーションを構築しており、プロセスのバースト時間を作成するための指数ジェネレーターを持っています。

通常、次のような値を返します: 3.14707,1.04998. しかし、10 分の 1 の確率で、2.64823e-307 という数値が得られます。

これはジェネレーターのコードです（srand ( time(NULL) );プログラムの最初に使用しています）：

double exponential(float u)
{
   double x,mean;
   mean = 10;
   // generate a U(0,1) random variate
   x = rand();
   u = x / RAND_MAX;
   return (-mean * log(u));
}

そして、これが私が値を割り当てる方法です。内部のwhile部分は、そのような値を取り除くための私の努力ですが、うまくいきませんでした:

for (int i = 0; i < nPages; i++)
{
    index[i] = i;
arrival[i]= poisson(r);
burst[i]=exponential(u);
while (burst[i]<1 || burst[i]>150)
    {
        cout<<"P"<<i<<endl;
        burst[i]=(burst[i-1]+burst[i+1])/2;
    }
}

Answer 1

なぜ C++ ライブラリではなく C ライブラリを使用するのですか??

std::random_device rd;
std::default_random_engine gen(rd());
std::exponential_distribution<double> dist(lambda);

double x = dist(gen);