ErikFabrizzi/NuSiF-Solver
1
0
Fork 0
forked from moebiusband/NuSiF-Solver
Code Pull Requests Activity

fix: corrected bd-> result still not okay

Browse Source
This commit is contained in:
Erik Fabrizzi 2025-04-28 19:03:44 +02:00
parent 765883ab88
commit 356c4fe4fa
3 changed files with 13 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
PoissonSolver/2D-mpi/src/main.c
View File

@ -35,7 +35,7 @@ int main(int argc, char** argv)
}
initSolver(&solver, &params, 2);
solve(&solver);
solveRB(&solver);
getResult(&solver);
MPI_Finalize();

17
PoissonSolver/2D-mpi/src/solver.c
View File

@ -310,12 +310,17 @@ int solveRB(Solver *solver) {
}
jsw = 3 - jsw;
}
for (int i = 1; i < imax + 1; i++) {
P(i, 0) = P(i, 1);
P(i, jmaxLocal + 1) = P(i, jmaxLocal);
}
if(solver->rank==0)
for(int i = 1 ; i < imax+1 ;i++)P(i,0)=P(i,1);
if(solver->rank== solver->size-1)
for(int i = 1; i < imax +1; i++)P(i, jmaxLocal + 1) = P(i, jmaxLocal);
// for (int i = 1; i < imax + 1; i++) {
// P(i, 0) = P(i, 1);
// P(i, jmaxLocal + 1) = P(i, jmaxLocal);
// }
for (int j = 1; j < jmaxLocal + 1; j++) {
P(0, j) = P(1, j);
P(imax + 1, j) = P(imax, j);

2
PoissonSolver/2D-mpi/src/solver.h
View File

@ -25,5 +25,5 @@ extern void debug(Solver*);
extern void initSolver(Solver*, Parameter*, int problem);
extern void getResult(Solver*);
extern void writeResult(Solver*, double*, char*);
extern int solve(Solver*);
extern int solveRB(Solver*);
#endif