Lines Matching refs:tia
41 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()
88 nz = tia[A->rmap->n / bs]; in MatGetRowIJ_SeqSBAIJ()
90 for (i = 0; i < A->rmap->n / bs + 1; i++) tia[i] = tia[i] + 1; in MatGetRowIJ_SeqSBAIJ()
91 *ia = tia; in MatGetRowIJ_SeqSBAIJ()
107 *ia = tia; in MatGetRowIJ_SeqSBAIJ()