Searched refs:rowlengths (Results 1 – 4 of 4) sorted by relevance
| /petsc/src/mat/impls/sbaij/seq/ |
| H A D | aijsbaij.c | 10 … a->i, *aj = a->j, m = A->rmap->N, n = A->cmap->n, i, j, k, *bi, *bj, *rowlengths, nz, *rowstart, … in MatConvert_SeqSBAIJ_SeqAIJ() local
17 PetscCall(PetscMalloc2(m, &rowlengths, m + 1, &rowstart)); in MatConvert_SeqSBAIJ_SeqAIJ()
19 for (i = 0; i < mbs; i++) rowlengths[i * bs] = 0; in MatConvert_SeqSBAIJ_SeqAIJ()
24 rowlengths[k] += nz; /* no. of upper triangular blocks */ in MatConvert_SeqSBAIJ_SeqAIJ()
31 rowlengths[(*aj) * bs]++; in MatConvert_SeqSBAIJ_SeqAIJ()
35 rowlengths[k] *= bs; in MatConvert_SeqSBAIJ_SeqAIJ()
36 for (j = 1; j < bs; j++) rowlengths[k + j] = rowlengths[k]; in MatConvert_SeqSBAIJ_SeqAIJ()
44 PetscCall(MatSeqAIJSetPreallocation(B, 0, rowlengths)); in MatConvert_SeqSBAIJ_SeqAIJ()
58 b->ilen[i * bs + j] = rowlengths[i * bs]; in MatConvert_SeqSBAIJ_SeqAIJ()
59 rowstart[i * bs + j + 1] = rowstart[i * bs + j] + rowlengths[i * bs]; in MatConvert_SeqSBAIJ_SeqAIJ()
[all …]
|
| /petsc/src/mat/impls/baij/seq/ |
| H A D | aijbaij.c | 10 PetscInt *rowlengths, *rows, *cols, maxlen = 0, ncols; in MatConvert_SeqBAIJ_SeqAIJ() local
18 PetscCall(PetscMalloc1(n * bs, &rowlengths)); in MatConvert_SeqBAIJ_SeqAIJ()
21 for (j = 0; j < bs; j++) rowlengths[i * bs + j] = bs * (ai[i + 1] - ai[i]); in MatConvert_SeqBAIJ_SeqAIJ()
27 PetscCall(MatSeqAIJSetPreallocation(B, 0, rowlengths)); in MatConvert_SeqBAIJ_SeqAIJ()
28 PetscCall(PetscFree(rowlengths)); in MatConvert_SeqBAIJ_SeqAIJ()
97 PetscInt m = A->rmap->N, n = A->cmap->N, *rowlengths, bs = A->rmap->bs; in MatConvert_SeqAIJ_SeqBAIJ() local
101 PetscCall(MatConvert_SeqAIJ_SeqBAIJ_Preallocate(A, &rowlengths)); in MatConvert_SeqAIJ_SeqBAIJ()
105 PetscCall(MatSeqBAIJSetPreallocation(B, bs, 0, rowlengths)); in MatConvert_SeqAIJ_SeqBAIJ()
106 PetscCall(PetscFree(rowlengths)); in MatConvert_SeqAIJ_SeqBAIJ()
|
| /petsc/src/mat/impls/sell/seq/ |
| H A D | sell.c | 277 PetscInt *ai = a->i, m = A->rmap->N, n = A->cmap->N, i, *rowlengths, row, ncols; in MatConvert_SeqAIJ_SeqSELL() local
286 PetscCall(PetscMalloc1(m, &rowlengths)); in MatConvert_SeqAIJ_SeqSELL()
287 for (i = 0; i < m; i++) rowlengths[i] = ai[i + 1] - ai[i]; in MatConvert_SeqAIJ_SeqSELL()
291 PetscCall(PetscArraycmp(rowlengths, a->ilen, m, &eq)); in MatConvert_SeqAIJ_SeqSELL()
293 PetscCall(PetscFree(rowlengths)); in MatConvert_SeqAIJ_SeqSELL()
294 rowlengths = a->ilen; in MatConvert_SeqAIJ_SeqSELL()
295 } else if (a->ilen) rowlengths = a->ilen; in MatConvert_SeqAIJ_SeqSELL()
299 PetscCall(MatSeqSELLSetPreallocation(B, 0, rowlengths)); in MatConvert_SeqAIJ_SeqSELL()
300 if (rowlengths != a->ilen) PetscCall(PetscFree(rowlengths)); in MatConvert_SeqAIJ_SeqSELL()
|
| /petsc/share/petsc/matlab/ |
| H A D | PetscBinaryRead.m | 111 str = sprintf('No-Nonzeros sum-rowlengths do not match %d %d',nz,sum_nz);
|