Home
last modified time | relevance | path

Searched refs:jmin (Results 1 – 17 of 17) sorted by relevance

/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
88 jmin = ai[rip[k]] + 1; /* exclude diag[k] */ in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
90 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
112 jmin = 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
288 jmin = il[prow] + 1; /* index of the 2nd nzero entry in U(prow,k:mbs-1) */ in MatCholeskyFactorSymbolic_SeqSBAIJ()
290 ncols = jmax - jmin; in MatCholeskyFactorSymbolic_SeqSBAIJ()
291 … uj_ptr = ui_ptr[prow] + jmin - ui[prow]; /* points to the 2nd nzero entry in U(prow,k:mbs-1) */ in MatCholeskyFactorSymbolic_SeqSBAIJ()
296 if (jmin < jmax) { in MatCholeskyFactorSymbolic_SeqSBAIJ()
[all …]
H A Dsbaijfact4.c11 PetscInt *ai, *aj, k, k1, jmin, jmax, *jl, *il, vj, nexti, ili; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering() local
33 jmin = ai[k]; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
35 if (jmin < jmax) { in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
36 ap = aa + jmin * 9; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
37 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
89 jmin = ili + 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()
111 il[i] = jmin; /* update il(i) in column k+1, ... mbs-1 */ 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
41 jmin = ai[i]; in MatCholeskyFactorNumeric_SeqSBAIJ_3()
43 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorNumeric_SeqSBAIJ_3()
71 jmin = ai[perm_ptr[k]]; in MatCholeskyFactorNumeric_SeqSBAIJ_3()
73 if (jmin < jmax) { in MatCholeskyFactorNumeric_SeqSBAIJ_3()
74 ap = aa + jmin * 9; in MatCholeskyFactorNumeric_SeqSBAIJ_3()
75 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorNumeric_SeqSBAIJ_3()
127 jmin = ili + 1; in MatCholeskyFactorNumeric_SeqSBAIJ_3()
129 if (jmin < jmax) { in MatCholeskyFactorNumeric_SeqSBAIJ_3()
130 for (j = jmin; j < jmax; j++) { 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
42 jmin = ai[i]; in MatCholeskyFactorNumeric_SeqSBAIJ_4()
44 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorNumeric_SeqSBAIJ_4()
72 jmin = ai[perm_ptr[k]]; in MatCholeskyFactorNumeric_SeqSBAIJ_4()
74 if (jmin < jmax) { in MatCholeskyFactorNumeric_SeqSBAIJ_4()
75 ap = aa + jmin * 16; in MatCholeskyFactorNumeric_SeqSBAIJ_4()
76 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorNumeric_SeqSBAIJ_4()
144 jmin = ili + 1; in MatCholeskyFactorNumeric_SeqSBAIJ_4()
146 if (jmin < jmax) { in MatCholeskyFactorNumeric_SeqSBAIJ_4()
147 for (j = jmin; j < jmax; j++) { 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
34 jmin = ai[k]; in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
36 if (jmin < jmax) { in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
37 ap = aa + jmin * 16; in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
38 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
106 jmin = ili + 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()
136 il[i] = jmin; /* update il(i) in column k+1, ... mbs-1 */ in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
[all …]
H A Dsbaijfact7.c10 PetscInt i, j, *a2anew, k, k1, jmin, jmax, *jl, *il, vj, nexti, ili, ipvt[5]; in MatCholeskyFactorNumeric_SeqSBAIJ_5() local
41 jmin = ai[i]; in MatCholeskyFactorNumeric_SeqSBAIJ_5()
43 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorNumeric_SeqSBAIJ_5()
71 jmin = ai[perm_ptr[k]]; in MatCholeskyFactorNumeric_SeqSBAIJ_5()
73 if (jmin < jmax) { in MatCholeskyFactorNumeric_SeqSBAIJ_5()
74 ap = aa + jmin * 25; in MatCholeskyFactorNumeric_SeqSBAIJ_5()
75 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorNumeric_SeqSBAIJ_5()
163 jmin = ili + 1; in MatCholeskyFactorNumeric_SeqSBAIJ_5()
165 if (jmin < jmax) { in MatCholeskyFactorNumeric_SeqSBAIJ_5()
166 for (j = jmin; j < jmax; j++) { 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
33 jmin = ai[k]; in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
35 if (jmin < jmax) { in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
36 ap = aa + jmin * 25; in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
37 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
125 jmin = ili + 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()
165 il[i] = jmin; /* update il(i) in column k+1, ... mbs-1 */ in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
[all …]
H A Dsbaij2.c293 PetscInt mbs = a->mbs, i, n, cval, j, jmin; in MatMult_SeqSBAIJ_2() local
311 jmin = 0; in MatMult_SeqSBAIJ_2()
317 jmin++; in MatMult_SeqSBAIJ_2()
319 …PetscPrefetchBlock(ib + jmin + n, n, 0, PETSC_PREFETCH_HINT_NTA); /* Indices for the next row (ass… in MatMult_SeqSBAIJ_2()
321 for (j = jmin; j < n; j++) { in MatMult_SeqSBAIJ_2()
347 PetscInt mbs = a->mbs, i, n, cval, j, jmin; in MatMult_SeqSBAIJ_3() local
366 jmin = 0; in MatMult_SeqSBAIJ_3()
373 jmin++; in MatMult_SeqSBAIJ_3()
375 …PetscPrefetchBlock(ib + jmin + n, n, 0, PETSC_PREFETCH_HINT_NTA); /* Indices for the next row (ass… in MatMult_SeqSBAIJ_3()
377 for (j = jmin; j < n; j++) { in MatMult_SeqSBAIJ_3()
[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()
120 for (j=jmin; j<jmax; j++) { in MatReorderingSeqSBAIJ()
H A Dsbaijfact9.c10 PetscInt i, j, *a2anew, k, k1, jmin, jmax, *jl, *il, vj, nexti, ili; in MatCholeskyFactorNumeric_SeqSBAIJ_6() local
43 jmin = ai[i]; in MatCholeskyFactorNumeric_SeqSBAIJ_6()
45 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorNumeric_SeqSBAIJ_6()
73 jmin = ai[perm_ptr[k]]; in MatCholeskyFactorNumeric_SeqSBAIJ_6()
75 if (jmin < jmax) { in MatCholeskyFactorNumeric_SeqSBAIJ_6()
76 ap = aa + jmin * 36; in MatCholeskyFactorNumeric_SeqSBAIJ_6()
77 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorNumeric_SeqSBAIJ_6()
226 jmin = ili + 1; in MatCholeskyFactorNumeric_SeqSBAIJ_6()
228 if (jmin < jmax) { in MatCholeskyFactorNumeric_SeqSBAIJ_6()
229 for (j = jmin; j < jmax; j++) { 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
41 jmin = ai[k]; in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
43 if (jmin < jmax) { in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
44 ap = aa + jmin * 36; in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
45 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
231 jmin = ili + 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()
321 il[i] = jmin; /* update il(i) in column k+1, ... mbs-1 */ 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
44 jmin = ai[i]; in MatCholeskyFactorNumeric_SeqSBAIJ_7()
46 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorNumeric_SeqSBAIJ_7()
74 jmin = ai[perm_ptr[k]]; in MatCholeskyFactorNumeric_SeqSBAIJ_7()
76 if (jmin < jmax) { in MatCholeskyFactorNumeric_SeqSBAIJ_7()
77 ap = aa + jmin * 49; in MatCholeskyFactorNumeric_SeqSBAIJ_7()
78 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorNumeric_SeqSBAIJ_7()
268 jmin = ili + 1; in MatCholeskyFactorNumeric_SeqSBAIJ_7()
270 if (jmin < jmax) { in MatCholeskyFactorNumeric_SeqSBAIJ_7()
271 for (j = jmin; j < jmax; j++) { 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
36 jmin = ai[k]; in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
38 if (jmin < jmax) { in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
39 ap = aa + jmin * 49; in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
40 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
231 jmin = ili + 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()
348 il[i] = jmin; /* update il(i) in column k+1, ... mbs-1 */ in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
[all …]
H A Drelax.h73 PetscInt nonzerorow = 0, jmin; in MatMult_SeqSBAIJ_1_ushort() local
87 jmin = 0; in MatMult_SeqSBAIJ_1_ushort()
91 jmin++; in MatMult_SeqSBAIJ_1_ushort()
96 for (j = jmin; j < nz; j++) { in MatMult_SeqSBAIJ_1_ushort()
103 for (j = jmin; j < nz; j++) { in MatMult_SeqSBAIJ_1_ushort()
H A Dsbaijfact2.c2285 PetscInt *prowl, *q, jmin, jmax, juidx, nzk, qm, *iu, *ju, k, j, vj, umax, maxadd; in MatICCFactorSymbolic_SeqSBAIJ_MSR() local
2348 jmin = 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
2559 jmin = il[prow] + 1; /* index of the 2nd nzero entry in U(prow,k:am-1) */ in MatICCFactorSymbolic_SeqSBAIJ()
2561 ncols = jmax - jmin; in MatICCFactorSymbolic_SeqSBAIJ()
2562 i = jmin - ui[prow]; in MatICCFactorSymbolic_SeqSBAIJ()
2571 if (jmin < jmax) { in MatICCFactorSymbolic_SeqSBAIJ()
2572 il[prow] = jmin; in MatICCFactorSymbolic_SeqSBAIJ()
2679 PetscInt jmin, jmax, nzk, k, j, *jl, prow, *il, nextprow; in MatICCFactorSymbolic_SeqSBAIJ_inplace() local
[all …]
/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
791 jmin = ai[rip[k]]; in MatCholeskyFactorNumeric_SeqBAIJ_N()
793 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorNumeric_SeqBAIJ_N()
818 jmin = ili + 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()
823 il[i] = jmin; in MatCholeskyFactorNumeric_SeqBAIJ_N()
824 j = bj[jmin]; in MatCholeskyFactorNumeric_SeqBAIJ_N()
834 jmin = bi[k] + 1; in MatCholeskyFactorNumeric_SeqBAIJ_N()
835 nz = bi[k + 1] - jmin; in MatCholeskyFactorNumeric_SeqBAIJ_N()
[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
1767 jmin = ai[rip[k]]; in MatCholeskyFactorNumeric_SeqAIJ()
1769 for (j = jmin; j < jmax; j++) { in MatCholeskyFactorNumeric_SeqAIJ()
1793 jmin = ili + 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()
1798 il[i] = jmin; in MatCholeskyFactorNumeric_SeqAIJ()
1799 j = bj[jmin]; in MatCholeskyFactorNumeric_SeqAIJ()
1808 jmin = bi[k]; in MatCholeskyFactorNumeric_SeqAIJ()
1810 if (jmin < jmax) { in MatCholeskyFactorNumeric_SeqAIJ()
[all …]