Compare commits

...

1 Commits

Author SHA1 Message Date
356c4fe4fa fix: corrected bd-> result still not okay 2025-04-28 19:03:44 +02:00
3 changed files with 13 additions and 8 deletions

View File

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

View File

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

View File

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