| /petsc/src/mat/impls/sbaij/seq/ |
| H A D | sbaijfact.c | 45 PetscInt *jl, *q, jmin, jmax, juidx, nzk, qm, *iu, *ju, k, j, vj, umax, maxadd; in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR() local
85 …nzk = 0; /* num. of nz blocks in k-th block row with diagonal block exc… in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
99 nzk++; in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
122 nzk++; in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
132 if (nzk > 0) { in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
137 iu[k + 1] = iu[k] + nzk; in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
145 if (maxadd < nzk) maxadd = (mbs - k) * (nzk + 1) / 2; in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
158 while (nzk--) { in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
228 PetscInt *jl, jmin, jmax, nzk, *ui, k, j, *il, nextprow; in MatCholeskyFactorSymbolic_SeqSBAIJ() local
272 nzk = 0; in MatCholeskyFactorSymbolic_SeqSBAIJ()
[all …]
|
| H A D | sbaijfact2.c | 2285 PetscInt *prowl, *q, jmin, jmax, juidx, nzk, qm, *iu, *ju, k, j, vj, umax, maxadd; in MatICCFactorSymbolic_SeqSBAIJ_MSR() local
2320 nzk = 0; in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2334 nzk++; in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2360 nzk++; in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2371 if (nzk > 1) { in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2376 iu[k + 1] = iu[k] + nzk; in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2384 if (maxadd < nzk) maxadd = (mbs - k) * (nzk + 1) / 2; in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2402 while (nzk--) { in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2477 PetscInt jmin, jmax, nzk, k, j, *jl, prow, *il, nextprow; in MatICCFactorSymbolic_SeqSBAIJ() local
2545 nzk = 0; in MatICCFactorSymbolic_SeqSBAIJ()
[all …]
|
| /petsc/src/mat/impls/baij/seq/ |
| H A D | baijfact.c | 1024 PetscInt jmin, jmax, nzk, k, j, *jl, prow, *il, nextprow; in MatICCFactorSymbolic_SeqBAIJ() local
1105 nzk = 0; in MatICCFactorSymbolic_SeqBAIJ()
1118 nzk += nlnk; in MatICCFactorSymbolic_SeqBAIJ()
1134 nzk += nlnk; in MatICCFactorSymbolic_SeqBAIJ()
1148 if (current_space->local_remaining < nzk) { in MatICCFactorSymbolic_SeqBAIJ()
1150 …i = PetscMin(PetscIntMultTruncate(i, nzk), PetscIntMultTruncate(i, i - 1)); /* i*nzk, i*(i-1): est… in MatICCFactorSymbolic_SeqBAIJ()
1157 …PetscCall(PetscIncompleteLLClean(am, am, nzk, lnk, lnk_lvl, current_space->array, current_space_lv… in MatICCFactorSymbolic_SeqBAIJ()
1160 if (nzk - 1 > 0) { in MatICCFactorSymbolic_SeqBAIJ()
1169 current_space->array += nzk; in MatICCFactorSymbolic_SeqBAIJ()
1170 current_space->local_used += nzk; in MatICCFactorSymbolic_SeqBAIJ()
[all …]
|
| /petsc/src/mat/impls/aij/seq/ |
| H A D | aijfact.c | 2058 PetscInt jmin, jmax, nzk, k, j, *jl, prow, *il, nextprow; in MatICCFactorSymbolic_SeqAIJ() local
2121 nzk = 0; in MatICCFactorSymbolic_SeqAIJ()
2133 nzk += nlnk; in MatICCFactorSymbolic_SeqAIJ()
2150 nzk += nlnk; in MatICCFactorSymbolic_SeqAIJ()
2163 if (current_space->local_remaining < nzk) { in MatICCFactorSymbolic_SeqAIJ()
2165 …i = PetscIntMultTruncate(i, PetscMin(nzk, i - 1)); /* i*nzk, i*(i-1): estimated and max additional… in MatICCFactorSymbolic_SeqAIJ()
2172 …PetscCheck(nzk != 0, PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Empty row %" PetscInt_FMT " in ICC mat… in MatICCFactorSymbolic_SeqAIJ()
2173 …PetscCall(PetscIncompleteLLClean(am, am, nzk, lnk, lnk_lvl, current_space->array, current_space_lv… in MatICCFactorSymbolic_SeqAIJ()
2176 if (nzk > 1) { in MatICCFactorSymbolic_SeqAIJ()
2185 current_space->array += nzk; in MatICCFactorSymbolic_SeqAIJ()
[all …]
|