From 1d99310b2e552472e64375da3d131df63c63e9c7 Mon Sep 17 00:00:00 2001 From: Jan Eitzinger Date: Sun, 18 Jun 2023 07:53:01 +0200 Subject: [PATCH] Fix bug in domain decomposition with MPI Cart --- BasicSolver/2D-mpi-v3/src/solver.c | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/BasicSolver/2D-mpi-v3/src/solver.c b/BasicSolver/2D-mpi-v3/src/solver.c index 3378e94..cf9baf5 100644 --- a/BasicSolver/2D-mpi-v3/src/solver.c +++ b/BasicSolver/2D-mpi-v3/src/solver.c @@ -333,8 +333,8 @@ void initSolver(Solver* solver, Parameter* params) &solver->jNeighbours[1]); MPI_Cart_get(solver->comm, NDIMS, solver->dims, periods, solver->coords); - solver->imaxLocal = sizeOfRank(solver->rank, dims[IDIM], solver->imax); - solver->jmaxLocal = sizeOfRank(solver->rank, dims[JDIM], solver->jmax); + solver->imaxLocal = sizeOfRank(solver->coords[IDIM], dims[IDIM], solver->imax); + solver->jmaxLocal = sizeOfRank(solver->coords[JDIM], dims[JDIM], solver->jmax); MPI_Datatype jBufferType; MPI_Type_contiguous(solver->imaxLocal, MPI_DOUBLE, &jBufferType);