[petsc-users] How subdomain talks to neighbors?
(Rebecca) Xuefei YUAN
xy2102 at columbia.edu
Wed Mar 10 10:43:23 CST 2010
Dear all,
I have a question about how subdomains "talk" to neighbor subdomains.
For example, I have a 6X5 mesh, the solution at each mesh f[j][i] is given
for (j = jFirst; j <= jLast; j++){
for (i = iFirst; i <= iLast; i++){
f[j][i] = i+j*info.mx;
}
}
i.e.,
0 1 2 3 4 5
6 7 8 9 10 11
12 13 14 15 16 17
18 19 20 21 22 23
24 25 26 27 28 29
When I read this solution (binary file) as an input of another set of
code with np=2, the array gives me
rank=0: xs=0,x=3,ys=0,y=5
0 1 2 3 4 5
3 4 5 6 7 8
6 7 8 9 10 11
9 10 11 12 13 14
12 13 14 X X X
rank=1: xs=3,x=3,ys=0,y=5
X X X 15 16 17
15 16 17 18 19 20
18 19 20 21 22 23
21 22 23 24 25 26
24 25 26 27 28 29
or np=3 gives
rank=0: xs=0,x=2,ys=0,y=5
0 1 2 3 4 5
2 3 4 5 6 7
4 5 6 7 8 9
6 7 8 9 X X
8 9 X X X X
rank=1: xs=2,x=2,ys=0,y=5
X X 10 11 12 13
10 11 12 13 14 15
12 13 14 15 16 17
14 15 16 17 18 19
16 17 18 19 X X
rank=2: xs=4,x=2,ys=0,y=5
X X X X 20 21
X X 20 21 22 23
20 21 22 23 24 25
22 23 24 25 26 27
24 25 26 27 28 29
, where X means this location is not visited. Still I do not know how
subdomains talk to each other. Any references on how to understand
those communication between subdomains?
Thanks a lot!
Rebecca
More information about the petsc-users
mailing list