Home
last modified time | relevance | path

Searched refs:jmax (Results 1 – 22 of 22) sorted by relevance

/petsc/src/dm/tutorials/
H A Dex13f90aux.F9021 function dfdt_vdp(t, dt, ib1, ibn, jb1, jbn, kb1, kbn, imax, jmax, kmax, n, f)
27 PetscInt, intent(in) :: ib1, ibn, jb1, jbn, kb1, kbn, imax, jmax, kmax, n
29 PetscReal, dimension(n, imax, jmax, kmax) :: dfdt_vdp
38 …recursive subroutine forw_euler(t, dt, ib1, ibn, jb1, jbn, kb1, kbn, imax, jmax, kmax, neq, y, dfd… argument
40 PetscInt, intent(in) :: ib1, ibn, jb1, jbn, kb1, kbn, imax, jmax, kmax, neq
46 function dfdt(t, dt, ib1, ibn, jb1, jbn, kb1, kbn, imax, jmax, kmax, n, f)
49 PetscInt, intent(in) :: ib1, ibn, jb1, jbn, kb1, kbn, imax, jmax, kmax, n
51 PetscReal, dimension(n, imax, jmax, kmax) :: dfdt
56 …y(:, 1:imax, 1:jmax, 1:kmax) = y(:, 1:imax, 1:jmax, 1:kmax) + dt*dfdt(t, dt, ib1, ibn, jb1, jbn, k…
H A Dex13f90.F9037 …PetscInt :: i, j, k, igmax, jgmax, kgmax, ib1, ibn, jb1, jbn, kb1, kbn, imax, jmax, kmax, itime, m…
83 PetscCallA(DMDAGetCorners(SolScal, i, j, k, imax, jmax, kmax, ierr))
88 jbn = jmax + stw
104 xmax = grid(1, imax, jmax, kmax)
105 ymax = grid(2, imax, jmax, kmax)
106 zmax = grid(3, imax, jmax, kmax)
140 call forw_euler(t, dt, ib1, ibn, jb1, jbn, kb1, kbn, imax, jmax, kmax, dof, f, dfdt_vdp)
/petsc/src/mat/impls/sbaij/seq/
H A Dsbaijfact.c45 PetscInt *jl, *q, jmin, jmax, juidx, nzk, qm, *iu, *ju, k, j, vj, umax, maxadd; in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR() local
89 jmax = ai[rip[k] + 1]; in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
90 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
113 jmax = iu[prow + 1]; in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
115 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
228 PetscInt *jl, jmin, jmax, nzk, *ui, k, j, *il, nextprow; in MatCholeskyFactorSymbolic_SeqSBAIJ() local
289 jmax = ui[prow + 1]; in MatCholeskyFactorSymbolic_SeqSBAIJ()
290 ncols = jmax - jmin; in MatCholeskyFactorSymbolic_SeqSBAIJ()
296 if (jmin < jmax) { in MatCholeskyFactorSymbolic_SeqSBAIJ()
393 PetscInt *jl, jmin, jmax, nzk, *ui, k, j, *il, nextprow; in MatCholeskyFactorSymbolic_SeqSBAIJ_inplace() local
[all …]
H A Dsbaijfact4.c11 PetscInt *ai, *aj, k, k1, jmin, jmax, *jl, *il, vj, nexti, ili; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering() local
34 jmax = ai[k + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
35 if (jmin < jmax) { in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
37 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
90 jmax = bi[i + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
91 if (jmin < jmax) { in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
92 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
108 PetscCall(PetscLogFlops(2.0 * 27.0 * (jmax - jmin))); in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
128 jmax = bi[k + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
129 if (jmin < jmax) { in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
[all …]
H A Dsbaijfact3.c10 PetscInt *a2anew, i, j, k, k1, jmin, jmax, *jl, *il, vj, nexti, ili; in MatCholeskyFactorNumeric_SeqSBAIJ_3() local
42 jmax = ai[i + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_3()
43 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorNumeric_SeqSBAIJ_3()
72 jmax = ai[perm_ptr[k] + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_3()
73 if (jmin < jmax) { in MatCholeskyFactorNumeric_SeqSBAIJ_3()
75 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorNumeric_SeqSBAIJ_3()
128 jmax = bi[i + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_3()
129 if (jmin < jmax) { in MatCholeskyFactorNumeric_SeqSBAIJ_3()
130 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorNumeric_SeqSBAIJ_3()
146 PetscCall(PetscLogFlops(2.0 * 27.0 * (jmax - jmin))); in MatCholeskyFactorNumeric_SeqSBAIJ_3()
[all …]
H A Dsbaijfact6.c10 PetscInt i, j, *a2anew, k, k1, jmin, jmax, *jl, *il, vj, nexti, ili; in MatCholeskyFactorNumeric_SeqSBAIJ_4() local
43 jmax = ai[i + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_4()
44 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorNumeric_SeqSBAIJ_4()
73 jmax = ai[perm_ptr[k] + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_4()
74 if (jmin < jmax) { in MatCholeskyFactorNumeric_SeqSBAIJ_4()
76 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorNumeric_SeqSBAIJ_4()
145 jmax = bi[i + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_4()
146 if (jmin < jmax) { in MatCholeskyFactorNumeric_SeqSBAIJ_4()
147 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorNumeric_SeqSBAIJ_4()
171 PetscCall(PetscLogFlops(2.0 * 64.0 * (jmax - jmin))); in MatCholeskyFactorNumeric_SeqSBAIJ_4()
[all …]
H A Dsbaijfact5.c11 PetscInt *ai, *aj, k, k1, jmin, jmax, *jl, *il, vj, nexti, ili; in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering() local
35 jmax = ai[k + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
36 if (jmin < jmax) { in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
38 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
107 jmax = bi[i + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
108 if (jmin < jmax) { in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
109 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
133 PetscCall(PetscLogFlops(2.0 * 64.0 * (jmax - jmin))); in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
157 jmax = bi[k + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
158 if (jmin < jmax) { in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
[all …]
H A Dsro.c42 PetscInt *nzr,nz,jmin,jmax,j,k,ajk,i; in MatReorderingSeqSBAIJ()
104 jmin = ai[0]; jmax = ai[mbs]; in MatReorderingSeqSBAIJ()
105 nz = jmax - jmin; in MatReorderingSeqSBAIJ()
106 j = jmax-1; in MatReorderingSeqSBAIJ()
120 for (j=jmin; j<jmax; j++) { in MatReorderingSeqSBAIJ()
H A Dsbaijfact7.c10 PetscInt i, j, *a2anew, k, k1, jmin, jmax, *jl, *il, vj, nexti, ili, ipvt[5]; in MatCholeskyFactorNumeric_SeqSBAIJ_5() local
42 jmax = ai[i + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_5()
43 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorNumeric_SeqSBAIJ_5()
72 jmax = ai[perm_ptr[k] + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_5()
73 if (jmin < jmax) { in MatCholeskyFactorNumeric_SeqSBAIJ_5()
75 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorNumeric_SeqSBAIJ_5()
164 jmax = bi[i + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_5()
165 if (jmin < jmax) { in MatCholeskyFactorNumeric_SeqSBAIJ_5()
166 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorNumeric_SeqSBAIJ_5()
200 PetscCall(PetscLogFlops(2.0 * 125.0 * (jmax - jmin))); in MatCholeskyFactorNumeric_SeqSBAIJ_5()
[all …]
H A Dsbaijfact8.c11 PetscInt *ai, *aj, k, k1, jmin, jmax, *jl, *il, vj, nexti, ili, ipvt[5]; in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering() local
34 jmax = ai[k + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
35 if (jmin < jmax) { in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
37 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
126 jmax = bi[i + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
127 if (jmin < jmax) { in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
128 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
162 PetscCall(PetscLogFlops(2.0 * 125.0 * (jmax - jmin))); in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
182 jmax = bi[k + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
183 if (jmin < jmax) { in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
[all …]
H A Dsbaijfact9.c10 PetscInt i, j, *a2anew, k, k1, jmin, jmax, *jl, *il, vj, nexti, ili; in MatCholeskyFactorNumeric_SeqSBAIJ_6() local
44 jmax = ai[i + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_6()
45 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorNumeric_SeqSBAIJ_6()
74 jmax = ai[perm_ptr[k] + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_6()
75 if (jmin < jmax) { in MatCholeskyFactorNumeric_SeqSBAIJ_6()
77 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorNumeric_SeqSBAIJ_6()
227 jmax = bi[i + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_6()
228 if (jmin < jmax) { in MatCholeskyFactorNumeric_SeqSBAIJ_6()
229 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorNumeric_SeqSBAIJ_6()
313 PetscCall(PetscLogFlops(2.0 * 216.0 * (jmax - jmin))); in MatCholeskyFactorNumeric_SeqSBAIJ_6()
[all …]
H A Dsbaijfact10.c11 PetscInt *ai, *aj, k, k1, jmin, jmax, *jl, *il, vj, nexti, ili; in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering() local
42 jmax = ai[k + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
43 if (jmin < jmax) { in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
45 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
232 jmax = bi[i + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
233 if (jmin < jmax) { in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
234 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
318 PetscCall(PetscLogFlops(2.0 * 216.0 * (jmax - jmin))); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
338 jmax = bi[k + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
339 if (jmin < jmax) { in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
[all …]
H A Dsbaijfact11.c10 PetscInt *a2anew, k, k1, jmin, jmax, *jl, *il, vj, nexti, ili, i, j, mbs = a->mbs; in MatCholeskyFactorNumeric_SeqSBAIJ_7() local
45 jmax = ai[i + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_7()
46 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorNumeric_SeqSBAIJ_7()
75 jmax = ai[perm_ptr[k] + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_7()
76 if (jmin < jmax) { in MatCholeskyFactorNumeric_SeqSBAIJ_7()
78 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorNumeric_SeqSBAIJ_7()
269 jmax = bi[i + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_7()
270 if (jmin < jmax) { in MatCholeskyFactorNumeric_SeqSBAIJ_7()
271 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorNumeric_SeqSBAIJ_7()
382 PetscCall(PetscLogFlops(2.0 * 343.0 * (jmax - jmin))); in MatCholeskyFactorNumeric_SeqSBAIJ_7()
[all …]
H A Dsbaijfact12.c11 PetscInt *ai, *aj, k, k1, jmin, jmax, *jl, *il, vj, nexti, ili; in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering() local
37 jmax = ai[k + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
38 if (jmin < jmax) { in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
40 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
232 jmax = bi[i + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
233 if (jmin < jmax) { in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
234 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
345 PetscCall(PetscLogFlops(2.0 * 343.0 * (jmax - jmin))); in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
365 jmax = bi[k + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
366 if (jmin < jmax) { in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
[all …]
H A Dsbaij2.c1279 …PetscInt i, j, k, bs = A->rmap->bs, bs2 = a->bs2, k1, mbs = a->mbs, jmin, jmax, nexti, ik,… in MatNorm_SeqSBAIJ() local
1290 jmax = a->i[k + 1]; in MatNorm_SeqSBAIJ()
1292 if (jmax - jmin > 0 && *col == k) { /* diagonal block */ in MatNorm_SeqSBAIJ()
1299 for (j = jmin; j < jmax; j++) { /* off-diagonal blocks */ in MatNorm_SeqSBAIJ()
1332 jmax = a->i[i + 1]; in MatNorm_SeqSBAIJ()
1333 if (jmin < jmax) { in MatNorm_SeqSBAIJ()
1343 jmax = a->i[k + 1]; in MatNorm_SeqSBAIJ()
1344 for (i = jmin; i < jmax; i++) { in MatNorm_SeqSBAIJ()
1355 if (jmax - jmin > 0 && *col == k) jmin++; in MatNorm_SeqSBAIJ()
1356 if (jmin < jmax) { in MatNorm_SeqSBAIJ()
H A Dsbaijfact2.c2285 PetscInt *prowl, *q, jmin, jmax, juidx, nzk, qm, *iu, *ju, k, j, vj, umax, maxadd; in MatICCFactorSymbolic_SeqSBAIJ_MSR() local
2349 jmax = iu[prow + 1]; in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2351 for (j = jmin; j < jmax; j++) { in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2477 PetscInt jmin, jmax, nzk, k, j, *jl, prow, *il, nextprow; in MatICCFactorSymbolic_SeqSBAIJ() local
2560 jmax = ui[prow + 1]; in MatICCFactorSymbolic_SeqSBAIJ()
2561 ncols = jmax - jmin; in MatICCFactorSymbolic_SeqSBAIJ()
2571 if (jmin < jmax) { in MatICCFactorSymbolic_SeqSBAIJ()
2679 PetscInt jmin, jmax, nzk, k, j, *jl, prow, *il, nextprow; in MatICCFactorSymbolic_SeqSBAIJ_inplace() local
2757 jmax = ui[prow + 1]; in MatICCFactorSymbolic_SeqSBAIJ_inplace()
2758 ncols = jmax - jmin; in MatICCFactorSymbolic_SeqSBAIJ_inplace()
[all …]
/petsc/src/mat/impls/baij/mpi/
H A Dbaijov.c277 PetscInt is_no, ct1, max, *rbuf2_i, isz_i, *data_i, jmax; in MatIncreaseOverlap_MPIBAIJ_Once() local
285 jmax = rbuf2[idex][0]; in MatIncreaseOverlap_MPIBAIJ_Once()
286 for (PetscInt j = 1; j <= jmax; j++) { in MatIncreaseOverlap_MPIBAIJ_Once()
625 PetscInt **lens, is_no, ncols, *cols, mat_i, *mat_j, tmp2, jmax; local
758 jmax = nrow[i];
761 PetscCall(PetscMalloc1(jmax, &row2proc_i));
765 for (PetscInt j = 0; j < jmax; j++) {
844 jmax = nrow[i];
845 for (PetscInt j = 0; j < jmax; j++) { /* parse the indices of each IS */
982 jmax = ncol[i];
[all …]
/petsc/src/mat/impls/dense/mpi/
H A Dmmdense.c61 PetscInt is_no, jmax, **rmap, *rmap_i; in MatCreateSubMatrices_MPIDense_Local() local
97 jmax = C->rmap->range[i + 1]; in MatCreateSubMatrices_MPIDense_Local()
98 for (; j < jmax; j++) rtable[j] = i; in MatCreateSubMatrices_MPIDense_Local()
108 jmax = nrow[i]; in MatCreateSubMatrices_MPIDense_Local()
110 for (PetscInt j = 0; j < jmax; j++) { in MatCreateSubMatrices_MPIDense_Local()
187 jmax = nrow[i]; in MatCreateSubMatrices_MPIDense_Local()
188 for (PetscInt j = 0; j < jmax; j++) { /* parse the indices of each IS */ in MatCreateSubMatrices_MPIDense_Local()
326 jmax = nrow[i]; in MatCreateSubMatrices_MPIDense_Local()
327 for (PetscInt j = 0; j < jmax; j++) rmap_i[irow_i[j]] = j; in MatCreateSubMatrices_MPIDense_Local()
/petsc/src/mat/impls/baij/seq/
H A Dbaijfact.c755 PetscInt k, jmin, jmax, *jl, *il, col, nexti, ili, nz; in MatCholeskyFactorNumeric_SeqBAIJ_N() local
792 jmax = ai[rip[k] + 1]; in MatCholeskyFactorNumeric_SeqBAIJ_N()
793 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorNumeric_SeqBAIJ_N()
819 jmax = bi[i + 1]; in MatCholeskyFactorNumeric_SeqBAIJ_N()
820 if (jmin < jmax) { in MatCholeskyFactorNumeric_SeqBAIJ_N()
821 for (j = jmin; j < jmax; j++) rtmp[bj[j]] += uikdi * ba[j]; in MatCholeskyFactorNumeric_SeqBAIJ_N()
853 jmax = bi[k + 1]; in MatCholeskyFactorNumeric_SeqBAIJ_N()
854 if (jmin < jmax) { in MatCholeskyFactorNumeric_SeqBAIJ_N()
855 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorNumeric_SeqBAIJ_N()
1024 PetscInt jmin, jmax, nzk, k, j, *jl, prow, *il, nextprow; in MatICCFactorSymbolic_SeqBAIJ() local
[all …]
/petsc/src/mat/impls/aij/seq/
H A Daijfact.c1712 PetscInt k, jmin, jmax, *c2r, *il, col, nexti, ili, nz; in MatCholeskyFactorNumeric_SeqAIJ() local
1768 jmax = ai[rip[k] + 1]; in MatCholeskyFactorNumeric_SeqAIJ()
1769 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorNumeric_SeqAIJ()
1794 jmax = bi[i + 1]; in MatCholeskyFactorNumeric_SeqAIJ()
1795 if (jmin < jmax) { in MatCholeskyFactorNumeric_SeqAIJ()
1796 for (j = jmin; j < jmax; j++) rtmp[bj[j]] += uikdi * ba[j]; in MatCholeskyFactorNumeric_SeqAIJ()
1809 jmax = bi[k + 1] - 1; in MatCholeskyFactorNumeric_SeqAIJ()
1810 if (jmin < jmax) { in MatCholeskyFactorNumeric_SeqAIJ()
1811 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorNumeric_SeqAIJ()
1879 PetscInt k, jmin, jmax, *jl, *il, col, nexti, ili, nz; in MatCholeskyFactorNumeric_SeqAIJ_inplace() local
[all …]
H A Dinode.c58 const PetscInt *j, *jmax, *ai = a->i, *aj = a->j; in MatGetRowIJ_SeqAIJ_Inode_Symmetric() local
89 jmax = aj + ai[row + 1] + ishift; in MatGetRowIJ_SeqAIJ_Inode_Symmetric()
90 if (j == jmax) continue; /* empty row */ in MatGetRowIJ_SeqAIJ_Inode_Symmetric()
93 while (i2 < i1 && j < jmax) { /* 1.[-xx-d-xx--] 2.[-xx-------],off-diagonal elements */ in MatGetRowIJ_SeqAIJ_Inode_Symmetric()
97 while ((j < jmax) && ((col = *j + ishift) < tns[i2])) ++j; in MatGetRowIJ_SeqAIJ_Inode_Symmetric()
119 jmax = aj + ai[row + 1] + ishift; in MatGetRowIJ_SeqAIJ_Inode_Symmetric()
120 if (j == jmax) continue; /* empty row */ in MatGetRowIJ_SeqAIJ_Inode_Symmetric()
123 while (i2 < i1 && j < jmax) { in MatGetRowIJ_SeqAIJ_Inode_Symmetric()
127 …while ((j < jmax) && ((col = *j + ishift) < tns[i2])) ++j; /* Skip rest col indices in this node */ in MatGetRowIJ_SeqAIJ_Inode_Symmetric()
/petsc/src/mat/impls/aij/mpi/
H A Dmpiov.c740 PetscInt is_no, ct1, max, *rbuf2_i, isz_i, jmax; in MatIncreaseOverlap_MPIAIJ_Once() local
749 jmax = rbuf2[idex][0]; in MatIncreaseOverlap_MPIAIJ_Once()
750 for (PetscInt j = 1; j <= jmax; j++) { in MatIncreaseOverlap_MPIAIJ_Once()
2046 PetscInt **lens, is_no, ncols, *cols, mat_i, *mat_j, tmp2, jmax; in MatCreateSubMatrices_MPIAIJ_Local() local
2161 jmax = nrow[i]; in MatCreateSubMatrices_MPIAIJ_Local()
2164 PetscCall(PetscMalloc1(jmax, &row2proc_i)); in MatCreateSubMatrices_MPIAIJ_Local()
2168 for (PetscInt j = 0; j < jmax; j++) { in MatCreateSubMatrices_MPIAIJ_Local()
2247 jmax = nrow[i]; in MatCreateSubMatrices_MPIAIJ_Local()
2248 for (PetscInt j = 0; j < jmax; j++) { /* parse the indices of each IS */ in MatCreateSubMatrices_MPIAIJ_Local()
2392 jmax = ncol[i]; in MatCreateSubMatrices_MPIAIJ_Local()
[all …]