| /petsc/src/mat/impls/sbaij/mpi/ |
| H A D | mpisbaij.c | 208 bcol = col / bs; \
215 if (rp[t] > bcol) high = t; \
219 if (rp[_i] > bcol) break; \
220 if (rp[_i] == bcol) { \
229 …MatSeqXAIJReallocateAIJ(A, a->mbs, bs2, nrow, brow, bcol, rmax, aa, ai, aj, rp, ap, aimax, a->none…
235 rp[_i] = bcol; \
248 bcol = col / bs; \
255 if (rp[t] > bcol) high = t; \
259 if (rp[_i] > bcol) break; \
260 if (rp[_i] == bcol) { \
[all …]
|
| /petsc/src/mat/impls/sbaij/seq/ |
| H A D | sbaijfact.c | 1194 PetscInt *a2anew, i, j, mbs = a->mbs, *bi = b->i, *bj = b->j, *bcol; in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace() local
1287 bcol = bj + jmin; in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace()
1289 rs += PetscAbsScalar(rtmp[*bcol]); in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace()
1290 bcol++; in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace()
1504 PetscInt k, jmin, *jl, *il, nexti, ili, *acol, *bcol, nz; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering_inplace() local
1563 bcol = bj + jmin; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering_inplace()
1566 while (nz--) rtmp[*bcol++] += uikdi * (*bval++); in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering_inplace()
1583 bcol = bj + jmin; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering_inplace()
1585 rs += PetscAbsScalar(rtmp[*bcol]); in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering_inplace()
1586 bcol++; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering_inplace()
[all …]
|
| H A D | sbaij.c | 515 PetscInt brow, bcol, ridx, cidx, bs = A->rmap->bs, bs2 = a->bs2; in MatGetValues_SeqSBAIJ() local
537 bcol = col / bs; in MatGetValues_SeqSBAIJ()
544 if (rp[t] > bcol) high = t; in MatGetValues_SeqSBAIJ()
548 if (rp[i] > bcol) break; in MatGetValues_SeqSBAIJ()
549 if (rp[i] == bcol) { in MatGetValues_SeqSBAIJ()
742 PetscInt *aj = a->j, nonew = a->nonew, bs = A->rmap->bs, brow, bcol; in MatSetValues_SeqSBAIJ() local
762 bcol = col / bs; /* block col number */ in MatSetValues_SeqSBAIJ()
764 if (brow > bcol) { in MatSetValues_SeqSBAIJ()
771 if ((brow == bcol && ridx <= cidx) || (brow < bcol)) { in MatSetValues_SeqSBAIJ()
783 if (rp[t] > bcol) high = t; in MatSetValues_SeqSBAIJ()
[all …]
|
| H A D | sbaij2.c | 11 PetscInt brow, i, j, k, l, mbs, n, *nidx, isz, bcol, bcol_max, start, end, *ai, *aj, bs; in MatIncreaseOverlap_SeqSBAIJ() local
58 bcol = aj[l]; in MatIncreaseOverlap_SeqSBAIJ()
59 if (!PetscBTLookupSet(table_out, bcol)) { in MatIncreaseOverlap_SeqSBAIJ()
60 nidx[isz++] = bcol; in MatIncreaseOverlap_SeqSBAIJ()
61 if (bcol_max < bcol) bcol_max = bcol; in MatIncreaseOverlap_SeqSBAIJ()
68 bcol = aj[l]; in MatIncreaseOverlap_SeqSBAIJ()
69 if (bcol > bcol_max) break; in MatIncreaseOverlap_SeqSBAIJ()
70 if (PetscBTLookup(table_in, bcol)) { in MatIncreaseOverlap_SeqSBAIJ()
1521 PetscInt ncols, brow, bcol, krow, kcol; in MatGetRowMaxAbs_SeqSBAIJ() local
1541 bcol = bs * (*aj); in MatGetRowMaxAbs_SeqSBAIJ()
[all …]
|
| /petsc/src/mat/impls/baij/seq/ |
| H A D | baijfact.c | 753 PetscInt i, j, mbs = a->mbs, bs = A->rmap->bs, *bi = b->i, *bj = b->j, *bcol; in MatCholeskyFactorNumeric_SeqBAIJ_N() local
837 bcol = bj + jmin; in MatCholeskyFactorNumeric_SeqBAIJ_N()
839 rs += PetscAbsScalar(rtmp[*bcol]); in MatCholeskyFactorNumeric_SeqBAIJ_N()
840 bcol++; in MatCholeskyFactorNumeric_SeqBAIJ_N()
892 PetscInt k, jmin, *jl, *il, nexti, ili, *acol, *bcol, nz; in MatCholeskyFactorNumeric_SeqBAIJ_N_NaturalOrdering() local
946 bcol = bj + jmin; in MatCholeskyFactorNumeric_SeqBAIJ_N_NaturalOrdering()
948 while (nz--) rtmp[*bcol++] += uikdi * (*bval++); in MatCholeskyFactorNumeric_SeqBAIJ_N_NaturalOrdering()
964 bcol = bj + jmin; in MatCholeskyFactorNumeric_SeqBAIJ_N_NaturalOrdering()
966 rs += PetscAbsScalar(rtmp[*bcol]); in MatCholeskyFactorNumeric_SeqBAIJ_N_NaturalOrdering()
967 bcol++; in MatCholeskyFactorNumeric_SeqBAIJ_N_NaturalOrdering()
[all …]
|
| H A D | baij.c | 1362 PetscInt *aj = a->j, brow, bcol; in matsetvalues4_() local
1377 bcol = col / 4; in matsetvalues4_()
1386 if (rp[t] > bcol) high = t; in matsetvalues4_()
1390 if (rp[i] > bcol) break; in matsetvalues4_()
1391 if (rp[i] == bcol) { in matsetvalues4_()
1403 rp[i] = bcol; in matsetvalues4_()
2059 PetscInt brow, bcol, ridx, cidx, bs = A->rmap->bs, bs2 = a->bs2; in MatGetValues_SeqBAIJ() local
2081 bcol = col / bs; in MatGetValues_SeqBAIJ()
2088 if (rp[t] > bcol) high = t; in MatGetValues_SeqBAIJ()
2092 if (rp[i] > bcol) break; in MatGetValues_SeqBAIJ()
[all …]
|
| H A D | baij2.c | 3180 PetscInt *bcol = a->j; in MatNorm_SeqBAIJ() local
3184 k1 = bs * (*bcol) + j; /* column index */ in MatNorm_SeqBAIJ()
3190 bcol++; in MatNorm_SeqBAIJ()
|
| /petsc/src/mat/impls/blockmat/seq/ |
| H A D | blockmat.c | 219 PetscInt *aj = a->j, nonew = a->nonew, bs = A->rmap->bs, brow, bcol; in MatSetValues_BlockMat() local
240 bcol = col / bs; in MatSetValues_BlockMat()
241 if (A->symmetric == PETSC_BOOL3_TRUE && brow > bcol) continue; in MatSetValues_BlockMat()
252 if (rp[t] > bcol) high = t; in MatSetValues_BlockMat()
256 if (rp[i] > bcol) break; in MatSetValues_BlockMat()
257 if (rp[i] == bcol) goto noinsert1; in MatSetValues_BlockMat()
262 …MatSeqXAIJReallocateAIJ(A, a->mbs, 1, nrow, brow, bcol, rmax, aa, ai, aj, rp, ap, imax, nonew, Mat… in MatSetValues_BlockMat()
271 rp[i] = bcol; in MatSetValues_BlockMat()
|
| /petsc/src/mat/impls/aij/seq/ |
| H A D | matmatmult.c | 725 PetscInt bcol = bj[bb[j]]; in MatMatMultSymbolic_SeqAIJ_SeqAIJ_BTHeap() local
726 if (!PetscBTLookupSet(bt, bcol)) { /* new entry */ in MatMatMultSymbolic_SeqAIJ_SeqAIJ_BTHeap()
727 PetscCall(PetscHeapAdd(h, j, bcol)); in MatMatMultSymbolic_SeqAIJ_SeqAIJ_BTHeap()
748 PetscInt bcol = bj[bb[j]]; in MatMatMultSymbolic_SeqAIJ_SeqAIJ_BTHeap() local
749 if (!PetscBTLookupSet(bt, bcol)) { /* new entry */ in MatMatMultSymbolic_SeqAIJ_SeqAIJ_BTHeap()
750 PetscCall(PetscHeapAdd(h, j, bcol)); in MatMatMultSymbolic_SeqAIJ_SeqAIJ_BTHeap()
1111 PetscInt *ci, *cj, bcol; in MatMatMultSymbolic_SeqAIJ_SeqAIJ_Sorted() local
1137 bcol = bj[k]; in MatMatMultSymbolic_SeqAIJ_SeqAIJ_Sorted()
1138 if (!seen[bcol]) { /* new entry */ in MatMatMultSymbolic_SeqAIJ_SeqAIJ_Sorted()
1141 *slot = bcol; in MatMatMultSymbolic_SeqAIJ_SeqAIJ_Sorted()
[all …]
|
| H A D | aijfact.c | 1877 PetscInt i, j, mbs = A->rmap->n, *bi = b->i, *bj = b->j, *bcol, *bjtmp; in MatCholeskyFactorNumeric_SeqAIJ_inplace() local
1975 bcol = bj + jmin; in MatCholeskyFactorNumeric_SeqAIJ_inplace()
1976 for (j = 0; j < nz; j++) rs += PetscAbsScalar(rtmp[bcol[j]]); in MatCholeskyFactorNumeric_SeqAIJ_inplace()
|
| /petsc/src/mat/impls/baij/mpi/ |
| H A D | mpibaij.c | 167 bcol = col / bs; \
174 if (rp[t] > bcol) high = t; \
178 if (rp[_i] > bcol) break; \
179 if (rp[_i] == bcol) { \
188 …MatSeqXAIJReallocateAIJ(A, a->mbs, bs2, nrow, brow, bcol, rmax, aa, ai, aj, rp, ap, aimax, a->none…
194 rp[_i] = bcol; \
207 bcol = col / bs; \
214 if (rp[t] > bcol) high = t; \
218 if (rp[_i] > bcol) break; \
219 if (rp[_i] == bcol) { \
[all …]
|
| /petsc/src/mat/impls/aij/mpi/ |
| H A D | mpiaij.c | 1835 Vec col, bcol; in MatNorm_MPIAIJ() local
1846 PetscCall(MatCreateVecs(aij->B, &bcol, NULL)); in MatNorm_MPIAIJ()
1847 PetscCall(VecSet(bcol, 0.0)); in MatNorm_MPIAIJ()
1848 PetscCall(VecGetArrayWrite(bcol, &array)); in MatNorm_MPIAIJ()
1853 PetscCall(VecRestoreArrayWrite(bcol, &array)); in MatNorm_MPIAIJ()
1854 PetscCall(VecDestroy(&bcol)); in MatNorm_MPIAIJ()
|