| /petsc/src/mat/impls/sbaij/seq/ |
| H A D | sbaij2.c | 294 const PetscInt *aj = a->j, *ai = a->i, *ib; in MatMult_SeqSBAIJ_2() local
310 ib = aj + *ai; in MatMult_SeqSBAIJ_2()
313 if (*ib == i) { /* (diag of A)*x */ in MatMult_SeqSBAIJ_2()
319 …PetscPrefetchBlock(ib + jmin + n, n, 0, PETSC_PREFETCH_HINT_NTA); /* Indices for the next row (ass… in MatMult_SeqSBAIJ_2()
323 cval = ib[j] * 2; in MatMult_SeqSBAIJ_2()
348 const PetscInt *aj = a->j, *ai = a->i, *ib; in MatMult_SeqSBAIJ_3() local
365 ib = aj + *ai; in MatMult_SeqSBAIJ_3()
368 if (*ib == i) { /* (diag of A)*x */ in MatMult_SeqSBAIJ_3()
375 …PetscPrefetchBlock(ib + jmin + n, n, 0, PETSC_PREFETCH_HINT_NTA); /* Indices for the next row (ass… in MatMult_SeqSBAIJ_3()
379 cval = ib[j] * 3; in MatMult_SeqSBAIJ_3()
[all …]
|
| H A D | relax.h | 67 const unsigned short *ib = a->jshort; in MatMult_SeqSBAIJ_1_ushort() local
70 const PetscInt *ib = a->j; in MatMult_SeqSBAIJ_1_ushort() local
89 if (ib[0] == i) { in MatMult_SeqSBAIJ_1_ushort()
93 …PetscPrefetchBlock(ib + nz, nz, 0, PETSC_PREFETCH_HINT_NTA); /* Indices for the next row (assumes … in MatMult_SeqSBAIJ_1_ushort()
97 ibt = ib[j]; in MatMult_SeqSBAIJ_1_ushort()
104 ibt = ib[j]; in MatMult_SeqSBAIJ_1_ushort()
112 ib += nz; in MatMult_SeqSBAIJ_1_ushort()
|
| /petsc/src/ksp/pc/impls/vpbjacobi/ |
| H A D | vpbjacobi.c | 9 PetscInt ib, jb, bs; in PCApply_VPBJacobi() local
92 for (ib = 0; ib < bs; ib++) { in PCApply_VPBJacobi()
94 for (jb = 0; jb < bs; jb++) rowsum += diag[ib + jb * bs] * xx[ncnt + jb]; in PCApply_VPBJacobi()
95 yy[ncnt + ib] = rowsum; in PCApply_VPBJacobi()
111 PetscInt ib, jb, bs; in PCApplyTranspose_VPBJacobi() local
194 for (ib = 0; ib < bs; ib++) { in PCApplyTranspose_VPBJacobi()
196 for (jb = 0; jb < bs; jb++) rowsum += diag[ib * bs + jb] * xx[ncnt + jb]; in PCApplyTranspose_VPBJacobi()
197 yy[ncnt + ib] = rowsum; in PCApplyTranspose_VPBJacobi()
|
| /petsc/src/ksp/pc/impls/pbjacobi/ |
| H A D | pbjacobi.c | 7 PetscInt i, ib, jb; in PCApply_PBJacobi() local
112 for (ib = 0; ib < bs; ib++) { in PCApply_PBJacobi()
114 for (jb = 0; jb < bs; jb++) rowsum += diag[ib + jb * bs] * xx[bs * i + jb]; in PCApply_PBJacobi()
115 yy[bs * i + ib] = rowsum; in PCApply_PBJacobi()
129 PetscInt i, ib, jb; in PCApplyTranspose_PBJacobi() local
234 for (ib = 0; ib < bs; ib++) { in PCApplyTranspose_PBJacobi()
236 for (jb = 0; jb < bs; jb++) rowsum += diag[ib * bs + jb] * xx[bs * i + jb]; in PCApplyTranspose_PBJacobi()
237 yy[bs * i + ib] = rowsum; in PCApplyTranspose_PBJacobi()
|
| /petsc/src/mat/impls/baij/mpi/ |
| H A D | mpibaij.c | 2239 PetscInt ib, jb, bs = A->rmap->bs; in MatGetColumnReductions_MPIBAIJ() local
2252 for (ib = 0; ib < bs; ib++) { in MatGetColumnReductions_MPIBAIJ()
2260 for (ib = 0; ib < bs; ib++) { in MatGetColumnReductions_MPIBAIJ()
2269 for (ib = 0; ib < bs; ib++) { in MatGetColumnReductions_MPIBAIJ()
2277 for (ib = 0; ib < bs; ib++) { in MatGetColumnReductions_MPIBAIJ()
2286 for (ib = 0; ib < bs; ib++) { in MatGetColumnReductions_MPIBAIJ()
2295 for (ib = 0; ib < bs; ib++) { in MatGetColumnReductions_MPIBAIJ()
2305 for (ib = 0; ib < bs; ib++) { in MatGetColumnReductions_MPIBAIJ()
2313 for (ib = 0; ib < bs; ib++) { in MatGetColumnReductions_MPIBAIJ()
2322 for (ib = 0; ib < bs; ib++) { in MatGetColumnReductions_MPIBAIJ()
[all …]
|
| /petsc/src/mat/impls/baij/seq/ |
| H A D | baij2.c | 2832 const PetscInt *idx, *ii, *ib, *ridx = NULL; in MatMultHermitianTransposeAdd_SeqBAIJ() local
2858 ib = idx + ii[0]; in MatMultHermitianTransposeAdd_SeqBAIJ()
2862 rval = ib[j]; in MatMultHermitianTransposeAdd_SeqBAIJ()
2874 ib = idx + ii[0]; in MatMultHermitianTransposeAdd_SeqBAIJ()
2878 rval = ib[j] * 2; in MatMultHermitianTransposeAdd_SeqBAIJ()
2892 ib = idx + ii[0]; in MatMultHermitianTransposeAdd_SeqBAIJ()
2896 rval = ib[j] * 3; in MatMultHermitianTransposeAdd_SeqBAIJ()
2912 ib = idx + ii[0]; in MatMultHermitianTransposeAdd_SeqBAIJ()
2916 rval = ib[j] * 4; in MatMultHermitianTransposeAdd_SeqBAIJ()
2934 ib = idx + ii[0]; in MatMultHermitianTransposeAdd_SeqBAIJ()
[all …]
|
| H A D | baij.c | 37 PetscInt m, n, ib, jb, bs = A->rmap->bs; in MatGetColumnReductions_SeqBAIJ() local
46 for (ib = 0; ib < bs; ib++) { in MatGetColumnReductions_SeqBAIJ()
55 for (ib = 0; ib < bs; ib++) { in MatGetColumnReductions_SeqBAIJ()
64 for (ib = 0; ib < bs; ib++) { in MatGetColumnReductions_SeqBAIJ()
74 for (ib = 0; ib < bs; ib++) { in MatGetColumnReductions_SeqBAIJ()
83 for (ib = 0; ib < bs; ib++) { in MatGetColumnReductions_SeqBAIJ()
|
| /petsc/src/ts/tutorials/advection-diffusion-reaction/ |
| H A D | shashi.F90 | 402 PetscInt jj, jb, ii3, id, ib, i
484 do ib = 1, 26
485 d_eq(ib, jb) = 0.0d0
|
| /petsc/src/dm/dt/interface/ |
| H A D | dt.c | 1572 PetscInt ia, ib; in PetscDTGaussJacobiQuadrature_Newton_Internal() local
1575 ib = (PetscInt)b; in PetscDTGaussJacobiQuadrature_Newton_Internal()
1579 } else if (b == b && ib >= 0) { /* compute ratio of rising factorials wrt b */ in PetscDTGaussJacobiQuadrature_Newton_Internal()
1580 for (k = 0; k < ib; k++) gf *= (npoints + 1. + k) / (npoints + a + 1. + k); in PetscDTGaussJacobiQuadrature_Newton_Internal()
1644 PetscInt ia = (PetscInt)a, ib = (PetscInt)b; in PetscDTGaussJacobiQuadrature_GolubWelsch_Internal() local
1646 …PetscCheck(ia == a && ib == b && ia + 1 > 0 && ib + 1 > 0 && ia + ib + 2 > 0, PETSC_COMM_SELF, PET… in PetscDTGaussJacobiQuadrature_GolubWelsch_Internal()
1649 PetscCall(PetscDTFactorial(ib, &gb)); in PetscDTGaussJacobiQuadrature_GolubWelsch_Internal()
1650 PetscCall(PetscDTFactorial(ia + ib + 1, &gab)); in PetscDTGaussJacobiQuadrature_GolubWelsch_Internal()
|
| /petsc/src/mat/impls/aij/mpi/ |
| H A D | mpiov.c | 1234 PetscInt *cworkB, lwrite, *subcols, ib, jb; in MatCreateSubMatrices_MPIAIJ_SingleIS_Local() local
1602 PetscCall(ISGetBlockSize(isrow[0], &ib)); in MatCreateSubMatrices_MPIAIJ_SingleIS_Local()
1604 if (ib > 1 || jb > 1) PetscCall(MatSetBlockSizes(submat, ib, jb)); in MatCreateSubMatrices_MPIAIJ_SingleIS_Local()
|
| H A D | mpiaij.c | 208 const PetscInt *ia, *ib; in MatFindNonzeroRows_MPIAIJ() local
217 ib = b->i; in MatFindNonzeroRows_MPIAIJ()
222 nb = ib[i + 1] - ib[i]; in MatFindNonzeroRows_MPIAIJ()
231 bb = PetscSafePointerPlusOffset(bav, ib[i]); in MatFindNonzeroRows_MPIAIJ()
248 nb = ib[i + 1] - ib[i]; in MatFindNonzeroRows_MPIAIJ()
257 bb = PetscSafePointerPlusOffset(bav, ib[i]); in MatFindNonzeroRows_MPIAIJ()
|