Lines Matching refs:dd
5 DM_DA *dd = (DM_DA *)da->data; in DMSetFromOptions_DA() local
12 …PetscCheck(dd->M >= 0, PetscObjectComm((PetscObject)da), PETSC_ERR_ARG_OUTOFRANGE, "Dimension must… in DMSetFromOptions_DA()
13 …PetscCheck(dd->N >= 0, PetscObjectComm((PetscObject)da), PETSC_ERR_ARG_OUTOFRANGE, "Dimension must… in DMSetFromOptions_DA()
14 …PetscCheck(dd->P >= 0, PetscObjectComm((PetscObject)da), PETSC_ERR_ARG_OUTOFRANGE, "Dimension must… in DMSetFromOptions_DA()
17 …Int("-da_grid_x", "Number of grid points in x direction", "DMDASetSizes", dd->M, &dd->M, NULL, 1)); in DMSetFromOptions_DA()
18 …Int("-da_grid_y", "Number of grid points in y direction", "DMDASetSizes", dd->N, &dd->N, NULL, 1)); in DMSetFromOptions_DA()
19 …Int("-da_grid_z", "Number of grid points in z direction", "DMDASetSizes", dd->P, &dd->P, NULL, 1)); in DMSetFromOptions_DA()
21 …verlap", "Decomposition overlap in all directions", "DMDASetOverlap", dd->xol, &dd->xol, &flg, 0)); in DMSetFromOptions_DA()
22 if (flg) PetscCall(DMDASetOverlap(da, dd->xol, dd->xol, dd->xol)); in DMSetFromOptions_DA()
23 …overlap_x", "Decomposition overlap in x direction", "DMDASetOverlap", dd->xol, &dd->xol, NULL, 0)); in DMSetFromOptions_DA()
24 …overlap_y", "Decomposition overlap in y direction", "DMDASetOverlap", dd->yol, &dd->yol, NULL, 0)); in DMSetFromOptions_DA()
25 …overlap_z", "Decomposition overlap in z direction", "DMDASetOverlap", dd->zol, &dd->zol, NULL, 0)); in DMSetFromOptions_DA()
27 …nsBoundedInt("-da_local_subdomains", "", "DMDASetNumLocalSubdomains", dd->Nsub, &dd->Nsub, &flg, P… in DMSetFromOptions_DA()
28 if (flg) PetscCall(DMDASetNumLocalSubDomains(da, dd->Nsub)); in DMSetFromOptions_DA()
31 …ocessors_x", "Number of processors in x direction", "DMDASetNumProcs", dd->m, &dd->m, NULL, PETSC_… in DMSetFromOptions_DA()
32 …ocessors_y", "Number of processors in y direction", "DMDASetNumProcs", dd->n, &dd->n, NULL, PETSC_… in DMSetFromOptions_DA()
33 …ocessors_z", "Number of processors in z direction", "DMDASetNumProcs", dd->p, &dd->p, NULL, PETSC_… in DMSetFromOptions_DA()
35 … direction", "DMDASetBoundaryType", DMBoundaryTypes, (PetscEnum)dd->bx, (PetscEnum *)&dd->bx, NULL… in DMSetFromOptions_DA()
36 … direction", "DMDASetBoundaryType", DMBoundaryTypes, (PetscEnum)dd->by, (PetscEnum *)&dd->by, NULL… in DMSetFromOptions_DA()
37 … direction", "DMDASetBoundaryType", DMBoundaryTypes, (PetscEnum)dd->bz, (PetscEnum *)&dd->bz, NULL… in DMSetFromOptions_DA()
41 …e_x", "Refinement ratio in x direction", "DMDASetRefinementFactor", dd->refine_x, &dd->refine_x, N… in DMSetFromOptions_DA()
42 …e_y", "Refinement ratio in y direction", "DMDASetRefinementFactor", dd->refine_y, &dd->refine_y, N… in DMSetFromOptions_DA()
43 …e_z", "Refinement ratio in z direction", "DMDASetRefinementFactor", dd->refine_z, &dd->refine_z, N… in DMSetFromOptions_DA()
44 dd->coarsen_x = dd->refine_x; in DMSetFromOptions_DA()
45 dd->coarsen_y = dd->refine_y; in DMSetFromOptions_DA()
46 dd->coarsen_z = dd->refine_z; in DMSetFromOptions_DA()
58 dd->refine_x = refx[0]; in DMSetFromOptions_DA()
59 dd->refine_x_hier_n = n; in DMSetFromOptions_DA()
60 PetscCall(PetscMalloc1(n, &dd->refine_x_hier)); in DMSetFromOptions_DA()
61 PetscCall(PetscArraycpy(dd->refine_x_hier, refx, n)); in DMSetFromOptions_DA()
67 dd->refine_y = refy[0]; in DMSetFromOptions_DA()
68 dd->refine_y_hier_n = n; in DMSetFromOptions_DA()
69 PetscCall(PetscMalloc1(n, &dd->refine_y_hier)); in DMSetFromOptions_DA()
70 PetscCall(PetscArraycpy(dd->refine_y_hier, refy, n)); in DMSetFromOptions_DA()
77 dd->refine_z = refz[0]; in DMSetFromOptions_DA()
78 dd->refine_z_hier_n = n; in DMSetFromOptions_DA()
79 PetscCall(PetscMalloc1(n, &dd->refine_z_hier)); in DMSetFromOptions_DA()
80 PetscCall(PetscArraycpy(dd->refine_z_hier, refz, n)); in DMSetFromOptions_DA()
88 if (dd->bx == DM_BOUNDARY_PERIODIC || dd->interptype == DMDA_Q0) { in DMSetFromOptions_DA()
89 PetscCall(PetscIntMultError(dd->refine_x, dd->M, &dd->M)); in DMSetFromOptions_DA()
91 PetscCall(PetscIntMultError(dd->refine_x, dd->M - 1, &dd->M)); in DMSetFromOptions_DA()
92 dd->M += 1; in DMSetFromOptions_DA()
94 if (dim > 1 && (dd->by == DM_BOUNDARY_PERIODIC || dd->interptype == DMDA_Q0)) { in DMSetFromOptions_DA()
95 PetscCall(PetscIntMultError(dd->refine_y, dd->N, &dd->N)); in DMSetFromOptions_DA()
97 PetscCall(PetscIntMultError(dd->refine_y, dd->N - 1, &dd->N)); in DMSetFromOptions_DA()
98 dd->N += 1; in DMSetFromOptions_DA()
100 if (dim > 2 && (dd->bz == DM_BOUNDARY_PERIODIC || dd->interptype == DMDA_Q0)) { in DMSetFromOptions_DA()
101 PetscCall(PetscIntMultError(dd->refine_z, dd->P, &dd->P)); in DMSetFromOptions_DA()
103 PetscCall(PetscIntMultError(dd->refine_z, dd->P - 1, &dd->P)); in DMSetFromOptions_DA()
104 dd->P += 1; in DMSetFromOptions_DA()
108 dd->refine_x = refx[da->levelup - da->leveldown]; in DMSetFromOptions_DA()
109 dd->refine_y = refy[da->levelup - da->leveldown]; in DMSetFromOptions_DA()
110 dd->refine_z = refz[da->levelup - da->leveldown]; in DMSetFromOptions_DA()
113 dd->coarsen_x = refx[da->levelup - da->leveldown - 1]; in DMSetFromOptions_DA()
114 dd->coarsen_y = refy[da->levelup - da->leveldown - 1]; in DMSetFromOptions_DA()
115 dd->coarsen_z = refz[da->levelup - da->leveldown - 1]; in DMSetFromOptions_DA()
210 DM_DA *dd = (DM_DA *)dm->data; in DMCreateFieldDecomposition_DA() local
211 PetscInt dof = dd->w; in DMCreateFieldDecomposition_DA()
228 … PetscCheck(dd->fieldname, PETSC_COMM_SELF, PETSC_ERR_SUP, "Currently DMDA must have fieldnames"); in DMCreateFieldDecomposition_DA()
229 for (i = 0; i < dof; i++) PetscCall(PetscStrallocpy(dd->fieldname[i], &(*namelist)[i])); in DMCreateFieldDecomposition_DA()
236 PetscCall(DMDASetSizes(da, dd->M, dd->N, dd->P)); in DMCreateFieldDecomposition_DA()
237 PetscCall(DMDASetNumProcs(da, dd->m, dd->n, dd->p)); in DMCreateFieldDecomposition_DA()
238 PetscCall(DMDASetBoundaryType(da, dd->bx, dd->by, dd->bz)); in DMCreateFieldDecomposition_DA()
240 PetscCall(DMDASetStencilType(da, dd->stencil_type)); in DMCreateFieldDecomposition_DA()
241 PetscCall(DMDASetStencilWidth(da, dd->s)); in DMCreateFieldDecomposition_DA()
335 DM_DA *dd; in DMCreate_DA() local
339 PetscCall(PetscNew(&dd)); in DMCreate_DA()
340 da->data = dd; in DMCreate_DA()
343 dd->interptype = DMDA_Q1; in DMCreate_DA()
344 dd->refine_x = 2; in DMCreate_DA()
345 dd->refine_y = 2; in DMCreate_DA()
346 dd->refine_z = 2; in DMCreate_DA()
347 dd->coarsen_x = 2; in DMCreate_DA()
348 dd->coarsen_y = 2; in DMCreate_DA()
349 dd->coarsen_z = 2; in DMCreate_DA()
350 dd->fieldname = NULL; in DMCreate_DA()
351 dd->nlocal = -1; in DMCreate_DA()
352 dd->Nlocal = -1; in DMCreate_DA()
353 dd->M = -1; in DMCreate_DA()
354 dd->N = -1; in DMCreate_DA()
355 dd->P = -1; in DMCreate_DA()
356 dd->m = -1; in DMCreate_DA()
357 dd->n = -1; in DMCreate_DA()
358 dd->p = -1; in DMCreate_DA()
359 dd->w = -1; in DMCreate_DA()
360 dd->s = -1; in DMCreate_DA()
362 dd->xs = -1; in DMCreate_DA()
363 dd->xe = -1; in DMCreate_DA()
364 dd->ys = -1; in DMCreate_DA()
365 dd->ye = -1; in DMCreate_DA()
366 dd->zs = -1; in DMCreate_DA()
367 dd->ze = -1; in DMCreate_DA()
368 dd->Xs = -1; in DMCreate_DA()
369 dd->Xe = -1; in DMCreate_DA()
370 dd->Ys = -1; in DMCreate_DA()
371 dd->Ye = -1; in DMCreate_DA()
372 dd->Zs = -1; in DMCreate_DA()
373 dd->Ze = -1; in DMCreate_DA()
375 dd->Nsub = 1; in DMCreate_DA()
376 dd->xol = 0; in DMCreate_DA()
377 dd->yol = 0; in DMCreate_DA()
378 dd->zol = 0; in DMCreate_DA()
379 dd->xo = 0; in DMCreate_DA()
380 dd->yo = 0; in DMCreate_DA()
381 dd->zo = 0; in DMCreate_DA()
382 dd->Mo = -1; in DMCreate_DA()
383 dd->No = -1; in DMCreate_DA()
384 dd->Po = -1; in DMCreate_DA()
386 dd->gtol = NULL; in DMCreate_DA()
387 dd->ltol = NULL; in DMCreate_DA()
388 dd->ao = NULL; in DMCreate_DA()
389 PetscCall(PetscStrallocpy(AOBASIC, (char **)&dd->aotype)); in DMCreate_DA()
390 dd->base = -1; in DMCreate_DA()
391 dd->bx = DM_BOUNDARY_NONE; in DMCreate_DA()
392 dd->by = DM_BOUNDARY_NONE; in DMCreate_DA()
393 dd->bz = DM_BOUNDARY_NONE; in DMCreate_DA()
394 dd->stencil_type = DMDA_STENCIL_BOX; in DMCreate_DA()
395 dd->interptype = DMDA_Q1; in DMCreate_DA()
396 dd->lx = NULL; in DMCreate_DA()
397 dd->ly = NULL; in DMCreate_DA()
398 dd->lz = NULL; in DMCreate_DA()
400 dd->elementtype = DMDA_ELEMENT_Q1; in DMCreate_DA()