Lines Matching refs:B
19 Mat B = sell->B, Bnew; in MatDisAssemble_MPISELL() local
20 Mat_SeqSELL *Bsell = (Mat_SeqSELL *)B->data; in MatDisAssemble_MPISELL()
37 PetscCall(MatAssemblyBegin(B, MAT_FINAL_ASSEMBLY)); in MatDisAssemble_MPISELL()
38 PetscCall(MatAssemblyEnd(B, MAT_FINAL_ASSEMBLY)); in MatDisAssemble_MPISELL()
42 PetscCall(MatSetSizes(Bnew, B->rmap->n, N, B->rmap->n, N)); in MatDisAssemble_MPISELL()
44 PetscCall(MatSetType(Bnew, ((PetscObject)B)->type_name)); in MatDisAssemble_MPISELL()
54 Bnew->nonzerostate = B->nonzerostate; in MatDisAssemble_MPISELL()
56 totalslices = PetscCeilInt(B->rmap->n, Bsell->sliceheight); in MatDisAssemble_MPISELL()
60 …Bnew, Bsell->sliceheight * i + row, sell->garray[Bsell->colidx[j]], Bsell->val[j], B->insertmode)); in MatDisAssemble_MPISELL()
65 PetscCall(MatDestroy(&B)); in MatDisAssemble_MPISELL()
67 sell->B = Bnew; in MatDisAssemble_MPISELL()
76 Mat_SeqSELL *B = (Mat_SeqSELL *)sell->B->data; in MatSetUpMultiply_MPISELL() local
77 PetscInt i, j, *bcolidx = B->colidx, ec = 0, *garray, totalslices; in MatSetUpMultiply_MPISELL()
90 totalslices = PetscCeilInt(sell->B->rmap->n, B->sliceheight); in MatSetUpMultiply_MPISELL()
95 PetscCall(PetscHMapICreateWithSize(sell->B->rmap->n, &gid1_lid1)); in MatSetUpMultiply_MPISELL()
97 for (j = B->sliidx[i]; j < B->sliidx[i + 1]; j++) { in MatSetUpMultiply_MPISELL()
98 …isnonzero = (PetscBool)((j - B->sliidx[i]) / B->sliceheight < B->rlen[i * B->sliceheight + j % B->… in MatSetUpMultiply_MPISELL()
126 for (j = B->sliidx[i]; j < B->sliidx[i + 1]; j++) { in MatSetUpMultiply_MPISELL()
127 …isnonzero = (PetscBool)((j - B->sliidx[i]) / B->sliceheight < B->rlen[i * B->sliceheight + j % B->… in MatSetUpMultiply_MPISELL()
136 PetscCall(PetscLayoutDestroy(&sell->B->cmap)); in MatSetUpMultiply_MPISELL()
137 …PetscCall(PetscLayoutCreateFromSizes(PetscObjectComm((PetscObject)sell->B), ec, ec, 1, &sell->B->c… in MatSetUpMultiply_MPISELL()
144 for (j = B->sliidx[i]; j < B->sliidx[i + 1]; j++) { in MatSetUpMultiply_MPISELL()
145 …isnonzero = (PetscBool)((j - B->sliidx[i]) / B->sliceheight < B->rlen[i * B->sliceheight + j % B->… in MatSetUpMultiply_MPISELL()
165 for (j = B->sliidx[i]; j < B->sliidx[i + 1]; j++) { in MatSetUpMultiply_MPISELL()
166 …isnonzero = (PetscBool)((j - B->sliidx[i]) / B->sliceheight < B->rlen[i * B->sliceheight + j % B->… in MatSetUpMultiply_MPISELL()
170 PetscCall(PetscLayoutDestroy(&sell->B->cmap)); in MatSetUpMultiply_MPISELL()
171 …PetscCall(PetscLayoutCreateFromSizes(PetscObjectComm((PetscObject)sell->B), ec, ec, 1, &sell->B->c… in MatSetUpMultiply_MPISELL()
225 for (i = 0; i < ina->B->cmap->n; i++) lindices[garray[i]] = i + 1; in MatMPISELLDiagonalScaleLocalSetUp()
268 PetscCall(MatDiagonalScale(a->B, NULL, auglyoo)); in MatDiagonalScaleLocal_MPISELL()