| /petsc/src/ksp/ksp/guess/impls/pod/ |
| H A D | pod.c | 28 PetscInt st; /* first eigenvector of correlation matrix to be retained */ member
159 …PetscCallBLAS("BLASgemv", BLASgemv_("T", &bN, &bNen, &one, pod->eigv + pod->st * pod->n, &bN, pod-… in KSPGuessFormGuess_POD()
196 …PetscCallBLAS("BLASgemv", BLASgemv_("N", &bN, &bNen, &one, pod->eigv + pod->st * pod->n, &bN, pod-… in KSPGuessFormGuess_POD()
321 pod->st = -1; in KSPGuessUpdate_POD()
325 if (!pod->eigs[i]) pod->st = i; in KSPGuessUpdate_POD()
328 for (i = pod->n - 1, parten = 0; i > pod->st && toten > 0; i--) { in KSPGuessUpdate_POD()
333 pod->st = pod->n - pod->nen; in KSPGuessUpdate_POD()
336 for (i = pod->st; i < pod->n; i++) { in KSPGuessUpdate_POD()
338 const PetscInt st = pod->n * i; in KSPGuessUpdate_POD() local
341 for (j = 0; j < pod->n; j++) pod->eigv[st + j] *= v; in KSPGuessUpdate_POD()
[all …]
|
| /petsc/src/ksp/pc/impls/telescope/ |
| H A D | telescope.c | 132 PetscInt m, M, bs, st, ed; in PCTelescopeSetUp_default() local
171 PetscCall(VecGetOwnershipRange(xred, &st, &ed)); in PCTelescopeSetUp_default()
172 PetscCall(ISCreateStride(comm, ed - st, st, 1, &isin)); in PCTelescopeSetUp_default()
174 PetscCall(VecGetOwnershipRange(x, &st, &ed)); in PCTelescopeSetUp_default()
175 PetscCall(ISCreateStride(comm, 0, st, 1, &isin)); in PCTelescopeSetUp_default()
249 PetscInt st, ed; in PCTelescopeSubNullSpaceCreate_Telescope() local
258 PetscCall(VecGetOwnershipRange(sub_vecs[k], &st, &ed)); in PCTelescopeSubNullSpaceCreate_Telescope()
260 for (i = 0; i < ed - st; i++) LA_sub_vec[i] = x_array[i]; in PCTelescopeSubNullSpaceCreate_Telescope()
571 PetscInt i, st, ed; in PCApply_Telescope() local
592 PetscCall(VecGetOwnershipRange(xred, &st, &ed)); in PCApply_Telescope()
[all …]
|
| H A D | telescope_coarsedm.c | 47 PetscInt m, bs, st, ed; in PCTelescopeSetUp_scatters_CoarseDM() local
62 PetscCall(VecGetOwnershipRange(xred, &st, &ed)); in PCTelescopeSetUp_scatters_CoarseDM()
63 PetscCall(ISCreateStride(comm, ed - st, st, 1, &isin)); in PCTelescopeSetUp_scatters_CoarseDM()
66 PetscCall(VecGetOwnershipRange(x, &st, &ed)); in PCTelescopeSetUp_scatters_CoarseDM()
67 PetscCall(ISCreateStride(comm, 0, st, 1, &isin)); in PCTelescopeSetUp_scatters_CoarseDM()
|
| H A D | telescope_dmda.c | 632 PetscInt m, bs, st, ed; in PCTelescopeSetUp_dmda_scatters() local
647 PetscCall(VecGetOwnershipRange(xred, &st, &ed)); in PCTelescopeSetUp_dmda_scatters()
648 PetscCall(ISCreateStride(comm, ed - st, st, 1, &isin)); in PCTelescopeSetUp_dmda_scatters()
651 PetscCall(VecGetOwnershipRange(x, &st, &ed)); in PCTelescopeSetUp_dmda_scatters()
652 PetscCall(ISCreateStride(comm, 0, st, 1, &isin)); in PCTelescopeSetUp_dmda_scatters()
804 PetscInt st, ed; in PCTelescopeSubNullSpaceCreate_dmda_Telescope() local
817 PetscCall(VecGetOwnershipRange(sub_vecs[k], &st, &ed)); in PCTelescopeSubNullSpaceCreate_dmda_Telescope()
819 for (i = 0; i < ed - st; i++) LA_sub_vec[i] = x_array[i]; in PCTelescopeSubNullSpaceCreate_dmda_Telescope()
874 PetscInt i, st, ed; in PCApply_Telescope_dmda() local
902 PetscCall(VecGetOwnershipRange(xred, &st, &ed)); in PCApply_Telescope_dmda()
[all …]
|
| /petsc/src/dm/tests/ |
| H A D | ex7.c | 8 … i, start, end, P = 5, N = 6, m = PETSC_DECIDE, n = PETSC_DECIDE, p = PETSC_DECIDE, pt = 0, st = 0; in main() local
30 PetscCall(PetscOptionsGetInt(NULL, NULL, "-stencil_type", &st, NULL)); in main()
32 stencil_type = (DMDAStencilType)st; in main()
|
| H A D | ex4.c | 20 DMDAStencilType st = DMDA_STENCIL_BOX; in main() local
49 if (flg) st = DMDA_STENCIL_STAR; in main()
52 if (flg) st = DMDA_STENCIL_BOX; in main()
73 PetscCall(DMDACreate2d(PETSC_COMM_WORLD, bx, by, st, M, N, m, n, w, s, lx, ly, &da)); in main()
|
| H A D | ex33.c | 11 PetscInt stencil_width = 1, pt = 0, st = 0; in main() local
43 PetscCall(PetscOptionsGetInt(NULL, NULL, "-stencil_type", &st, NULL)); in main()
44 stencil_type = (DMDAStencilType)st; in main()
|
| H A D | ex23.c | 9 …CIDE, n = PETSC_DECIDE, p = PETSC_DECIDE, M = 8, dof = 1, stencil_width = 1, P = 5, pt = 0, st = 0; in main() local
37 PetscCall(PetscOptionsGetInt(NULL, NULL, "-stencil_type", &st, NULL)); in main()
38 stencil_type = (DMDAStencilType)st; in main()
|
| H A D | ex10k.kokkos.cxx | 33 DMDAStencilType st = DMDA_STENCIL_STAR; in main() local
49 …PetscCall(DMDACreate2d(PETSC_COMM_WORLD, bx, by, st, M, N, PETSC_DECIDE, PETSC_DECIDE, dof, sw, NU… in main()
122 …PetscCall(DMDACreate2d(PETSC_COMM_WORLD, bx, by, st, M, N, PETSC_DECIDE, PETSC_DECIDE, sizeof(Node… in main()
|
| /petsc/src/mat/tests/ |
| H A D | ex195.c | 73 PetscInt st; in main() local
76 PetscCall(MatGetOwnershipRange(A1, &st, NULL)); in main()
77 PetscCall(ISCreateStride(comm, m, st, 2, &is[0])); in main()
78 PetscCall(ISCreateStride(comm, m, st + 1, 2, &is[1])); in main()
|
| /petsc/src/mat/impls/h2opus/cuda/ |
| H A D | math2opusutils.cu | 95 PetscInt st, en; in VecSetDelta() local
98 PetscCall(VecGetOwnershipRange(x, &st, &en)); in VecSetDelta()
105 StandardBasis_Functor delta(ax, i - st); in VecSetDelta()
106 …vice, thrust::counting_iterator<PetscInt>(0), thrust::counting_iterator<PetscInt>(en - st), delta); in VecSetDelta()
112 if (st <= i && i < en) PetscCall(VecSetValue(x, i, 1.0, INSERT_VALUES)); in VecSetDelta()
|
| /petsc/src/dm/impls/da/ |
| H A D | fdda.c | 265 DMDAStencilType st; in DMCreateColoring_DA_2d_MPIAIJ() local
275 PetscCall(DMDAGetInfo(da, &dim, &m, &n, NULL, &M, &N, NULL, &nc, &s, &bx, &by, NULL, &st)); in DMCreateColoring_DA_2d_MPIAIJ()
282 if (st == DMDA_STENCIL_STAR && s == 1) { in DMCreateColoring_DA_2d_MPIAIJ()
329 DMDAStencilType st; in DMCreateColoring_DA_3d_MPIAIJ() local
338 PetscCall(DMDAGetInfo(da, &dim, &m, &n, &p, &M, &N, &P, &nc, &s, &bx, &by, &bz, &st)); in DMCreateColoring_DA_3d_MPIAIJ()
830 DMDAStencilType st; in DMCreateMatrix_DA_2d_MPISELL() local
837 PetscCall(DMDAGetInfo(da, &dim, &m, &n, NULL, NULL, NULL, NULL, &nc, &s, &bx, &by, NULL, &st)); in DMCreateMatrix_DA_2d_MPISELL()
863 … if (st == DMDA_STENCIL_BOX || !l || !p) { /* entries on star have either l = 0 or p = 0 */ in DMCreateMatrix_DA_2d_MPISELL()
901 … if (st == DMDA_STENCIL_BOX || !l || !p) { /* entries on star have either l = 0 or p = 0 */ in DMCreateMatrix_DA_2d_MPISELL()
932 DMDAStencilType st; in DMCreateMatrix_DA_3d_MPISELL() local
[all …]
|
| H A D | daview.c | 108 …eOp DMBoundaryType *bx, PeOp DMBoundaryType *by, PeOp DMBoundaryType *bz, PeOp DMDAStencilType *st) in DMDAGetInfo() argument
135 if (st) *st = dd->stencil_type; in DMDAGetInfo()
178 info->st = dd->stencil_type; in DMDAGetLocalInfo()
|
| /petsc/src/sys/tests/output/ |
| H A D | ex7_1.out | 1 1st test
|
| H A D | ex35_1.out | 1 1st test
|
| /petsc/src/mat/tests/output/ |
| H A D | ex116_2.out | 4 1st 10 of 36 singular values:
|
| H A D | ex116_1.out | 4 1st 10 of 36 singular values:
|
| /petsc/include/ |
| H A D | petscdmdatypes.h | 67 DMDAStencilType st; member
|
| /petsc/src/dm/partitioner/impls/simple/ |
| H A D | partsimple.c | 165 PetscInt st, j; in PetscPartitionerPartition_Simple() local
172 for (np = 0, st = 0; np < nparts; ++np) { in PetscPartitionerPartition_Simple()
174 for (j = 0; j < tpwgts[np]; j++) PetscCall(VecSetValue(v, st + j, np, INSERT_VALUES)); in PetscPartitionerPartition_Simple()
176 st += tpwgts[np]; in PetscPartitionerPartition_Simple()
|
| /petsc/src/ksp/ksp/tutorials/ |
| H A D | ex73.c | 282 PetscInt m, bs, st, ed; in PCTelescopeSetUp_dmda_scatters() local
311 PetscCall(VecGetOwnershipRange(xred, &st, &ed)); in PCTelescopeSetUp_dmda_scatters()
312 PetscCall(ISCreateStride(comm, ed - st, st, 1, &isin)); in PCTelescopeSetUp_dmda_scatters()
315 PetscCall(VecGetOwnershipRange(x, &st, &ed)); in PCTelescopeSetUp_dmda_scatters()
316 PetscCall(ISCreateStride(comm, 0, st, 1, &isin)); in PCTelescopeSetUp_dmda_scatters()
438 PetscInt i, st, ed; in DMShellDAFieldScatter_Forward() local
460 PetscCall(VecGetOwnershipRange(xc, &st, &ed)); in DMShellDAFieldScatter_Forward()
463 for (i = 0; i < ed - st; i++) LA_xred[i] = x_array[i]; in DMShellDAFieldScatter_Forward()
476 PetscInt i, st, ed; in DMShellDAFieldScatter_Reverse() local
493 PetscCall(VecGetOwnershipRange(yc, &st, &ed)); in DMShellDAFieldScatter_Reverse()
[all …]
|
| /petsc/src/snes/tutorials/ |
| H A D | ex30.c | 310 PetscScalar st, ct, th, c = param->c, d = param->d; in HorizVelocity() local
316 st = z / r; in HorizVelocity()
319 return ct * (c * th * st + d * (st + th * ct)) + st * (c * (st - th * ct) + d * th * st); in HorizVelocity()
327 PetscScalar st, ct, th, c = param->c, d = param->d; in VertVelocity() local
333 st = z / r; in VertVelocity()
336 return st * (c * th * st + d * (st + th * ct)) - ct * (c * (st - th * ct) + d * th * st); in VertVelocity()
344 PetscScalar x, z, r, st, ct, c = param->c, d = param->d; in Pressure() local
349 st = z / r; in Pressure()
351 return -2.0 * (c * ct - d * st) / r; in Pressure()
|
| /petsc/src/dm/impls/swarm/ |
| H A D | data_ex.c | 261 PetscMPIInt st, rt; in _get_tags() local
263 st = (N * r0 + r1 + N * N * counter) % maxtag; in _get_tags()
265 *_st = st; in _get_tags()
333 PetscMPIInt symm_nn = 0, *symm_procs = NULL, r0, n, st, rt, size, *maxtag, iflg; in DMSwarmDataExTopologyFinalize() local
362 _get_tags(d->instance, size, r0, r1, *maxtag, &st, &rt); in DMSwarmDataExTopologyFinalize()
363 d->send_tags[n] = st; in DMSwarmDataExTopologyFinalize()
|
| /petsc/src/ts/tests/ |
| H A D | ex21.c | 141 DMDAStencilType st = DMDA_STENCIL_STAR; in main() local
163 PetscCall(DMDACreate2d(PETSC_COMM_WORLD, bt, bt, st, N, N, n, n, 1, sw, NULL, NULL, &da)); in main()
|
| /petsc/src/dm/impls/plex/transform/impls/extrude/ |
| H A D | plextrextrude.c | 539 for (PetscInt st = 0; st < starSize * 2; st += 2) { in DMPlexTransformSetUp_Extrude() local
540 const PetscInt face = star[st]; in DMPlexTransformSetUp_Extrude()
831 for (PetscInt st = 0; st < starSize * 2; st += 2) { in DMPlexTransformMapCoordinates_Extrude() local
832 if ((star[st] >= cStart) && (star[st] < cEnd)) { in DMPlexTransformMapCoordinates_Extrude()
835 PetscCall(DMPlexComputeCellGeometryFVM(dm, star[st], NULL, NULL, cnormal)); in DMPlexTransformMapCoordinates_Extrude()
|
| /petsc/src/vec/is/is/impls/general/ |
| H A D | general.c | 495 PetscInt st = 0; in ISView_General() local
497 if (fmt == PETSC_VIEWER_ASCII_INDEX) st = is->map->rstart; in ISView_General()
499 …CIISynchronizedPrintf(viewer, "[%d] %" PetscInt_FMT " %" PetscInt_FMT "\n", rank, i + st, idx[i])); in ISView_General()
510 PetscInt st = 0; in ISView_General() local
512 if (fmt == PETSC_VIEWER_ASCII_INDEX) st = is->map->rstart; in ISView_General()
514 …tscViewerASCIISynchronizedPrintf(viewer, "%" PetscInt_FMT " %" PetscInt_FMT "\n", i + st, idx[i])); in ISView_General()
|