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

.

Browse Source
This commit is contained in:
WickedJack99
2025-06-29 22:34:59 +02:00
parent c68e74cd52
commit f9a1248740
1 changed files with 6 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

10
lab12/exc4/jacobi.cpp
View File

@@ -240,7 +240,7 @@ Jacobi::Result Jacobi::run(const Matrix &init, double eps, int maxNumIter)
int t0 = 0; // array index of current timestep
int t1 = 1; // array index of next timestep
#pragma omp parallel firstprivate(t0, t1) lastprivate(nIter)
#pragma omp parallel
while (dist > eps && nIter < maxNumIter)
{
dist = 0;
@@ -262,10 +262,12 @@ Jacobi::Result Jacobi::run(const Matrix &init, double eps, int maxNumIter)
}
#pragma omp single
MPI_Allreduce(MPI_IN_PLACE, &dist, 1, MPI_DOUBLE, MPI_MAX, MPI_COMM_WORLD);
{
MPI_Allreduce(MPI_IN_PLACE, &dist, 1, MPI_DOUBLE, MPI_MAX, MPI_COMM_WORLD);
nIter++;
std::swap(t0, t1);
nIter++;
std::swap(t0, t1);
}
}
return {removeHaloLayers(phi[t0]), dist, nIter};