Bug Summary

File:dm/impls/da/dadd.c
Warning:line 75, column 18
Array access (from variable 'lx') results in a null pointer dereference

Annotated Source Code

[?] Use j/k keys for keyboard navigation

1#include <petsc/private/dmdaimpl.h> /*I "petscdmda.h" I*/
2
3/*@
4 DMDACreatePatchIS - Creates an index set corresponding to a patch of the DA.
5
6 Not Collective
7
8 Input Parameters:
9+ da - the DMDA
10. lower - a matstencil with i, j and k corresponding to the lower corner of the patch
11- upper - a matstencil with i, j and k corresponding to the upper corner of the patch
12
13 Output Parameters:
14. is - the IS corresponding to the patch
15
16 Level: developer
17
18.seealso: DMDACreateDomainDecomposition(), DMDACreateDomainDecompositionScatters()
19@*/
20PetscErrorCode DMDACreatePatchIS(DM da,MatStencil *lower,MatStencil *upper,IS *is)
21{
22 PetscInt ms=0,ns=0,ps=0;
23 PetscInt me=1,ne=1,pe=1;
24 PetscInt mr=0,nr=0,pr=0;
25 PetscInt ii,jj,kk;
26 PetscInt si,sj,sk;
27 PetscInt i,j,k,l,idx;
28 PetscInt base;
29 PetscInt xm=1,ym=1,zm=1;
30 const PetscInt *lx,*ly,*lz;
31 PetscInt ox,oy,oz;
32 PetscInt m,n,p,M,N,P,dof;
33 PetscInt nindices;
34 PetscInt *indices;
35 DM_DA *dd = (DM_DA*)da->data;
36 PetscErrorCode ierr;
37
38 PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"; petscstack
->line[petscstack->currentsize] = 38; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);
39 /* need to get the sizes of the actual DM rather than the "global" space of a subdomain DM */
40 M = dd->M;N = dd->N;P=dd->P;
41 m = dd->m;n = dd->n;p=dd->p;
42 dof = dd->w;
43 ierr = DMDAGetOffset(da,&ox,&oy,&oz,NULL((void*)0),NULL((void*)0),NULL((void*)0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),43,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
44 ierr = DMDAGetOwnershipRanges(da,&lx,&ly,&lz);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),44,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
52
←
Value assigned to 'lx'
→
45 nindices = (upper->i - lower->i)*(upper->j - lower->j)*(upper->k - lower->k)*dof;
46 ierr = PetscMalloc1(nindices,&indices)PetscMallocA(1,PETSC_FALSE,46,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,(size_t)(nindices)*sizeof(**(&indices)),(&indices));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),46,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
47 /* start at index 0 on processor 0 */
48 mr = 0;
49 nr = 0;
50 pr = 0;
51 ms = 0;
52 ns = 0;
53 ps = 0;
54 if (lx) me = lx[0];
13
←
Assuming 'lx' is null
→
14
←
Taking false branch
→
26
←
Assuming 'lx' is null
→
27
←
Taking false branch
→
39
←
Assuming 'lx' is null
→
40
←
Taking false branch
→
53
←
Assuming 'lx' is null
→
54
←
Taking false branch
→
55 if (ly) ne = ly[0];
15
←
Assuming 'ly' is null
→
16
←
Taking false branch
→
28
←
Assuming 'ly' is null
→
29
←
Taking false branch
→
41
←
Assuming 'ly' is null
→
42
←
Taking false branch
→
55
←
Assuming 'ly' is null
→
56
←
Taking false branch
→
56 if (lz) pe = lz[0];
17
←
Assuming 'lz' is null
→
18
←
Taking false branch
→
30
←
Assuming 'lz' is null
→
31
←
Taking false branch
→
43
←
Assuming 'lz' is null
→
44
←
Taking false branch
→
57
←
Assuming 'lz' is null
→
58
←
Taking false branch
→
57 idx = 0;
58 for (k=lower->k-oz;k<upper->k-oz;k++) {
19
←
Loop condition is false. Execution continues on line 123
→
32
←
Loop condition is false. Execution continues on line 123
→
45
←
Loop condition is false. Execution continues on line 123
→
59
←
Loop condition is true. Entering loop body
→
59 for (j=lower->j-oy;j < upper->j-oy;j++) {
60
←
Loop condition is true. Entering loop body
→
62
←
Loop condition is true. Entering loop body
→
64
←
Loop condition is true. Entering loop body
→
60 for (i=lower->i-ox;i < upper->i-ox;i++) {
61
←
Loop condition is false. Execution continues on line 59
→
63
←
Loop condition is false. Execution continues on line 59
→
65
←
Loop condition is true. Entering loop body
→
61 /* "actual" indices rather than ones outside of the domain */
62 ii = i;
63 jj = j;
64 kk = k;
65 if (ii < 0) ii = M + ii;
66
←
Assuming 'ii' is >= 0
→
67
←
Taking false branch
→
66 if (jj < 0) jj = N + jj;
68
←
Assuming 'jj' is >= 0
→
69
←
Taking false branch
→
67 if (kk < 0) kk = P + kk;
70
←
Assuming 'kk' is >= 0
→
71
←
Taking false branch
→
68 if (ii > M-1) ii = ii - M;
72
←
Taking false branch
→
69 if (jj > N-1) jj = jj - N;
73
←
Taking false branch
→
70 if (kk > P-1) kk = kk - P;
74
←
Taking false branch
→
71 /* gone out of processor range on x axis */
72 while(ii > me-1 || ii < ms) {
75
←
Assuming the condition is true
→
73 if (mr == m-1) {
76
←
Assuming the condition is true
→
77
←
Taking true branch
→
74 ms = 0;
75 me = lx[0];
78
←
Array access (from variable 'lx') results in a null pointer dereference
76 mr = 0;
77 } else {
78 mr++;
79 ms = me;
80 me += lx[mr];
81 }
82 }
83 /* gone out of processor range on y axis */
84 while(jj > ne-1 || jj < ns) {
85 if (nr == n-1) {
86 ns = 0;
87 ne = ly[0];
88 nr = 0;
89 } else {
90 nr++;
91 ns = ne;
92 ne += ly[nr];
93 }
94 }
95 /* gone out of processor range on z axis */
96 while(kk > pe-1 || kk < ps) {
97 if (pr == p-1) {
98 ps = 0;
99 pe = lz[0];
100 pr = 0;
101 } else {
102 pr++;
103 ps = pe;
104 pe += lz[pr];
105 }
106 }
107 /* compute the vector base on owning processor */
108 xm = me - ms;
109 ym = ne - ns;
110 zm = pe - ps;
111 base = ms*ym*zm + ns*M*zm + ps*M*N;
112 /* compute the local coordinates on owning processor */
113 si = ii - ms;
114 sj = jj - ns;
115 sk = kk - ps;
116 for (l=0;l<dof;l++) {
117 indices[idx] = l + dof*(base + si + xm*sj + xm*ym*sk);
118 idx++;
119 }
120 }
121 }
122 }
123 ISCreateGeneral(PETSC_COMM_SELF((MPI_Comm)0x44000001),idx,indices,PETSC_OWN_POINTER,is);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),123,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
124 PetscFunctionReturn(0)do { do { ; if (petscstack && petscstack->currentsize
> 0) { petscstack->currentsize--; petscstack->function
[petscstack->currentsize] = 0; petscstack->file[petscstack
->currentsize] = 0; petscstack->line[petscstack->currentsize
] = 0; petscstack->petscroutine[petscstack->currentsize
] = PETSC_FALSE; } if (petscstack) { petscstack->hotdepth =
(((petscstack->hotdepth-1)<(0)) ? (0) : (petscstack->
hotdepth-1)); } ; } while (0); return(0);} while (0);
125}
126
127PetscErrorCode DMDASubDomainDA_Private(DM dm, PetscInt *nlocal, DM **sdm)
128{
129 DM *da;
130 PetscInt dim,size,i,j,k,idx;
131 PetscErrorCode ierr;
132 DMDALocalInfo info;
133 PetscInt xsize,ysize,zsize;
134 PetscInt xo,yo,zo;
135 PetscInt xs,ys,zs;
136 PetscInt xm=1,ym=1,zm=1;
137 PetscInt xol,yol,zol;
138 PetscInt m=1,n=1,p=1;
139 PetscInt M,N,P;
140 PetscInt pm,mtmp;
141
142 PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"; petscstack
->line[petscstack->currentsize] = 142; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);
143 ierr = DMDAGetLocalInfo(dm,&info);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),143,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
144 ierr = DMDAGetOverlap(dm,&xol,&yol,&zol);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),144,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
145 ierr = DMDAGetNumLocalSubDomains(dm,&size);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),145,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
146 ierr = PetscMalloc1(size,&da)PetscMallocA(1,PETSC_FALSE,146,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,(size_t)(size)*sizeof(**(&da)),(&da));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),146,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
147
148 if (nlocal) *nlocal = size;
149
150 dim = info.dim;
151
152 M = info.xm;
153 N = info.ym;
154 P = info.zm;
155
156 if (dim == 1) {
157 m = size;
158 } else if (dim == 2) {
159 m = (PetscInt)(0.5 + PetscSqrtReal(((PetscReal)M)*((PetscReal)size)/((PetscReal)N))sqrt(((PetscReal)M)*((PetscReal)size)/((PetscReal)N)));
160 while (m > 0) {
161 n = size/m;
162 if (m*n*p == size) break;
163 m--;
164 }
165 } else if (dim == 3) {
166 n = (PetscInt)(0.5 + PetscPowReal(((PetscReal)N*N)*((PetscReal)size)/((PetscReal)P*M),(PetscReal)(1./3.))pow(((PetscReal)N*N)*((PetscReal)size)/((PetscReal)P*M),(PetscReal
)(1./3.))); if (!n) n = 1;
167 while (n > 0) {
168 pm = size/n;
169 if (n*pm == size) break;
170 n--;
171 }
172 if (!n) n = 1;
173 m = (PetscInt)(0.5 + PetscSqrtReal(((PetscReal)M)*((PetscReal)size)/((PetscReal)P*n))sqrt(((PetscReal)M)*((PetscReal)size)/((PetscReal)P*n)));
174 if (!m) m = 1;
175 while (m > 0) {
176 p = size/(m*n);
177 if (m*n*p == size) break;
178 m--;
179 }
180 if (M > P && m < p) {mtmp = m; m = p; p = mtmp;}
181 }
182
183 zs = info.zs;
184 idx = 0;
185 for (k = 0; k < p; k++) {
186 ys = info.ys;
187 for (j = 0; j < n; j++) {
188 xs = info.xs;
189 for (i = 0; i < m; i++) {
190 if (dim == 1) {
191 xm = M/m + ((M % m) > i);
192 } else if (dim == 2) {
193 xm = M/m + ((M % m) > i);
194 ym = N/n + ((N % n) > j);
195 } else if (dim == 3) {
196 xm = M/m + ((M % m) > i);
197 ym = N/n + ((N % n) > j);
198 zm = P/p + ((P % p) > k);
199 }
200
201 xsize = xm;
202 ysize = ym;
203 zsize = zm;
204 xo = xs;
205 yo = ys;
206 zo = zs;
207
208 ierr = DMDACreate(PETSC_COMM_SELF((MPI_Comm)0x44000001),&(da[idx]));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),208,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
209 ierr = DMSetOptionsPrefix(da[idx],"sub_");CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),209,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
210 ierr = DMSetDimension(da[idx], info.dim);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),210,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
211 ierr = DMDASetDof(da[idx], info.dof);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),211,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
212
213 ierr = DMDASetStencilType(da[idx],info.st);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),213,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
214 ierr = DMDASetStencilWidth(da[idx],info.sw);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),214,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
215
216 if (info.bx == DM_BOUNDARY_PERIODIC || (xs != 0)) {
217 xsize += xol;
218 xo -= xol;
219 }
220 if (info.by == DM_BOUNDARY_PERIODIC || (ys != 0)) {
221 ysize += yol;
222 yo -= yol;
223 }
224 if (info.bz == DM_BOUNDARY_PERIODIC || (zs != 0)) {
225 zsize += zol;
226 zo -= zol;
227 }
228
229 if (info.bx == DM_BOUNDARY_PERIODIC || (xs+xm != info.mx)) xsize += xol;
230 if (info.by == DM_BOUNDARY_PERIODIC || (ys+ym != info.my)) ysize += yol;
231 if (info.bz == DM_BOUNDARY_PERIODIC || (zs+zm != info.mz)) zsize += zol;
232
233 if (info.bx != DM_BOUNDARY_PERIODIC) {
234 if (xo < 0) {
235 xsize += xo;
236 xo = 0;
237 }
238 if (xo+xsize > info.mx-1) {
239 xsize -= xo+xsize - info.mx;
240 }
241 }
242 if (info.by != DM_BOUNDARY_PERIODIC) {
243 if (yo < 0) {
244 ysize += yo;
245 yo = 0;
246 }
247 if (yo+ysize > info.my-1) {
248 ysize -= yo+ysize - info.my;
249 }
250 }
251 if (info.bz != DM_BOUNDARY_PERIODIC) {
252 if (zo < 0) {
253 zsize += zo;
254 zo = 0;
255 }
256 if (zo+zsize > info.mz-1) {
257 zsize -= zo+zsize - info.mz;
258 }
259 }
260
261 ierr = DMDASetSizes(da[idx], xsize, ysize, zsize);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),261,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
262 ierr = DMDASetNumProcs(da[idx], 1, 1, 1);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),262,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
263 ierr = DMDASetBoundaryType(da[idx], DM_BOUNDARY_GHOSTED, DM_BOUNDARY_GHOSTED, DM_BOUNDARY_GHOSTED);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),263,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
264
265 /* set up as a block instead */
266 ierr = DMSetUp(da[idx]);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),266,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
267
268 /* nonoverlapping region */
269 ierr = DMDASetNonOverlappingRegion(da[idx],xs,ys,zs,xm,ym,zm);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),269,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
270
271 /* this alters the behavior of DMDAGetInfo, DMDAGetLocalInfo, DMDAGetCorners, and DMDAGetGhostedCorners and should be used with care */
272 ierr = DMDASetOffset(da[idx],xo,yo,zo,info.mx,info.my,info.mz);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),272,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
273 xs += xm;
274 idx++;
275 }
276 ys += ym;
277 }
278 zs += zm;
279 }
280 *sdm = da;
281 PetscFunctionReturn(0)do { do { ; if (petscstack && petscstack->currentsize
> 0) { petscstack->currentsize--; petscstack->function
[petscstack->currentsize] = 0; petscstack->file[petscstack
->currentsize] = 0; petscstack->line[petscstack->currentsize
] = 0; petscstack->petscroutine[petscstack->currentsize
] = PETSC_FALSE; } if (petscstack) { petscstack->hotdepth =
(((petscstack->hotdepth-1)<(0)) ? (0) : (petscstack->
hotdepth-1)); } ; } while (0); return(0);} while (0);
282}
283
284/*
285 Fills the local vector problem on the subdomain from the global problem.
286
287 Right now this assumes one subdomain per processor.
288
289 */
290PetscErrorCode DMCreateDomainDecompositionScatters_DA(DM dm,PetscInt nsubdms,DM *subdms,VecScatter **iscat,VecScatter **oscat, VecScatter **lscat)
291{
292 PetscErrorCode ierr;
293 DMDALocalInfo info,subinfo;
294 DM subdm;
295 MatStencil upper,lower;
296 IS idis,isis,odis,osis,gdis;
297 Vec svec,dvec,slvec;
298 PetscInt xm,ym,zm,xs,ys,zs;
299 PetscInt i;
300
301 PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"; petscstack
->line[petscstack->currentsize] = 301; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);
302
303 /* allocate the arrays of scatters */
304 if (iscat) {ierr = PetscMalloc1(nsubdms,iscat)PetscMallocA(1,PETSC_FALSE,304,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,(size_t)(nsubdms)*sizeof(**(iscat)),(iscat));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),304,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}
305 if (oscat) {ierr = PetscMalloc1(nsubdms,oscat)PetscMallocA(1,PETSC_FALSE,305,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,(size_t)(nsubdms)*sizeof(**(oscat)),(oscat));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),305,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}
306 if (lscat) {ierr = PetscMalloc1(nsubdms,lscat)PetscMallocA(1,PETSC_FALSE,306,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,(size_t)(nsubdms)*sizeof(**(lscat)),(lscat));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),306,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}
307
308 ierr = DMDAGetLocalInfo(dm,&info);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),308,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
309 for (i = 0; i < nsubdms; i++) {
310 subdm = subdms[i];
311 ierr = DMDAGetLocalInfo(subdm,&subinfo);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),311,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
312 ierr = DMDAGetNonOverlappingRegion(subdm,&xs,&ys,&zs,&xm,&ym,&zm);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),312,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
313
314 /* create the global and subdomain index sets for the inner domain */
315 lower.i = xs;
316 lower.j = ys;
317 lower.k = zs;
318 upper.i = xs+xm;
319 upper.j = ys+ym;
320 upper.k = zs+zm;
321 ierr = DMDACreatePatchIS(dm,&lower,&upper,&idis);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),321,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
322 ierr = DMDACreatePatchIS(subdm,&lower,&upper,&isis);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),322,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
323
324 /* create the global and subdomain index sets for the outer subdomain */
325 lower.i = subinfo.xs;
326 lower.j = subinfo.ys;
327 lower.k = subinfo.zs;
328 upper.i = subinfo.xs+subinfo.xm;
329 upper.j = subinfo.ys+subinfo.ym;
330 upper.k = subinfo.zs+subinfo.zm;
331 ierr = DMDACreatePatchIS(dm,&lower,&upper,&odis);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),331,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
332 ierr = DMDACreatePatchIS(subdm,&lower,&upper,&osis);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),332,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
333
334 /* global and subdomain ISes for the local indices of the subdomain */
335 /* todo - make this not loop over at nonperiodic boundaries, which will be more involved */
336 lower.i = subinfo.gxs;
337 lower.j = subinfo.gys;
338 lower.k = subinfo.gzs;
339 upper.i = subinfo.gxs+subinfo.gxm;
340 upper.j = subinfo.gys+subinfo.gym;
341 upper.k = subinfo.gzs+subinfo.gzm;
342
343 ierr = DMDACreatePatchIS(dm,&lower,&upper,&gdis);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),343,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
344
345 /* form the scatter */
346 ierr = DMGetGlobalVector(dm,&dvec);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),346,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
347 ierr = DMGetGlobalVector(subdm,&svec);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),347,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
348 ierr = DMGetLocalVector(subdm,&slvec);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),348,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
349
350 if (iscat) {ierr = VecScatterCreate(dvec,idis,svec,isis,&(*iscat)[i]);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),350,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}
351 if (oscat) {ierr = VecScatterCreate(dvec,odis,svec,osis,&(*oscat)[i]);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),351,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}
352 if (lscat) {ierr = VecScatterCreate(dvec,gdis,slvec,NULL((void*)0),&(*lscat)[i]);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),352,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}
353
354 ierr = DMRestoreGlobalVector(dm,&dvec);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),354,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
355 ierr = DMRestoreGlobalVector(subdm,&svec);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),355,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
356 ierr = DMRestoreLocalVector(subdm,&slvec);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),356,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
357
358 ierr = ISDestroy(&idis);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),358,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
359 ierr = ISDestroy(&isis);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),359,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
360
361 ierr = ISDestroy(&odis);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),361,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
362 ierr = ISDestroy(&osis);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),362,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
363
364 ierr = ISDestroy(&gdis);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),364,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
365 }
366 PetscFunctionReturn(0)do { do { ; if (petscstack && petscstack->currentsize
> 0) { petscstack->currentsize--; petscstack->function
[petscstack->currentsize] = 0; petscstack->file[petscstack
->currentsize] = 0; petscstack->line[petscstack->currentsize
] = 0; petscstack->petscroutine[petscstack->currentsize
] = PETSC_FALSE; } if (petscstack) { petscstack->hotdepth =
(((petscstack->hotdepth-1)<(0)) ? (0) : (petscstack->
hotdepth-1)); } ; } while (0); return(0);} while (0);
367}
368
369PetscErrorCode DMDASubDomainIS_Private(DM dm,PetscInt n,DM *subdm,IS **iis,IS **ois)
370{
371 PetscErrorCode ierr;
372 PetscInt i;
373 DMDALocalInfo info,subinfo;
374 MatStencil lower,upper;
375
376 PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"; petscstack
->line[petscstack->currentsize] = 376; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);
377 ierr = DMDAGetLocalInfo(dm,&info);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),377,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
378 if (iis) {ierr = PetscMalloc1(n,iis)PetscMallocA(1,PETSC_FALSE,378,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,(size_t)(n)*sizeof(**(iis)),(iis));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),378,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}
4
←
Assuming 'iis' is null
→
5
←
Taking false branch
→
379 if (ois) {ierr = PetscMalloc1(n,ois)PetscMallocA(1,PETSC_FALSE,379,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,(size_t)(n)*sizeof(**(ois)),(ois));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),379,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}
6
←
Assuming 'ois' is non-null
→
7
←
Taking true branch
→
380
381 for (i = 0;i < n; i++) {
8
←
Assuming 'i' is < 'n'
→
9
←
Loop condition is true. Entering loop body
→
21
←
Assuming 'i' is < 'n'
→
22
←
Loop condition is true. Entering loop body
→
34
←
Assuming 'i' is < 'n'
→
35
←
Loop condition is true. Entering loop body
→
47
←
Assuming 'i' is < 'n'
→
48
←
Loop condition is true. Entering loop body
→
382 ierr = DMDAGetLocalInfo(subdm[i],&subinfo);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),382,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
383 if (iis) {
10
←
Taking false branch
→
23
←
Taking false branch
→
36
←
Taking false branch
→
49
←
Taking false branch
→
384 /* create the inner IS */
385 lower.i = info.xs;
386 lower.j = info.ys;
387 lower.k = info.zs;
388 upper.i = info.xs+info.xm;
389 upper.j = info.ys+info.ym;
390 upper.k = info.zs+info.zm;
391 ierr = DMDACreatePatchIS(dm,&lower,&upper,&(*iis)[i]);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),391,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
392 }
393
394 if (ois) {
11
←
Taking true branch
→
24
←
Taking true branch
→
37
←
Taking true branch
→
50
←
Taking true branch
→
395 /* create the outer IS */
396 lower.i = subinfo.xs;
397 lower.j = subinfo.ys;
398 lower.k = subinfo.zs;
399 upper.i = subinfo.xs+subinfo.xm;
400 upper.j = subinfo.ys+subinfo.ym;
401 upper.k = subinfo.zs+subinfo.zm;
402 ierr = DMDACreatePatchIS(dm,&lower,&upper,&(*ois)[i]);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),402,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
12
←
Calling 'DMDACreatePatchIS'
→
20
←
Returning from 'DMDACreatePatchIS'
→
25
←
Calling 'DMDACreatePatchIS'
→
33
←
Returning from 'DMDACreatePatchIS'
→
38
←
Calling 'DMDACreatePatchIS'
→
46
←
Returning from 'DMDACreatePatchIS'
→
51
←
Calling 'DMDACreatePatchIS'
→
403 }
404 }
405 PetscFunctionReturn(0)do { do { ; if (petscstack && petscstack->currentsize
> 0) { petscstack->currentsize--; petscstack->function
[petscstack->currentsize] = 0; petscstack->file[petscstack
->currentsize] = 0; petscstack->line[petscstack->currentsize
] = 0; petscstack->petscroutine[petscstack->currentsize
] = PETSC_FALSE; } if (petscstack) { petscstack->hotdepth =
(((petscstack->hotdepth-1)<(0)) ? (0) : (petscstack->
hotdepth-1)); } ; } while (0); return(0);} while (0);
406}
407
408PetscErrorCode DMCreateDomainDecomposition_DA(DM dm,PetscInt *len,char ***names,IS **iis,IS **ois,DM **subdm)
409{
410 PetscErrorCode ierr;
411 DM *sdm;
412 PetscInt n,i;
413
414 PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"; petscstack
->line[petscstack->currentsize] = 414; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);
415 ierr = DMDASubDomainDA_Private(dm,&n,&sdm);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),415,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
416 if (names) {
1
Assuming 'names' is null
→
2
←
Taking false branch
→
417 ierr = PetscMalloc1(n,names)PetscMallocA(1,PETSC_FALSE,417,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,(size_t)(n)*sizeof(**(names)),(names));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),417,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
418 for (i=0;i<n;i++) (*names)[i] = 0;
419 }
420 ierr = DMDASubDomainIS_Private(dm,n,sdm,iis,ois);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),420,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
3
←
Calling 'DMDASubDomainIS_Private'
→
421 if (subdm) *subdm = sdm;
422 else {
423 for (i=0;i<n;i++) {
424 ierr = DMDestroy(&sdm[i]);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),424,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dadd.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
425 }
426 }
427 if (len) *len = n;
428 PetscFunctionReturn(0)do { do { ; if (petscstack && petscstack->currentsize
> 0) { petscstack->currentsize--; petscstack->function
[petscstack->currentsize] = 0; petscstack->file[petscstack
->currentsize] = 0; petscstack->line[petscstack->currentsize
] = 0; petscstack->petscroutine[petscstack->currentsize
] = PETSC_FALSE; } if (petscstack) { petscstack->hotdepth =
(((petscstack->hotdepth-1)<(0)) ? (0) : (petscstack->
hotdepth-1)); } ; } while (0); return(0);} while (0);
429}