| /petsc/src/mat/impls/baij/seq/ |
| H A D | baijfact4.c | 13 PetscInt *ajtmpold, *ajtmp, nz, row, *ai = a->i, *aj = a->j, k, flg; in MatILUFactorNumeric_SeqBAIJ_N_inplace() local
34 ajtmp = bj + bi[i]; in MatILUFactorNumeric_SeqBAIJ_N_inplace()
35 for (j = 0; j < nz; j++) PetscCall(PetscArrayzero(rtmp + bs2 * ajtmp[j], bs2)); in MatILUFactorNumeric_SeqBAIJ_N_inplace()
41 row = *ajtmp++; in MatILUFactorNumeric_SeqBAIJ_N_inplace()
60 row = *ajtmp++; in MatILUFactorNumeric_SeqBAIJ_N_inplace()
|
| H A D | baijfact13.c | 16 PetscInt *ajtmpold, *ajtmp, nz, row, *ai = a->i, *aj = a->j; in MatILUFactorNumeric_SeqBAIJ_3_inplace() local
38 ajtmp = bj + bi[i]; in MatILUFactorNumeric_SeqBAIJ_3_inplace()
40 x = rtmp + 9 * ajtmp[j]; in MatILUFactorNumeric_SeqBAIJ_3_inplace()
61 row = *ajtmp++; in MatILUFactorNumeric_SeqBAIJ_3_inplace()
124 row = *ajtmp++; in MatILUFactorNumeric_SeqBAIJ_3_inplace()
175 const PetscInt *ajtmp, *bjtmp, *bdiag = b->diag, *pj, bs2 = a->bs2; in MatLUFactorNumeric_SeqBAIJ_3() local
204 ajtmp = aj + ai[r[i]]; in MatLUFactorNumeric_SeqBAIJ_3()
206 for (j = 0; j < nz; j++) PetscCall(PetscArraycpy(rtmp + bs2 * ic[ajtmp[j]], v + bs2 * j, bs2)); in MatLUFactorNumeric_SeqBAIJ_3()
276 PetscInt *ajtmpold, *ajtmp, nz, row; in MatILUFactorNumeric_SeqBAIJ_3_NaturalOrdering_inplace() local
296 ajtmp = bj + bi[i]; in MatILUFactorNumeric_SeqBAIJ_3_NaturalOrdering_inplace()
[all …]
|
| H A D | baijfact11.c | 16 PetscInt *ajtmpold, *ajtmp, nz, row; in MatILUFactorNumeric_SeqBAIJ_4_inplace() local
41 ajtmp = bj + bi[i]; in MatILUFactorNumeric_SeqBAIJ_4_inplace()
43 x = rtmp + 16 * ajtmp[j]; in MatILUFactorNumeric_SeqBAIJ_4_inplace()
72 row = *ajtmp++; in MatILUFactorNumeric_SeqBAIJ_4_inplace()
174 row = *ajtmp++; in MatILUFactorNumeric_SeqBAIJ_4_inplace()
237 const PetscInt *ajtmp, *bjtmp, *bdiag = b->diag, *pj, bs2 = a->bs2; in MatLUFactorNumeric_SeqBAIJ_4() local
272 ajtmp = aj + ai[r[i]]; in MatLUFactorNumeric_SeqBAIJ_4()
274 for (j = 0; j < nz; j++) PetscCall(PetscArraycpy(rtmp + bs2 * ic[ajtmp[j]], v + bs2 * j, bs2)); in MatLUFactorNumeric_SeqBAIJ_4()
347 PetscInt *ajtmpold, *ajtmp, nz, row; in MatILUFactorNumeric_SeqBAIJ_4_NaturalOrdering_inplace() local
370 ajtmp = bj + bi[i]; in MatILUFactorNumeric_SeqBAIJ_4_NaturalOrdering_inplace()
[all …]
|
| H A D | baijfact3.c | 158 PetscInt *bi, *bj, *ajtmp; in MatLUFactorSymbolic_SeqBAIJ() local
200 ajtmp = aj + ai[r[i]]; in MatLUFactorSymbolic_SeqBAIJ()
201 PetscCall(PetscLLAddPerm(nnz, ajtmp, ic, n, &nlnk, lnk, lnkbt)); in MatLUFactorSymbolic_SeqBAIJ()
208 ajtmp = bi_ptr[row] + nzbd; /* points to the entry next to the diagonal */ in MatLUFactorSymbolic_SeqBAIJ()
209 PetscCall(PetscLLAddSortedLU(ajtmp, row, &nlnk, lnk, lnkbt, i, nzbd, im)); in MatLUFactorSymbolic_SeqBAIJ()
|
| H A D | baijfact.c | 14 const PetscInt *ajtmp, *bjtmp, *bdiag = b->diag; in MatLUFactorNumeric_SeqBAIJ_2() local
43 ajtmp = aj + ai[r[i]]; in MatLUFactorNumeric_SeqBAIJ_2()
45 …for (PetscInt j = 0; j < nz; j++) PetscCall(PetscArraycpy(rtmp + bs2 * ic[ajtmp[j]], v + bs2 * j, … in MatLUFactorNumeric_SeqBAIJ_2()
123 const PetscInt *ajtmp, *bjtmp, *bdiag = b->diag; in MatLUFactorNumeric_SeqBAIJ_2_NaturalOrdering() local
151 ajtmp = aj + ai[i]; in MatLUFactorNumeric_SeqBAIJ_2_NaturalOrdering()
153 for (j = 0; j < nz; j++) PetscCall(PetscArraycpy(rtmp + bs2 * ajtmp[j], v + bs2 * j, bs2)); in MatLUFactorNumeric_SeqBAIJ_2_NaturalOrdering()
230 PetscInt *ajtmpold, *ajtmp, nz, row; in MatILUFactorNumeric_SeqBAIJ_2_inplace() local
251 ajtmp = bj + bi[i]; in MatILUFactorNumeric_SeqBAIJ_2_inplace()
253 x = rtmp + 4 * ajtmp[j]; in MatILUFactorNumeric_SeqBAIJ_2_inplace()
269 row = *ajtmp++; in MatILUFactorNumeric_SeqBAIJ_2_inplace()
[all …]
|
| H A D | baijfact9.c | 15 PetscInt *bi = b->i, *bj = b->j, *ajtmpold, *ajtmp; in MatILUFactorNumeric_SeqBAIJ_5_inplace() local
45 ajtmp = bj + bi[i]; in MatILUFactorNumeric_SeqBAIJ_5_inplace()
48 PetscCall(PetscArrayzero(rtmp + 25 * ajtmp[j], 25)); in MatILUFactorNumeric_SeqBAIJ_5_inplace()
50 x = rtmp + 25 * ajtmp[j]; in MatILUFactorNumeric_SeqBAIJ_5_inplace()
94 row = *ajtmp++; in MatILUFactorNumeric_SeqBAIJ_5_inplace()
243 row = *ajtmp++; in MatILUFactorNumeric_SeqBAIJ_5_inplace()
315 const PetscInt *ajtmp, *bjtmp, *bdiag = b->diag, *pj, bs2 = a->bs2; in MatLUFactorNumeric_SeqBAIJ_5() local
344 ajtmp = aj + ai[r[i]]; in MatLUFactorNumeric_SeqBAIJ_5()
346 for (j = 0; j < nz; j++) PetscCall(PetscArraycpy(rtmp + bs2 * ic[ajtmp[j]], v + bs2 * j, bs2)); in MatLUFactorNumeric_SeqBAIJ_5()
419 PetscInt *ajtmpold, *ajtmp, nz, row; in MatILUFactorNumeric_SeqBAIJ_5_NaturalOrdering_inplace() local
[all …]
|
| H A D | baijfact7.c | 14 …const PetscInt *ajtmpold, *ajtmp, *diag_offset = b->diag, *r, *ic, *bi = b->i, *bj = b->j, *ai = a… in MatILUFactorNumeric_SeqBAIJ_6_inplace() local
41 ajtmp = bj + bi[i]; in MatILUFactorNumeric_SeqBAIJ_6_inplace()
43 x = rtmp + 36 * ajtmp[j]; in MatILUFactorNumeric_SeqBAIJ_6_inplace()
95 row = *ajtmp++; in MatILUFactorNumeric_SeqBAIJ_6_inplace()
301 row = *ajtmp++; in MatILUFactorNumeric_SeqBAIJ_6_inplace()
373 const PetscInt *ajtmp, *bjtmp, *bdiag = b->diag, *pj, bs2 = a->bs2; in MatLUFactorNumeric_SeqBAIJ_6() local
402 ajtmp = aj + ai[r[i]]; in MatLUFactorNumeric_SeqBAIJ_6()
404 for (j = 0; j < nz; j++) PetscCall(PetscArraycpy(rtmp + bs2 * ic[ajtmp[j]], v + bs2 * j, bs2)); in MatLUFactorNumeric_SeqBAIJ_6()
474 PetscInt *ajtmpold, *ajtmp, nz, row; in MatILUFactorNumeric_SeqBAIJ_6_NaturalOrdering_inplace() local
500 ajtmp = bj + bi[i]; in MatILUFactorNumeric_SeqBAIJ_6_NaturalOrdering_inplace()
[all …]
|
| 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
50 ajtmp = aj + ai[i]; in MatLUFactorNumeric_SeqBAIJ_15_NaturalOrdering()
52 for (j = 0; j < nz; j++) PetscCall(PetscArraycpy(rtmp + bs2 * ajtmp[j], v + bs2 * j, bs2)); in MatLUFactorNumeric_SeqBAIJ_15_NaturalOrdering()
121 PetscInt *ajtmp, *bjtmp, nz, nzL, row, *bdiag = b->diag, *pj; in MatLUFactorNumeric_SeqBAIJ_N() local
152 ajtmp = aj + ai[r[i]]; in MatLUFactorNumeric_SeqBAIJ_N()
154 for (j = 0; j < nz; j++) PetscCall(PetscArraycpy(rtmp + bs2 * ic[ajtmp[j]], v + bs2 * j, bs2)); in MatLUFactorNumeric_SeqBAIJ_N()
|
| H A D | baijfact5.c | 13 const PetscInt *r, *ic, *bi = b->i, *bj = b->j, *ajtmp, *ai = a->i, *aj = a->j, *pj, *ajtmpold; in MatILUFactorNumeric_SeqBAIJ_7_inplace() local
44 ajtmp = bj + bi[i]; in MatILUFactorNumeric_SeqBAIJ_7_inplace()
46 x = rtmp + 49 * ajtmp[j]; in MatILUFactorNumeric_SeqBAIJ_7_inplace()
112 row = *ajtmp++; in MatILUFactorNumeric_SeqBAIJ_7_inplace()
384 row = *ajtmp++; in MatILUFactorNumeric_SeqBAIJ_7_inplace()
469 const PetscInt *ajtmp, *bjtmp, *bdiag = b->diag, *pj, bs2 = a->bs2; in MatLUFactorNumeric_SeqBAIJ_7() local
498 ajtmp = aj + ai[r[i]]; in MatLUFactorNumeric_SeqBAIJ_7()
500 for (j = 0; j < nz; j++) PetscCall(PetscArraycpy(rtmp + bs2 * ic[ajtmp[j]], v + bs2 * j, bs2)); in MatLUFactorNumeric_SeqBAIJ_7()
570 PetscInt *ajtmpold, *ajtmp, nz, row; in MatILUFactorNumeric_SeqBAIJ_7_NaturalOrdering_inplace() local
599 ajtmp = bj + bi[i]; in MatILUFactorNumeric_SeqBAIJ_7_NaturalOrdering_inplace()
[all …]
|
| /petsc/src/mat/impls/aij/seq/ |
| H A D | aijfact.c | 58 const PetscInt *r, *ic, *ai = a->i, *aj = a->j, *ajtmp; in MatLUFactorSymbolic_SeqAIJ() local
98 ajtmp = aj + ai[r[i]]; in MatLUFactorSymbolic_SeqAIJ()
99 PetscCall(PetscLLAddPerm(nnz, ajtmp, ic, n, &nlnk, lnk, lnkbt)); in MatLUFactorSymbolic_SeqAIJ()
106 ajtmp = bi_ptr[row] + nzbd; /* points to the entry next to the diagonal */ in MatLUFactorSymbolic_SeqAIJ()
107 PetscCall(PetscLLAddSortedLU(ajtmp, row, &nlnk, lnk, lnkbt, i, nzbd, im)); in MatLUFactorSymbolic_SeqAIJ()
226 const PetscInt *ajtmp, *bjtmp; in MatLUFactorNumeric_SeqAIJ() local
281 ajtmp = aj + ai[r[i]]; in MatLUFactorNumeric_SeqAIJ()
283 for (j = 0; j < nz; j++) rtmp[ics[ajtmp[j]]] = v[j]; in MatLUFactorNumeric_SeqAIJ()
405 const PetscInt *ajtmp, *bjtmp, *ddiag, *pj; in MatLUFactorNumeric_SeqAIJ_inplace() local
453 ajtmp = aj + ai[r[i]]; in MatLUFactorNumeric_SeqAIJ_inplace()
[all …]
|
| H A D | matmatmult.c | 1536 const PetscInt ajtmp = aj[j]; in MatMatMultNumericAdd_SeqAIJ_SeqDense() local
1537 r1 += aatmp * b1[ajtmp]; in MatMatMultNumericAdd_SeqAIJ_SeqDense()
1538 r2 += aatmp * b2[ajtmp]; in MatMatMultNumericAdd_SeqAIJ_SeqDense()
1539 r3 += aatmp * b3[ajtmp]; in MatMatMultNumericAdd_SeqAIJ_SeqDense()
1540 r4 += aatmp * b4[ajtmp]; in MatMatMultNumericAdd_SeqAIJ_SeqDense()
1587 const PetscInt ajtmp = aj[j]; in MatMatMultNumericAdd_SeqAIJ_SeqDense() local
1588 r1 += aatmp * b1[ajtmp]; in MatMatMultNumericAdd_SeqAIJ_SeqDense()
1589 r2 += aatmp * b2[ajtmp]; in MatMatMultNumericAdd_SeqAIJ_SeqDense()
1607 const PetscInt ajtmp = aj[j]; in MatMatMultNumericAdd_SeqAIJ_SeqDense() local
1608 r1 += aatmp * b1[ajtmp]; in MatMatMultNumericAdd_SeqAIJ_SeqDense()
[all …]
|
| H A D | inode.c | 1206 const PetscInt *ajtmp, *bjtmp; in MatLUFactorNumeric_SeqAIJ_Inode() local
1307 ajtmp = aj + ai[r[i]]; in MatLUFactorNumeric_SeqAIJ_Inode()
1309 for (j = 0; j < nz; j++) rtmp1[ics[ajtmp[j]]] = v[j]; in MatLUFactorNumeric_SeqAIJ_Inode()
1386 ajtmp = aj + ai[r[i]]; in MatLUFactorNumeric_SeqAIJ_Inode()
1390 col = ics[ajtmp[j]]; in MatLUFactorNumeric_SeqAIJ_Inode()
1520 ajtmp = aj + ai[r[i]]; in MatLUFactorNumeric_SeqAIJ_Inode()
1525 col = ics[ajtmp[j]]; in MatLUFactorNumeric_SeqAIJ_Inode()
1708 ajtmp = aj + ai[r[i]]; in MatLUFactorNumeric_SeqAIJ_Inode()
1714 col = ics[ajtmp[j]]; in MatLUFactorNumeric_SeqAIJ_Inode()
|
| /petsc/src/mat/impls/sbaij/seq/ |
| H A D | sbaijfact.c | 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()
1380 rtmp[ajtmp[j]] += PetscAbsScalar(v[j]); in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
|