| /petsc/src/ts/tutorials/ |
| H A D | ex22f_mf.F90 | 53 PetscInt xs, xe, gxs, gxe, dof, gdof 105 PetscCallA(GetLayout(da, mx, xs, xe, gxs, gxe, ierr)) 106 dof = i2*(xe - xs + 1) 147 PetscCallA(GetLayout(da, mx, xs, xe, gxs, gxe, ierr)) 148 dof = i2*(xe - xs + 1) 150 call SaveSolutionToDisk(da, X, gdof, xs, xe) 165 subroutine GetLayout(da, mx, xs, xe, gxs, gxe, ierr) argument 170 PetscInt mx, xs, xe, gxs, gxe 178 xe = xs + xm - 1 182 subroutine FormIFunctionLocal(mx, xs, xe, gxs, gxe, x, xdot, f, a, k, s, ierr) argument [all …]
|
| H A D | ex22f.F90 | 137 subroutine GetLayout(da, mx, xs, xe, gxs, gxe, ierr) argument 142 PetscInt mx, xs, xe, gxs, gxe 150 xe = xs + xm - 1 154 subroutine FormIFunctionLocal(mx, xs, xe, gxs, gxe, x, xdot, f, a, k, s, ierr) argument 156 PetscInt mx, xs, xe, gxs, gxe 157 PetscScalar x(2, xs:xe) 158 PetscScalar xdot(2, xs:xe) 159 PetscScalar f(2, xs:xe) 163 do i = xs, xe 180 PetscInt mx, xs, xe, gxs, gxe [all …]
|
| /petsc/src/dm/impls/da/ |
| H A D | dagetelem.c | 6 PetscInt i, xs, xe, Xs, Xe; in DMDAGetElements_1D() local 14 PetscCall(DMDAGetCorners(dm, &xs, NULL, NULL, &xe, NULL, NULL)); in DMDAGetElements_1D() 16 xe += xs; in DMDAGetElements_1D() 19 da->ne = 1 * (xe - xs - 1); in DMDAGetElements_1D() 21 for (i = xs; i < xe - 1; i++) { in DMDAGetElements_1D() 28 corners[1] = (xe - 1 - Xs); in DMDAGetElements_1D() 40 PetscInt i, xs, xe, Xs, Xe; in DMDAGetElements_2D() local 65 PetscCall(DMDAGetCorners(dm, &xs, &ys, NULL, &xe, &ye, NULL)); in DMDAGetElements_2D() 67 xe += xs; in DMDAGetElements_2D() 73 da->ne = ns * (xe - xs - 1) * (ye - ys - 1); in DMDAGetElements_2D() [all …]
|
| H A D | dageometry.c | 41 PetscInt xs, xe, Xs, Xe; in DMGetLocalBoundingBox_DA() local 52 PetscCall(DMDAGetCorners(da, &xs, &ys, &zs, &xe, &ye, &ze)); in DMGetLocalBoundingBox_DA() 56 xe += xs; in DMGetLocalBoundingBox_DA() 80 c1 = (xe - 2 - Xs + 1); in DMGetLocalBoundingBox_DA() 84 c1 = (xe - 2 - Xs + 1) + (ye - 2 - Ys + 1) * (Xe - Xs); in DMGetLocalBoundingBox_DA() 88 …c1 = (xe - 2 - Xs + 1) + (ye - 2 - Ys + 1) * (Xe - Xs) + (ze - 2 - Zs + 1) * (Xe - Xs) * (Ye - Ys); in DMGetLocalBoundingBox_DA() 105 if (bx == DM_BOUNDARY_PERIODIC && xe == M) { in DMGetLocalBoundingBox_DA() 107 ++xe; in DMGetLocalBoundingBox_DA() 123 ce[0] = xe; in DMGetLocalBoundingBox_DA() 134 PetscInt xs, xe, mxlocal; in private_DMDALocatePointsIS_2D_Regular() local [all …]
|
| H A D | da1.c | 96 xmax = (dd->xe / dd->w) - 1; in DMView_DA_1d() 138 PetscInt i, *idx, nn, left, xs, xe, x, Xs, Xe, start, m, IXs, IXe; in DMSetUp_DA_1D() local 196 xe = xs + x; in DMSetUp_DA_1D() 207 if (xe + sDist <= M) { in DMSetUp_DA_1D() 208 Xe = xe + sDist; in DMSetUp_DA_1D() 209 IXe = xe + sDist; in DMSetUp_DA_1D() 211 if (bx) Xe = xe + sDist; in DMSetUp_DA_1D() 218 Xe = xe + sDist; in DMSetUp_DA_1D() 220 IXe = xe + sDist; in DMSetUp_DA_1D() 249 if ((xe + i) < M) idx[nn++] = xe + i; in DMSetUp_DA_1D() [all …]
|
| H A D | daltol.c | 33 PetscCall(PetscMalloc1(dd->xe - dd->xs, &idx)); in DMLocalToLocalCreate_DA() 34 for (j = 0; j < dd->xe - dd->xs; j++) idx[j] = left + j; in DMLocalToLocalCreate_DA() 39 PetscCall(PetscMalloc1((dd->xe - dd->xs) * (up - down), &idx)); in DMLocalToLocalCreate_DA() 42 for (j = 0; j < dd->xe - dd->xs; j++) idx[count++] = left + i * (dd->Xe - dd->Xs) + j; in DMLocalToLocalCreate_DA() 50 count = (dd->xe - dd->xs) * (top - bottom) * (up - down); in DMLocalToLocalCreate_DA() 55 …for (k = 0; k < dd->xe - dd->xs; k++) idx[count++] = (left + j * (dd->Xe - dd->Xs)) + i * (dd->Xe … in DMLocalToLocalCreate_DA()
|
| H A D | daindex.c | 18 Nlocal = (dd->xe - dd->xs); in DMDAGetNatural_Private() 25 for (i = dd->xs; i < dd->xe; i++) { in DMDAGetNatural_Private() 31 for (i = dd->xs; i < dd->xe; i++) { 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()
|
| H A D | da2.c | 90 xmax = (dd->xe - 1) / dd->w; in DMView_DA_2d() 202 PetscInt xs, xe, ys, ye, x, y, Xs, Xe, Ys, Ye, IXs, IXe, IYs, IYe; in DMSetUp_DA_2D() local 303 xe = xs + x; in DMSetUp_DA_2D() 318 if (xe + s <= M) { in DMSetUp_DA_2D() 319 Xe = xe + s; in DMSetUp_DA_2D() 320 IXe = xe + s; in DMSetUp_DA_2D() 324 Xe = xe + s; in DMSetUp_DA_2D() 333 IXe = xe + s; in DMSetUp_DA_2D() 335 Xe = xe + s; in DMSetUp_DA_2D() 649 } else if (Xe - xe > 0 && ys - Ys > 0) { in DMSetUp_DA_2D() [all …]
|
| H A D | da3.c | 115 xmax = (dd->xe - 1) / dd->w + (dd->M + 1) * k; in DMView_DA_3d() 123 xmax = (dd->xe - 1) / dd->w; 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() 208 PetscInt xs = 0, xe, ys = 0, ye, zs = 0, ze, x = 0, y = 0, z = 0; in DMSetUp_DA_3D() local 376 xe = xs + x; in DMSetUp_DA_3D() 388 if (xe + s <= M) { in DMSetUp_DA_3D() 389 Xe = xe + s; in DMSetUp_DA_3D() 390 IXe = xe + s; in DMSetUp_DA_3D() 394 Xe = xe + s; in DMSetUp_DA_3D() 401 IXe = xe + s; in DMSetUp_DA_3D() [all …]
|
| H A D | daview.c | 184 info->xm = (dd->xe - dd->xs) / w; in DMDAGetLocalInfo()
|
| /petsc/src/binding/petsc4py/demo/python_types/ |
| H A D | mat.py | 46 (xs, xe), (ys, ye) = self.da.getRanges() 57 for i in range(xs, xe): 76 (xs, xe), (ys, ye) = self.da.getRanges() 82 for i in range(xs, xe): 90 (xs, xe), (ys, ye) = self.da.getRanges() 92 for i in range(xs, xe):
|
| /petsc/src/tao/unconstrained/tutorials/ |
| H A D | eptorsion2.c | 159 PetscInt xs, ys, xm, ym, gxm, gym, gxs, gys, xe, ye; in FormInitialGuess() local 168 xe = xs + xm; in FormInitialGuess() 172 for (i = xs; i < xe; i++) { /* for (i=0; i<mx; i++) */ in FormInitialGuess() 200 PetscInt xe, ye, xsm, ysm, xep, yep; in FormFunctionGradient() local 232 xe = xs + xm; in FormFunctionGradient() 238 if (xe == mx) xep = xe + 1; in FormFunctionGradient() 239 else xep = xe; in FormFunctionGradient() 245 for (i = xsm; i < xe; i++) { /* for (i=-1; i<mx; i++) */ in FormFunctionGradient() 322 PetscCall(PetscLogFlops((ye - ysm) * (xe - xsm) * 20 + (xep - xs) * (yep - ys) * 16)); in FormFunctionGradient()
|
| H A D | eptorsion2f.F90 | 55 PetscInt i, j, k, xe, ye 70 xe = xs + xm 74 do i = xs, xe - 1 121 PetscInt xe, ye, xsm, ysm 152 xe = xs + xm 164 if (xe == mx) then 165 xep = xe + 1 167 xep = xe 178 do i = xsm, xe - 1 255 PetscCall(PetscLogFlops(20.0d0*(ye - ysm)*(xe - xsm) + 16.0d0*(xep - xs)*(yep - ys), ierr))
|
| /petsc/src/snes/tutorials/ |
| H A D | ex5f.F90 | 39 PetscInt xs, xe, xm, gxs, gxe, gxm 114 PetscScalar x(xs:xe, ys:ye) 131 do i = xs, xe 165 PetscScalar f(xs:xe, ys:ye) 175 xe = xs + info%XM - 1 192 do i = xs, xe 289 do i = xs, xe 458 xe = xs + xm - 1
|
| H A D | ex5f90t.F90 | 40 PetscInt xs, xe, xm, gxs, gxe, gxm 177 PetscScalar x(ctx%xs:ctx%xe, ctx%ys:ctx%ye) 195 do i = ctx%xs, ctx%xe 225 PetscScalar f(ctx%xs:ctx%xe, ctx%ys:ctx%ye) 244 do i = ctx%xs, ctx%xe 424 do i = ctx%xs, ctx%xe 535 ctx%xe = ctx%xs + ctx%xm - 1
|
| H A D | ex5f90.F90 | 40 PetscInt xs, xe, xm, gxs, gxe, gxm 187 PetscScalar x(ctx%xs:ctx%xe, ctx%ys:ctx%ye) 205 do i = ctx%xs, ctx%xe 235 PetscScalar f(ctx%xs:ctx%xe, ctx%ys:ctx%ye) 254 do i = ctx%xs, ctx%xe 435 do i = ctx%xs, ctx%xe 547 ctx%xe = ctx%xs + ctx%xm - 1
|
| /petsc/src/binding/petsc4py/demo/legacy/poisson2d/ |
| H A D | poisson2d.py | 39 (xs, xe), (ys, ye) = self.da.getRanges() 41 for i in range(xs, xe): 52 (xs, xe), (ys, ye) = self.da.getRanges() 54 for i in range(xs, xe):
|
| /petsc/src/binding/petsc4py/demo/legacy/wrap-f2py/ |
| H A D | Bratu2D.F90 | 28 PetscInt mx, xs, xe, xm, gxs, gxe, gxm 50 grd%xe = grd%xs + grd%xm - 1 59 PetscScalar x(grd%xs:grd%xe, grd%ys:grd%ye) 73 do i = grd%xs, grd%xe 91 PetscScalar f(grd%xs:grd%xe, grd%ys:grd%ye) 108 do i = grd%xs, grd%xe 149 do i = grd%xs, grd%xe
|
| /petsc/src/ksp/pc/impls/kaczmarz/ |
| H A D | kaczmarz.c | 18 PetscInt xs, xe, ys, ye, ncols, i, j; in PCApply_Kaczmarz() local 27 PetscCall(MatGetOwnershipRange(pc->pmat, &xs, &xe)); in PCApply_Kaczmarz() 32 for (i = xs; i < xe; i++) { in PCApply_Kaczmarz() 49 for (i = xe - 1; i >= xs; i--) { in PCApply_Kaczmarz()
|
| /petsc/src/binding/petsc4py/demo/legacy/bratu3d/ |
| H A D | bratu3d.py | 46 (xs, xe), (ys, ye), (zs, ze) = self.da.getRanges() 51 for i in range(xs, xe): 76 (xs, xe), (ys, ye), (zs, ze) = self.da.getRanges() 79 for i in range(xs, xe): 114 (xs, xe), (ys, ye), (zs, ze) = self.da.getRanges() 117 for i in range(xs, xe):
|
| /petsc/src/sys/classes/draw/impls/image/ |
| H A D | drawimage.h | 23 int x, xs = img->clip[0], xe = img->clip[2]; in PetscImageClear() local 26 for (x = xs; x < xe; x++) img->buffer[y * img->w + x] = 0; in PetscImageClear() 77 int xs = PetscMax(x, img->clip[0]), xe = PetscMin(x + w, img->clip[2]); in PetscImageDrawRectangle() local 79 if (xs >= xe || ys >= ye) return; in PetscImageDrawRectangle() 81 for (x = xs; x < xe; x++) img->buffer[y * img->w + x] = (unsigned char)c; in PetscImageDrawRectangle()
|
| /petsc/src/dm/impls/plex/exodusii/ |
| H A D | plexexodusii2.c | 1327 PetscInt xs, xe, bs; in VecViewPlex_ExodusII_Nodal_Internal() local 1347 PetscCall(VecGetOwnershipRange(vNatural, &xs, &xe)); in VecViewPlex_ExodusII_Nodal_Internal() 1351 … PetscCallExternal(ex_put_partial_var, exoid, step, EX_NODAL, offset, 1, xs + 1, xe - xs, varray); in VecViewPlex_ExodusII_Nodal_Internal() 1357 PetscCall(ISCreateStride(comm, (xe - xs) / bs, xs, bs, &compIS)); in VecViewPlex_ExodusII_Nodal_Internal() 1359 PetscCall(ISStrideSetStride(compIS, (xe - xs) / bs, xs + c, bs)); in VecViewPlex_ExodusII_Nodal_Internal() 1362 …PetscCallExternal(ex_put_partial_var, exoid, step, EX_NODAL, offset + c, 1, xs / bs + 1, (xe - xs)… in VecViewPlex_ExodusII_Nodal_Internal() 1379 PetscInt xs, xe, bs; in VecLoadPlex_ExodusII_Nodal_Internal() local 1392 PetscCall(VecGetOwnershipRange(vNatural, &xs, &xe)); in VecLoadPlex_ExodusII_Nodal_Internal() 1396 … PetscCallExternal(ex_get_partial_var, exoid, step, EX_NODAL, offset, 1, xs + 1, xe - xs, varray); in VecLoadPlex_ExodusII_Nodal_Internal() 1402 PetscCall(ISCreateStride(comm, (xe - xs) / bs, xs, bs, &compIS)); in VecLoadPlex_ExodusII_Nodal_Internal() [all …]
|
| /petsc/src/mat/impls/elemental/ |
| H A D | matelem.cxx | 226 El::DistMatrix<PetscElemScalar, El::VC, El::STAR> xe, ye; in MatMult_Elemental() local 227 xe.LockedAttach(A->cmap->N, 1, *a->grid, 0, 0, x, A->cmap->n); in MatMult_Elemental() 229 El::Gemv(El::NORMAL, one, *a->emat, xe, zero, ye); in MatMult_Elemental() 247 El::DistMatrix<PetscElemScalar, El::VC, El::STAR> xe, ye; in MatMultTranspose_Elemental() local 248 xe.LockedAttach(A->rmap->N, 1, *a->grid, 0, 0, x, A->rmap->n); in MatMultTranspose_Elemental() 250 El::Gemv(El::TRANSPOSE, one, *a->emat, xe, zero, ye); in MatMultTranspose_Elemental() 269 El::DistMatrix<PetscElemScalar, El::VC, El::STAR> xe, ze; in MatMultAdd_Elemental() local 270 xe.LockedAttach(A->cmap->N, 1, *a->grid, 0, 0, x, A->cmap->n); in MatMultAdd_Elemental() 272 El::Gemv(El::NORMAL, one, *a->emat, xe, one, ze); in MatMultAdd_Elemental() 291 El::DistMatrix<PetscElemScalar, El::VC, El::STAR> xe, ze; in MatMultTransposeAdd_Elemental() local [all …]
|
| /petsc/src/dm/impls/swarm/tests/ |
| H A D | ex5.c | 455 const PetscReal xe[3] = {r0 * ct, r0 * st, 0.0}; in ComputeError() local 460 e[(p * 2 + 0) * dim + d] = x[d] - xe[d]; in ComputeError() 466 …]), (double)x[0], (double)x[1], (double)v[0], (double)v[1], (double)xe[0], (double)xe[1], (double)… in ComputeError()
|
| /petsc/include/petsc/private/ |
| H A D | dmdaimpl.h | 15 PetscInt xs, xe, ys, ye, zs, ze; /* range of local values */ member
|