Lines Matching refs:dd
13 DM_DA *dd = (DM_DA *)da->data; in DMView_DA_3d() local
56 … w %" PetscInt_FMT " s %" PetscInt_FMT "\n", rank, dd->M, dd->N, dd->P, dd->m, dd->n, dd->p, dd->w… in DMView_DA_3d()
76 PetscReal ymin = -1.0, ymax = (PetscReal)dd->N; in DMView_DA_3d()
77 PetscReal xmin = -1.0, xmax = (PetscReal)((dd->M + 2) * dd->P), x, y, ycoord, xcoord; in DMView_DA_3d()
94 for (k = 0; k < dd->P; k++) { in DMView_DA_3d()
96 ymax = (PetscReal)(dd->N - 1); in DMView_DA_3d()
97 …for (xmin = (PetscReal)(k * (dd->M + 1)); xmin < (PetscReal)(dd->M + (k * (dd->M + 1))); xmin++) P… in DMView_DA_3d()
98 xmin = (PetscReal)(k * (dd->M + 1)); in DMView_DA_3d()
99 xmax = xmin + (PetscReal)(dd->M - 1); in DMView_DA_3d()
100 …for (ymin = 0; ymin < (PetscReal)dd->N; ymin++) PetscCall(PetscDrawLine(draw, xmin, ymin, xmax, ym… in DMView_DA_3d()
109 for (k = 0; k < dd->P; k++) { in DMView_DA_3d()
110 if ((k >= dd->zs) && (k < dd->ze)) { in DMView_DA_3d()
112 ymin = dd->ys; in DMView_DA_3d()
113 ymax = dd->ye - 1; in DMView_DA_3d()
114 xmin = dd->xs / dd->w + (dd->M + 1) * k; in DMView_DA_3d()
115 xmax = (dd->xe - 1) / dd->w + (dd->M + 1) * k; in DMView_DA_3d()
122 xmin = dd->xs / dd->w; in DMView_DA_3d()
123 xmax = (dd->xe - 1) / dd->w; in DMView_DA_3d()
127 … PetscCall(PetscDrawString(draw, xmin + (dd->M + 1) * k + .2, ymin + .3, PETSC_DRAW_RED, node)); in DMView_DA_3d()
129 base = (dd->base + (dd->xe - dd->xs) * (dd->ye - dd->ys) * (k - dd->zs)) / dd->w; in DMView_DA_3d()
131 for (x = xmin + (dd->M + 1) * k; x <= xmax + (dd->M + 1) * k; x++) { in DMView_DA_3d()
143 for (k = 0 - dd->s; k < dd->P + dd->s; k++) { in DMView_DA_3d()
145 if ((k >= dd->Zs) && (k < dd->Ze)) { in DMView_DA_3d()
147 base = (dd->Xe - dd->Xs) * (dd->Ye - dd->Ys) * (k - dd->Zs) / dd->w; in DMView_DA_3d()
151 if (k < 0) plane = dd->P + k; in DMView_DA_3d()
152 if (k >= dd->P) plane = k - dd->P; in DMView_DA_3d()
153 ymin = dd->Ys; in DMView_DA_3d()
154 ymax = dd->Ye; in DMView_DA_3d()
155 xmin = (dd->M + 1) * plane * dd->w; in DMView_DA_3d()
156 xmax = (dd->M + 1) * plane * dd->w + dd->M * dd->w; in DMView_DA_3d()
158 for (x = xmin + dd->Xs; x < xmin + dd->Xe; x += dd->w) { in DMView_DA_3d()
162 if (y < 0) ycoord = dd->N + y; in DMView_DA_3d()
163 if (y >= dd->N) ycoord = y - dd->N; in DMView_DA_3d()
167 PetscCall(PetscDrawString(draw, xcoord / dd->w, ycoord, PETSC_DRAW_BLUE, node)); in DMView_DA_3d()
192 DM_DA *dd = (DM_DA *)da->data; in DMSetUp_DA_3D() local
193 const PetscInt M = dd->M; in DMSetUp_DA_3D()
194 const PetscInt N = dd->N; in DMSetUp_DA_3D()
195 const PetscInt P = dd->P; in DMSetUp_DA_3D()
197 const PetscInt dof = dd->w; in DMSetUp_DA_3D()
198 const PetscInt s = dd->s; in DMSetUp_DA_3D()
199 DMBoundaryType bx = dd->bx; in DMSetUp_DA_3D()
200 DMBoundaryType by = dd->by; in DMSetUp_DA_3D()
201 DMBoundaryType bz = dd->bz; in DMSetUp_DA_3D()
202 DMDAStencilType stencil_type = dd->stencil_type; in DMSetUp_DA_3D()
203 PetscInt *lx = dd->lx; in DMSetUp_DA_3D()
204 PetscInt *ly = dd->ly; in DMSetUp_DA_3D()
205 PetscInt *lz = dd->lz; in DMSetUp_DA_3D()
228 PetscCall(PetscMPIIntCast(dd->m, &m)); in DMSetUp_DA_3D()
229 PetscCall(PetscMPIIntCast(dd->n, &n)); in DMSetUp_DA_3D()
230 PetscCall(PetscMPIIntCast(dd->p, &p)); in DMSetUp_DA_3D()
336 PetscCall(PetscMalloc1(m, &dd->lx)); in DMSetUp_DA_3D()
337 lx = dd->lx; in DMSetUp_DA_3D()
346 PetscCall(PetscMalloc1(n, &dd->ly)); in DMSetUp_DA_3D()
347 ly = dd->ly; in DMSetUp_DA_3D()
357 PetscCall(PetscMalloc1(p, &dd->lz)); in DMSetUp_DA_3D()
358 lz = dd->lz; in DMSetUp_DA_3D()
469 dd->Nlocal = x * y * z * dof; in DMSetUp_DA_3D()
470 PetscCall(VecCreateMPIWithArray(comm, dof, dd->Nlocal, PETSC_DECIDE, NULL, &global)); in DMSetUp_DA_3D()
471 dd->nlocal = (Xe - Xs) * (Ye - Ys) * (Ze - Zs) * dof; in DMSetUp_DA_3D()
472 PetscCall(VecCreateSeqWithArray(PETSC_COMM_SELF, dof, dd->nlocal, NULL, &local)); in DMSetUp_DA_3D()
752 PetscCall(PetscMalloc1(27, &dd->neighbors)); in DMSetUp_DA_3D()
754 dd->neighbors[0] = n0; in DMSetUp_DA_3D()
755 dd->neighbors[1] = n1; in DMSetUp_DA_3D()
756 dd->neighbors[2] = n2; in DMSetUp_DA_3D()
757 dd->neighbors[3] = n3; in DMSetUp_DA_3D()
758 dd->neighbors[4] = n4; in DMSetUp_DA_3D()
759 dd->neighbors[5] = n5; in DMSetUp_DA_3D()
760 dd->neighbors[6] = n6; in DMSetUp_DA_3D()
761 dd->neighbors[7] = n7; in DMSetUp_DA_3D()
762 dd->neighbors[8] = n8; in DMSetUp_DA_3D()
763 dd->neighbors[9] = n9; in DMSetUp_DA_3D()
764 dd->neighbors[10] = n10; in DMSetUp_DA_3D()
765 dd->neighbors[11] = n11; in DMSetUp_DA_3D()
766 dd->neighbors[12] = n12; in DMSetUp_DA_3D()
767 dd->neighbors[13] = rank; in DMSetUp_DA_3D()
768 dd->neighbors[14] = n14; in DMSetUp_DA_3D()
769 dd->neighbors[15] = n15; in DMSetUp_DA_3D()
770 dd->neighbors[16] = n16; in DMSetUp_DA_3D()
771 dd->neighbors[17] = n17; in DMSetUp_DA_3D()
772 dd->neighbors[18] = n18; in DMSetUp_DA_3D()
773 dd->neighbors[19] = n19; in DMSetUp_DA_3D()
774 dd->neighbors[20] = n20; in DMSetUp_DA_3D()
775 dd->neighbors[21] = n21; in DMSetUp_DA_3D()
776 dd->neighbors[22] = n22; in DMSetUp_DA_3D()
777 dd->neighbors[23] = n23; in DMSetUp_DA_3D()
778 dd->neighbors[24] = n24; in DMSetUp_DA_3D()
779 dd->neighbors[25] = n25; in DMSetUp_DA_3D()
780 dd->neighbors[26] = n26; in DMSetUp_DA_3D()
1412 dd->m = m; in DMSetUp_DA_3D()
1413 dd->n = n; in DMSetUp_DA_3D()
1414 dd->p = p; in DMSetUp_DA_3D()
1416 dd->xs = xs * dof; in DMSetUp_DA_3D()
1417 dd->xe = xe * dof; in DMSetUp_DA_3D()
1418 dd->ys = ys; in DMSetUp_DA_3D()
1419 dd->ye = ye; in DMSetUp_DA_3D()
1420 dd->zs = zs; in DMSetUp_DA_3D()
1421 dd->ze = ze; in DMSetUp_DA_3D()
1422 dd->Xs = Xs * dof; in DMSetUp_DA_3D()
1423 dd->Xe = Xe * dof; in DMSetUp_DA_3D()
1424 dd->Ys = Ys; in DMSetUp_DA_3D()
1425 dd->Ye = Ye; in DMSetUp_DA_3D()
1426 dd->Zs = Zs; in DMSetUp_DA_3D()
1427 dd->Ze = Ze; in DMSetUp_DA_3D()
1432 dd->gtol = gtol; in DMSetUp_DA_3D()
1433 dd->base = base; in DMSetUp_DA_3D()
1435 dd->ltol = NULL; in DMSetUp_DA_3D()
1436 dd->ao = NULL; in DMSetUp_DA_3D()