Lines Matching refs:ai
42 const PetscInt *rip, *ai, *aj; in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR() local
57 ai = a->i; in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
64 ai = a->inew; in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
70 umax = (PetscInt)(f * ai[mbs] + 1); in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
88 jmin = ai[rip[k]] + 1; /* exclude diag[k] */ in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
89 jmax = ai[rip[k] + 1]; in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
165 if (ai[mbs] != 0) { in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
166 PetscReal af = ((PetscReal)iu[mbs]) / ((PetscReal)ai[mbs]); in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
206 if (ai[mbs] != 0) { in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
207 F->info.fill_ratio_needed = ((PetscReal)iu[mbs]) / ((PetscReal)ai[mbs]); in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
226 const PetscInt *rip, *ai = a->i, *aj = a->j; in MatCholeskyFactorSymbolic_SeqSBAIJ() local
267 PetscCall(PetscFreeSpaceGet(PetscRealIntMultTruncate(fill, ai[mbs] + 1), &free_space)); in MatCholeskyFactorSymbolic_SeqSBAIJ()
273 ncols = ai[k + 1] - ai[k]; in MatCholeskyFactorSymbolic_SeqSBAIJ()
276 i = *(aj + ai[k] + j); in MatCholeskyFactorSymbolic_SeqSBAIJ()
366 if (ai[mbs] != 0) { in MatCholeskyFactorSymbolic_SeqSBAIJ()
367 fact->info.fill_ratio_needed = ((PetscReal)ui[mbs]) / ai[mbs]; in MatCholeskyFactorSymbolic_SeqSBAIJ()
372 if (ai[mbs] != 0) { in MatCholeskyFactorSymbolic_SeqSBAIJ()
391 const PetscInt *rip, *ai, *aj; in MatCholeskyFactorSymbolic_SeqSBAIJ_inplace() local
421 ai = a->i; in MatCholeskyFactorSymbolic_SeqSBAIJ_inplace()
442 PetscCall(PetscFreeSpaceGet(PetscRealIntMultTruncate(fill, ai[mbs] + 1), &free_space)); in MatCholeskyFactorSymbolic_SeqSBAIJ_inplace()
448 ncols = ai[rip[k] + 1] - ai[rip[k]]; in MatCholeskyFactorSymbolic_SeqSBAIJ_inplace()
450 i = *(aj + ai[rip[k]] + j); in MatCholeskyFactorSymbolic_SeqSBAIJ_inplace()
539 if (ai[mbs] != 0) { in MatCholeskyFactorSymbolic_SeqSBAIJ_inplace()
540 fact->info.fill_ratio_needed = ((PetscReal)ui[mbs]) / ai[mbs]; in MatCholeskyFactorSymbolic_SeqSBAIJ_inplace()
545 if (ai[mbs] != 0) { in MatCholeskyFactorSymbolic_SeqSBAIJ_inplace()
562 const PetscInt *ai, *aj, *perm_ptr, mbs = a->mbs, *bi = b->i, *bj = b->j; in MatCholeskyFactorNumeric_SeqSBAIJ_N() local
588 ai = a->i; in MatCholeskyFactorNumeric_SeqSBAIJ_N()
592 ai = a->inew; in MatCholeskyFactorNumeric_SeqSBAIJ_N()
594 PetscCall(PetscMalloc1(bs2 * ai[mbs], &aa)); in MatCholeskyFactorNumeric_SeqSBAIJ_N()
595 PetscCall(PetscArraycpy(aa, a->a, bs2 * ai[mbs])); in MatCholeskyFactorNumeric_SeqSBAIJ_N()
596 PetscCall(PetscMalloc1(ai[mbs], &a2anew)); in MatCholeskyFactorNumeric_SeqSBAIJ_N()
597 PetscCall(PetscArraycpy(a2anew, a->a2anew, ai[mbs])); in MatCholeskyFactorNumeric_SeqSBAIJ_N()
600 jmin = ai[i]; in MatCholeskyFactorNumeric_SeqSBAIJ_N()
601 jmax = ai[i + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_N()
629 jmin = ai[perm_ptr[k]]; in MatCholeskyFactorNumeric_SeqSBAIJ_N()
630 jmax = ai[perm_ptr[k] + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_N()
737 PetscInt *ai, *aj, k, k1, jmin, jmax, *jl, *il, vj, nexti, ili; in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering() local
755 ai = a->i; in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering()
762 jmin = ai[k]; in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering()
763 jmax = ai[k + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering()
869 const PetscInt *ai, *aj, *perm_ptr; in MatCholeskyFactorNumeric_SeqSBAIJ_2() local
898 ai = a->i; in MatCholeskyFactorNumeric_SeqSBAIJ_2()
902 ai = a->inew; in MatCholeskyFactorNumeric_SeqSBAIJ_2()
904 PetscCall(PetscMalloc1(4 * ai[mbs], &aa)); in MatCholeskyFactorNumeric_SeqSBAIJ_2()
905 PetscCall(PetscArraycpy(aa, a->a, 4 * ai[mbs])); in MatCholeskyFactorNumeric_SeqSBAIJ_2()
906 PetscCall(PetscMalloc1(ai[mbs], &a2anew)); in MatCholeskyFactorNumeric_SeqSBAIJ_2()
907 PetscCall(PetscArraycpy(a2anew, a->a2anew, ai[mbs])); in MatCholeskyFactorNumeric_SeqSBAIJ_2()
910 jmin = ai[i]; in MatCholeskyFactorNumeric_SeqSBAIJ_2()
911 jmax = ai[i + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_2()
938 jmin = ai[perm_ptr[k]]; in MatCholeskyFactorNumeric_SeqSBAIJ_2()
939 jmax = ai[perm_ptr[k] + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_2()
1050 PetscInt *ai, *aj, k, k1, jmin, jmax, *jl, *il, vj, nexti, ili; in MatCholeskyFactorNumeric_SeqSBAIJ_2_NaturalOrdering() local
1073 ai = a->i; in MatCholeskyFactorNumeric_SeqSBAIJ_2_NaturalOrdering()
1080 jmin = ai[k]; in MatCholeskyFactorNumeric_SeqSBAIJ_2_NaturalOrdering()
1081 jmax = ai[k + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_2_NaturalOrdering()
1193 const PetscInt *ai, *aj, *rip; in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace() local
1206 ai = a->i; in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace()
1210 ai = a->inew; in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace()
1212 nz = ai[mbs]; in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace()
1240 jmin = ai[rip[k]]; in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace()
1241 jmax = ai[rip[k] + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace()
1351 PetscInt *ai = a->i, *aj = a->j, *ajtmp; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering() local
1375 ajtmp = aj + ai[i] + 1; /* exclude diagonal */ in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1376 v = aa + ai[i] + 1; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1377 nz = ai[i + 1] - ai[i] - 1; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1408 jmin = ai[k]; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1409 jmax = ai[k + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1503 PetscInt *ai = a->i, *aj = a->j, *bi = b->i, *bj = b->j; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering_inplace() local
1535 nz = ai[k + 1] - ai[k]; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering_inplace()
1536 acol = aj + ai[k]; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering_inplace()
1537 aval = aa + ai[k]; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering_inplace()