Studium_Angewandte_Informatik_Bachelor/hpc_laboratory_exercises
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

fix

Browse Source
This commit is contained in:
kai
2025-04-26 13:14:42 +02:00
parent 9872088a59
commit b28774b428
1 changed files with 13 additions and 10 deletions
Download Patch File Download Diff File Expand all files Collapse all files

23
lab06/multi_monte.cpp
View File

@@ -13,14 +13,17 @@ struct bench_result {
double monteCarloPi(int64_t N) {
int localCount = 0;
std::mt19937 generator(std::random_device{}());
std::uniform_real_distribution<double> dist(0.0, 1.0);
#pragma omp parallel for reduction(+ : localCount)
for (int i = 0; i < N; ++i) {
double x = dist(generator);
double y = dist(generator);
if (x * x + y * y <= 1.0) {
localCount++;
#pragma omp parallel
{
std::mt19937 generator(std::random_device{}() + omp_get_thread_num());
std::uniform_real_distribution<double> dist(0.0, 1.0);
#pragma omp for reduction(+ : localCount)
for (int i = 0; i < N; ++i) {
double x = dist(generator);
double y = dist(generator);
if (x * x + y * y <= 1.0) {
localCount++;
}
}
}
@@ -35,7 +38,7 @@ bench_result benchmark(int numThreads, int64_t N, double function(int64_t)) {
auto result = bench_result{0.0, 0.0};
// perform benchmark
const int maxNumExec = 10;
const int maxNumExec = 4;
for (int nExec = 0; nExec < maxNumExec; ++nExec) {
result.pi = function(N);
}
@@ -54,7 +57,7 @@ int main() {
fout << "pi" << ",";
fout << "time" << "\n";
for (int64_t i = 1000; i <= max_n; i += 1000) {
for (int64_t i = 1000; i <= max_n; i += 5000) {
const bench_result res = benchmark(numThreads, i, monteCarloPi);
fout << res.pi << ",";