| /petsc/src/mat/impls/aij/seq/cholmod/ |
| H A D | aijcholmod.c | 10 PetscInt m = A->rmap->n, i, j, k, nz, *ci, *cj; in MatWrapCholmod_seqaij() local
16 PetscCall(PetscMalloc2(m + 1, &ci, nz, &cj)); in MatWrapCholmod_seqaij()
25 cj[k] = aj[j]; in MatWrapCholmod_seqaij()
40 C->i = cj; in MatWrapCholmod_seqaij()
|
| /petsc/src/mat/impls/aij/seq/ |
| H A D | fdaij.c | 179 const PetscInt *is, *row, *ci, *cj; in MatFDColoringSetUp_SeqXAIJ() local
222 …PetscCall(MatGetColumnIJ_SeqBAIJ_Color(mat, 0, PETSC_FALSE, PETSC_FALSE, &ncols, &ci, &cj, &spidx,… in MatFDColoringSetUp_SeqXAIJ()
224 …PetscCall(MatGetColumnIJ_SeqSELL_Color(mat, 0, PETSC_FALSE, PETSC_FALSE, &ncols, &ci, &cj, &spidx,… in MatFDColoringSetUp_SeqXAIJ()
226 …PetscCall(MatGetColumnIJ_SeqAIJ_Color(mat, 0, PETSC_FALSE, PETSC_FALSE, &ncols, &ci, &cj, &spidx, … in MatFDColoringSetUp_SeqXAIJ()
248 row = cj + tmp; in MatFDColoringSetUp_SeqXAIJ()
286 …PetscCall(MatRestoreColumnIJ_SeqBAIJ_Color(mat, 0, PETSC_FALSE, PETSC_FALSE, &ncols, &ci, &cj, &sp… in MatFDColoringSetUp_SeqXAIJ()
289 …PetscCall(MatRestoreColumnIJ_SeqSELL_Color(mat, 0, PETSC_FALSE, PETSC_FALSE, &ncols, &ci, &cj, &sp… in MatFDColoringSetUp_SeqXAIJ()
291 …PetscCall(MatRestoreColumnIJ_SeqAIJ_Color(mat, 0, PETSC_FALSE, PETSC_FALSE, &ncols, &ci, &cj, &spi… in MatFDColoringSetUp_SeqXAIJ()
|
| H A D | matmatmult.c | 147 PetscInt *ai = a->i, *bi = b->i, *ci, *cj; in MatMatMultSymbolic_SeqAIJ_SeqAIJ_LLCondensed() local
210 PetscCall(PetscMalloc1(ci[am] + 1, &cj)); in MatMatMultSymbolic_SeqAIJ_SeqAIJ_LLCondensed()
211 PetscCall(PetscFreeSpaceContiguous(&free_space, cj)); in MatMatMultSymbolic_SeqAIJ_SeqAIJ_LLCondensed()
215 …PetscCall(MatSetSeqAIJWithArrays_private(PetscObjectComm((PetscObject)A), am, bn, ci, cj, NULL, ((… in MatMatMultSymbolic_SeqAIJ_SeqAIJ_LLCondensed()
252 PetscInt *ai = a->i, *aj = a->j, *bi = b->i, *bj = b->j, *bjj, *ci = c->i, *cj = c->j; in MatMatMultNumeric_SeqAIJ_SeqAIJ_Sorted() local
301 ca[k] += ab_dense[cj[k]]; in MatMatMultNumeric_SeqAIJ_SeqAIJ_Sorted()
302 ab_dense[cj[k]] = 0.0; /* zero ab_dense */ in MatMatMultNumeric_SeqAIJ_SeqAIJ_Sorted()
305 cj = PetscSafePointerPlusOffset(cj, cnzi); in MatMatMultNumeric_SeqAIJ_SeqAIJ_Sorted()
325 PetscInt *ai = a->i, *aj = a->j, *bi = b->i, *bj = b->j, *bjj, *ci = c->i, *cj = c->j; in MatMatMultNumeric_SeqAIJ_SeqAIJ_Scalable() local
358 if (cj[k] == bjj[nextb]) { /* ccol == bcol */ in MatMatMultNumeric_SeqAIJ_SeqAIJ_Scalable()
[all …]
|
| H A D | matptap.c | 69 PetscInt *ci, *cj, *ptadenserow, *ptasparserow, *ptaj, nspacedouble = 0; in MatPtAPSymbolic_SeqAIJ_SeqAIJ_SparseAxpy() local
148 PetscCall(PetscMalloc1(ci[pn], &cj)); in MatPtAPSymbolic_SeqAIJ_SeqAIJ_SparseAxpy()
149 PetscCall(PetscFreeSpaceContiguous(&free_space, cj)); in MatPtAPSymbolic_SeqAIJ_SeqAIJ_SparseAxpy()
156 …PetscCall(MatSetSeqAIJWithArrays_private(PetscObjectComm((PetscObject)A), pn, pn, ci, cj, ca, ((Pe… in MatPtAPSymbolic_SeqAIJ_SeqAIJ_SparseAxpy()
194 PetscInt *ci = c->i, *cj = c->j, *cjj; in MatPtAPNumeric_SeqAIJ_SeqAIJ_SparseAxpy() local
251 cjj = cj + ci[crow]; in MatPtAPNumeric_SeqAIJ_SeqAIJ_SparseAxpy()
|
| /petsc/src/dm/partitioner/impls/simple/ |
| H A D | partsimple.c | 50 PetscInt Np = 1, Nr, np, nk, nj, ni, pk, pj, pi, ck, cj, ci, i; in PetscPartitionerPartition_Simple_Grid() local
97 for (cj = 0; cj < pcells[1]; ++cj) { in PetscPartitionerPartition_Simple_Grid()
99 …] + pk) * pcells[2] + ck) * cells[1] + ((nj * procs[1] + pj) * pcells[1] + cj)) * cells[0] + (ni *… in PetscPartitionerPartition_Simple_Grid()
|
| /petsc/src/mat/impls/aij/seq/spqr/ |
| H A D | aijspqr.c | 22 SuiteSparse_long *ci, *cj; /* SuiteSparse_long is the only choice for SPQR */ in MatWrapCholmod_SPQR_seqaij() local
58 PetscCall(PetscMalloc2(n + 1, &cj, nz, &ci)); in MatWrapCholmod_SPQR_seqaij()
65 cj[j] = k; in MatWrapCholmod_SPQR_seqaij()
74 cj[j] = k; in MatWrapCholmod_SPQR_seqaij()
87 C->p = cj; in MatWrapCholmod_SPQR_seqaij()
|
| /petsc/src/ksp/pc/impls/ml/ |
| H A D | ml.c | 176 PetscInt *ci, *cj, ncols; in MatConvert_MPIAIJ_ML() local
186 PetscCall(PetscMalloc1(1 + ci[am], &cj)); in MatConvert_MPIAIJ_ML()
194 cj[k] = *aj++; in MatConvert_MPIAIJ_ML()
200 cj[k] = an + (*bj); in MatConvert_MPIAIJ_ML()
209 PetscCall(MatCreateSeqAIJWithArrays(PETSC_COMM_SELF, am, an, ci, cj, ca, Aloc)); in MatConvert_MPIAIJ_ML()
221 cj = mat->j; in MatConvert_MPIAIJ_ML()
|
| /petsc/src/mat/impls/maij/ |
| H A D | maij.c | 592 const PetscInt *ci = c->i, *cj = c->j, *cjj; in MatPtAPNumeric_SeqAIJ_SeqMAIJ() local
643 cjj = cj + ci[crow]; in MatPtAPNumeric_SeqAIJ_SeqMAIJ()
675 PetscInt *ci, *cj, *ptadenserow, *ptasparserow, *denserow, *sparserow, *ptaj; in MatPtAPSymbolic_SeqAIJ_SeqMAIJ() local
768 PetscCall(PetscMalloc1(ci[cn], &cj)); in MatPtAPSymbolic_SeqAIJ_SeqMAIJ()
769 PetscCall(PetscFreeSpaceContiguous(&free_space, cj)); in MatPtAPSymbolic_SeqAIJ_SeqMAIJ()
776 …PetscCall(MatSetSeqAIJWithArrays_private(PetscObjectComm((PetscObject)A), cn, cn, ci, cj, ca, NULL… in MatPtAPSymbolic_SeqAIJ_SeqMAIJ()
|
| /petsc/src/ksp/ksp/tutorials/ |
| H A D | ex49.c | 914 PetscInt ii, jj, ci, cj; in solve_elasticity_2d() local
931 cj = j - jj * (opts_t + opts_w + opts_t); in solve_elasticity_2d()
938 if ((cj >= opts_t) && (cj < opts_t + opts_w)) { in solve_elasticity_2d()
|
| /petsc/src/mat/impls/aij/mpi/ |
| H A D | mpimatmatmult.c | 1596 PetscInt i, j, k, anz, pnz, row, *cj, nexta; in MatTransposeMatMultNumeric_MPIAIJ_MPIAIJ() local
1657 cj = coj + coi[row]; in MatTransposeMatMultNumeric_MPIAIJ_MPIAIJ()
1663 if (cj[k] == adj[nexta]) { in MatTransposeMatMultNumeric_MPIAIJ_MPIAIJ()
1677 cj = bj + bi[row]; in MatTransposeMatMultNumeric_MPIAIJ_MPIAIJ()
1683 if (cj[k] == adj[nexta]) { in MatTransposeMatMultNumeric_MPIAIJ_MPIAIJ()
1732 cj = buf_rj[k] + *nextci[k]; in MatTransposeMatMultNumeric_MPIAIJ_MPIAIJ()
1736 if (bj_i[j] == cj[nextcj]) { /* bcol == ccol */ in MatTransposeMatMultNumeric_MPIAIJ_MPIAIJ()
|
| H A D | mpiaij.c | 5129 PetscInt *ci, *cj, col, ncols_d, ncols_o, jo; in MatMPIAIJGetLocalMat() local
5161 PetscCall(PetscMalloc1(1 + ci[am], &cj)); in MatMPIAIJGetLocalMat()
5171 cj[k] = col; in MatMPIAIJGetLocalMat()
5177 cj[k] = cstart + *aj++; in MatMPIAIJGetLocalMat()
5182 cj[k] = cmap[*bj++]; in MatMPIAIJGetLocalMat()
5187 PetscCall(MatCreateSeqAIJWithArrays(PETSC_COMM_SELF, am, A->cmap->N, ci, cj, ca, A_loc)); in MatMPIAIJGetLocalMat()
5197 cj = mat->j; in MatMPIAIJGetLocalMat()
5277 PetscInt *ci, *cj; in MatMPIAIJGetLocalMatMerge() local
5289 PetscCall(PetscMalloc1(ai[am] + bi[am], &cj)); in MatMPIAIJGetLocalMatMerge()
5298 cj[k] = *aj++; in MatMPIAIJGetLocalMatMerge()
[all …]
|
| /petsc/src/mat/impls/aij/seq/kokkos/ |
| H A D | aijkok.kokkos.cxx | 717 MatColIdxKokkosView aj, bj, cj; in MatSeqAIJKokkosMergeMats() local
748 auto cj = MatColIdxKokkosView("j", aj.extent(0) + bj.extent(0)); in MatSeqAIJKokkosMergeMats() local
764 cj(coffset + k) = aj(ai(i) + k); in MatSeqAIJKokkosMergeMats()
767 … cj(coffset + k) = bj(bi(i) + k - alen) + aN; /* Entries in B get new column indices in C */ in MatSeqAIJKokkosMergeMats()
771 PetscCallCXX(ckok = new Mat_SeqAIJKokkos(m, n, nnz, ci, cj, ca)); in MatSeqAIJKokkosMergeMats()
|
| /petsc/src/snes/interface/ |
| H A D | snes.c | 3168 const PetscInt *bj, *cj; in SNESComputeJacobian() local
3172 PetscCall(MatGetRow(Bfd, i, &cn, &cj, &ca)); in SNESComputeJacobian()
3199 PetscCall(MatRestoreRow(Bfd, i, &cn, &cj, &ca)); in SNESComputeJacobian()
|