forked from moebiusband/NuSiF-Solver
Compare commits
1 Commits
main
...
poisson/2D
Author | SHA1 | Date | |
---|---|---|---|
356c4fe4fa |
@ -35,7 +35,7 @@ int main(int argc, char** argv)
|
|||||||
}
|
}
|
||||||
|
|
||||||
initSolver(&solver, ¶ms, 2);
|
initSolver(&solver, ¶ms, 2);
|
||||||
solve(&solver);
|
solveRB(&solver);
|
||||||
getResult(&solver);
|
getResult(&solver);
|
||||||
|
|
||||||
MPI_Finalize();
|
MPI_Finalize();
|
||||||
|
@ -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);
|
||||||
|
@ -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
|
||||||
|
Loading…
x
Reference in New Issue
Block a user