| /petsc/src/dm/impls/stag/tests/ |
| H A D | ex7.c | 10 …tz, nx, ny, nz, i, j, k, d, is, js, ks, dof0, dof1, dof2, dof3, dofTotal, stencilWidth, Nx, Ny, Nz; in main() local 21 stencilWidth = 2; in main() 22 …TSC_DECIDE, PETSC_DECIDE, dof0, dof1, dof2, dof3, DMSTAG_STENCIL_BOX, stencilWidth, NULL, NULL, NU… in main() 27 PetscCall(DMStagGetStencilWidth(dm, &stencilWidth)); in main() 47 for (ks = -stencilWidth; ks <= stencilWidth; ++ks) { in main() 48 for (js = -stencilWidth; js <= stencilWidth; ++js) { in main() 49 …for (is = -stencilWidth; is <= stencilWidth; ++is) a2[k][j][i][d] += a1[k + ks][j + js][i + is][d]; in main() 67 for (d = 0; d < 3; ++d) expected *= (2 * stencilWidth + 1); in main() 75 …PetscCheck(stencilWidth <= 1, PETSC_COMM_WORLD, PETSC_ERR_SUP, "Check implemented assuming stencil… in main() 89 …for (dd = 0; dd < 3; ++dd) expected *= (bnd[dd] ? stencilWidth + 1 + extra[dd] : 2 * stencilWidth … in main() [all …]
|
| H A D | ex9.c | 10 … nx, ny, nz, i, j, k, d, is, js, ks, dof0, dof1, dof2, dof3, dofTotal, stencilWidth, ngx, ngy, ngz; in main() local 21 stencilWidth = 2; in main() 22 …SC_DECIDE, PETSC_DECIDE, dof0, dof1, dof2, dof3, DMSTAG_STENCIL_STAR, stencilWidth, NULL, NULL, NU… in main() 27 PetscCall(DMStagGetStencilWidth(dm, &stencilWidth)); in main() 47 for (ks = -stencilWidth; ks <= stencilWidth; ++ks) a2[k][j][i][d] += a1[k + ks][j][i][d]; in main() 48 for (js = -stencilWidth; js <= stencilWidth; ++js) a2[k][j][i][d] += a1[k][j + js][i][d]; in main() 49 for (is = -stencilWidth; is <= stencilWidth; ++is) a2[k][j][i][d] += a1[k][j][i + is][d]; in main() 65 …expected = (ngx * ngy * ngz - 8 * stencilWidth * stencilWidth * stencilWidth - 4 * stencilWidth * … in main() 70 expected = 1 + 6 * stencilWidth; in main()
|
| H A D | ex10.c | 10 …PetscInt startx, starty, nx, ny, i, j, d, is, js, dof0, dof1, dof2, dofTotal, stencilWidth, … in main() local 20 stencilWidth = 2; in main() 21 … 4, PETSC_DECIDE, PETSC_DECIDE, dof0, dof1, dof2, DMSTAG_STENCIL_BOX, stencilWidth, NULL, NULL, &d… in main() 26 PetscCall(DMStagGetStencilWidth(dm, &stencilWidth)); in main() 45 for (js = -stencilWidth; js <= stencilWidth; ++js) { in main() 46 for (is = -stencilWidth; is <= stencilWidth; ++is) a2[j][i][d] += a1[j + js][i + is][d]; in main() 62 for (d = 0; d < 2; ++d) expected *= (2 * stencilWidth + 1); in main() 70 …PetscCheck(stencilWidth <= 1, PETSC_COMM_WORLD, PETSC_ERR_SUP, "Non-periodic check implemented ass… in main() 81 …for (dd = 0; dd < 2; ++dd) expected *= (bnd[dd] ? stencilWidth + 1 + extra[dd] : 2 * stencilWidth … in main() 89 PetscScalar expected = (bnd[1] ? stencilWidth + 1 + extra[1] : 2 * stencilWidth + 1); in main() [all …]
|
| H A D | ex12.c | 10 …PetscInt startx, starty, nx, ny, i, j, d, is, js, dof0, dof1, dof2, dofTotal, stencilWidth, … in main() local 20 stencilWidth = 2; in main() 21 …4, PETSC_DECIDE, PETSC_DECIDE, dof0, dof1, dof2, DMSTAG_STENCIL_STAR, stencilWidth, NULL, NULL, &d… in main() 26 PetscCall(DMStagGetStencilWidth(dm, &stencilWidth)); in main() 45 for (js = -stencilWidth; js <= stencilWidth; ++js) a2[j][i][d] += a1[j + js][i][d]; in main() 46 for (is = -stencilWidth; is <= stencilWidth; ++is) a2[j][i][d] += a1[j][i + is][d]; in main() 61 expected = (ngx * ngy - 4 * stencilWidth * stencilWidth) * dofTotal; in main() 66 expected = 1 + 4 * stencilWidth; in main()
|
| H A D | ex18.c | 106 const PetscInt stencilWidth = 1; in main() local 107 … 8, PETSC_DECIDE, PETSC_DECIDE, dof0, dof1, dof2, DMSTAG_STENCIL_BOX, stencilWidth, NULL, NULL, &d… in main()
|
| H A D | ex30.c | 48 const PetscInt stencilWidth = 1; in main() local 49 …TSC_DECIDE, PETSC_DECIDE, dof0, dof1, dof2, dof3, DMSTAG_STENCIL_BOX, stencilWidth, NULL, NULL, NU… in main()
|
| /petsc/src/dm/impls/stag/ |
| H A D | stagintern.c | 47 stag->stencilWidth, stag->l[0], stag->l[1], stag->l[2], *newdm)); in DMStagDuplicateWithoutSetup() 63 …PetscInt dof2, PetscInt dof3, DMStagStencilType stencilType, PetscInt stencilWidth, const PetscInt… in DMStagInitialize() argument 71 PetscCall(DMStagSetStencilWidth(dm, stencilWidth)); in DMStagInitialize()
|
| H A D | stag1d.c | 39 …PetscInt dof0, PetscInt dof1, DMStagStencilType stencilType, PetscInt stencilWidth, const PetscInt… in DMStagCreate1d() argument 47 …DM_BOUNDARY_NONE, M, 0, 0, size, 0, 0, dof0, dof1, 0, 0, stencilType, stencilWidth, lx, NULL, NULL… in DMStagCreate1d() 196 … stag->startGhost[0] = stag->firstRank[0] ? stag->start[0] : stag->start[0] - stag->stencilWidth; in DMSetUp_Stag_1d() 198 stag->nGhost[0] += stag->firstRank[0] ? 0 : stag->stencilWidth; in DMSetUp_Stag_1d() 199 stag->nGhost[0] += stag->lastRank[0] ? 1 : stag->stencilWidth; in DMSetUp_Stag_1d() 213 stag->startGhost[0] = stag->start[0] - stag->stencilWidth; /* This value may be negative */ in DMSetUp_Stag_1d() 214 …stag->nGhost[0] = stag->n[0] + 2 * stag->stencilWidth + (stag->lastRank[0] && stag->stencilWid… in DMSetUp_Stag_1d() 228 stag->startGhost[0] = stag->start[0] - stag->stencilWidth; /* This value may be negative */ in DMSetUp_Stag_1d() 229 stag->nGhost[0] = stag->n[0] + 2 * stag->stencilWidth; in DMSetUp_Stag_1d() 276 …stencilWidth, PETSC_COMM_SELF, PETSC_ERR_SUP, "DMStag 1d setup does not support local sizes (%" Pe… in DMSetUp_Stag_1d() [all …]
|
| H A D | stagda.c | 10 PetscInt dim, i, j, stencilWidth, dof, N[DMSTAG_MAX_DIM]; in DMStagCreateCompatibleDMDA() local 102 stencilWidth = stag->stencilWidth; in DMStagCreateCompatibleDMDA() 106 stencilWidth = stag->stencilWidth; in DMStagCreateCompatibleDMDA() 115 …e1d(PetscObjectComm((PetscObject)dm), stag->boundaryType[0], N[0], dof, stencilWidth, l[0], dmda)); in DMStagCreateCompatibleDMDA() 118 …e[1], stencilType, N[0], N[1], stag->nRanks[0], stag->nRanks[1], dof, stencilWidth, l[0], l[1], dm… in DMStagCreateCompatibleDMDA() 121 …, N[1], N[2], stag->nRanks[0], stag->nRanks[1], stag->nRanks[2], dof, stencilWidth, l[0], l[1], l[… in DMStagCreateCompatibleDMDA()
|
| H A D | stagutils.c | 636 PetscErrorCode DMStagGetStencilWidth(DM dm, PetscInt *stencilWidth) in DMStagGetStencilWidth() argument 642 if (stencilWidth) *stencilWidth = stag->stencilWidth; in DMStagGetStencilWidth() 1242 PetscErrorCode DMStagSetStencilWidth(DM dm, PetscInt stencilWidth) in DMStagSetStencilWidth() argument 1248 PetscValidLogicalCollectiveInt(dm, stencilWidth, 2); in DMStagSetStencilWidth() 1250 …PetscCheck(stencilWidth >= 0, PetscObjectComm((PetscObject)dm), PETSC_ERR_ARG_OUTOFRANGE, "Stencil… in DMStagSetStencilWidth() 1251 stag->stencilWidth = stencilWidth; in DMStagSetStencilWidth() 1542 …ag->boundaryType[d], stag->N[d], dof0, dof1, stag->stencilType, stag->stencilWidth, stag->l[d], &s… in DMStagSetUniformCoordinatesProduct()
|
| H A D | stag2d.c | 48 …PetscInt dof1, PetscInt dof2, DMStagStencilType stencilType, PetscInt stencilWidth, const PetscInt… in DMStagCreate2d() argument 53 …DM_BOUNDARY_NONE, M, N, 0, m, n, 0, dof0, dof1, dof2, 0, stencilType, stencilWidth, lx, ly, NULL, … in DMStagCreate2d() 285 …lType == DMSTAG_STENCIL_NONE || !(stag->n[0] < stag->stencilWidth || stag->n[1] < stag->stencilWid… in DMSetUp_Stag_2d() 286 stag->n[0], stag->n[1], stag->stencilWidth); in DMSetUp_Stag_2d() 302 stag->nGhost[d] += stag->stencilWidth; /* add interior ghost elements */ in DMSetUp_Stag_2d() 303 stag->startGhost[d] -= stag->stencilWidth; in DMSetUp_Stag_2d() 306 stag->nGhost[d] += stag->stencilWidth; /* add interior ghost elements */ in DMSetUp_Stag_2d() 323 stag->startGhost[d] = stag->start[d] - stag->stencilWidth; /* This value may be negative */ in DMSetUp_Stag_2d() 324 …stag->nGhost[d] = stag->n[d] + 2 * stag->stencilWidth + (stag->lastRank[d] && stag->stencilWid… in DMSetUp_Stag_2d() 338 stag->nGhost[d] = stag->n[d] + 2 * stag->stencilWidth; in DMSetUp_Stag_2d() [all …]
|
| H A D | stag.c | 636 if (stag->stencilWidth != stag2->stencilWidth) { in DMGetCompatibility_Stag() 637 …tencil widths: %" PetscInt_FMT " != %" PetscInt_FMT "\n", stag->stencilWidth, stag->stencilWidth)); in DMGetCompatibility_Stag() 840 PetscCall(PetscViewerASCIIPrintf(viewer, ", width %" PetscInt_FMT "\n", stag->stencilWidth)); in DMView_Stag() 894 …"Elementwise stencil width", "DMStagSetStencilWidth", stag->stencilWidth, &stag->stencilWidth, NUL… in DMSetFromOptions_Stag() 968 stag->stencilWidth = 0; in DMCreate_Stag()
|
| H A D | stag3d.c | 55 …PetscInt dof2, PetscInt dof3, DMStagStencilType stencilType, PetscInt stencilWidth, const PetscInt… in DMStagCreate3d() argument 60 …x, bndy, bndz, M, N, P, m, n, p, dof0, dof1, dof2, dof3, stencilType, stencilWidth, lx, ly, lz, *d… in DMStagCreate3d() 391 stag->nGhost[d] += stag->stencilWidth; /* add interior ghost elements */ in DMSetUp_Stag_3d() 392 stag->startGhost[d] -= stag->stencilWidth; in DMSetUp_Stag_3d() 395 stag->nGhost[d] += stag->stencilWidth; /* add interior ghost elements */ in DMSetUp_Stag_3d() 412 stag->startGhost[d] = stag->start[d] - stag->stencilWidth; /* This value may be negative */ in DMSetUp_Stag_3d() 413 …stag->nGhost[d] = stag->n[d] + 2 * stag->stencilWidth + (stag->lastRank[d] && stag->stencilWid… in DMSetUp_Stag_3d() 427 stag->nGhost[d] = stag->n[d] + 2 * stag->stencilWidth; in DMSetUp_Stag_3d() 428 stag->startGhost[d] = stag->start[d] - stag->stencilWidth; in DMSetUp_Stag_3d() 972 …TENCIL_NONE && (stag->n[0] < stag->stencilWidth || stag->n[1] < stag->stencilWidth || stag->n[2] <… in DMStagSetUpBuildScatter_3d() [all …]
|
| /petsc/include/petsc/private/ |
| H A D | dmstagimpl.h | 20 PetscInt stencilWidth; /* Elementwise ghost width */ member
|
| /petsc/src/dm/impls/stag/tutorials/ |
| H A D | ex2.c | 94 const PetscInt stencilWidth = 1; in main() local 95 … 9, PETSC_DECIDE, PETSC_DECIDE, dof0, dof1, dof2, DMSTAG_STENCIL_BOX, stencilWidth, NULL, NULL, &d… in main()
|
| H A D | ex6.c | 93 const PetscInt stencilWidth = 1; in main() local 100 …00, PETSC_DECIDE, PETSC_DECIDE, dof0, dof1, dof2, DMSTAG_STENCIL_BOX, stencilWidth, NULL, NULL, &c… in main() 107 …TSC_DECIDE, PETSC_DECIDE, dof0, dof1, dof2, dof3, DMSTAG_STENCIL_BOX, stencilWidth, NULL, NULL, NU… in main()
|
| H A D | ex3.c | 123 const PetscInt stencilWidth = 1; in main() local 124 …TSC_DECIDE, PETSC_DECIDE, dof0, dof1, dof2, dof3, DMSTAG_STENCIL_BOX, stencilWidth, NULL, NULL, NU… in main()
|