Searched refs:uikdi (Results 1 – 3 of 3) sorted by relevance
| /petsc/src/mat/impls/sbaij/seq/ |
| H A D | sbaijfact.c | 1196 MatScalar *rtmp, *ba = b->a, *bval, *aa, dk, uikdi; in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace() local
1261 uikdi = -ba[ili] * ba[bi[i]]; /* diagonal(k) */ in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace()
1262 dk += uikdi * ba[ili]; in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace()
1263 ba[ili] = uikdi; /* -U(i,k) */ in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace()
1269 for (j = jmin; j < jmax; j++) rtmp[bj[j]] += uikdi * ba[j]; in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace()
1353 MatScalar *rtmp, *ba = b->a, *bval, *aa = a->a, dk, uikdi; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering() local
1427 uikdi = -ba[ili] * ba[bdiag[i]]; /* diagonal(k) */ in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1428 dk += uikdi * ba[ili]; /* update diag[k] */ in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1429 ba[ili] = uikdi; /* -U(i,k) */ in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1435 for (j = jmin; j < jmax; j++) rtmp[bj[j]] += uikdi * ba[j]; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
[all …]
|
| /petsc/src/mat/impls/baij/seq/ |
| H A D | baijfact.c | 756 MatScalar *rtmp, *ba = b->a, *bval, *aa = a->a, dk, uikdi; in MatCholeskyFactorNumeric_SeqBAIJ_N() local
813 uikdi = -ba[ili] * ba[bi[i]]; /* diagonal(k) */ in MatCholeskyFactorNumeric_SeqBAIJ_N()
814 dk += uikdi * ba[ili]; in MatCholeskyFactorNumeric_SeqBAIJ_N()
815 ba[ili] = uikdi; /* -U(i,k) */ in MatCholeskyFactorNumeric_SeqBAIJ_N()
821 for (j = jmin; j < jmax; j++) rtmp[bj[j]] += uikdi * ba[j]; in MatCholeskyFactorNumeric_SeqBAIJ_N()
893 MatScalar *rtmp, *ba = b->a, *aa = a->a, dk, uikdi, *aval, *bval; in MatCholeskyFactorNumeric_SeqBAIJ_N_NaturalOrdering() local
938 uikdi = -ba[ili] * ba[bi[i]]; in MatCholeskyFactorNumeric_SeqBAIJ_N_NaturalOrdering()
939 dk += uikdi * ba[ili]; in MatCholeskyFactorNumeric_SeqBAIJ_N_NaturalOrdering()
940 ba[ili] = uikdi; /* -U(i,k) */ in MatCholeskyFactorNumeric_SeqBAIJ_N_NaturalOrdering()
948 while (nz--) rtmp[*bcol++] += uikdi * (*bval++); in MatCholeskyFactorNumeric_SeqBAIJ_N_NaturalOrdering()
|
| /petsc/src/mat/impls/aij/seq/ |
| H A D | aijfact.c | 1713 MatScalar *rtmp, *ba = b->a, *bval, dk, uikdi; in MatCholeskyFactorNumeric_SeqAIJ() local
1788 uikdi = -ba[ili] * ba[bdiag[i]]; /* diagonal(k) */ in MatCholeskyFactorNumeric_SeqAIJ()
1789 dk += uikdi * ba[ili]; /* update diag[k] */ in MatCholeskyFactorNumeric_SeqAIJ()
1790 ba[ili] = uikdi; /* -U(i,k) */ in MatCholeskyFactorNumeric_SeqAIJ()
1796 for (j = jmin; j < jmax; j++) rtmp[bj[j]] += uikdi * ba[j]; in MatCholeskyFactorNumeric_SeqAIJ()
1880 MatScalar *rtmp, *ba = b->a, *bval, dk, uikdi; in MatCholeskyFactorNumeric_SeqAIJ_inplace() local
1952 uikdi = -ba[ili] * ba[bi[i]]; /* diagonal(k) */ in MatCholeskyFactorNumeric_SeqAIJ_inplace()
1953 dk += uikdi * ba[ili]; in MatCholeskyFactorNumeric_SeqAIJ_inplace()
1954 ba[ili] = uikdi; /* -U(i,k) */ in MatCholeskyFactorNumeric_SeqAIJ_inplace()
1960 for (j = jmin; j < jmax; j++) rtmp[bj[j]] += uikdi * ba[j]; in MatCholeskyFactorNumeric_SeqAIJ_inplace()
|