78 PetscErrorCode MatSeqSELLSetPreallocation(Mat B, PetscInt rlenmax, const PetscInt rlen[])  in MatSeqSELLSetPreallocation()  argument
83 …ryMethod(B, "MatSeqSELLSetPreallocation_C", (Mat, PetscInt, const PetscInt[]), (B, rlenmax, rlen));  in MatSeqSELLSetPreallocation()
87 PetscErrorCode MatSeqSELLSetPreallocation_SeqSELL(Mat B, PetscInt maxallocrow, const PetscInt rlen[…  in MatSeqSELLSetPreallocation_SeqSELL()  argument
97   if (maxallocrow >= 0 || rlen) realalloc = PETSC_TRUE;  in MatSeqSELLSetPreallocation_SeqSELL()
109   if (rlen) {  in MatSeqSELLSetPreallocation_SeqSELL()
111 …eck(rlen[i] >= 0, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "rlen cannot be less than 0: local ro…  in MatSeqSELLSetPreallocation_SeqSELL()
112 …rlen[i] <= B->cmap->n, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "rlen cannot be greater than row…  in MatSeqSELLSetPreallocation_SeqSELL()
135     if (!rlen) { /* if rlen is not provided, allocate same space for all the slices */  in MatSeqSELLSetPreallocation_SeqSELL()
152 …for (j = 0; j < b->sliceheight; j++) b->sliidx[i] = PetscMax(b->sliidx[i], rlen[b->sliceheight * (…  in MatSeqSELLSetPreallocation_SeqSELL()
164 …ices - 1); j < B->rmap->n; j++) b->sliidx[totalslices] = PetscMax(b->sliidx[totalslices], rlen[j]);  in MatSeqSELLSetPreallocation_SeqSELL()
181     PetscCall(PetscCalloc1(b->sliceheight * totalslices, &b->rlen));  in MatSeqSELLSetPreallocation_SeqSELL()
211   if (nz) *nz = a->rlen[row];  in MatGetRow_SeqSELL()
216     for (j = 0; j < a->rlen[row]; j++) a->getrowcols[j] = a->colidx[shift + a->sliceheight * j];  in MatGetRow_SeqSELL()
221     for (j = 0; j < a->rlen[row]; j++) a->getrowvals[j] = a->val[shift + a->sliceheight * j];  in MatGetRow_SeqSELL()
247     PetscCall(MatSeqAIJSetPreallocation(B, 0, a->rlen));  in MatConvert_SeqSELL_SeqAIJ()
436         nnz_in_row = a->rlen[row];  in MatMult_SeqSELL()
476         nnz_in_row = a->rlen[row];  in MatMult_SeqSELL()
653         nnz_in_row = a->rlen[row];  in MatMultAdd_SeqSELL()
740           nnz_in_row = a->rlen[row];  in MatMultTransposeAdd_SeqSELL()
789         for (PetscInt j = 0; j < a->rlen[i]; j++) {  in MatGetDiagonalMarkers_SeqSELL()
811         for (PetscInt j = 0; j < a->rlen[i]; j++) {  in MatGetDiagonalMarkers_SeqSELL()
897   PetscCall(PetscFree(a->rlen));  in MatDestroy_SeqSELL()
985     for (j = 0; j < a->rlen[i]; j++) {  in MatGetDiagonal_SeqSELL()
1062       high = a->rlen[row];  in MatGetValues_SeqSELL()
1109       for (j = 0; j < a->rlen[i]; j++) {  in MatView_SeqSELL_ASCII()
1136       for (j = 0; j < a->rlen[i]; j++) {  in MatView_SeqSELL_ASCII()
1170         if (jcnt < a->rlen[i] && j == a->colidx[shift + a->sliceheight * j]) {  in MatView_SeqSELL_ASCII()
1200       for (j = 0; j < a->rlen[i]; j++) {  in MatView_SeqSELL_ASCII()
1262         for (j = a->diag[i] + 1; j < shift + a->sliceheight * a->rlen[i]; j += a->sliceheight) {  in MatView_SeqSELL_ASCII()
1281         for (j = 0; j < a->rlen[i]; j++) {  in MatView_SeqSELL_ASCII()
1329       for (j = 0; j < a->rlen[i]; j++) {  in MatView_SeqSELL_Draw_Zoom()
1341       for (j = 0; j < a->rlen[i]; j++) {  in MatView_SeqSELL_Draw_Zoom()
1353       for (j = 0; j < a->rlen[i]; j++) {  in MatView_SeqSELL_Draw_Zoom()
1379       for (j = 0; j < a->rlen[i]; j++) {  in MatView_SeqSELL_Draw_Zoom()
1459       nrow = a->rlen[row]; /* number of nonzeros in row */  in MatAssemblyEnd_SeqSELL()
1470           if (a->rlen[a->sliceheight * i + j]) {  in MatAssemblyEnd_SeqSELL()
1552     nrow  = a->rlen[row];  in MatSetValues_SeqSELL()
1602       a->rlen[row]++;  in MatSetValues_SeqSELL()
1614     a->rlen[row] = nrow;  in MatSetValues_SeqSELL()
1729         n     = a->rlen[i] - (diag[i] - shift) / a->sliceheight - 1;  in MatSOR_SeqSELL()
1751         n = a->rlen[i] - (diag[i] - shift) / a->sliceheight - 1;  in MatSOR_SeqSELL()
1764           n = a->rlen[i];  in MatSOR_SeqSELL()
1768           n = a->rlen[i] - (diag[i] - shift) / a->sliceheight - 1;  in MatSOR_SeqSELL()
2238   PetscCall(PetscMalloc1(c->sliceheight * totalslices, &c->rlen));  in MatDuplicateNoCreate_SeqSELL()
2241   for (i = 0; i < m; i++) c->rlen[i] = a->rlen[i];  in MatDuplicateNoCreate_SeqSELL()
2386 …eateSeqSELL(MPI_Comm comm, PetscInt m, PetscInt n, PetscInt rlenmax, const PetscInt rlen[], Mat *A)  in MatCreateSeqSELL()  argument
2392   PetscCall(MatSeqSELLSetPreallocation_SeqSELL(*A, rlenmax, rlen));  in MatCreateSeqSELL()