Marsaglia C++の通常のランダム変数

Question

私のクエリはかなりシンプルです。私はJavaで数値メソッドをコード化することを好むが、しばしばC++でいくつかのことを行う必要がある。私はMarsagliaアルゴリズムを使用しており、通常の確率変数を両方保持しているので、Javaのガウス確率変数が気に入っています。最初の呼び出しで1を返し、2番目の呼び出しで2番目の呼び出しを返し、3番目の呼び出しまで高価な計算を再度実行しません。下のoracleリンクを使用して（プログラムのコメント内で）、このコードをC++で実装しようとしましたが、C++バージョンの「Synchronized」パブリック・メソッドを作成して、両方の正規乱数変数を使用できるようにする方法はわかりません。私はプロのプログラマーではないので、どんな指導も大変ありがとうございます。

私は維持したいと思います要するに

：

V2 *乗数

// This function is Similar to the GNU 
// Java Implementation as seen on 
// http://docs.oracle.com/javase/1.4.2/docs/api/java/util/Random.html#nextGaussian%28%29 
double nextGaussian() { 

    double v1, v2, s, nextNextGaussian; 
    do { 
     v1 = 2 * nextUniform() - 1; // between -1.0 and 1.0 
     v2 = 2 * nextUniform() - 1; // between -1.0 and 1.0 
     s = v1 * v1 + v2 * v2; 
    } while (s >= 1 || s == 0); 
    double multiplier = sqrt(-2 * log(s)/s); 
    nextNextGaussian = v2 * multiplier; 
    return v1 * multiplier; 

} 

Answer 1

だけnextGaussianValのその後値すなわち

static double nextGaussianVal; 

、静的としてnextGaussianValを宣言はnexで利用可能になりますメソッドが呼び出された時刻。あなたも現在のカウントに追いつくために、別の静的変数が必要な場合がありますので、同じよう：

double nextGaussian() 
{ 
    static int count = 0; 
    static double nextGaussianVal; 
    double firstGaussianVal, v1, v2, s; 

    if (count == 0) { 
     do { 
      v1 = 2 * nextUniform() - 1; // between -1.0 and 1.0 
      v2 = 2 * nextUniform() - 1; // between -1.0 and 1.0 
      s = v1 * v1 + v2 * v2; 
     } while (s >= 1 || s == 0); 
     double multiplier = sqrt(-2 * log(s)/s); 
     nextGaussianVal = v2 * multiplier; 
     firstGaussianVal = v1 * multiplier; 
     count = 1; 
     return firstGaussianVal; 
    } 

    count = 0; 
    return nextGaussianVal; 
} 

編集：より詳細な説明 - 関数が呼び出される初めて、カウントがゼロに初期化されます。 if文に基づき、問題の計算が行われ、数は1の値が割り当てられ、firstGaussianValとnextGaussianValに値が割り当てられていることを想定して、firstGaussianValが返されます。 カウントには、前に割り当てられた値が1つあり、nextGaussianValには、最初の呼び出しで割り当てられた値が入ります。つまり、カウントが1になり、関数は、ステートメントの場合に基づいて、カウントに0を割り当て、を返します.GaussianVal。リンス、繰り返し...

Answer 2

オブジェクト指向の問題では、「乱数ジェネレータ」オブジェクトにそのようなものを保持する必要があります。たとえば、次のコードを見てください： https://code.cor-lab.org/projects/nemomath/repository/entry/trunk/nemomath/src/nemo/Random.h クラス「ガウス」は、アルゴリズム的な方法で実装したいものです。

Answer 3

上記の提案の問題は、COUNT変数が整数ではなく、ブール値であることです。さらに、格納されたガウス関数もまたSTATICである必要があります。

私が望む正しい解決策に私が手伝ってくれたことに皆様のおかげで感謝したいと思います。私はこのスニペットがオブジェクトに属していることを認識しています。私は今、1つのファイルに保存する理由があります。

double nextGaussian() { 

    // Static variables allow the function to make use of 
    // both Gaussian Random variables. Generated by the polar Marsaglia 
    // method This makes the function much more efficient with will 
    // pay off for simulations 

    static bool hasNextGaussian = false; 
    static double nextNextGaussian; 
    double v1, v2, s; 

    if (!hasNextGaussian) { 

     do { 

      v1 = 2 * nextUniform() - 1; // between -1.0 and 1.0 
      v2 = 2 * nextUniform() - 1; // between -1.0 and 1.0 
      s = v1 * v1 + v2 * v2; 

     } while (s >= 1 || s == 0); 

     double multiplier = sqrt(-2 * log(s)/s); 
     nextNextGaussian = v2 * multiplier; 
     hasNextGaussian = true; 
     return v1 * multiplier; 

    } else { 

     hasNextGaussian = false; 
     return nextNextGaussian; 

    } 

} 

ダブルnextUniform（）{

double Uniform = rand()/double(RAND_MAX); 
return Uniform; 

}