Lines Matching refs:ai
58 const PetscInt *r, *ic, *ai = a->i, *aj = a->j, *ajtmp; in MatLUFactorSymbolic_SeqAIJ() local
91 PetscCall(PetscFreeSpaceGet(PetscRealIntMultTruncate(f, ai[n] + 1), &free_space)); in MatLUFactorSymbolic_SeqAIJ()
97 nnz = ai[r[i] + 1] - ai[r[i]]; in MatLUFactorSymbolic_SeqAIJ()
98 ajtmp = aj + ai[r[i]]; in MatLUFactorSymbolic_SeqAIJ()
175 if (ai[n]) { in MatLUFactorSymbolic_SeqAIJ()
176 B->info.fill_ratio_needed = ((PetscReal)(bdiag[0] + 1)) / ((PetscReal)ai[n]); in MatLUFactorSymbolic_SeqAIJ()
181 if (ai[n] != 0) { in MatLUFactorSymbolic_SeqAIJ()
224 const PetscInt n = A->rmap->n, *ai = a->i, *aj = a->j, *bi = b->i, *bj = b->j, *bdiag = b->diag; in MatLUFactorNumeric_SeqAIJ() local
249 v = aa + ai[i]; in MatLUFactorNumeric_SeqAIJ()
250 nz = ai[i + 1] - ai[i]; in MatLUFactorNumeric_SeqAIJ()
280 nz = ai[r[i] + 1] - ai[r[i]]; in MatLUFactorNumeric_SeqAIJ()
281 ajtmp = aj + ai[r[i]]; in MatLUFactorNumeric_SeqAIJ()
282 v = aa + ai[r[i]]; in MatLUFactorNumeric_SeqAIJ()
404 const PetscInt *ai = a->i, *aj = a->j, *bi = b->i, *bj = b->j; in MatLUFactorNumeric_SeqAIJ_inplace() local
428 v = aa + ai[i]; in MatLUFactorNumeric_SeqAIJ_inplace()
429 nz = ai[i + 1] - ai[i]; in MatLUFactorNumeric_SeqAIJ_inplace()
452 nz = ai[r[i] + 1] - ai[r[i]]; in MatLUFactorNumeric_SeqAIJ_inplace()
453 ajtmp = aj + ai[r[i]]; in MatLUFactorNumeric_SeqAIJ_inplace()
454 v = aa + ai[r[i]]; in MatLUFactorNumeric_SeqAIJ_inplace()
560 PetscInt i, j, n = A->rmap->n, *ai = a->i, *aj = a->j; in MatLUFactorNumeric_SeqAIJ_InplaceWithPerm() local
587 vtmp = aa + ai[i]; in MatLUFactorNumeric_SeqAIJ_InplaceWithPerm()
588 nz = ai[i + 1] - ai[i]; in MatLUFactorNumeric_SeqAIJ_InplaceWithPerm()
617 v = aa + ai[i]; in MatLUFactorNumeric_SeqAIJ_InplaceWithPerm()
618 nz = ai[i + 1] - ai[i]; in MatLUFactorNumeric_SeqAIJ_InplaceWithPerm()
637 nz = ai[r[i] + 1] - ai[r[i]]; in MatLUFactorNumeric_SeqAIJ_InplaceWithPerm()
638 ajtmp = aj + ai[r[i]]; in MatLUFactorNumeric_SeqAIJ_InplaceWithPerm()
639 v = aa + ai[r[i]]; in MatLUFactorNumeric_SeqAIJ_InplaceWithPerm()
644 diag[r[i]] = ai[r[i]]; in MatLUFactorNumeric_SeqAIJ_InplaceWithPerm()
660 nz = ai[r[row] + 1] - diag[r[row]] - 1; in MatLUFactorNumeric_SeqAIJ_InplaceWithPerm()
667 pv = aa + ai[r[i]]; in MatLUFactorNumeric_SeqAIJ_InplaceWithPerm()
668 pj = aj + ai[r[i]]; in MatLUFactorNumeric_SeqAIJ_InplaceWithPerm()
669 nz = ai[r[i] + 1] - ai[r[i]]; in MatLUFactorNumeric_SeqAIJ_InplaceWithPerm()
670 nbdiag = diag[r[i]] - ai[r[i]]; /* num of entries before the diagonal */ in MatLUFactorNumeric_SeqAIJ_InplaceWithPerm()
745 PetscInt i, n = A->rmap->n, *vi, *ai = a->i, *aj = a->j; in MatSolve_SeqAIJ_inplace() local
771 v = aa + ai[i]; in MatSolve_SeqAIJ_inplace()
772 vi = aj + ai[i]; in MatSolve_SeqAIJ_inplace()
773 nz = adiag[i] - ai[i]; in MatSolve_SeqAIJ_inplace()
783 nz = ai[i + 1] - adiag[i] - 1; in MatSolve_SeqAIJ_inplace()
801 PetscInt i, n = A->rmap->n, *vi, *ai = a->i, *aj = a->j; in MatMatSolve_SeqAIJ_inplace() local
832 v = aa + ai[i]; in MatMatSolve_SeqAIJ_inplace()
833 vi = aj + ai[i]; in MatMatSolve_SeqAIJ_inplace()
834 nz = adiag[i] - ai[i]; in MatMatSolve_SeqAIJ_inplace()
843 nz = ai[i + 1] - adiag[i] - 1; in MatMatSolve_SeqAIJ_inplace()
864 PetscInt i, n = A->rmap->n, *vi, *ai = a->i, *aj = a->j; in MatMatSolve_SeqAIJ() local
896 nz = ai[i + 1] - ai[i]; in MatMatSolve_SeqAIJ()
928 PetscInt i, n = A->rmap->n, *vi, *ai = a->i, *aj = a->j, j; in MatMatSolveTranspose_SeqAIJ() local
971 v = aa + ai[i]; in MatMatSolveTranspose_SeqAIJ()
972 vi = aj + ai[i]; in MatMatSolveTranspose_SeqAIJ()
973 nz = ai[i + 1] - ai[i]; in MatMatSolveTranspose_SeqAIJ()
997 PetscInt i, n = A->rmap->n, *vi, *ai = a->i, *aj = a->j; in MatSolve_SeqAIJ_InplaceWithPerm() local
1022 v = aa + ai[i]; in MatSolve_SeqAIJ_InplaceWithPerm()
1023 vi = aj + ai[i]; in MatSolve_SeqAIJ_InplaceWithPerm()
1024 nz = adiag[i] - ai[i]; in MatSolve_SeqAIJ_InplaceWithPerm()
1035 nz = ai[i + 1] - adiag[i] - 1; in MatSolve_SeqAIJ_InplaceWithPerm()
1054 const PetscInt *ai = a->i, *aj = a->j, *adiag; in MatSolve_SeqAIJ_NaturalOrdering_inplace() local
1073 fortransolveaij_(&n, x, ai, aj, adiag, aa, b); in MatSolve_SeqAIJ_NaturalOrdering_inplace()
1078 ai_i = ai[i]; in MatSolve_SeqAIJ_NaturalOrdering_inplace()
1092 nz = ai[i + 1] - adiag_i - 1; in MatSolve_SeqAIJ_NaturalOrdering_inplace()
1111 const PetscInt *rout, *cout, *r, *c, *vi, *ai = a->i, *aj = a->j, *adiag; in MatSolveAdd_SeqAIJ_inplace() local
1133 v = aa + ai[i]; in MatSolveAdd_SeqAIJ_inplace()
1134 vi = aj + ai[i]; in MatSolveAdd_SeqAIJ_inplace()
1135 nz = adiag[i] - ai[i]; in MatSolveAdd_SeqAIJ_inplace()
1145 nz = ai[i + 1] - adiag[i] - 1; in MatSolveAdd_SeqAIJ_inplace()
1167 const PetscInt *rout, *cout, *r, *c, *vi, *ai = a->i, *aj = a->j, *adiag; in MatSolveAdd_SeqAIJ() local
1191 nz = ai[i + 1] - ai[i]; in MatSolveAdd_SeqAIJ()
1225 const PetscInt *rout, *cout, *r, *c, *diag = a->diag, *ai = a->i, *aj = a->j, *vi; in MatSolveTranspose_SeqAIJ_inplace() local
1250 nz = ai[i + 1] - diag[i] - 1; in MatSolveTranspose_SeqAIJ_inplace()
1261 nz = diag[i] - ai[i]; in MatSolveTranspose_SeqAIJ_inplace()
1283 const PetscInt *rout, *cout, *r, *c, *adiag, *ai = a->i, *aj = a->j, *vi; in MatSolveTranspose_SeqAIJ() local
1318 v = aa + ai[i]; in MatSolveTranspose_SeqAIJ()
1319 vi = aj + ai[i]; in MatSolveTranspose_SeqAIJ()
1320 nz = ai[i + 1] - ai[i]; in MatSolveTranspose_SeqAIJ()
1342 const PetscInt *rout, *cout, *r, *c, *diag = a->diag, *ai = a->i, *aj = a->j, *vi; in MatSolveTransposeAdd_SeqAIJ_inplace() local
1368 nz = ai[i + 1] - diag[i] - 1; in MatSolveTransposeAdd_SeqAIJ_inplace()
1379 nz = diag[i] - ai[i]; in MatSolveTransposeAdd_SeqAIJ_inplace()
1401 const PetscInt *rout, *cout, *r, *c, *adiag, *ai = a->i, *aj = a->j, *vi; in MatSolveTransposeAdd_SeqAIJ() local
1437 v = aa + ai[i]; in MatSolveTransposeAdd_SeqAIJ()
1438 vi = aj + ai[i]; in MatSolveTransposeAdd_SeqAIJ()
1439 nz = ai[i + 1] - ai[i]; in MatSolveTransposeAdd_SeqAIJ()
1476 const PetscInt n = A->rmap->n, *ai = a->i, *aj, *adiag; in MatILUFactorSymbolic_SeqAIJ_ilu0() local
1487 PetscCall(PetscShmgetAllocateArray(ai[n], sizeof(PetscScalar), (void **)&b->a)); in MatILUFactorSymbolic_SeqAIJ_ilu0()
1488 PetscCall(PetscShmgetAllocateArray(ai[n], sizeof(PetscInt), (void **)&b->j)); in MatILUFactorSymbolic_SeqAIJ_ilu0()
1490 if (n > 0) PetscCall(PetscArrayzero(b->a, ai[n])); in MatILUFactorSymbolic_SeqAIJ_ilu0()
1504 nz = adiag[i] - ai[i]; in MatILUFactorSymbolic_SeqAIJ_ilu0()
1506 aj = a->j + ai[i]; in MatILUFactorSymbolic_SeqAIJ_ilu0()
1513 nz = ai[i + 1] - adiag[i] - 1; in MatILUFactorSymbolic_SeqAIJ_ilu0()
1543 PetscInt n = A->rmap->n, *ai = a->i, *aj = a->j; in MatILUFactorSymbolic_SeqAIJ() local
1588 PetscCall(PetscFreeSpaceGet(PetscRealIntMultTruncate(f, ai[n] + 1), &free_space)); in MatILUFactorSymbolic_SeqAIJ()
1590 PetscCall(PetscFreeSpaceGet(PetscRealIntMultTruncate(f, ai[n] + 1), &free_space_lvl)); in MatILUFactorSymbolic_SeqAIJ()
1595 nnz = ai[r[i] + 1] - ai[r[i]]; in MatILUFactorSymbolic_SeqAIJ()
1597 cols = aj + ai[r[i]]; in MatILUFactorSymbolic_SeqAIJ()
1664 PetscReal af = ((PetscReal)(bdiag[0] + 1)) / ((PetscReal)ai[n]); in MatILUFactorSymbolic_SeqAIJ()
1696 fact->info.fill_ratio_needed = ((PetscReal)(bdiag[0] + 1)) / ((PetscReal)ai[n]); in MatILUFactorSymbolic_SeqAIJ()
1711 PetscInt *ai = a->i, *aj = a->j; in MatCholeskyFactorNumeric_SeqAIJ() local
1733 v = aa + ai[i]; in MatCholeskyFactorNumeric_SeqAIJ()
1734 nz = ai[i + 1] - ai[i]; in MatCholeskyFactorNumeric_SeqAIJ()
1767 jmin = ai[rip[k]]; in MatCholeskyFactorNumeric_SeqAIJ()
1768 jmax = ai[rip[k] + 1]; in MatCholeskyFactorNumeric_SeqAIJ()
1878 PetscInt *ai = a->i, *aj = a->j; in MatCholeskyFactorNumeric_SeqAIJ_inplace() local
1900 v = aa + ai[i]; in MatCholeskyFactorNumeric_SeqAIJ_inplace()
1901 nz = ai[i + 1] - ai[i]; in MatCholeskyFactorNumeric_SeqAIJ_inplace()
1931 jmin = ai[rip[k]]; in MatCholeskyFactorNumeric_SeqAIJ_inplace()
1932 jmax = ai[rip[k] + 1]; in MatCholeskyFactorNumeric_SeqAIJ_inplace()
2056 PetscInt reallocs = 0, i, *ai = a->i, *aj = a->j, am = A->rmap->n, *ui, *udiag; in MatICCFactorSymbolic_SeqAIJ() local
2082 ncols = ai[i + 1] - adiag[i]; in MatICCFactorSymbolic_SeqAIJ()
2090 ncols = ai[i + 1] - adiag[i] - 1; in MatICCFactorSymbolic_SeqAIJ()
2114 PetscCall(PetscFreeSpaceGet(PetscRealIntMultTruncate(fill, (ai[am] + am) / 2), &free_space)); in MatICCFactorSymbolic_SeqAIJ()
2116 … PetscCall(PetscFreeSpaceGet(PetscRealIntMultTruncate(fill, (ai[am] + am) / 2), &free_space_lvl)); in MatICCFactorSymbolic_SeqAIJ()
2122 ncols = ai[rip[k] + 1] - ai[rip[k]]; in MatICCFactorSymbolic_SeqAIJ()
2126 i = *(aj + ai[rip[k]] + j); /* unpermuted column index */ in MatICCFactorSymbolic_SeqAIJ()
2232 if (ai[am] != 0) { in MatICCFactorSymbolic_SeqAIJ()
2234 fact->info.fill_ratio_needed = ((PetscReal)2 * ui[am]) / (ai[am] + am); in MatICCFactorSymbolic_SeqAIJ()
2239 if (ai[am] != 0) { in MatICCFactorSymbolic_SeqAIJ()
2259 PetscInt i, am = A->rmap->n, *ai = a->i, *aj = a->j, reallocs = 0, prow; in MatCholeskyFactorSymbolic_SeqAIJ() local
2296 PetscCall(PetscFreeSpaceGet(PetscRealIntMultTruncate(fill, (ai[am] + am) / 2), &free_space)); in MatCholeskyFactorSymbolic_SeqAIJ()
2302 ncols = ai[rip[k] + 1] - ai[rip[k]]; in MatCholeskyFactorSymbolic_SeqAIJ()
2306 i = riip[*(aj + ai[rip[k]] + j)]; in MatCholeskyFactorSymbolic_SeqAIJ()
2399 if (ai[am] != 0) { in MatCholeskyFactorSymbolic_SeqAIJ()
2401 fact->info.fill_ratio_needed = ((PetscReal)2 * ui[am]) / (ai[am] + am); in MatCholeskyFactorSymbolic_SeqAIJ()
2406 if (ai[am] != 0) { in MatCholeskyFactorSymbolic_SeqAIJ()
2423 const PetscInt *ai = a->i, *aj = a->j, *adiag = a->diag, *vi; in MatSolve_SeqAIJ_NaturalOrdering() local
2441 nz = ai[i + 1] - ai[i]; in MatSolve_SeqAIJ_NaturalOrdering()
2470 PetscInt i, n = A->rmap->n, *vi, *ai = a->i, *aj = a->j, *adiag = a->diag, nz; in MatSolve_SeqAIJ() local
2494 nz = ai[i + 1] - ai[i]; in MatSolve_SeqAIJ()