Lines Matching refs:mbs
11 PetscInt brow, i, j, k, l, mbs, n, *nidx, isz, bcol, bcol_max, start, end, *ai, *aj, bs; in MatIncreaseOverlap_SeqSBAIJ() local
17 mbs = a->mbs; in MatIncreaseOverlap_SeqSBAIJ()
21 PetscCall(PetscBTCreate(mbs, &table_out)); in MatIncreaseOverlap_SeqSBAIJ()
22 PetscCall(PetscMalloc1(mbs + 1, &nidx)); in MatIncreaseOverlap_SeqSBAIJ()
23 PetscCall(PetscBTCreate(mbs, &table_in)); in MatIncreaseOverlap_SeqSBAIJ()
27 PetscCall(PetscBTMemzero(mbs, table_out)); in MatIncreaseOverlap_SeqSBAIJ()
37 … PetscCheck(brow < mbs, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "index greater than mat-dim"); in MatIncreaseOverlap_SeqSBAIJ()
49 PetscCall(PetscBTMemzero(mbs, table_in)); in MatIncreaseOverlap_SeqSBAIJ()
53 for (brow = 0; brow < mbs; brow++) { in MatIncreaseOverlap_SeqSBAIJ()
142 …PetscCheck(c->mbs == nrows && c->nbs == ncols && (*B)->rmap->bs == bs, PETSC_COMM_SELF, PETSC_ERR_… in MatCreateSubMatrix_SeqSBAIJ_Private()
143 PetscCall(PetscArraycmp(c->ilen, lens, c->mbs, &flag)); in MatCreateSubMatrix_SeqSBAIJ_Private()
145 PetscCall(PetscArrayzero(c->ilen, c->mbs)); in MatCreateSubMatrix_SeqSBAIJ_Private()
149 …PetscCheck(d->mbs == nrows && d->nbs == ncols && (*B)->rmap->bs == bs, PETSC_COMM_SELF, PETSC_ERR_… in MatCreateSubMatrix_SeqSBAIJ_Private()
150 PetscCall(PetscArraycmp(d->ilen, lens, d->mbs, &flag)); in MatCreateSubMatrix_SeqSBAIJ_Private()
152 PetscCall(PetscArrayzero(d->ilen, d->mbs)); in MatCreateSubMatrix_SeqSBAIJ_Private()
293 PetscInt mbs = a->mbs, i, n, cval, j, jmin; in MatMult_SeqSBAIJ_2() local
306 for (i = 0; i < mbs; i++) { in MatMult_SeqSBAIJ_2()
347 PetscInt mbs = a->mbs, i, n, cval, j, jmin; in MatMult_SeqSBAIJ_3() local
360 for (i = 0; i < mbs; i++) { in MatMult_SeqSBAIJ_3()
405 PetscInt mbs = a->mbs, i, n, cval, j, jmin; in MatMult_SeqSBAIJ_4() local
418 for (i = 0; i < mbs; i++) { in MatMult_SeqSBAIJ_4()
467 PetscInt mbs = a->mbs, i, n, cval, j, jmin; in MatMult_SeqSBAIJ_5() local
480 for (i = 0; i < mbs; i++) { in MatMult_SeqSBAIJ_5()
533 PetscInt mbs = a->mbs, i, n, cval, j, jmin; in MatMult_SeqSBAIJ_6() local
546 for (i = 0; i < mbs; i++) { in MatMult_SeqSBAIJ_6()
603 PetscInt mbs = a->mbs, i, n, cval, j, jmin; in MatMult_SeqSBAIJ_7() local
616 for (i = 0; i < mbs; i++) { in MatMult_SeqSBAIJ_7()
679 PetscInt mbs = a->mbs, i, bs = A->rmap->bs, j, n, bs2 = a->bs2, ncols, k; in MatMult_SeqSBAIJ_N() local
699 for (i = 0; i < mbs; i++) { in MatMult_SeqSBAIJ_N()
752 PetscInt mbs = a->mbs, i, n, cval, j, jmin; in MatMultAdd_SeqSBAIJ_1() local
765 for (i = 0; i < mbs; i++) { in MatMultAdd_SeqSBAIJ_1()
805 PetscInt mbs = a->mbs, i, n, cval, j, jmin; in MatMultAdd_SeqSBAIJ_2() local
818 for (i = 0; i < mbs; i++) { in MatMultAdd_SeqSBAIJ_2()
859 PetscInt mbs = a->mbs, i, n, cval, j, jmin; in MatMultAdd_SeqSBAIJ_3() local
872 for (i = 0; i < mbs; i++) { in MatMultAdd_SeqSBAIJ_3()
918 PetscInt mbs = a->mbs, i, n, cval, j, jmin; in MatMultAdd_SeqSBAIJ_4() local
931 for (i = 0; i < mbs; i++) { in MatMultAdd_SeqSBAIJ_4()
981 PetscInt mbs = a->mbs, i, n, cval, j, jmin; in MatMultAdd_SeqSBAIJ_5() local
994 for (i = 0; i < mbs; i++) { in MatMultAdd_SeqSBAIJ_5()
1048 PetscInt mbs = a->mbs, i, n, cval, j, jmin; in MatMultAdd_SeqSBAIJ_6() local
1061 for (i = 0; i < mbs; i++) { in MatMultAdd_SeqSBAIJ_6()
1119 PetscInt mbs = a->mbs, i, n, cval, j, jmin; in MatMultAdd_SeqSBAIJ_7() local
1132 for (i = 0; i < mbs; i++) { in MatMultAdd_SeqSBAIJ_7()
1194 PetscInt mbs = a->mbs, i, bs = A->rmap->bs, j, n, bs2 = a->bs2, ncols, k; in MatMultAdd_SeqSBAIJ_N() local
1213 for (i = 0; i < mbs; i++) { in MatMultAdd_SeqSBAIJ_N()
1279 …PetscInt i, j, k, bs = A->rmap->bs, bs2 = a->bs2, k1, mbs = a->mbs, jmin, jmax, nexti, ik,… in MatNorm_SeqSBAIJ() local
1288 for (k = 0; k < mbs; k++) { in MatNorm_SeqSBAIJ()
1309 PetscCall(PetscMalloc3(bs, &sum, mbs, &il, mbs, &jl)); in MatNorm_SeqSBAIJ()
1310 for (i = 0; i < mbs; i++) jl[i] = mbs; in MatNorm_SeqSBAIJ()
1314 for (k = 0; k < mbs; k++) { /* k_th block row */ in MatNorm_SeqSBAIJ()
1320 while (i < mbs) { in MatNorm_SeqSBAIJ()
1367 PetscCall(PetscLogFlops(PetscMax(mbs * a->nz - 1, 0))); in MatNorm_SeqSBAIJ()
1384 PetscCall(PetscArraycmp(a->i, b->i, a->mbs + 1, flg)); in MatEqual_SeqSBAIJ()
1400 const PetscInt bs = A->rmap->bs, ambs = a->mbs, bs2 = a->bs2; in MatGetDiagonal_SeqSBAIJ()
1447 PetscInt i, j, k, lm, M, m, mbs, tmp, bs, bs2; in MatDiagonalScale_SeqSBAIJ() local
1457 mbs = a->mbs; in MatDiagonalScale_SeqSBAIJ()
1463 for (i = 0; i < mbs; i++) { /* for each block row */ in MatDiagonalScale_SeqSBAIJ()
1509 PetscCall(PetscArrayzero(a->a, a->bs2 * a->i[a->mbs])); in MatZeroEntries_SeqSBAIJ()
1516 PetscInt i, j, n, row, col, bs, mbs; in MatGetRowMaxAbs_SeqSBAIJ() local
1530 mbs = a->mbs; in MatGetRowMaxAbs_SeqSBAIJ()
1536 for (i = 0; i < mbs; i++) { in MatGetRowMaxAbs_SeqSBAIJ()
1574 PetscInt mbs = a->mbs, i, *idx = a->j, *ii = a->i, j, *jj, n, k; in MatMatMult_SeqSBAIJ_1_Private() local
1578 for (i = 0; i < mbs; i++) { in MatMatMult_SeqSBAIJ_1_Private()
1609 PetscInt mbs = a->mbs, i, *idx = a->j, *ii = a->i, j, *jj, n, k; in MatMatMult_SeqSBAIJ_2_Private() local
1612 for (i = 0; i < mbs; i++) { in MatMatMult_SeqSBAIJ_2_Private()
1648 PetscInt mbs = a->mbs, i, *idx = a->j, *ii = a->i, j, *jj, n, k; in MatMatMult_SeqSBAIJ_3_Private() local
1651 for (i = 0; i < mbs; i++) { in MatMatMult_SeqSBAIJ_3_Private()
1690 PetscInt mbs = a->mbs, i, *idx = a->j, *ii = a->i, j, *jj, n, k; in MatMatMult_SeqSBAIJ_4_Private() local
1693 for (i = 0; i < mbs; i++) { in MatMatMult_SeqSBAIJ_4_Private()
1735 PetscInt mbs = a->mbs, i, *idx = a->j, *ii = a->i, j, *jj, n, k; in MatMatMult_SeqSBAIJ_5_Private() local
1738 for (i = 0; i < mbs; i++) { in MatMatMult_SeqSBAIJ_5_Private()
1782 PetscInt mbs, i, bs = A->rmap->bs, j, n, bs2 = a->bs2; in MatMatMultNumeric_SeqSBAIJ_SeqDense() local
1821 mbs = a->mbs; in MatMatMultNumeric_SeqSBAIJ_SeqDense()
1824 for (i = 0; i < mbs; i++) { in MatMatMultNumeric_SeqSBAIJ_SeqDense()