Skip to main content

次のタグが付いている話題の回答:

3

Open MPで並列化したら精度が落ちるモンテカルロ法

疑似乱数生成(メルセンヌ・ツイスター)の並列化とほぼ同じ問題を抱えています。 OpenMPの並列ループ内で乱数engineを参照しています。マルチスレッド安全ではない気がします。最善でもスレッド間で同じ乱数値を取得することになりその時点で精度は1スレッド分程度に落ちます。つまり8並列であればN/8回の施行しか意味を持たなくなる。最悪な場合、std::...
sayuri's user avatar
  • 4.3万
2
承認済み

疑似乱数生成(メルセンヌ・ツイスター)の並列化

本質的でないuniform_realとvariate_generatorを取り除いて説明します。 void test1() { printf("%s\n", __func__); std::mt19937 gen{ 123 }; #pragma omp parallel for private(gen) for (int i = 0; i < 16; ++i)...
sayuri's user avatar
  • 4.3万

無効化に投票