| /petsc/src/ksp/pc/impls/wb/ |
| H A D | wb.c | 19 …nt, Nface, Nwire, Nsurf, *Iint, *Isurf, cint = 0, csurf = 0, istart, jstart, kstart, *II, N, c = 0; in DMDAGetWireBasketInterpolation() local 39 PetscCall(DMDAGetGhostCorners(da, &istart, &jstart, &kstart, &mwidth, &nwidth, &pwidth)); in DMDAGetWireBasketInterpolation() 42 kstart = kstart ? -1 : 0; in DMDAGetWireBasketInterpolation() 57 N = (m - istart) * (n - jstart) * (p - kstart); in DMDAGetWireBasketInterpolation() 58 Nint = (m - 2 - istart) * (n - 2 - jstart) * (p - 2 - kstart); in DMDAGetWireBasketInterpolation() 59 …art) * (n - 2 - jstart) + (m - 2 - istart) * (p - 2 - kstart) + (n - 2 - jstart) * (p - 2 - kstart… in DMDAGetWireBasketInterpolation() 60 Nwire = 4 * ((m - 2 - istart) + (n - 2 - jstart) + (p - 2 - kstart)) + 8; in DMDAGetWireBasketInterpolation() 72 …PetscCheck(p - kstart >= 3, PETSC_COMM_SELF, PETSC_ERR_SUP, "Number of grid points per process in … in DMDAGetWireBasketInterpolation() 90 for (k = 1; k < p - 1 - kstart; k++) { in DMDAGetWireBasketInterpolation() 141 for (k = 0; k < p - kstart; k++) { in DMDAGetWireBasketInterpolation() [all …]
|
| /petsc/src/tao/constrained/impls/ipm/ |
| H A D | ipm.c | 810 PetscInt hstart, hend, kstart, kend; in IPMUpdateK() local 832 PetscCall(VecGetOwnershipRange(ipmP->bigrhs, &kstart, &kend)); in IPMUpdateK() 834 klocalsize = kend - kstart; in IPMUpdateK() 837 PetscCall(PetscMalloc1(kend - kstart, &nonzeros)); in IPMUpdateK() 845 nonzeros[i - kstart] = ipmP->n; in IPMUpdateK() 847 nonzeros[i - kstart] = ipmP->n + 1; in IPMUpdateK() 849 nonzeros[i - kstart] = 2; in IPMUpdateK() 859 PetscCall(PetscMalloc1(kend - kstart, &d_nonzeros)); in IPMUpdateK() 860 PetscCall(PetscMalloc1(kend - kstart, &o_nonzeros)); in IPMUpdateK() 861 for (i = kstart; i < kend; i++) { in IPMUpdateK() [all …]
|
| /petsc/src/dm/impls/da/ |
| H A D | fdda.c | 927 …PetscInt istart, iend, jstart, jend, kstart, kend, zs, nz, gzs, gnz, ii, jj, kk, M, … in DMCreateMatrix_DA_3d_MPISELL() local 958 kstart = (bz == DM_BOUNDARY_PERIODIC) ? -s : (PetscMax(-s, -k)); in DMCreateMatrix_DA_3d_MPISELL() 967 for (kk = kstart; kk < kend + 1; kk++) { in DMCreateMatrix_DA_3d_MPISELL() 1000 kstart = (bz == DM_BOUNDARY_PERIODIC) ? -s : (PetscMax(-s, -k)); in DMCreateMatrix_DA_3d_MPISELL() 1009 for (kk = kstart; kk < kend + 1; kk++) { in DMCreateMatrix_DA_3d_MPISELL() 1273 …PetscInt istart, iend, jstart, jend, kstart, kend, zs, nz, gzs, gnz, ii, jj, kk, M, … in DMCreateMatrix_DA_3d_MPIAIJ() local 1314 kstart = (bz == DM_BOUNDARY_PERIODIC) ? -s : (PetscMax(-s, -k)); in DMCreateMatrix_DA_3d_MPIAIJ() 1323 for (kk = kstart; kk < kend + 1; kk++) { in DMCreateMatrix_DA_3d_MPIAIJ() 1357 kstart = (bz == DM_BOUNDARY_PERIODIC) ? -s : (PetscMax(-s, -k)); in DMCreateMatrix_DA_3d_MPIAIJ() 1363 for (kk = kstart; kk < kend + 1; kk++) { in DMCreateMatrix_DA_3d_MPIAIJ() [all …]
|
| H A D | gr1.c | 34 PetscInt i, j, k, M, N, P, istart, isize, jstart, jsize, kstart, ksize, dim, cnt; in DMDASetUniformCoordinates() local 44 PetscCall(DMDAGetCorners(da, &istart, &jstart, &kstart, &isize, &jsize, &ksize)); in DMDASetUniformCoordinates() 81 coors[cnt++] = zmin + hz_ * (k + kstart); in DMDASetUniformCoordinates()
|
| /petsc/src/mat/impls/sbaij/seq/ |
| H A D | sbaij2.c | 109 PetscInt *smap, i, k, kstart, kend, oldcols = a->nbs, *lens; in MatCreateSubMatrix_SeqSBAIJ_Private() local 130 kstart = ai[irow[i]]; in MatCreateSubMatrix_SeqSBAIJ_Private() 131 kend = kstart + a->ilen[irow[i]]; in MatCreateSubMatrix_SeqSBAIJ_Private() 133 for (k = kstart; k < kend; k++) { in MatCreateSubMatrix_SeqSBAIJ_Private() 170 kstart = ai[row]; in MatCreateSubMatrix_SeqSBAIJ_Private() 171 kend = kstart + a->ilen[row]; in MatCreateSubMatrix_SeqSBAIJ_Private() 183 for (k = kstart; k < kend; k++) { in MatCreateSubMatrix_SeqSBAIJ_Private()
|
| /petsc/src/mat/impls/aij/seq/ |
| H A D | aij.c | 2391 PetscInt *smap, i, k, kstart, kend, oldcols = A->cmap->n, *lens; in MatCreateSubMatrix_SeqAIJ() local 2418 kstart = ai[irow[i]]; in MatCreateSubMatrix_SeqAIJ() 2419 kend = kstart + ailen[irow[i]]; in MatCreateSubMatrix_SeqAIJ() 2420 starts[i] = kstart; in MatCreateSubMatrix_SeqAIJ() 2421 for (k = kstart; k < kend; k++) { in MatCreateSubMatrix_SeqAIJ() 2483 kstart = ai[irow[i]]; in MatCreateSubMatrix_SeqAIJ() 2484 kend = kstart + a->ilen[irow[i]]; in MatCreateSubMatrix_SeqAIJ() 2486 for (k = kstart; k < kend; k++) { in MatCreateSubMatrix_SeqAIJ() 2516 kstart = ai[row]; in MatCreateSubMatrix_SeqAIJ() 2517 kend = kstart + a->ilen[row]; in MatCreateSubMatrix_SeqAIJ() [all …]
|
| /petsc/src/mat/impls/baij/seq/ |
| H A D | baij2.c | 73 PetscInt *smap, i, k, kstart, kend, oldcols = a->nbs, *lens; in MatCreateSubMatrix_SeqBAIJ_Private() local 94 kstart = ai[irow[i]]; in MatCreateSubMatrix_SeqBAIJ_Private() 95 kend = kstart + a->ilen[irow[i]]; in MatCreateSubMatrix_SeqBAIJ_Private() 97 for (k = kstart; k < kend; k++) { in MatCreateSubMatrix_SeqBAIJ_Private() 119 kstart = ai[row]; in MatCreateSubMatrix_SeqBAIJ_Private() 120 kend = kstart + a->ilen[row]; in MatCreateSubMatrix_SeqBAIJ_Private() 125 for (k = kstart; k < kend; k++) { in MatCreateSubMatrix_SeqBAIJ_Private()
|