Lines Matching refs:mbs

41   PetscInt      i, j, n = a->mbs, nz = a->i[n], *tia, *tja, bs = A->rmap->bs, k, l, cnt;  in MatGetRowIJ_SeqSBAIJ()
322 for (i = 0; i < a->mbs; i++) { /* for row block i */ in MatView_SeqSBAIJ_ASCII()
354 for (i = 0; i < a->mbs; i++) { /* for row block i */ in MatView_SeqSBAIJ_ASCII()
387 PetscInt row, i, j, k, l, mbs = a->mbs, bs = A->rmap->bs, bs2 = a->bs2; in MatView_SeqSBAIJ_Draw_Zoom() local
403 for (i = 0, row = 0; i < mbs; i++, row += bs) { in MatView_SeqSBAIJ_Draw_Zoom()
419 for (i = 0, row = 0; i < mbs; i++, row += bs) { in MatView_SeqSBAIJ_Draw_Zoom()
435 for (i = 0, row = 0; i < mbs; i++, row += bs) { in MatView_SeqSBAIJ_Draw_Zoom()
591 …heck(row < a->mbs, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Block index row too large %" PetscI… in MatSetValuesBlocked_SeqSBAIJ()
648 …MatSeqXAIJReallocateAIJ(A, a->mbs, bs2, nrow, row, col, rmax, aa, ai, aj, rp, ap, imax, nonew, Mat… in MatSetValuesBlocked_SeqSBAIJ()
679 PetscInt mbs = a->mbs, bs2 = a->bs2, rmax = 0; in MatAssemblyEnd_SeqSBAIJ() local
686 for (i = 1; i < mbs; i++) { in MatAssemblyEnd_SeqSBAIJ()
699 if (mbs) { in MatAssemblyEnd_SeqSBAIJ()
700 fshift += imax[mbs - 1] - ailen[mbs - 1]; in MatAssemblyEnd_SeqSBAIJ()
701 ai[mbs] = ai[mbs - 1] + ailen[mbs - 1]; in MatAssemblyEnd_SeqSBAIJ()
704 for (i = 0; i < mbs; i++) ailen[i] = imax[i] = ai[i + 1] - ai[i]; in MatAssemblyEnd_SeqSBAIJ()
705 a->nz = ai[mbs]; in MatAssemblyEnd_SeqSBAIJ()
804 …MatSeqXAIJReallocateAIJ(A, a->mbs, bs2, nrow, brow, bcol, rmax, aa, ai, aj, rp, ap, imax, nonew, M… in MatSetValues_SeqSBAIJ()
920 …PetscCheck(a->i[a->mbs] == b->i[b->mbs], PETSC_COMM_SELF, PETSC_ERR_ARG_INCOMP, "Number of nonzero… in MatCopy_SeqSBAIJ()
921 …PetscCheck(a->mbs == b->mbs, PETSC_COMM_SELF, PETSC_ERR_ARG_INCOMP, "Number of rows in two matrice… in MatCopy_SeqSBAIJ()
923 PetscCall(PetscArraycpy(b->a, a->a, a->bs2 * a->i[a->mbs])); in MatCopy_SeqSBAIJ()
951 PetscInt bs = Y->rmap->bs, mbs = Y->rmap->N / bs; in MatAXPYGetPreallocation_SeqSBAIJ() local
957 PetscCall(MatAXPYGetPreallocation_SeqX_private(mbs, x->i, x->j, y->i, y->j, nnz)); in MatAXPYGetPreallocation_SeqSBAIJ()
969 PetscBool e = x->nz == y->nz && x->mbs == y->mbs ? PETSC_TRUE : PETSC_FALSE; in MatAXPY_SeqSBAIJ()
971 PetscCall(PetscArraycmp(x->i, y->i, x->mbs + 1, &e)); in MatAXPY_SeqSBAIJ()
973 PetscCall(PetscArraycmp(x->j, y->j, x->i[x->mbs], &e)); in MatAXPY_SeqSBAIJ()
1024 PetscInt i, nz = a->bs2 * a->i[a->mbs]; in MatConjugate_SeqSBAIJ()
1035 PetscInt i, nz = a->bs2 * a->i[a->mbs]; in MatRealPart_SeqSBAIJ()
1046 PetscInt i, nz = a->bs2 * a->i[a->mbs]; in MatImaginaryPart_SeqSBAIJ()
1145 PetscInt mbs = a->mbs, bs2 = a->bs2, rmax = 0; in MatEliminateZeros_SeqSBAIJ() local
1152 for (i = 1; i <= mbs; i++) { in MatEliminateZeros_SeqSBAIJ()
1173 if (mbs) { in MatEliminateZeros_SeqSBAIJ()
1174 ai[mbs] -= fshift; in MatEliminateZeros_SeqSBAIJ()
1175 a->nz = ai[mbs]; in MatEliminateZeros_SeqSBAIJ()
1365 PetscInt i, mbs, nbs, bs2; in MatSeqSBAIJSetPreallocation_SeqSBAIJ() local
1389 mbs = B->rmap->N / bs; in MatSeqSBAIJSetPreallocation_SeqSBAIJ()
1393 …PetscCheck(mbs * bs == B->rmap->N && nbs * bs == B->cmap->n, PETSC_COMM_SELF, PETSC_ERR_ARG_SIZ, "… in MatSeqSBAIJSetPreallocation_SeqSBAIJ()
1403 for (i = 0; i < mbs; i++) { in MatSeqSBAIJSetPreallocation_SeqSBAIJ()
1462 b->mbs = mbs; in MatSeqSBAIJSetPreallocation_SeqSBAIJ()
1466 PetscCall(PetscMalloc2(mbs, &b->imax, mbs, &b->ilen)); in MatSeqSBAIJSetPreallocation_SeqSBAIJ()
1474 for (i = 0; i < mbs; i++) b->imax[i] = nz; in MatSeqSBAIJSetPreallocation_SeqSBAIJ()
1475 PetscCall(PetscIntMultError(nz, mbs, &nz)); in MatSeqSBAIJSetPreallocation_SeqSBAIJ()
1478 for (i = 0; i < mbs; i++) { in MatSeqSBAIJSetPreallocation_SeqSBAIJ()
1485 for (i = 0; i < mbs; i++) b->ilen[i] = 0; in MatSeqSBAIJSetPreallocation_SeqSBAIJ()
1502 for (i = 1; i < mbs + 1; i++) b->i[i] = b->i[i - 1] + b->imax[i - 1]; in MatSeqSBAIJSetPreallocation_SeqSBAIJ()
1983 PetscInt i, mbs = a->mbs, nz = a->nz, bs2 = a->bs2; in MatDuplicate_SeqSBAIJ() local
1987 PetscCheck(a->i[mbs] == nz, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Corrupt matrix"); in MatDuplicate_SeqSBAIJ()
2007 c->mbs = a->mbs; in MatDuplicate_SeqSBAIJ()
2015 PetscCall(PetscMalloc2(mbs + 1, &c->imax, mbs + 1, &c->ilen)); in MatDuplicate_SeqSBAIJ()
2016 for (i = 0; i < mbs; i++) { in MatDuplicate_SeqSBAIJ()
2037 PetscCall(PetscShmgetAllocateArray(mbs + 1, sizeof(PetscInt), (void **)&c->i)); in MatDuplicate_SeqSBAIJ()
2038 PetscCall(PetscArraycpy(c->i, a->i, mbs + 1)); in MatDuplicate_SeqSBAIJ()
2041 if (mbs > 0) { in MatDuplicate_SeqSBAIJ()