Home
last modified time | relevance | path

Searched refs:jstart (Results 1 – 12 of 12) sorted by relevance

/petsc/src/mat/impls/aij/seq/ftn-kernels/
H A Dfmult.F9014 PetscInt :: i, jstart, jend
18 jstart = jend
20 y(jj(jstart:jend - 1)) = y(jj(jstart:jend - 1)) + x(i - 1)*a(jstart:jend - 1)
31 PetscInt :: i, jstart, jend
37 jstart = ii(i)
39 y(i) = sum(a(jstart:jend - 1)*x(jj(jstart:jend - 1)))
H A Dfsolve.F9016 PetscInt i, jstart, jend
22 jstart = ai(i)
24 x(i) = b(i) - sum(aa(jstart:jend)*x(aj(jstart:jend)))
30 jstart = adiag(i) + 1
32 x(i) = x(i) - sum(aa(jstart:jend)*x(aj(jstart:jend)))*aa(adiag(i))
H A Dfmultadd.F9014 PetscInt :: i, jstart, jend
18 jstart = jend
20 z(i) = y(i) + sum(a(jstart:jend - 1)*x(jj(jstart:jend - 1)))
/petsc/src/mat/impls/baij/seq/ftn-kernels/
H A Dfsolvebaij.F9019 PetscInt :: i, j, jstart, jend
36 jstart = ai(i)
38 ax = 16*jstart
41 do j = jstart, jend
57 jstart = adiag(i) + 1
59 ax = 16*jstart
61 do j = jstart, jend
90 PetscInt :: jstart, jend, idx, ax, jdx, kdx, nn
110 jstart = ai(i)
113 if (jend - jstart >= 500) error stop 'Overflowing vector FortranSolveBAIJ4()'
[all …]
/petsc/src/ksp/pc/impls/wb/
H A Dwb.c19 …, Nint, Nface, Nwire, Nsurf, *Iint, *Isurf, cint = 0, csurf = 0, istart, jstart, kstart, *II, N, c… in DMDAGetWireBasketInterpolation() local
39 PetscCall(DMDAGetGhostCorners(da, &istart, &jstart, &kstart, &mwidth, &nwidth, &pwidth)); in DMDAGetWireBasketInterpolation()
41 jstart = jstart ? -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 …Nface = 2 * ((m - 2 - istart) * (n - 2 - jstart) + (m - 2 - istart) * (p - 2 - kstart) + (n - 2 - … in DMDAGetWireBasketInterpolation()
60 Nwire = 4 * ((m - 2 - istart) + (n - 2 - jstart) + (p - 2 - kstart)) + 8; in DMDAGetWireBasketInterpolation()
71 …PetscCheck(n - jstart >= 3, PETSC_COMM_SELF, PETSC_ERR_SUP, "Number of grid points per process in … in DMDAGetWireBasketInterpolation()
80 for (j = 1; j < n - 1 - jstart; j++) { in DMDAGetWireBasketInterpolation()
95 for (j = 1; j < n - 1 - jstart; j++) { in DMDAGetWireBasketInterpolation()
[all …]
/petsc/src/mat/impls/aij/seq/aijperm/
H A Daijperm.c299 PetscInt jstart, jend; in MatMult_SeqAIJPERM() local
327 jstart = xgroup[igroup]; in MatMult_SeqAIJPERM()
334 for (i = jstart; i <= jend; i++) y[iperm[i]] = 0.0; in MatMult_SeqAIJPERM()
336 for (i = jstart; i <= jend; i++) { in MatMult_SeqAIJPERM()
345 for (istart = jstart; istart <= jend; istart += NDIM) { in MatMult_SeqAIJPERM()
476 PetscInt jstart, jend; in MatMultAdd_SeqAIJPERM() local
508 jstart = xgroup[igroup]; in MatMultAdd_SeqAIJPERM()
516 for (i = jstart; i <= jend; i++) { in MatMultAdd_SeqAIJPERM()
521 for (i = jstart; i <= jend; i++) { in MatMultAdd_SeqAIJPERM()
532 for (istart = jstart; istart <= jend; istart += NDIM) { in MatMultAdd_SeqAIJPERM()
/petsc/src/mat/impls/aij/seq/bas/
H A Dspbas.c848 PetscInt jstart; in spbas_keep_upper() local
853 …for (jstart = 0; (jstart < inout_matrix->row_nnz[i]) && (inout_matrix->icols[i][jstart] < 0); jsta… in spbas_keep_upper()
854 if (jstart > 0) { in spbas_keep_upper()
855 …(j = 0; j < inout_matrix->row_nnz[i] - jstart; j++) inout_matrix->icols[i][j] = inout_matrix->icol… in spbas_keep_upper()
858 …j = 0; j < inout_matrix->row_nnz[i] - jstart; j++) inout_matrix->values[i][j] = inout_matrix->valu… in spbas_keep_upper()
861 inout_matrix->row_nnz[i] -= jstart; in spbas_keep_upper()
866 inout_matrix->nnz -= jstart; in spbas_keep_upper()
/petsc/src/snes/tutorials/ex10d/
H A Dex10.c89 PetscInt i, j, jstart, inode, nb, nbrs, Nvneighborstotal = 0; in main() local
237 jstart = 0; in main()
243 tmp[j + jstart] = user.AdjM[i][j]; in main()
245 jstart += user.itot[i]; in main()
260 jstart = 0; in main()
264 user.AdjM[i][j] = tmp[j + jstart]; in main()
268 jstart += user.itot[i]; in main()
/petsc/src/dm/impls/da/
H A Dfdda.c927 …PetscInt istart, iend, jstart, jend, kstart, kend, zs, nz, gzs, gnz, ii, jj, kk, M, … in DMCreateMatrix_DA_3d_MPISELL() local
955 jstart = (by == DM_BOUNDARY_PERIODIC) ? -s : (PetscMax(-s, -j)); in DMCreateMatrix_DA_3d_MPISELL()
966 for (jj = jstart; jj < jend + 1; jj++) { in DMCreateMatrix_DA_3d_MPISELL()
997 jstart = (by == DM_BOUNDARY_PERIODIC) ? -s : (PetscMax(-s, -j)); in DMCreateMatrix_DA_3d_MPISELL()
1008 for (jj = jstart; jj < jend + 1; jj++) { 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
1311 jstart = (by == DM_BOUNDARY_PERIODIC) ? -s : (PetscMax(-s, -j)); in DMCreateMatrix_DA_3d_MPIAIJ()
1322 for (jj = jstart; jj < jend + 1; jj++) { in DMCreateMatrix_DA_3d_MPIAIJ()
1354 jstart = (by == DM_BOUNDARY_PERIODIC) ? -s : (PetscMax(-s, -j)); in DMCreateMatrix_DA_3d_MPIAIJ()
1364 for (jj = jstart; jj < jend + 1; jj++) { in DMCreateMatrix_DA_3d_MPIAIJ()
[all …]
H A Dgr1.c34 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()
63 coors[cnt++] = ymin + hy * (j + jstart); in DMDASetUniformCoordinates()
80 coors[cnt++] = ymin + hy * (j + jstart); in DMDASetUniformCoordinates()
/petsc/src/binding/petsc4py/test/
H A Dtest_dmda.py297 istart, iend, jstart, jend = da.getRanges()
301 for j in range(jstart, jend):
/petsc/src/tao/constrained/impls/ipm/
H A Dipm.c675 PetscInt astart, aend, jstart, jend; in IPMUpdateAi() local
720 PetscCall(MatGetOwnershipRange(tao->jacobian_inequality, &jstart, &jend)); in IPMUpdateAi()
721 for (i = jstart; i < jend; i++) { in IPMUpdateAi()