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

add tags

Browse Source
This commit is contained in:
kai
2025-06-01 12:49:51 +02:00
parent c4277be067
commit 5b0bcec2ac
1 changed files with 16 additions and 16 deletions
Download Patch File Download Diff File Expand all files Collapse all files

32
lab10/jacobi/jacobi_mpi.cpp
View File

@@ -53,12 +53,12 @@ Jacobi::Result JacobiMPI::run(const Matrix &init, double epsilon,
if (rank == 0) {
MPI_Request send_lower;
std::vector<double> send_vec_low = phi[t0].get_row(numRows - 1);
MPI_Isend(send_vec_low.data(), numCols, MPI_DOUBLE, neighborLower, 0,
MPI_COMM_WORLD, &send_lower);
MPI_Isend(send_vec_low.data(), numCols, MPI_DOUBLE, neighborLower,
indexRowGlobalEnd, MPI_COMM_WORLD, &send_lower);
MPI_Request requestLower;
MPI_Irecv(haloLower.data(), numCols, MPI_DOUBLE, neighborLower, 0,
MPI_COMM_WORLD, &requestLower);
MPI_Irecv(haloLower.data(), numCols, MPI_DOUBLE, neighborLower,
rowHaloLowerIndex, MPI_COMM_WORLD, &requestLower);
MPI_Wait(&send_lower, MPI_STATUS_IGNORE);
MPI_Wait(&requestLower, MPI_STATUS_IGNORE);
@@ -79,12 +79,12 @@ Jacobi::Result JacobiMPI::run(const Matrix &init, double epsilon,
else if (rank == (numProc - 1)) {
MPI_Request send_upper;
std::vector<double> send_vec_up = phi[t0].get_row(rank * numRows);
MPI_Isend(send_vec_up.data(), numCols, MPI_DOUBLE, neighborUpper, 0,
MPI_COMM_WORLD, &send_upper);
MPI_Isend(send_vec_up.data(), numCols, MPI_DOUBLE, neighborUpper,
indexRowGlobalStart, MPI_COMM_WORLD, &send_upper);
MPI_Request request_upper;
MPI_Irecv(haloUpper.data(), numCols, MPI_DOUBLE, neighborUpper, 0,
MPI_COMM_WORLD, &request_upper);
MPI_Irecv(haloUpper.data(), numCols, MPI_DOUBLE, neighborUpper,
rowHaloUpperIndex, MPI_COMM_WORLD, &request_upper);
MPI_Wait(&send_upper, MPI_STATUS_IGNORE);
MPI_Wait(&request_upper, MPI_STATUS_IGNORE);
@@ -106,20 +106,20 @@ Jacobi::Result JacobiMPI::run(const Matrix &init, double epsilon,
MPI_Request send_lower;
std::vector<double> send_vec_up = phi[t0].get_row(rank * numRows);
MPI_Isend(send_vec_up.data(), numCols, MPI_DOUBLE, neighborUpper, 0,
MPI_COMM_WORLD, &send_upper);
MPI_Isend(send_vec_up.data(), numCols, MPI_DOUBLE, neighborUpper,
indexRowGlobalStart, MPI_COMM_WORLD, &send_upper);
std::vector<double> send_vec_low =
phi[t0].get_row(rank * numRows + numRows);
MPI_Isend(send_vec_low.data(), numCols, MPI_DOUBLE, neighborLower, 0,
MPI_COMM_WORLD, &send_lower);
MPI_Isend(send_vec_low.data(), numCols, MPI_DOUBLE, neighborLower,
indexRowGlobalEnd, MPI_COMM_WORLD, &send_lower);
MPI_Request request_upper;
MPI_Request request_lower;
MPI_Irecv(haloUpper.data(), numCols, MPI_DOUBLE, neighborUpper, 0,
MPI_COMM_WORLD, &request_upper);
MPI_Irecv(haloLower.data(), numCols, MPI_DOUBLE, neighborLower, 0,
MPI_COMM_WORLD, &request_lower);
MPI_Irecv(haloUpper.data(), numCols, MPI_DOUBLE, neighborUpper,
rowHaloUpperIndex, MPI_COMM_WORLD, &request_upper);
MPI_Irecv(haloLower.data(), numCols, MPI_DOUBLE, neighborLower,
rowHaloLowerIndex, MPI_COMM_WORLD, &request_lower);
MPI_Wait(&send_upper, MPI_STATUS_IGNORE);
MPI_Wait(&send_lower, MPI_STATUS_IGNORE);