| /petsc/src/dm/impls/da/ |
| H A D | dagetarray.c | 45 PetscInt xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, N, dim, dof; in DMDAVecGetArray() local 51 PetscCall(DMDAGetCorners(da, &xs, &ys, &zs, &xm, &ym, &zm)); in DMDAVecGetArray() 63 gzs = zs; in DMDAVecGetArray() 94 PetscInt xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, N, dim, dof; in DMDAVecRestoreArray() local 100 PetscCall(DMDAGetCorners(da, &xs, &ys, &zs, &xm, &ym, &zm)); in DMDAVecRestoreArray() 112 gzs = zs; in DMDAVecRestoreArray() 169 PetscInt xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, N, dim, dof; in DMDAVecGetArrayWrite() local 179 PetscCall(DMDAGetCorners(da, &xs, &ys, &zs, &xm, &ym, &zm)); in DMDAVecGetArrayWrite() 191 gzs = zs; in DMDAVecGetArrayWrite() 221 PetscInt xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, N, dim, dof; in DMDAVecRestoreArrayWrite() local [all …]
|
| H A D | dagetelem.c | 107 PetscInt k, zs, ze, Zs, Ze; in DMDAGetElements_3D() local 131 PetscCall(DMDAGetCorners(dm, &xs, &ys, &zs, &xe, &ye, &ze)); in DMDAGetElements_3D() 139 ze += zs; in DMDAGetElements_3D() 141 if (zs != Zs) zs -= 1; in DMDAGetElements_3D() 142 da->ne = ns * (xe - xs - 1) * (ye - ys - 1) * (ze - zs - 1); in DMDAGetElements_3D() 144 for (k = zs; k < ze - 1; k++) { in DMDAGetElements_3D() 165 corners[0] = (xs - Xs) + (ys - Ys) * (Xe - Xs) + (zs - Zs) * (Xe - Xs) * (Ye - Ys); in DMDAGetElements_3D() 166 corners[1] = (xe - 1 - Xs) + (ys - Ys) * (Xe - Xs) + (zs - Zs) * (Xe - Xs) * (Ye - Ys); in DMDAGetElements_3D() 167 corners[2] = (xs - Xs) + (ye - 1 - Ys) * (Xe - Xs) + (zs - Zs) * (Xe - Xs) * (Ye - Ys); in DMDAGetElements_3D() 168 corners[3] = (xe - 1 - Xs) + (ye - 1 - Ys) * (Xe - Xs) + (zs - Zs) * (Xe - Xs) * (Ye - Ys); in DMDAGetElements_3D() [all …]
|
| H A D | dadd.c | 272 PetscInt xs, ys, zs; in DMDASubDomainDA_Private() local 325 zs = info.zs; in DMDASubDomainDA_Private() 348 zo = zs; in DMDASubDomainDA_Private() 366 if (info.bz == DM_BOUNDARY_PERIODIC || (zs != 0)) { in DMDASubDomainDA_Private() 373 if (info.bz == DM_BOUNDARY_PERIODIC || (zs + zm != info.mz)) zsize += zol; in DMDASubDomainDA_Private() 405 PetscCall(DMDASetNonOverlappingRegion(da[idx], xs, ys, zs, xm, ym, zm)); in DMDASubDomainDA_Private() 414 zs += zm; in DMDASubDomainDA_Private() 433 PetscInt xm, ym, zm, xs, ys, zs; in DMCreateDomainDecompositionScatters_DA() local 447 PetscCall(DMDAGetNonOverlappingRegion(subdm, &xs, &ys, &zs, &xm, &ym, &zm)); in DMCreateDomainDecompositionScatters_DA() 452 lower.k = zs; in DMCreateDomainDecompositionScatters_DA() [all …]
|
| H A D | dageometry.c | 43 PetscInt zs, ze, Zs, Ze; in DMGetLocalBoundingBox_DA() local 52 PetscCall(DMDAGetCorners(da, &xs, &ys, &zs, &xe, &ye, &ze)); in DMGetLocalBoundingBox_DA() 60 ze += zs; in DMGetLocalBoundingBox_DA() 65 if (zs != Zs && Zs >= 0) zs -= 1; in DMGetLocalBoundingBox_DA() 87 c0 = (xs - Xs) + (ys - Ys) * (Xe - Xs) + (zs - Zs) * (Xe - Xs) * (Ye - Ys); in DMGetLocalBoundingBox_DA() 120 if (dim > 2) cs[2] = zs; in DMGetLocalBoundingBox_DA() 198 PetscInt zs, ze, mzlocal; in private_DMDALocatePointsIS_3D_Regular() local 208 zs = cs[2]; in private_DMDALocatePointsIS_3D_Regular() 216 mzlocal = ze - zs - 1; in private_DMDALocatePointsIS_3D_Regular() 252 if (mi[2] < zs) continue; in private_DMDALocatePointsIS_3D_Regular() [all …]
|
| H A D | dalocal.c | 368 ind[2] = k + da->zs; in DMDASetVertexCoordinates() 410 PetscInt j, i, xs, ys, xm, ym, zs, zm; in DMDAGetArray() local 430 zs = dd->Zs; in DMDAGetArray() 447 zs = dd->zs; in DMDAGetArray() 450 zm = dd->ze - dd->zs; in DMDAGetArray() 478 ptr = (void ***)(iarray_start + xm * ym * zm * sizeof(PetscScalar) - zs * sizeof(void *)); in DMDAGetArray() 480 for (i = zs; i < zs + zm; i++) ptr[i] = bptr + ((i - zs) * ym - ys); in DMDAGetArray() 481 for (i = zs; i < zs + zm; i++) { in DMDAGetArray() 482 … ym; j++) ptr[i][j] = iarray_start + sizeof(PetscScalar) * (xm * ym * (i - zs) + xm * (j - ys) - x… in DMDAGetArray()
|
| H A D | grvtk.c | 79 rloc[4] = info.zs; in DMDAVTKWriteAll_VTS() 88 PetscInt xs = -1, xm = -1, ys = -1, ym = -1, zs = -1, zm = -1, nnodes = 0; in DMDAVTKWriteAll_VTS() local 95 zs = grloc[r][4]; in DMDAVTKWriteAll_VTS() 100 …T " %" PetscInt_FMT " %" PetscInt_FMT "\">\n", xs, xs + xm - 1, ys, ys + ym - 1, zs, zs + zm - 1)); in DMDAVTKWriteAll_VTS() 151 PetscInt xs = -1, xm = -1, ys = -1, ym = -1, zs = -1, zm = -1, nnodes = 0; in DMDAVTKWriteAll_VTS() local 158 zs = grloc[r][4]; in DMDAVTKWriteAll_VTS() 200 array2[Iloc * 3 + 2] = zs + k; in DMDAVTKWriteAll_VTS() 293 rloc[4] = info.zs; in DMDAVTKWriteAll_VTR() 302 PetscInt xs = -1, xm = -1, ys = -1, ym = -1, zs = -1, zm = -1, nnodes = 0; in DMDAVTKWriteAll_VTR() local 309 zs = grloc[r][4]; in DMDAVTKWriteAll_VTR() [all …]
|
| H A D | dasub.c | 203 PetscInt flag, *owners, xs, xm, ys, ym, zs, zm; in DMDAGetProcessorSubset() local 210 PetscCall(DMDAGetCorners(da, &xs, &ys, &zs, &xm, &ym, &zm)); in DMDAGetProcessorSubset() 215 if (gp >= zs && gp < zs + zm) flag = 1; in DMDAGetProcessorSubset() 276 PetscInt xs, xm, ys, ym, zs, zm, firstPoint; in DMDAGetProcessorSubsets() local 281 PetscCall(DMDAGetCorners(da, &xs, &ys, &zs, &xm, &ym, &zm)); in DMDAGetProcessorSubsets() 285 firstPoint = zs; in DMDAGetProcessorSubsets()
|
| H A D | da3.c | 58 … info.xs + info.xm, info.ys, info.ys + info.ym, info.zs, info.zs + info.zm)); in DMView_DA_3d() 110 if ((k >= dd->zs) && (k < dd->ze)) { 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 373 zs = 0; in DMSetUp_DA_3D() 374 for (i = 0; i < (rank / (m * n)); i++) zs += lz[i]; in DMSetUp_DA_3D() 377 ze = zs + z; in DMSetUp_DA_3D() 430 if (zs - s > 0) { in DMSetUp_DA_3D() 431 Zs = zs - s; in DMSetUp_DA_3D() 432 IZs = zs - s; in DMSetUp_DA_3D() [all …]
|
| H A D | dacorn.c | 252 if (z) *z = dd->zs + dd->zo; in DMDAGetCorners() 255 if (p) *p = (dd->ze - dd->zs); in DMDAGetCorners() 267 lmin[2] = info.zs; in DMGetLocalBoundingIndices_DMDA() 270 lmax[2] = info.zs + info.zm - 1; in DMGetLocalBoundingIndices_DMDA()
|
| /petsc/src/dm/tests/ |
| H A D | ex2k.kokkos.cxx | 36 PetscInt xm, ym, zm, xs, ys, zs, gxm, gym, gzm, gxs, gys, gzs; in Update1() local 39 PetscCall(DMDAGetCorners(da, &xs, &ys, &zs, &xm, &ym, &zm)); in Update1() 43 for (k = zs; k < zs + zm; k++) { in Update1() 59 PetscInt xm, ym, zm, xs, ys, zs, gxm, gym, gzm, gxs, gys, gzs; in Update2() local 62 PetscCall(DMDAGetCorners(da, &xs, &ys, &zs, &xm, &ym, &zm)); in Update2() 65 #define Y2(k, j, i) y2[(k - zs) * ym * xm + (j - ys) * xm + (i - xs)] in Update2() 68 for (k = zs; k < zs + zm; k++) { in Update2() 84 PetscInt xm, ym, zm, xs, ys, zs, gxm, gym, gzm, gxs, gys, gzs; in main() local 111 PetscCall(DMDAGetCorners(da, &xs, &ys, &zs, &xm, &ym, &zm)); in main() 146 …efaultHostExecutionSpace, Rank<3, Iterate::Right, Iterate::Right>>({zs, ys, xs}, {zs + zm, ys + ym… in main()
|
| H A D | ex47.c | 30 for (k = info.zs; k < info.zs + info.zm; k++) { in test_3d() 149 for (k = info.zs; k < info.zs + info.zm; k++) { in test_3d_nocoord()
|
| H A D | ex42.c | 32 for (k = info.zs; k < info.zs + info.zm; k++) { in test_3d() 151 for (k = info.zs; k < info.zs + info.zm; k++) { in test_3d_nocoord()
|
| /petsc/src/dm/impls/da/ftn-custom/ |
| H A D | zda1f90.c | 68 PetscInt xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, N, dim, dof; in dmdavecgetarray1_() local 71 *ierr = DMDAGetCorners(*da, &xs, &ys, &zs, &xm, &ym, &zm); in dmdavecgetarray1_() 87 gzs = zs; in dmdavecgetarray1_() 109 PetscInt xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, N, dim, dof; in dmdavecgetarray2_() local 112 *ierr = DMDAGetCorners(*da, &xs, &ys, &zs, &xm, &ym, &zm); in dmdavecgetarray2_() 128 gzs = zs; in dmdavecgetarray2_() 156 PetscInt xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, N, dim, dof; in dmdavecgetarray3_() local 159 *ierr = DMDAGetCorners(*da, &xs, &ys, &zs, &xm, &ym, &zm); in dmdavecgetarray3_() 175 gzs = zs; in dmdavecgetarray3_() 205 PetscInt xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, N, dim, dof, zero = 0; in dmdavecgetarray4_() local [all …]
|
| /petsc/src/snes/tutorials/ |
| H A D | ex16.c | 601 PetscInt xs = info->xs, ys = info->ys, zs = info->zs; in FormJacobianLocal() local 620 zes = zs; in FormJacobianLocal() 623 zee = zs + zm; in FormJacobianLocal() 626 if (zs > 0) zes = zs - 1; in FormJacobianLocal() 629 if (zs + zm == mz) zee = zs + zm - 1; in FormJacobianLocal() 646 …= xs && i + ii < xm + xs && j + jj >= ys && j + jj < ys + ym && k + kk >= zs && k + kk < zs + zm) { in FormJacobianLocal() 676 for (k = zs; k < zs + zm; k++) { in FormJacobianLocal() 710 PetscInt xs = info->xs, ys = info->ys, zs = info->zs; in FormFunctionLocal() local 723 PetscCall(DMDAGetCorners(info->da, &xs, &ys, &zs, &xm, &ym, &zm)); in FormFunctionLocal() 726 for (k = zs; k < zs + zm; k++) { in FormFunctionLocal() [all …]
|
| H A D | ex14.c | 204 PetscInt i, j, k, Mx, My, Mz, xs, ys, zs, xm, ym, zm; in FormInitialGuess() local 232 PetscCall(DMDAGetCorners(user->da, &xs, &ys, &zs, &xm, &ym, &zm)); in FormInitialGuess() 237 for (k = zs; k < zs + zm; k++) { in FormInitialGuess() 273 PetscInt i, j, k, Mx, My, Mz, xs, ys, zs, xm, ym, zm; in FormFunctionLocal() local 300 PetscCall(DMDAGetCorners(da, &xs, &ys, &zs, &xm, &ym, &zm)); in FormFunctionLocal() 305 for (k = zs; k < zs + zm; k++) { in FormFunctionLocal() 389 PetscInt xs, ys, zs, xm, ym, zm; in FormJacobian() local 424 PetscCall(DMDAGetCorners(da, &xs, &ys, &zs, &xm, &ym, &zm)); in FormJacobian() 437 for (k = zs; k < zs + zm; k++) { in FormJacobian()
|
| /petsc/src/binding/petsc4py/demo/legacy/wrap-swig/ |
| H A D | Bratu3D.c | 24 PetscInt i,j,k,Mx,My,Mz,xs,ys,zs,xm,ym,zm; in FormInitGuess() local 55 PetscCall(DMDAGetCorners(da,&xs,&ys,&zs,&xm,&ym,&zm)); in FormInitGuess() 60 for (k=zs; k<zs+zm; k++) { in FormInitGuess() 84 PetscInt i,j,k,Mx,My,Mz,xs,ys,zs,xm,ym,zm; in FormFunction() local 120 PetscCall(DMDAGetCorners(da,&xs,&ys,&zs,&xm,&ym,&zm)); in FormFunction() 125 for (k=zs; k<zs+zm; k++) { in FormFunction() 161 PetscInt i,j,k,Mx,My,Mz,xs,ys,zs,xm,ym,zm; in FormJacobian() local 197 PetscCall(DMDAGetCorners(da,&xs,&ys,&zs,&xm,&ym,&zm)); in FormJacobian() 214 for (k=zs; k<zs+zm; k++) { in FormJacobian()
|
| /petsc/src/binding/petsc4py/demo/legacy/wrap-cython/ |
| H A D | Bratu3Dimpl.c | 24 PetscInt i,j,k,Mx,My,Mz,xs,ys,zs,xm,ym,zm; in FormInitGuess() local 56 PetscCall(DMDAGetCorners(da,&xs,&ys,&zs,&xm,&ym,&zm)); in FormInitGuess() 61 for (k=zs; k<zs+zm; k++) { in FormInitGuess() 85 PetscInt i,j,k,Mx,My,Mz,xs,ys,zs,xm,ym,zm; in FormFunction() local 124 PetscCall(DMDAGetCorners(da,&xs,&ys,&zs,&xm,&ym,&zm)); in FormFunction() 129 for (k=zs; k<zs+zm; k++) { in FormFunction() 165 PetscInt i,j,k,Mx,My,Mz,xs,ys,zs,xm,ym,zm; in FormJacobian() local 204 PetscCall(DMDAGetCorners(da,&xs,&ys,&zs,&xm,&ym,&zm)); in FormJacobian() 221 for (k=zs; k<zs+zm; k++) { in FormJacobian()
|
| /petsc/src/ksp/ksp/tutorials/ |
| H A D | ex45.c | 76 PetscInt i, j, k, mx, my, mz, xm, ym, zm, xs, ys, zs; in ComputeRHS() local 91 PetscCall(DMDAGetCorners(dm, &xs, &ys, &zs, &xm, &ym, &zm)); in ComputeRHS() 94 for (k = zs; k < zs + zm; k++) { in ComputeRHS() 120 PetscInt i, j, k, mx, my, mz, xm, ym, zm, xs, ys, zs; in ComputeMatrix() local 134 PetscCall(DMDAGetCorners(da, &xs, &ys, &zs, &xm, &ym, &zm)); in ComputeMatrix() 136 for (k = zs; k < zs + zm; k++) { in ComputeMatrix()
|
| H A D | ex34.c | 31 PetscInt i, j, k, mx, my, mz, xm, ym, zm, xs, ys, zs, d, dof; in main() local 67 PetscCall(DMDAGetCorners(da, &xs, &ys, &zs, &xm, &ym, &zm)); in main() 70 for (k = zs; k < zs + zm; k++) { in main() 97 PetscInt d, dof, i, j, k, mx, my, mz, xm, ym, zm, xs, ys, zs; in ComputeRHS() local 109 PetscCall(DMDAGetCorners(da, &xs, &ys, &zs, &xm, &ym, &zm)); in ComputeRHS() 111 for (k = zs; k < zs + zm; k++) { in ComputeRHS() 135 PetscInt dof, i, j, k, d, mx, my, mz, xm, ym, zm, xs, ys, zs, num, numi, numj, numk; in ComputeMatrix() local 151 PetscCall(DMDAGetCorners(da, &xs, &ys, &zs, &xm, &ym, &zm)); in ComputeMatrix() 152 for (k = zs; k < zs + zm; k++) { in ComputeMatrix()
|
| H A D | ex22f.F90 | 43 PetscInt ym, zm, xs, ys, zs, i1, i7 60 PetscCall(DMDAGetCorners(da, xs, ys, zs, xm, ym, zm, ierr)) 62 do k = zs, zs + zm - 1
|
| /petsc/src/dm/impls/da/kokkos/ |
| H A D | dagetov.kokkos.cxx | 8 #define DMDA_VEC_GET_SHAPE(da, vec, xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, N, dim, d… argument 10 PetscCall(DMDAGetCorners(da, &xs, &ys, &zs, &xm, &ym, &zm)); \ 21 gzs = zs; \ 29 …PetscInt xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, N, di… in DMDAVecGetKokkosOffsetView_Private() local 36 DMDA_VEC_GET_SHAPE(da, vec, xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, N, dim, dof); in DMDAVecGetKokkosOffsetView_Private() 63 …PetscInt xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, … in DMDAVecGetKokkosOffsetView() local 70 DMDA_VEC_GET_SHAPE(da, vec, xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, N, dim, dof); in DMDAVecGetKokkosOffsetView() 95 …PetscInt xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, N, di… in DMDAVecGetKokkosOffsetView_Private() local 102 DMDA_VEC_GET_SHAPE(da, vec, xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, N, dim, dof); in DMDAVecGetKokkosOffsetView_Private() 129 …PetscInt xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, … in DMDAVecGetKokkosOffsetView() local [all …]
|
| /petsc/src/dm/tutorials/ |
| H A D | ex11f90.F90 | 15 PetscInt zs, zl, sw 76 PetscCallA(DMDAGetCorners(ada, xs, ys, zs, xl, yl, zl, ierr)) 80 do k = zs, zs + zl - 1 145 PetscCallA(DMDAGetCorners(ada, xs, ys, zs, xl, yl, zl, ierr)) 149 do k = zs, zs + zl - 1
|
| H A D | ex14.c | 28 for (k = info.zs; k < info.zs + info.zm; k++) { in FillLocalSubdomain() 104 lower.k = info.zs; in main() 107 upper.k = info.zs + info.zm; in main()
|
| /petsc/src/dm/impls/da/hypre/ |
| H A D | mhyp.h | 24 PetscInt rstart, gnx, gnxgny, xs, ys, zs, nx, ny, nxny; member 45 PetscInt rstart, gnx, gnxgny, gnxgnygnz, xs, ys, zs, nx, ny, nz, nxny, nxnynz; member
|
| /petsc/src/binding/petsc4py/demo/legacy/bratu3d/ |
| H A D | bratu3d.py | 46 (xs, xe), (ys, ye), (zs, ze) = self.da.getRanges() 47 for k in range(zs, ze): 76 (xs, xe), (ys, ye), (zs, ze) = self.da.getRanges() 77 for k in range(zs, ze): 114 (xs, xe), (ys, ye), (zs, ze) = self.da.getRanges() 115 for k in range(zs, ze):
|