3
票
Open MPで並列化したら精度が落ちるモンテカルロ法
疑似乱数生成(メルセンヌ・ツイスター)の並列化とほぼ同じ問題を抱えています。
OpenMPの並列ループ内で乱数engineを参照しています。マルチスレッド安全ではない気がします。最善でもスレッド間で同じ乱数値を取得することになりその時点で精度は1スレッド分程度に落ちます。つまり8並列であればN/8回の施行しか意味を持たなくなる。最悪な場合、std::...
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)...
無効化に投票