| /petsc/src/mat/impls/sbaij/seq/ |
| H A D | sbaijfact.c | 564 PetscInt *a2anew, k, k1, jmin, jmax, *jl, *il, vj, nexti, ili; in MatCholeskyFactorNumeric_SeqSBAIJ_N() local
647 ili = il[i]; /* index of first nonzero element in U(i,k:bms-1) */ in MatCholeskyFactorNumeric_SeqSBAIJ_N()
651 u = ba + ili * bs2; in MatCholeskyFactorNumeric_SeqSBAIJ_N()
660 PetscCall(PetscArraycpy(ba + ili * bs2, uik, bs2)); in MatCholeskyFactorNumeric_SeqSBAIJ_N()
663 jmin = ili + 1; in MatCholeskyFactorNumeric_SeqSBAIJ_N()
737 PetscInt *ai, *aj, k, k1, jmin, jmax, *jl, *il, vj, nexti, ili; in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering() local
779 ili = il[i]; /* index of first nonzero element in U(i,k:bms-1) */ in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering()
783 u = ba + ili * bs2; in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering()
792 PetscCall(PetscArraycpy(ba + ili * bs2, uik, bs2)); in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering()
795 jmin = ili + 1; in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering()
[all …]
|
| H A D | sbaijfact4.c | 11 PetscInt *ai, *aj, k, k1, jmin, jmax, *jl, *il, vj, nexti, ili; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering() local
52 ili = il[i]; /* index of first nonzero element in U(i,k:bms-1) */ in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
56 u = ba + ili * 9; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
86 PetscCall(PetscArraycpy(ba + ili * 9, uik, 9)); in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
89 jmin = ili + 1; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
|
| H A D | sbaijfact3.c | 10 PetscInt *a2anew, i, j, k, k1, jmin, jmax, *jl, *il, vj, nexti, ili; in MatCholeskyFactorNumeric_SeqSBAIJ_3() local
90 ili = il[i]; /* index of first nonzero element in U(i,k:bms-1) */ in MatCholeskyFactorNumeric_SeqSBAIJ_3()
94 u = ba + ili * 9; in MatCholeskyFactorNumeric_SeqSBAIJ_3()
124 PetscCall(PetscArraycpy(ba + ili * 9, uik, 9)); in MatCholeskyFactorNumeric_SeqSBAIJ_3()
127 jmin = ili + 1; in MatCholeskyFactorNumeric_SeqSBAIJ_3()
|
| H A D | sbaijfact5.c | 11 PetscInt *ai, *aj, k, k1, jmin, jmax, *jl, *il, vj, nexti, ili; in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering() local
53 ili = il[i]; /* index of first nonzero element in U(i,k:bms-1) */ in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
57 u = ba + ili * 16; in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
103 PetscCall(PetscArraycpy(ba + ili * 16, uik, 16)); in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
106 jmin = ili + 1; in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
|
| H A D | sbaijfact6.c | 10 PetscInt i, j, *a2anew, k, k1, jmin, jmax, *jl, *il, vj, nexti, ili; in MatCholeskyFactorNumeric_SeqSBAIJ_4() local
91 ili = il[i]; /* index of first nonzero element in U(i,k:bms-1) */ in MatCholeskyFactorNumeric_SeqSBAIJ_4()
95 u = ba + ili * 16; in MatCholeskyFactorNumeric_SeqSBAIJ_4()
141 PetscCall(PetscArraycpy(ba + ili * 16, uik, 16)); in MatCholeskyFactorNumeric_SeqSBAIJ_4()
144 jmin = ili + 1; in MatCholeskyFactorNumeric_SeqSBAIJ_4()
|
| H A D | sbaijfact8.c | 11 PetscInt *ai, *aj, k, k1, jmin, jmax, *jl, *il, vj, nexti, ili, ipvt[5]; in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering() local
52 ili = il[i]; /* index of first nonzero element in U(i,k:bms-1) */ in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
56 u = ba + ili * 25; in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
122 PetscCall(PetscArraycpy(ba + ili * 25, uik, 25)); in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
125 jmin = ili + 1; in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
|
| H A D | sbaijfact7.c | 10 PetscInt i, j, *a2anew, k, k1, jmin, jmax, *jl, *il, vj, nexti, ili, ipvt[5]; in MatCholeskyFactorNumeric_SeqSBAIJ_5() local
90 ili = il[i]; /* index of first nonzero element in U(i,k:bms-1) */ in MatCholeskyFactorNumeric_SeqSBAIJ_5()
94 u = ba + ili * 25; in MatCholeskyFactorNumeric_SeqSBAIJ_5()
160 PetscCall(PetscArraycpy(ba + ili * 25, uik, 25)); in MatCholeskyFactorNumeric_SeqSBAIJ_5()
163 jmin = ili + 1; in MatCholeskyFactorNumeric_SeqSBAIJ_5()
|
| H A D | sbaijfact10.c | 11 PetscInt *ai, *aj, k, k1, jmin, jmax, *jl, *il, vj, nexti, ili; in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering() local
60 ili = il[i]; /* index of first nonzero element in U(i,k:bms-1) */ in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
64 u = ba + ili * 36; in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
228 PetscCall(PetscArraycpy(ba + ili * 36, uik, 36)); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
231 jmin = ili + 1; in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
|
| H A D | sbaijfact9.c | 10 PetscInt i, j, *a2anew, k, k1, jmin, jmax, *jl, *il, vj, nexti, ili; in MatCholeskyFactorNumeric_SeqSBAIJ_6() local
92 ili = il[i]; /* index of first nonzero element in U(i,k:bms-1) */ in MatCholeskyFactorNumeric_SeqSBAIJ_6()
96 u = ba + ili * 36; in MatCholeskyFactorNumeric_SeqSBAIJ_6()
223 PetscCall(PetscArraycpy(ba + ili * 36, uik, 36)); in MatCholeskyFactorNumeric_SeqSBAIJ_6()
226 jmin = ili + 1; in MatCholeskyFactorNumeric_SeqSBAIJ_6()
|
| H A D | sbaijfact12.c | 11 PetscInt *ai, *aj, k, k1, jmin, jmax, *jl, *il, vj, nexti, ili; in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering() local
56 ili = il[i]; /* index of first nonzero element in U(i,k:bms-1) */ in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
60 u = ba + ili * 49; in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
228 PetscCall(PetscArraycpy(ba + ili * 49, uik, 49)); in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
231 jmin = ili + 1; in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
|
| H A D | sbaijfact11.c | 10 PetscInt *a2anew, k, k1, jmin, jmax, *jl, *il, vj, nexti, ili, i, j, mbs = a->mbs; in MatCholeskyFactorNumeric_SeqSBAIJ_7() local
93 ili = il[i]; /* index of first nonzero element in U(i,k:bms-1) */ in MatCholeskyFactorNumeric_SeqSBAIJ_7()
97 u = ba + ili * 49; in MatCholeskyFactorNumeric_SeqSBAIJ_7()
265 PetscCall(PetscArraycpy(ba + ili * 49, uik, 49)); in MatCholeskyFactorNumeric_SeqSBAIJ_7()
268 jmin = ili + 1; in MatCholeskyFactorNumeric_SeqSBAIJ_7()
|
| /petsc/src/mat/impls/baij/seq/ |
| H A D | baijfact.c | 755 PetscInt k, jmin, jmax, *jl, *il, col, nexti, ili, nz; in MatCholeskyFactorNumeric_SeqBAIJ_N() local
812 ili = il[i]; /* index of first nonzero element in U(i,k:bms-1) */ in MatCholeskyFactorNumeric_SeqBAIJ_N()
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()
818 jmin = ili + 1; in MatCholeskyFactorNumeric_SeqBAIJ_N()
892 PetscInt k, jmin, *jl, *il, nexti, ili, *acol, *bcol, nz; in MatCholeskyFactorNumeric_SeqBAIJ_N_NaturalOrdering() local
937 ili = il[i]; /* index of first nonzero element in U(i,k:bms-1) */ in MatCholeskyFactorNumeric_SeqBAIJ_N_NaturalOrdering()
938 uikdi = -ba[ili] * ba[bi[i]]; in MatCholeskyFactorNumeric_SeqBAIJ_N_NaturalOrdering()
939 dk += uikdi * ba[ili]; in MatCholeskyFactorNumeric_SeqBAIJ_N_NaturalOrdering()
[all …]
|
| /petsc/src/mat/impls/aij/seq/ |
| H A D | aijfact.c | 1712 PetscInt k, jmin, jmax, *c2r, *il, col, nexti, ili, nz; in MatCholeskyFactorNumeric_SeqAIJ() local
1787 ili = il[i]; /* index of first nonzero element in U(i,k:bms-1) */ in MatCholeskyFactorNumeric_SeqAIJ()
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()
1793 jmin = ili + 1; in MatCholeskyFactorNumeric_SeqAIJ()
1879 PetscInt k, jmin, jmax, *jl, *il, col, nexti, ili, nz; in MatCholeskyFactorNumeric_SeqAIJ_inplace() local
1951 ili = il[i]; /* index of first nonzero element in U(i,k:bms-1) */ in MatCholeskyFactorNumeric_SeqAIJ_inplace()
1952 uikdi = -ba[ili] * ba[bi[i]]; /* diagonal(k) */ in MatCholeskyFactorNumeric_SeqAIJ_inplace()
1953 dk += uikdi * ba[ili]; in MatCholeskyFactorNumeric_SeqAIJ_inplace()
[all …]
|