| /petsc/src/mat/impls/baij/seq/ |
| H A D | baijfact2.c | 20 …const PetscInt n = a->mbs, *ai = a->i, *aj = a->j, *bi = b->i, *bj = b->j, *ajtmp, *bjtmp, *bdia… in MatLUFactorNumeric_SeqBAIJ_15_NaturalOrdering() local
40 bjtmp = bj + bi[i]; in MatLUFactorNumeric_SeqBAIJ_15_NaturalOrdering()
41 for (j = 0; j < nz; j++) PetscCall(PetscArrayzero(rtmp + bs2 * bjtmp[j], bs2)); in MatLUFactorNumeric_SeqBAIJ_15_NaturalOrdering()
45 bjtmp = bj + bdiag[i + 1] + 1; in MatLUFactorNumeric_SeqBAIJ_15_NaturalOrdering()
46 for (j = 0; j < nz; j++) PetscCall(PetscArrayzero(rtmp + bs2 * bjtmp[j], bs2)); in MatLUFactorNumeric_SeqBAIJ_15_NaturalOrdering()
55 bjtmp = bj + bi[i]; in MatLUFactorNumeric_SeqBAIJ_15_NaturalOrdering()
58 row = bjtmp[k]; in MatLUFactorNumeric_SeqBAIJ_15_NaturalOrdering()
121 PetscInt *ajtmp, *bjtmp, nz, nzL, row, *bdiag = b->diag, *pj; in MatLUFactorNumeric_SeqBAIJ_N() local
142 bjtmp = bj + bi[i]; in MatLUFactorNumeric_SeqBAIJ_N()
143 for (j = 0; j < nz; j++) PetscCall(PetscArrayzero(rtmp + bs2 * bjtmp[j], bs2)); in MatLUFactorNumeric_SeqBAIJ_N()
[all …]
|
| H A D | baijfact13.c | 175 const PetscInt *ajtmp, *bjtmp, *bdiag = b->diag, *pj, bs2 = a->bs2; in MatLUFactorNumeric_SeqBAIJ_3() local
194 bjtmp = bj + bi[i]; in MatLUFactorNumeric_SeqBAIJ_3()
195 for (j = 0; j < nz; j++) PetscCall(PetscArrayzero(rtmp + bs2 * bjtmp[j], bs2)); in MatLUFactorNumeric_SeqBAIJ_3()
199 bjtmp = bj + bdiag[i + 1] + 1; in MatLUFactorNumeric_SeqBAIJ_3()
200 for (j = 0; j < nz; j++) PetscCall(PetscArrayzero(rtmp + bs2 * bjtmp[j], bs2)); in MatLUFactorNumeric_SeqBAIJ_3()
209 bjtmp = bj + bi[i]; in MatLUFactorNumeric_SeqBAIJ_3()
212 row = bjtmp[k]; in MatLUFactorNumeric_SeqBAIJ_3()
427 const PetscInt *ajtmp, *bjtmp, *bdiag = b->diag, *pj, bs2 = a->bs2; in MatLUFactorNumeric_SeqBAIJ_3_NaturalOrdering() local
444 bjtmp = bj + bi[i]; in MatLUFactorNumeric_SeqBAIJ_3_NaturalOrdering()
445 for (j = 0; j < nz; j++) PetscCall(PetscArrayzero(rtmp + bs2 * bjtmp[j], bs2)); in MatLUFactorNumeric_SeqBAIJ_3_NaturalOrdering()
[all …]
|
| H A D | baijfact11.c | 237 const PetscInt *ajtmp, *bjtmp, *bdiag = b->diag, *pj, bs2 = a->bs2; in MatLUFactorNumeric_SeqBAIJ_4() local
262 bjtmp = bj + bi[i]; in MatLUFactorNumeric_SeqBAIJ_4()
263 for (j = 0; j < nz; j++) PetscCall(PetscArrayzero(rtmp + bs2 * bjtmp[j], bs2)); in MatLUFactorNumeric_SeqBAIJ_4()
267 bjtmp = bj + bdiag[i + 1] + 1; in MatLUFactorNumeric_SeqBAIJ_4()
268 for (j = 0; j < nz; j++) PetscCall(PetscArrayzero(rtmp + bs2 * bjtmp[j], bs2)); in MatLUFactorNumeric_SeqBAIJ_4()
277 bjtmp = bj + bi[i]; in MatLUFactorNumeric_SeqBAIJ_4()
280 row = bjtmp[k]; in MatLUFactorNumeric_SeqBAIJ_4()
557 const PetscInt *ajtmp, *bjtmp, *bdiag = b->diag, *pj, bs2 = a->bs2; in MatLUFactorNumeric_SeqBAIJ_4_NaturalOrdering() local
580 bjtmp = bj + bi[i]; in MatLUFactorNumeric_SeqBAIJ_4_NaturalOrdering()
581 for (j = 0; j < nz; j++) PetscCall(PetscArrayzero(rtmp + bs2 * bjtmp[j], bs2)); in MatLUFactorNumeric_SeqBAIJ_4_NaturalOrdering()
[all …]
|
| H A D | baijfact.c | 14 const PetscInt *ajtmp, *bjtmp, *bdiag = b->diag; in MatLUFactorNumeric_SeqBAIJ_2() local
33 bjtmp = bj + bi[i]; in MatLUFactorNumeric_SeqBAIJ_2()
34 for (PetscInt j = 0; j < nz; j++) PetscCall(PetscArrayzero(rtmp + bs2 * bjtmp[j], bs2)); in MatLUFactorNumeric_SeqBAIJ_2()
38 bjtmp = bj + bdiag[i + 1] + 1; in MatLUFactorNumeric_SeqBAIJ_2()
39 for (PetscInt j = 0; j < nz; j++) PetscCall(PetscArrayzero(rtmp + bs2 * bjtmp[j], bs2)); in MatLUFactorNumeric_SeqBAIJ_2()
48 bjtmp = bj + bi[i]; in MatLUFactorNumeric_SeqBAIJ_2()
52 PetscInt row = bjtmp[k]; in MatLUFactorNumeric_SeqBAIJ_2()
123 const PetscInt *ajtmp, *bjtmp, *bdiag = b->diag; in MatLUFactorNumeric_SeqBAIJ_2_NaturalOrdering() local
141 bjtmp = bj + bi[i]; in MatLUFactorNumeric_SeqBAIJ_2_NaturalOrdering()
142 for (j = 0; j < nz; j++) PetscCall(PetscArrayzero(rtmp + bs2 * bjtmp[j], bs2)); in MatLUFactorNumeric_SeqBAIJ_2_NaturalOrdering()
[all …]
|
| H A D | baijfact9.c | 315 const PetscInt *ajtmp, *bjtmp, *bdiag = b->diag, *pj, bs2 = a->bs2; in MatLUFactorNumeric_SeqBAIJ_5() local
334 bjtmp = bj + bi[i]; in MatLUFactorNumeric_SeqBAIJ_5()
335 for (j = 0; j < nz; j++) PetscCall(PetscArrayzero(rtmp + bs2 * bjtmp[j], bs2)); in MatLUFactorNumeric_SeqBAIJ_5()
339 bjtmp = bj + bdiag[i + 1] + 1; in MatLUFactorNumeric_SeqBAIJ_5()
340 for (j = 0; j < nz; j++) PetscCall(PetscArrayzero(rtmp + bs2 * bjtmp[j], bs2)); in MatLUFactorNumeric_SeqBAIJ_5()
349 bjtmp = bj + bi[i]; in MatLUFactorNumeric_SeqBAIJ_5()
352 row = bjtmp[k]; in MatLUFactorNumeric_SeqBAIJ_5()
682 const PetscInt *ajtmp, *bjtmp, *bdiag = b->diag, *pj, bs2 = a->bs2; in MatLUFactorNumeric_SeqBAIJ_5_NaturalOrdering() local
699 bjtmp = bj + bi[i]; in MatLUFactorNumeric_SeqBAIJ_5_NaturalOrdering()
700 for (j = 0; j < nz; j++) PetscCall(PetscArrayzero(rtmp + bs2 * bjtmp[j], bs2)); in MatLUFactorNumeric_SeqBAIJ_5_NaturalOrdering()
[all …]
|
| H A D | baijfact7.c | 373 const PetscInt *ajtmp, *bjtmp, *bdiag = b->diag, *pj, bs2 = a->bs2; in MatLUFactorNumeric_SeqBAIJ_6() local
392 bjtmp = bj + bi[i]; in MatLUFactorNumeric_SeqBAIJ_6()
393 for (j = 0; j < nz; j++) PetscCall(PetscArrayzero(rtmp + bs2 * bjtmp[j], bs2)); in MatLUFactorNumeric_SeqBAIJ_6()
397 bjtmp = bj + bdiag[i + 1] + 1; in MatLUFactorNumeric_SeqBAIJ_6()
398 for (j = 0; j < nz; j++) PetscCall(PetscArrayzero(rtmp + bs2 * bjtmp[j], bs2)); in MatLUFactorNumeric_SeqBAIJ_6()
407 bjtmp = bj + bi[i]; in MatLUFactorNumeric_SeqBAIJ_6()
410 row = bjtmp[k]; in MatLUFactorNumeric_SeqBAIJ_6()
827 const PetscInt *ajtmp, *bjtmp, *bdiag = b->diag, *pj, bs2 = a->bs2; in MatLUFactorNumeric_SeqBAIJ_6_NaturalOrdering() local
844 bjtmp = bj + bi[i]; in MatLUFactorNumeric_SeqBAIJ_6_NaturalOrdering()
845 for (j = 0; j < nz; j++) PetscCall(PetscArrayzero(rtmp + bs2 * bjtmp[j], bs2)); in MatLUFactorNumeric_SeqBAIJ_6_NaturalOrdering()
[all …]
|
| H A D | baijfact5.c | 469 const PetscInt *ajtmp, *bjtmp, *bdiag = b->diag, *pj, bs2 = a->bs2; in MatLUFactorNumeric_SeqBAIJ_7() local
488 bjtmp = bj + bi[i]; in MatLUFactorNumeric_SeqBAIJ_7()
489 for (j = 0; j < nz; j++) PetscCall(PetscArrayzero(rtmp + bs2 * bjtmp[j], bs2)); in MatLUFactorNumeric_SeqBAIJ_7()
493 bjtmp = bj + bdiag[i + 1] + 1; in MatLUFactorNumeric_SeqBAIJ_7()
494 for (j = 0; j < nz; j++) PetscCall(PetscArrayzero(rtmp + bs2 * bjtmp[j], bs2)); in MatLUFactorNumeric_SeqBAIJ_7()
503 bjtmp = bj + bi[i]; in MatLUFactorNumeric_SeqBAIJ_7()
506 row = bjtmp[k]; in MatLUFactorNumeric_SeqBAIJ_7()
1019 const PetscInt *ajtmp, *bjtmp, *bdiag = b->diag, *pj, bs2 = a->bs2; in MatLUFactorNumeric_SeqBAIJ_7_NaturalOrdering() local
1036 bjtmp = bj + bi[i]; in MatLUFactorNumeric_SeqBAIJ_7_NaturalOrdering()
1037 for (j = 0; j < nz; j++) PetscCall(PetscArrayzero(rtmp + bs2 * bjtmp[j], bs2)); in MatLUFactorNumeric_SeqBAIJ_7_NaturalOrdering()
[all …]
|
| /petsc/src/mat/impls/aij/seq/ |
| H A D | aijfact.c | 226 const PetscInt *ajtmp, *bjtmp; in MatLUFactorNumeric_SeqAIJ() local
271 bjtmp = bj + bi[i]; in MatLUFactorNumeric_SeqAIJ()
272 for (j = 0; j < nz; j++) rtmp[bjtmp[j]] = 0.0; in MatLUFactorNumeric_SeqAIJ()
276 bjtmp = bj + bdiag[i + 1] + 1; in MatLUFactorNumeric_SeqAIJ()
277 for (j = 0; j < nz; j++) rtmp[bjtmp[j]] = 0.0; in MatLUFactorNumeric_SeqAIJ()
288 bjtmp = bj + bi[i]; in MatLUFactorNumeric_SeqAIJ()
289 row = *bjtmp++; in MatLUFactorNumeric_SeqAIJ()
305 row = *bjtmp++; in MatLUFactorNumeric_SeqAIJ()
405 const PetscInt *ajtmp, *bjtmp, *ddiag, *pj; in MatLUFactorNumeric_SeqAIJ_inplace() local
448 bjtmp = bj + bi[i]; in MatLUFactorNumeric_SeqAIJ_inplace()
[all …]
|
| H A D | inode.c | 1206 const PetscInt *ajtmp, *bjtmp; in MatLUFactorNumeric_SeqAIJ_Inode() local
1297 bjtmp = bj + bi[i]; in MatLUFactorNumeric_SeqAIJ_Inode()
1298 for (j = 0; j < nz; j++) rtmp1[bjtmp[j]] = 0.0; in MatLUFactorNumeric_SeqAIJ_Inode()
1302 bjtmp = bj + bdiag[i + 1] + 1; in MatLUFactorNumeric_SeqAIJ_Inode()
1303 for (j = 0; j < nz; j++) rtmp1[bjtmp[j]] = 0.0; in MatLUFactorNumeric_SeqAIJ_Inode()
1315 bjtmp = bj + bi[i]; in MatLUFactorNumeric_SeqAIJ_Inode()
1316 row = *bjtmp++; in MatLUFactorNumeric_SeqAIJ_Inode()
1330 row = *bjtmp++; in MatLUFactorNumeric_SeqAIJ_Inode()
1368 bjtmp = bj + bi[i]; in MatLUFactorNumeric_SeqAIJ_Inode()
1370 col = bjtmp[j]; in MatLUFactorNumeric_SeqAIJ_Inode()
[all …]
|
| /petsc/src/mat/impls/sbaij/seq/ |
| H A D | sbaijfact.c | 1350 PetscInt i, j, mbs = A->rmap->n, *bi = b->i, *bj = b->j, *bdiag = b->diag, *bjtmp; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering() local
1403 bjtmp = bj + bi[k]; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1404 for (j = 0; j < nz; j++) rtmp[bjtmp[j]] = 0.0; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
|