| /petsc/src/ksp/pc/impls/vpbjacobi/ |
| H A D | vpbjacobi.c | 9 PetscInt ib, jb, bs; in PCApply_VPBJacobi() local 94 for (jb = 0; jb < bs; jb++) rowsum += diag[ib + jb * bs] * xx[ncnt + jb]; in PCApply_VPBJacobi() 111 PetscInt ib, jb, bs; in PCApplyTranspose_VPBJacobi() local 196 for (jb = 0; jb < bs; jb++) rowsum += diag[ib * bs + jb] * xx[ncnt + jb]; in PCApplyTranspose_VPBJacobi()
|
| /petsc/src/ksp/pc/impls/pbjacobi/ |
| H A D | pbjacobi.c | 7 PetscInt i, ib, jb; in PCApply_PBJacobi() local 114 for (jb = 0; jb < bs; jb++) rowsum += diag[ib + jb * bs] * xx[bs * i + jb]; in PCApply_PBJacobi() 129 PetscInt i, ib, jb; in PCApplyTranspose_PBJacobi() local 236 for (jb = 0; jb < bs; jb++) rowsum += diag[ib * bs + jb] * xx[bs * i + jb]; in PCApplyTranspose_PBJacobi()
|
| /petsc/src/mat/impls/baij/mpi/ |
| H A D | mpibaij.c | 1131 PetscInt header[4], M, N, m, rs, cs, bs, cnt, i, j, ja, jb, k, l; in MatView_MPIBAIJ_Binary() local 1168 for (jb = B->i[i]; jb < B->i[i + 1]; jb++) { in MatView_MPIBAIJ_Binary() 1169 if (garray[B->j[jb]] > cs / bs) break; in MatView_MPIBAIJ_Binary() 1170 for (l = 0; l < bs; l++) colidxs[cnt++] = bs * garray[B->j[jb]] + l; in MatView_MPIBAIJ_Binary() 1174 for (; jb < B->i[i + 1]; jb++) in MatView_MPIBAIJ_Binary() 1175 for (l = 0; l < bs; l++) colidxs[cnt++] = bs * garray[B->j[jb]] + l; in MatView_MPIBAIJ_Binary() 1186 for (jb = B->i[i]; jb < B->i[i + 1]; jb++) { in MatView_MPIBAIJ_Binary() 1187 if (garray[B->j[jb]] > cs / bs) break; in MatView_MPIBAIJ_Binary() 1188 for (l = 0; l < bs; l++) matvals[cnt++] = B->a[bs * (bs * jb + l) + k]; in MatView_MPIBAIJ_Binary() 1192 for (; jb < B->i[i + 1]; jb++) in MatView_MPIBAIJ_Binary() [all …]
|
| /petsc/src/vec/vec/impls/mpi/ |
| H A D | pbvec.c | 207 PetscInt j, jb, bs; in VecAssemblyBegin_MPI_BTS() local 238 for (i = 0, j = 0, jb = 0; i < x->nsendranks; i++) { in VecAssemblyBegin_MPI_BTS() 253 x->sendptrs[i].intb = &X->bstash.idx[jb]; in VecAssemblyBegin_MPI_BTS() 254 x->sendptrs[i].scalarb = &X->bstash.array[jb * bs]; in VecAssemblyBegin_MPI_BTS() 255 …for (; jb < X->bstash.n && X->bstash.idx[jb] * bs < X->map->range[rank + 1]; jb++) x->sendhdr[i].b… in VecAssemblyBegin_MPI_BTS()
|
| /petsc/src/mat/impls/baij/seq/ |
| H A D | baij.c | 37 PetscInt m, n, ib, jb, bs = A->rmap->bs; in MatGetColumnReductions_SeqBAIJ() local 45 for (jb = 0; jb < bs; jb++) { in MatGetColumnReductions_SeqBAIJ() 47 reductions[A->cmap->rstart + a_aij->j[i] * bs + jb] += PetscAbsScalar(*a_val * *a_val); in MatGetColumnReductions_SeqBAIJ() 54 for (jb = 0; jb < bs; jb++) { in MatGetColumnReductions_SeqBAIJ() 56 reductions[A->cmap->rstart + a_aij->j[i] * bs + jb] += PetscAbsScalar(*a_val); in MatGetColumnReductions_SeqBAIJ() 63 for (jb = 0; jb < bs; jb++) { in MatGetColumnReductions_SeqBAIJ() 65 PetscInt col = A->cmap->rstart + a_aij->j[i] * bs + jb; in MatGetColumnReductions_SeqBAIJ() 73 for (jb = 0; jb < bs; jb++) { in MatGetColumnReductions_SeqBAIJ() 75 reductions[A->cmap->rstart + a_aij->j[i] * bs + jb] += PetscRealPart(*a_val); in MatGetColumnReductions_SeqBAIJ() 82 for (jb = 0; jb < bs; jb++) { in MatGetColumnReductions_SeqBAIJ() [all …]
|
| /petsc/src/ts/tutorials/advection-diffusion-reaction/ |
| H A D | shashi.F90 | 402 PetscInt jj, jb, ii3, id, ib, i 483 do jb = 1, 26 485 d_eq(ib, jb) = 0.0d0
|
| /petsc/src/mat/impls/aij/mpi/ |
| H A D | mpiaij.c | 1191 PetscInt header[4], M, N, m, rs, cs, cnt, i, ja, jb; in MatView_MPIAIJ_Binary() local 1226 for (jb = B->i[i]; jb < B->i[i + 1]; jb++) { in MatView_MPIAIJ_Binary() 1227 if (garray[B->j[jb]] > cs) break; in MatView_MPIAIJ_Binary() 1228 colidxs[cnt++] = garray[B->j[jb]]; in MatView_MPIAIJ_Binary() 1231 for (; jb < B->i[i + 1]; jb++) colidxs[cnt++] = garray[B->j[jb]]; in MatView_MPIAIJ_Binary() 1242 for (jb = B->i[i]; jb < B->i[i + 1]; jb++) { in MatView_MPIAIJ_Binary() 1243 if (garray[B->j[jb]] > cs) break; in MatView_MPIAIJ_Binary() 1244 matvals[cnt++] = ba[jb]; in MatView_MPIAIJ_Binary() 1247 for (; jb < B->i[i + 1]; jb++) matvals[cnt++] = ba[jb]; in MatView_MPIAIJ_Binary()
|
| H A D | mpiov.c | 1234 PetscInt *cworkB, lwrite, *subcols, ib, jb; in MatCreateSubMatrices_MPIAIJ_SingleIS_Local() local 1603 PetscCall(ISGetBlockSize(iscol[0], &jb)); in MatCreateSubMatrices_MPIAIJ_SingleIS_Local() 1604 if (ib > 1 || jb > 1) PetscCall(MatSetBlockSizes(submat, ib, jb)); in MatCreateSubMatrices_MPIAIJ_SingleIS_Local()
|