15   DM_DA   *dd = (DM_DA *)da->data;  in DMDAGetNatural_Private()  local
18   Nlocal = (dd->xe - dd->xs);  in DMDAGetNatural_Private()
19   if (dim > 1) Nlocal *= (dd->ye - dd->ys);  in DMDAGetNatural_Private()
20   if (dim > 2) Nlocal *= (dd->ze - dd->zs);  in DMDAGetNatural_Private()
25     for (i = dd->xs; i < dd->xe; i++) {  in DMDAGetNatural_Private()
30     for (j = dd->ys; j < dd->ye; j++) {  in DMDAGetNatural_Private()
31       for (i = dd->xs; i < dd->xe; i++) {  in DMDAGetNatural_Private()
33         lidx[lict++] = i + j * dd->M * dd->w;  in DMDAGetNatural_Private()
37     for (k = dd->zs; k < dd->ze; k++) {  in DMDAGetNatural_Private()
38       for (j = dd->ys; j < dd->ye; j++) {  in DMDAGetNatural_Private()
39 …for (i = dd->xs; i < dd->xe; i++) lidx[lict++] = i + j * dd->M * dd->w + k * dd->M * dd->N * dd->w;  in DMDAGetNatural_Private()
68   DM_DA    *dd;  in DMDASetAOType()  local
76   dd = (DM_DA *)da->data;  in DMDASetAOType()
77   if (dd->ao) { /* check if the already computed AO has the same type as requested */  in DMDASetAOType()
79     PetscCall(PetscObjectTypeCompare((PetscObject)dd->ao, aotype, &match));  in DMDASetAOType()
83   PetscCall(PetscFree(dd->aotype));  in DMDASetAOType()
84   PetscCall(PetscStrallocpy(aotype, (char **)&dd->aotype));  in DMDASetAOType()
116   DM_DA    *dd;  in DMDAGetAO()  local
125   dd = (DM_DA *)da->data;  in DMDAGetAO()
130   if (!dd->ao) {  in DMDAGetAO()
135     PetscCall(ISCreateStride(PetscObjectComm((PetscObject)da), Nlocal, dd->base, 1, &ispetsc));  in DMDAGetAO()
136     PetscCall(AOCreate(PetscObjectComm((PetscObject)da), &dd->ao));  in DMDAGetAO()
137     PetscCall(AOSetIS(dd->ao, isnatural, ispetsc));  in DMDAGetAO()
138     PetscCall(AOSetType(dd->ao, dd->aotype));  in DMDAGetAO()
142   *ao = dd->ao;  in DMDAGetAO()