Home
last modified time | relevance | path

Searched refs:tia (Results 1 – 3 of 3) sorted by relevance

/petsc/src/mat/impls/sbaij/seq/
H A Dsbaij.c41 PetscInt i, j, n = a->mbs, nz = a->i[n], *tia, *tja, bs = A->rmap->bs, k, l, cnt; in MatGetRowIJ_SeqSBAIJ() local
48 PetscCall(MatToSymmetricIJ_SeqAIJ(n, a->i, a->j, PETSC_FALSE, 0, 0, &tia, &tja)); in MatGetRowIJ_SeqSBAIJ()
49 nz = tia[n]; in MatGetRowIJ_SeqSBAIJ()
51 tia = a->i; in MatGetRowIJ_SeqSBAIJ()
61 for (j = 1; j < bs; j++) (*ia)[j] = (tia[1] - tia[0]) * bs + (*ia)[j - 1]; in MatGetRowIJ_SeqSBAIJ()
65 (*ia)[i * bs] = (tia[i] - tia[i - 1]) * bs + (*ia)[i * bs - 1]; in MatGetRowIJ_SeqSBAIJ()
66 … for (j = 1; j < bs; j++) (*ia)[i * bs + j] = (tia[i + 1] - tia[i]) * bs + (*ia)[i * bs + j - 1]; in MatGetRowIJ_SeqSBAIJ()
68 if (n) (*ia)[n * bs] = (tia[n] - tia[n - 1]) * bs + (*ia)[n * bs - 1]; in MatGetRowIJ_SeqSBAIJ()
75 for (k = tia[i]; k < tia[i + 1]; k++) { in MatGetRowIJ_SeqSBAIJ()
83 PetscCall(PetscFree(tia)); in MatGetRowIJ_SeqSBAIJ()
[all …]
/petsc/src/mat/impls/baij/seq/
H A Dbaij.c1416 PetscInt i, j, n = a->mbs, nz = a->i[n], *tia, *tja, bs = A->rmap->bs, k, l, cnt; in MatGetRowIJ_SeqBAIJ() local
1423 PetscCall(MatToSymmetricIJ_SeqAIJ(n, a->i, a->j, PETSC_TRUE, 0, 0, &tia, &tja)); in MatGetRowIJ_SeqBAIJ()
1424 nz = tia[n]; in MatGetRowIJ_SeqBAIJ()
1426 tia = a->i; in MatGetRowIJ_SeqBAIJ()
1436 for (j = 1; j < bs; j++) (*ia)[j] = (tia[1] - tia[0]) * bs + (*ia)[j - 1]; in MatGetRowIJ_SeqBAIJ()
1440 (*ia)[i * bs] = (tia[i] - tia[i - 1]) * bs + (*ia)[i * bs - 1]; in MatGetRowIJ_SeqBAIJ()
1441 … for (j = 1; j < bs; j++) (*ia)[i * bs + j] = (tia[i + 1] - tia[i]) * bs + (*ia)[i * bs + j - 1]; in MatGetRowIJ_SeqBAIJ()
1443 if (n) (*ia)[n * bs] = (tia[n] - tia[n - 1]) * bs + (*ia)[n * bs - 1]; in MatGetRowIJ_SeqBAIJ()
1450 for (k = tia[i]; k < tia[i + 1]; k++) { in MatGetRowIJ_SeqBAIJ()
1458 PetscCall(PetscFree(tia)); in MatGetRowIJ_SeqBAIJ()
[all …]
/petsc/src/mat/impls/aij/seq/
H A Daij.c206 PetscInt *tia; in MatGetRowIJ_SeqAIJ() local
210 PetscCall(PetscMalloc1(A->rmap->n + 1, &tia)); in MatGetRowIJ_SeqAIJ()
211 for (i = 0; i < A->rmap->n + 1; i++) tia[i] = a->i[i] + 1; in MatGetRowIJ_SeqAIJ()
212 *ia = tia; in MatGetRowIJ_SeqAIJ()