Home
last modified time | relevance | path

Searched refs:aj (Results 1 – 25 of 113) sorted by relevance

12345

/petsc/src/mat/impls/sbaij/seq/
H A Dsro.c41 PetscInt *ai,*aj,*r; in MatReorderingSeqSBAIJ()
54 PetscCall(PetscMalloc2(mbs+1,&ai, 2*a->i[mbs],&aj)); in MatReorderingSeqSBAIJ()
56 ai = a->inew; aj = a->jnew; in MatReorderingSeqSBAIJ()
59 PetscCall(PetscArraycpy(aj,a->j,a->i[mbs])); in MatReorderingSeqSBAIJ()
80 k = aj[j]; /* col. index */ in MatReorderingSeqSBAIJ()
83 if (rip[k] < rip[i]) aj[j] = i; in MatReorderingSeqSBAIJ()
109 if (aj[j] == i) r[j] = ai[i]; /* put diagonal nonzero at beginning of row */ in MatReorderingSeqSBAIJ()
116 a->a2anew = aj + ai[mbs]; in MatReorderingSeqSBAIJ()
123 ajk = aj[k]; aj[k] = aj[j]; aj[j] = ajk; in MatReorderingSeqSBAIJ()
130 a->jnew = aj; in MatReorderingSeqSBAIJ()
H A Daijsbaij.c10 …PetscInt *ai = a->i, *aj = a->j, m = A->rmap->N, n = A->cmap->n, i, j, k, *bi, *bj, *rowlength… in MatConvert_SeqSBAIJ_SeqAIJ() local
25 if (*aj == i) { in MatConvert_SeqSBAIJ_SeqAIJ()
26 aj++; in MatConvert_SeqSBAIJ_SeqAIJ()
31 rowlengths[(*aj) * bs]++; in MatConvert_SeqSBAIJ_SeqAIJ()
32 aj++; in MatConvert_SeqSBAIJ_SeqAIJ()
66 aj = a->j; in MatConvert_SeqSBAIJ_SeqAIJ()
71 if (nz && *aj == i) { in MatConvert_SeqSBAIJ_SeqAIJ()
76 *(bj + rowstart[itmp]) = (*aj) * bs + k; in MatConvert_SeqSBAIJ_SeqAIJ()
81 aj++; in MatConvert_SeqSBAIJ_SeqAIJ()
88 itmp = (*aj) * bs + j; in MatConvert_SeqSBAIJ_SeqAIJ()
[all …]
H A Dsbaijfact2.c13 PetscInt mbs = a->mbs, *ai = a->i, *aj = a->j; in MatSolve_SeqSBAIJ_N_inplace() local
40 vj = aj + ai[k]; in MatSolve_SeqSBAIJ_N_inplace()
59 vj = aj + ai[k]; in MatSolve_SeqSBAIJ_N_inplace()
92 …orwardSolve_SeqSBAIJ_N_NaturalOrdering(const PetscInt *ai, const PetscInt *aj, const MatScalar *aa… in MatForwardSolve_SeqSBAIJ_N_NaturalOrdering() argument
107 vj = aj + ai[k]; in MatForwardSolve_SeqSBAIJ_N_NaturalOrdering()
124 …ckwardSolve_SeqSBAIJ_N_NaturalOrdering(const PetscInt *ai, const PetscInt *aj, const MatScalar *aa… in MatBackwardSolve_SeqSBAIJ_N_NaturalOrdering() argument
136 vj = aj + ai[k]; in MatBackwardSolve_SeqSBAIJ_N_NaturalOrdering()
152 const PetscInt mbs = a->mbs, *ai = a->i, *aj = a->j; in MatSolve_SeqSBAIJ_N_NaturalOrdering_inplace() local
164 PetscCall(MatForwardSolve_SeqSBAIJ_N_NaturalOrdering(ai, aj, aa, mbs, bs, x)); in MatSolve_SeqSBAIJ_N_NaturalOrdering_inplace()
167 PetscCall(MatBackwardSolve_SeqSBAIJ_N_NaturalOrdering(ai, aj, aa, mbs, bs, x)); in MatSolve_SeqSBAIJ_N_NaturalOrdering_inplace()
[all …]
H A Drelax.h135 const unsigned short *aj = a->jshort, *vj, *vj1; in MatSOR_SeqSBAIJ_ushort() local
137 const PetscInt *aj = a->j, *vj, *vj1; in MatSOR_SeqSBAIJ_ushort() local
170 vj = aj + ai[i] + 1; in MatSOR_SeqSBAIJ_ushort()
188 vj = aj + 1; in MatSOR_SeqSBAIJ_ushort()
203 vj = aj + ai[m] - 1; in MatSOR_SeqSBAIJ_ushort()
218 vj = aj + ai[m - 1] + 1; in MatSOR_SeqSBAIJ_ushort()
252 vj = aj + ai[i] + 1; in MatSOR_SeqSBAIJ_ushort()
277 vj = aj + ai[i] + 1; in MatSOR_SeqSBAIJ_ushort()
284 vj = aj + ai[i] + 1; in MatSOR_SeqSBAIJ_ushort()
/petsc/src/binding/petsc4py/test/
H A Dtest_mat_aij.py82 ai, aj, av = self._set_values()
84 self._chk_aij(self.A, ai, aj)
87 ai, aj, av = self._set_values_ijv()
89 self._chk_aij(self.A, ai, aj)
99 ai, aj, av = self._set_values()
101 self._chk_aij(self.A, ai, aj)
104 ai, aj, av = self._set_values_ijv()
106 self._chk_aij(self.A, ai, aj)
111 _, ai, aj, _ = self._get_aijv()
112 csr = [ai, aj]
[all …]
/petsc/src/mat/impls/aij/seq/
H A Dsymtranspose.c16 PetscInt *ati, *atj, *atfill, *ai = a->i, *aj = a->j; in MatTransposeSymbolic_SeqAIJ() local
25 for (i = 0; i < ai[am]; i++) ati[aj[i] + 1] += 1; in MatTransposeSymbolic_SeqAIJ()
37 atj[atfill[*aj]] = i; in MatTransposeSymbolic_SeqAIJ()
38 atfill[*aj++] += 1; in MatTransposeSymbolic_SeqAIJ()
61 PetscInt *ati, *atj, *atfill, *ai = a->i, *aj = a->j; in MatTranspose_SeqAIJ() local
84 for (i = 0; i < ai[am]; i++) ati[aj[i] + 1] += 1; in MatTranspose_SeqAIJ()
105 atj[atfill[*aj]] = i; in MatTranspose_SeqAIJ()
106 ata[atfill[*aj]] = *aa++; in MatTranspose_SeqAIJ()
107 atfill[*aj++] += 1; in MatTranspose_SeqAIJ()
114 atj[atfill[*aj]] = i; in MatTranspose_SeqAIJ()
[all …]
H A Dij.c31 PetscErrorCode MatToSymmetricIJ_SeqAIJ(PetscInt m, PetscInt *ai, PetscInt *aj, PetscBool lower_tria… in MatToSymmetricIJ_SeqAIJ() argument
45 j = aj + ai[row] + shiftin; in MatToSymmetricIJ_SeqAIJ()
73 j = aj + ai[row] + shiftin; in MatToSymmetricIJ_SeqAIJ()
H A Daijfact.c58 const PetscInt *r, *ic, *ai = a->i, *aj = a->j, *ajtmp; in MatLUFactorSymbolic_SeqAIJ() local
98 ajtmp = aj + ai[r[i]]; in MatLUFactorSymbolic_SeqAIJ()
224 const PetscInt n = A->rmap->n, *ai = a->i, *aj = a->j, *bi = b->i, *bj = b->j, *bdiag = b->diag; in MatLUFactorNumeric_SeqAIJ() local
281 ajtmp = aj + ai[r[i]]; in MatLUFactorNumeric_SeqAIJ()
404 const PetscInt *ai = a->i, *aj = a->j, *bi = b->i, *bj = b->j; in MatLUFactorNumeric_SeqAIJ_inplace() local
453 ajtmp = aj + ai[r[i]]; in MatLUFactorNumeric_SeqAIJ_inplace()
560 PetscInt i, j, n = A->rmap->n, *ai = a->i, *aj = a->j; in MatLUFactorNumeric_SeqAIJ_InplaceWithPerm() local
638 ajtmp = aj + ai[r[i]]; in MatLUFactorNumeric_SeqAIJ_InplaceWithPerm()
656 pj = aj + diag[r[row]] + 1; in MatLUFactorNumeric_SeqAIJ_InplaceWithPerm()
668 pj = aj + ai[r[i]]; in MatLUFactorNumeric_SeqAIJ_InplaceWithPerm()
[all …]
/petsc/src/mat/impls/baij/seq/
H A Dbaijsolvnat1.c11 const PetscInt n = a->mbs, *vi, *ai = a->i, *aj = a->j, *diag = a->diag; in MatSolve_SeqBAIJ_1_NaturalOrdering_inplace() local
27 vi = aj + ai[i]; in MatSolve_SeqBAIJ_1_NaturalOrdering_inplace()
42 vi = aj + diag[i] + 1; in MatSolve_SeqBAIJ_1_NaturalOrdering_inplace()
64 const PetscInt n = a->mbs, *ai = a->i, *aj = a->j, *vi; in MatForwardSolve_SeqBAIJ_1_NaturalOrdering() local
79 vi = aj; in MatForwardSolve_SeqBAIJ_1_NaturalOrdering()
97 const PetscInt n = a->mbs, *aj = a->j, *adiag = a->diag, *vi; in MatBackwardSolve_SeqBAIJ_1_NaturalOrdering() local
112 vi = aj + adiag[i + 1] + 1; in MatBackwardSolve_SeqBAIJ_1_NaturalOrdering()
128 const PetscInt n = a->mbs, *ai = a->i, *aj = a->j, *adiag = a->diag, *vi; in MatSolve_SeqBAIJ_1_NaturalOrdering() local
143 vi = aj; in MatSolve_SeqBAIJ_1_NaturalOrdering()
156 vi = aj + adiag[i + 1] + 1; in MatSolve_SeqBAIJ_1_NaturalOrdering()
H A Dbaijsolvtrannat1.c6 const PetscInt *adiag = a->diag, *ai = a->i, *aj = a->j, *vi; in MatSolveTranspose_SeqBAIJ_1_NaturalOrdering() local
24 vi = aj + adiag[i + 1] + 1; in MatSolveTranspose_SeqBAIJ_1_NaturalOrdering()
35 vi = aj + ai[i]; in MatSolveTranspose_SeqBAIJ_1_NaturalOrdering()
55 const PetscInt *diag = a->diag, n = a->mbs, *vi, *ai = a->i, *aj = a->j; in MatSolveTranspose_SeqBAIJ_1_NaturalOrdering_inplace() local
68 vi = aj + diag[i] + 1; in MatSolveTranspose_SeqBAIJ_1_NaturalOrdering_inplace()
76 vi = aj + diag[i] - 1; in MatSolveTranspose_SeqBAIJ_1_NaturalOrdering_inplace()
H A Dbaijsolvnat2.c15 const PetscInt n = a->mbs, *vi, *ai = a->i, *aj = a->j, *diag = a->diag; in MatSolve_SeqBAIJ_2_NaturalOrdering_inplace() local
31 vi = aj + ai[i]; in MatSolve_SeqBAIJ_2_NaturalOrdering_inplace()
50 vi = aj + diag[i] + 1; in MatSolve_SeqBAIJ_2_NaturalOrdering_inplace()
77 const PetscInt n = a->mbs, *vi, *ai = a->i, *aj = a->j, *adiag = a->diag; in MatSolve_SeqBAIJ_2_NaturalOrdering() local
92 vi = aj + ai[i]; in MatSolve_SeqBAIJ_2_NaturalOrdering()
114 vi = aj + adiag[i + 1] + 1; in MatSolve_SeqBAIJ_2_NaturalOrdering()
143 const PetscInt n = a->mbs, *vi, *ai = a->i, *aj = a->j; in MatForwardSolve_SeqBAIJ_2_NaturalOrdering() local
158 vi = aj + ai[i]; in MatForwardSolve_SeqBAIJ_2_NaturalOrdering()
186 const PetscInt n = a->mbs, *vi, *aj = a->j, *adiag = a->diag; in MatBackwardSolve_SeqBAIJ_2_NaturalOrdering() local
199 vi = aj + adiag[i + 1] + 1; in MatBackwardSolve_SeqBAIJ_2_NaturalOrdering()
H A Dbaijsolvtran1.c8 const PetscInt *rout, *cout, *r, *c, *adiag = a->diag, *ai = a->i, *aj = a->j, *vi; in MatSolveTranspose_SeqBAIJ_1() local
31 vi = aj + adiag[i + 1] + 1; in MatSolveTranspose_SeqBAIJ_1()
42 vi = aj + ai[i]; in MatSolveTranspose_SeqBAIJ_1()
65 const PetscInt *diag = a->diag, n = a->mbs, *vi, *ai = a->i, *aj = a->j; in MatSolveTranspose_SeqBAIJ_1_inplace() local
89 vi = aj + diag[i] + 1; in MatSolveTranspose_SeqBAIJ_1_inplace()
97 vi = aj + diag[i] - 1; in MatSolveTranspose_SeqBAIJ_1_inplace()
H A Dbaijsolvtrannat2.c7 const PetscInt *diag = a->diag, *vi, n = a->mbs, *ai = a->i, *aj = a->j; in MatSolveTranspose_SeqBAIJ_2_NaturalOrdering_inplace() local
26 vi = aj + diag[i] + 1; in MatSolveTranspose_SeqBAIJ_2_NaturalOrdering_inplace()
41 vi = aj + diag[i] - 1; in MatSolveTranspose_SeqBAIJ_2_NaturalOrdering_inplace()
61 const PetscInt n = a->mbs, *vi, *ai = a->i, *aj = a->j, *diag = a->diag; in MatSolveTranspose_SeqBAIJ_2_NaturalOrdering() local
82 vi = aj + diag[i] - 1; in MatSolveTranspose_SeqBAIJ_2_NaturalOrdering()
97 vi = aj + ai[i]; in MatSolveTranspose_SeqBAIJ_2_NaturalOrdering()
H A Dbaijsolvtrannat3.c6 const PetscInt n = a->mbs, *vi, *ai = a->i, *aj = a->j, *diag = a->diag; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering_inplace() local
28 vi = aj + diag[i] + 1; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering_inplace()
45 vi = aj + diag[i] - 1; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering_inplace()
67 const PetscInt n = a->mbs, *vi, *ai = a->i, *aj = a->j, *diag = a->diag; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering() local
90 vi = aj + diag[i] - 1; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering()
107 vi = aj + ai[i]; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering()
H A Dbaijsolvnat3.c11 const PetscInt n = a->mbs, *ai = a->i, *aj = a->j; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace() local
29 vi = aj + ai[i]; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
52 vi = aj + diag[i] + 1; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
83 const PetscInt n = a->mbs, *vi, *ai = a->i, *aj = a->j, *adiag = a->diag; in MatSolve_SeqBAIJ_3_NaturalOrdering() local
101 vi = aj + ai[i]; in MatSolve_SeqBAIJ_3_NaturalOrdering()
127 vi = aj + adiag[i + 1] + 1; in MatSolve_SeqBAIJ_3_NaturalOrdering()
160 const PetscInt n = a->mbs, *vi, *ai = a->i, *aj = a->j; in MatForwardSolve_SeqBAIJ_3_NaturalOrdering() local
178 vi = aj + ai[i]; in MatForwardSolve_SeqBAIJ_3_NaturalOrdering()
210 const PetscInt n = a->mbs, *vi, *aj = a->j, *adiag = a->diag; in MatBackwardSolve_SeqBAIJ_3_NaturalOrdering() local
225 vi = aj + adiag[i + 1] + 1; in MatBackwardSolve_SeqBAIJ_3_NaturalOrdering()
H A Dbaijsolvnat4.c12 const PetscInt *ai = a->i, *aj = a->j; in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace() local
25 fortransolvebaij4_(&n, x, ai, aj, diag, aa, b, w); in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace()
28 fortransolvebaij4unroll_(&n, x, ai, aj, diag, aa, b); in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace()
44 vi = aj + ai[i]; in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace()
73 vi = aj + diag[i] + 1; in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace()
110 const PetscInt n = a->mbs, *vi, *ai = a->i, *aj = a->j, *adiag = a->diag; in MatSolve_SeqBAIJ_4_NaturalOrdering() local
129 vi = aj + ai[i]; in MatSolve_SeqBAIJ_4_NaturalOrdering()
159 vi = aj + adiag[i + 1] + 1; in MatSolve_SeqBAIJ_4_NaturalOrdering()
H A Dbaijsolvtrann.c8 const PetscInt *r, *c, *rout, *cout, *ai = a->i, *aj = a->j, *vi; in MatSolveTranspose_SeqBAIJ_N_inplace() local
37 vi = aj + a->diag[i] + 1; in MatSolveTranspose_SeqBAIJ_N_inplace()
48 vi = aj + ai[i]; in MatSolveTranspose_SeqBAIJ_N_inplace()
74 const PetscInt n = a->mbs, *ai = a->i, *aj = a->j, *vi, *diag = a->diag; in MatSolveTranspose_SeqBAIJ_N() local
103 vi = aj + diag[i] - 1; in MatSolveTranspose_SeqBAIJ_N()
114 vi = aj + ai[i]; in MatSolveTranspose_SeqBAIJ_N()
H A Dbaijsolvtran2.c9 const PetscInt *diag = a->diag, n = a->mbs, *vi, *ai = a->i, *aj = a->j; in MatSolveTranspose_SeqBAIJ_2_inplace() local
45 vi = aj + diag[i] + 1; in MatSolveTranspose_SeqBAIJ_2_inplace()
60 vi = aj + diag[i] - 1; in MatSolveTranspose_SeqBAIJ_2_inplace()
94 const PetscInt n = a->mbs, *vi, *ai = a->i, *aj = a->j, *diag = a->diag; in MatSolveTranspose_SeqBAIJ_2() local
131 vi = aj + diag[i] - 1; in MatSolveTranspose_SeqBAIJ_2()
146 vi = aj + ai[i]; in MatSolveTranspose_SeqBAIJ_2()
/petsc/src/mat/impls/aij/seq/ftn-kernels/
H A Dfsolve.F9013 PetscInt, intent(in) :: n, ai(0:*), aj(0:*), adiag(0:*)
24 x(i) = b(i) - sum(aa(jstart:jend)*x(aj(jstart:jend)))
32 x(i) = x(i) - sum(aa(jstart:jend)*x(aj(jstart:jend)))*aa(adiag(i))
/petsc/src/mat/impls/aij/seq/aijperm/
H A Daijperm.c282 const PetscInt *aj, *ai; in MatMult_SeqAIJPERM() local
316 aj = a->j; /* aj[k] gives column index for element aa[k]. */ in MatMult_SeqAIJPERM()
339 y[iold] = aa[ipos] * x[aj[ipos]]; in MatMult_SeqAIJPERM()
381 vec_idx = _mm256_loadu_si256((__m256i const *)&aj[ipos]); in MatMult_SeqAIJPERM()
389 vec_idx = _mm256_loadu_si256((__m256i const *)&aj[ipos]); in MatMult_SeqAIJPERM()
394 yp[i] += aa[ipos] * x[aj[ipos]]; in MatMult_SeqAIJPERM()
395 yp[i] += aa[ipos + 1] * x[aj[ipos + 1]]; in MatMult_SeqAIJPERM()
397 yp[i] += aa[ipos] * x[aj[ipos]]; in MatMult_SeqAIJPERM()
403 yp[i] += aa[ipos] * x[aj[ipos]]; in MatMult_SeqAIJPERM()
418 vec_idx = _mm256_i32gather_epi32(aj, vec_ipos, _MM_SCALE_4); in MatMult_SeqAIJPERM()
[all …]
/petsc/src/mat/impls/baij/seq/ftn-kernels/
H A Dfsolvebaij.F9017 PetscInt, intent(in) :: ai(0:*), aj(0:*), adiag(0:*)
27 PETSC_AssertAlignx(16, aj(1))
42 jdx = 4*aj(j)
62 jdx = 4*aj(j)
87 PetscInt, intent(in) :: ai(0:*), aj(0:*), adiag(0:*)
98 PETSC_AssertAlignx(16, aj(1))
116 jdx = 4*aj(j)
151 jdx = 4*aj(j)
/petsc/src/mat/tests/
H A Dex108.c8 const PetscInt *ai, *aj; in main() local
43 PetscCall(MatGetRowIJ(A, 0, PETSC_FALSE, PETSC_FALSE, &n, &ai, &aj, &flg)); in main()
46 …PetscCall(MatCreateSeqBAIJWithArrays(PETSC_COMM_SELF, 1, 5, 5, (PetscInt *)ai, (PetscInt *)aj, aa,… in main()
48 PetscCall(MatRestoreRowIJ(A, 0, PETSC_FALSE, PETSC_FALSE, &n, &ai, &aj, &flg)); in main()
/petsc/src/tao/constrained/impls/ipm/
H A Dpdipm.c279 const PetscInt *aj, *ranges, *Jranges, *rranges, *cranges; in TaoSNESJacobian_PDIPM() local
319 PetscCall(MatGetRow(tao->jacobian_equality, i + rjstart, &nc, &aj, &aa)); in TaoSNESJacobian_PDIPM()
322 while (aj[j] >= cranges[proc + 1]) proc++; in TaoSNESJacobian_PDIPM()
323 cols[0] = aj[j] - cranges[proc] + Jranges[proc]; in TaoSNESJacobian_PDIPM()
326 PetscCall(MatRestoreRow(tao->jacobian_equality, i + rjstart, &nc, &aj, &aa)); in TaoSNESJacobian_PDIPM()
339 PetscCall(MatGetRow(tao->jacobian_inequality, i + rjstart, &nc, &aj, &aa)); in TaoSNESJacobian_PDIPM()
342 while (aj[j] >= cranges[proc + 1]) proc++; in TaoSNESJacobian_PDIPM()
343 cols[0] = aj[j] - cranges[proc] + Jranges[proc]; in TaoSNESJacobian_PDIPM()
346 PetscCall(MatRestoreRow(tao->jacobian_inequality, i + rjstart, &nc, &aj, &aa)); in TaoSNESJacobian_PDIPM()
371 PetscCall(MatGetRow(tao->hessian, i + rjstart, &nc, &aj, &aa)); in TaoSNESJacobian_PDIPM()
[all …]
/petsc/src/mat/impls/aij/seq/aijmkl/
H A Daijmkl.c118 PetscInt *aj, *ai; in MatSeqAIJMKL_create_mkl_handle()
142 aj = a->j; /* aj[k] gives column index for element aa[k]. */ in MatSeqAIJMKL_create_mkl_handle()
148 …SE_INDEX_BASE_ZERO, (MKL_INT)m, (MKL_INT)n, (MKL_INT *)ai, (MKL_INT *)(ai + 1), (MKL_INT *)aj, aa); in MatSeqAIJMKL_create_mkl_handle()
167 PetscInt *aj, *ai, *unused; in MatSeqAIJMKL_setup_structure_from_mkl_handle() local
174 …ndexing, (MKL_INT *)&m, (MKL_INT *)&n, (MKL_INT **)&ai, (MKL_INT **)&unused, (MKL_INT **)&aj, &aa); in MatSeqAIJMKL_setup_structure_from_mkl_handle()
177 aj = ai = NULL; in MatSeqAIJMKL_setup_structure_from_mkl_handle()
188 PetscCall(MatSeqAIJSetPreallocationCSR(A, ai, aj, NULL)); in MatSeqAIJMKL_setup_structure_from_mkl_handle()
227 PetscInt *ai, *aj, *unused; in MatSeqAIJMKL_update_from_mkl_handle() local
237 … (MKL_INT *)&nrows, (MKL_INT *)&ncols, (MKL_INT **)&ai, (MKL_INT **)&unused, (MKL_INT **)&aj, &aa); in MatSeqAIJMKL_update_from_mkl_handle()
243 PetscCall(MatSetValues_SeqAIJ(A, 1, &i, nz, aj + ai[i], aa + ai[i], INSERT_VALUES)); in MatSeqAIJMKL_update_from_mkl_handle()
[all …]
/petsc/src/mat/tutorials/
H A Dex12.c19 const PetscInt *aj; in PadMatrix() local
42 PetscCall(MatGetRow(A, i, &nc, &aj, &aa)); in PadMatrix()
43 PetscCall(MatSetValues(*B, 1, &i, nc, aj, aa, INSERT_VALUES)); in PadMatrix()
44 PetscCall(MatRestoreRow(A, i, &nc, &aj, &aa)); in PadMatrix()

12345