| /petsc/src/mat/impls/sbaij/seq/ |
| H A D | sbaijfact.c | 45 PetscInt *jl, *q, jmin, jmax, juidx, nzk, qm, *iu, *ju, k, j, vj, umax, maxadd; in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR() local
76 PetscCall(PetscMalloc2(mbs, &jl, mbs, &q)); in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
78 jl[i] = mbs; in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
108 prow = jl[prow]; /* next pivot row (== mbs for symbolic factorization) */ in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
128 prow = jl[prow]; /* next pivot row */ in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
134 jl[k] = jl[i]; in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
135 jl[i] = k; in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
177 PetscCall(PetscFree2(jl, q)); in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
228 PetscInt *jl, jmin, jmax, nzk, *ui, k, j, *il, nextprow; in MatCholeskyFactorSymbolic_SeqSBAIJ() local
256 PetscCall(PetscMalloc4(mbs, &ui_ptr, mbs, &il, mbs, &jl, mbs, &cols)); in MatCholeskyFactorSymbolic_SeqSBAIJ()
[all …]
|
| H A D | sbaijfact4.c | 11 PetscInt *ai, *aj, k, k1, jmin, jmax, *jl, *il, vj, nexti, ili; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering() local
21 PetscCall(PetscMalloc2(mbs, &il, mbs, &jl)); in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
23 for (i = 0; i < mbs; i++) jl[i] = mbs; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
46 i = jl[k]; /* first row to be added to k_th row */ in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
49 nexti = jl[i]; /* next row to be added to k_th row */ in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
113 jl[i] = jl[j]; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
114 jl[j] = i; /* update jl */ in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
143 jl[k] = jl[i]; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
144 jl[i] = k; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
149 PetscCall(PetscFree2(il, jl)); in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
|
| H A D | sbaijfact3.c | 10 PetscInt *a2anew, i, j, k, k1, jmin, jmax, *jl, *il, vj, nexti, ili; in MatCholeskyFactorNumeric_SeqSBAIJ_3() local
20 PetscCall(PetscMalloc2(mbs, &il, mbs, &jl)); in MatCholeskyFactorNumeric_SeqSBAIJ_3()
22 for (i = 0; i < mbs; i++) jl[i] = mbs; in MatCholeskyFactorNumeric_SeqSBAIJ_3()
84 i = jl[k]; /* first row to be added to k_th row */ in MatCholeskyFactorNumeric_SeqSBAIJ_3()
87 nexti = jl[i]; /* next row to be added to k_th row */ in MatCholeskyFactorNumeric_SeqSBAIJ_3()
151 jl[i] = jl[j]; in MatCholeskyFactorNumeric_SeqSBAIJ_3()
152 jl[j] = i; /* update jl */ in MatCholeskyFactorNumeric_SeqSBAIJ_3()
181 jl[k] = jl[i]; in MatCholeskyFactorNumeric_SeqSBAIJ_3()
182 jl[i] = k; in MatCholeskyFactorNumeric_SeqSBAIJ_3()
187 PetscCall(PetscFree2(il, jl)); in MatCholeskyFactorNumeric_SeqSBAIJ_3()
|
| H A D | sbaijfact5.c | 11 PetscInt *ai, *aj, k, k1, jmin, jmax, *jl, *il, vj, nexti, ili; in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering() local
22 PetscCall(PetscMalloc2(mbs, &il, mbs, &jl)); in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
24 for (i = 0; i < mbs; i++) jl[i] = mbs; in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
47 i = jl[k]; /* first row to be added to k_th row */ in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
50 nexti = jl[i]; /* next row to be added to k_th row */ in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
138 jl[i] = jl[j]; in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
139 jl[j] = i; /* update jl */ in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
172 jl[k] = jl[i]; in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
173 jl[i] = k; in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
178 PetscCall(PetscFree2(il, jl)); in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
|
| H A D | sbaijfact6.c | 10 PetscInt i, j, *a2anew, k, k1, jmin, jmax, *jl, *il, vj, nexti, ili; in MatCholeskyFactorNumeric_SeqSBAIJ_4() local
21 PetscCall(PetscMalloc2(mbs, &il, mbs, &jl)); in MatCholeskyFactorNumeric_SeqSBAIJ_4()
23 for (i = 0; i < mbs; i++) jl[i] = mbs; in MatCholeskyFactorNumeric_SeqSBAIJ_4()
85 i = jl[k]; /* first row to be added to k_th row */ in MatCholeskyFactorNumeric_SeqSBAIJ_4()
88 nexti = jl[i]; /* next row to be added to k_th row */ in MatCholeskyFactorNumeric_SeqSBAIJ_4()
176 jl[i] = jl[j]; in MatCholeskyFactorNumeric_SeqSBAIJ_4()
177 jl[j] = i; /* update jl */ in MatCholeskyFactorNumeric_SeqSBAIJ_4()
211 jl[k] = jl[i]; in MatCholeskyFactorNumeric_SeqSBAIJ_4()
212 jl[i] = k; in MatCholeskyFactorNumeric_SeqSBAIJ_4()
217 PetscCall(PetscFree2(il, jl)); in MatCholeskyFactorNumeric_SeqSBAIJ_4()
|
| H A D | sbaijfact8.c | 11 PetscInt *ai, *aj, k, k1, jmin, jmax, *jl, *il, vj, nexti, ili, ipvt[5]; in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering() local
21 PetscCall(PetscMalloc2(mbs, &il, mbs, &jl)); in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
23 for (i = 0; i < mbs; i++) jl[i] = mbs; in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
46 i = jl[k]; /* first row to be added to k_th row */ in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
49 nexti = jl[i]; /* next row to be added to k_th row */ in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
167 jl[i] = jl[j]; in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
168 jl[j] = i; /* update jl */ in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
197 jl[k] = jl[i]; in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
198 jl[i] = k; in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
203 PetscCall(PetscFree2(il, jl)); in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
|
| H A D | sbaijfact7.c | 10 PetscInt i, j, *a2anew, k, k1, jmin, jmax, *jl, *il, vj, nexti, ili, ipvt[5]; in MatCholeskyFactorNumeric_SeqSBAIJ_5() local
20 PetscCall(PetscMalloc2(mbs, &il, mbs, &jl)); in MatCholeskyFactorNumeric_SeqSBAIJ_5()
22 for (i = 0; i < mbs; i++) jl[i] = mbs; in MatCholeskyFactorNumeric_SeqSBAIJ_5()
84 i = jl[k]; /* first row to be added to k_th row */ in MatCholeskyFactorNumeric_SeqSBAIJ_5()
87 nexti = jl[i]; /* next row to be added to k_th row */ in MatCholeskyFactorNumeric_SeqSBAIJ_5()
205 jl[i] = jl[j]; in MatCholeskyFactorNumeric_SeqSBAIJ_5()
206 jl[j] = i; /* update jl */ in MatCholeskyFactorNumeric_SeqSBAIJ_5()
235 jl[k] = jl[i]; in MatCholeskyFactorNumeric_SeqSBAIJ_5()
236 jl[i] = k; in MatCholeskyFactorNumeric_SeqSBAIJ_5()
241 PetscCall(PetscFree2(il, jl)); in MatCholeskyFactorNumeric_SeqSBAIJ_5()
|
| H A D | sbaijfact9.c | 10 PetscInt i, j, *a2anew, k, k1, jmin, jmax, *jl, *il, vj, nexti, ili; in MatCholeskyFactorNumeric_SeqSBAIJ_6() local
22 PetscCall(PetscMalloc2(mbs, &il, mbs, &jl)); in MatCholeskyFactorNumeric_SeqSBAIJ_6()
24 for (i = 0; i < mbs; i++) jl[i] = mbs; in MatCholeskyFactorNumeric_SeqSBAIJ_6()
86 i = jl[k]; /* first row to be added to k_th row */ in MatCholeskyFactorNumeric_SeqSBAIJ_6()
89 nexti = jl[i]; /* next row to be added to k_th row */ in MatCholeskyFactorNumeric_SeqSBAIJ_6()
318 jl[i] = jl[j]; in MatCholeskyFactorNumeric_SeqSBAIJ_6()
319 jl[j] = i; /* update jl */ in MatCholeskyFactorNumeric_SeqSBAIJ_6()
348 jl[k] = jl[i]; in MatCholeskyFactorNumeric_SeqSBAIJ_6()
349 jl[i] = k; in MatCholeskyFactorNumeric_SeqSBAIJ_6()
354 PetscCall(PetscFree2(il, jl)); in MatCholeskyFactorNumeric_SeqSBAIJ_6()
|
| H A D | sbaijfact10.c | 11 PetscInt *ai, *aj, k, k1, jmin, jmax, *jl, *il, vj, nexti, ili; in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering() local
29 PetscCall(PetscMalloc2(mbs, &il, mbs, &jl)); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
31 for (i = 0; i < mbs; i++) jl[i] = mbs; in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
54 i = jl[k]; /* first row to be added to k_th row */ in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
57 nexti = jl[i]; /* next row to be added to k_th row */ in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
323 jl[i] = jl[j]; in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
324 jl[j] = i; /* update jl */ in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
353 jl[k] = jl[i]; in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
354 jl[i] = k; in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
359 PetscCall(PetscFree2(il, jl)); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
|
| H A D | sbaijfact12.c | 11 PetscInt *ai, *aj, k, k1, jmin, jmax, *jl, *il, vj, nexti, ili; in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering() local
24 PetscCall(PetscMalloc2(mbs, &il, mbs, &jl)); in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
26 for (i = 0; i < mbs; i++) jl[i] = mbs; in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
50 i = jl[k]; /* first row to be added to k_th row */ in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
53 nexti = jl[i]; /* next row to be added to k_th row */ in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
350 jl[i] = jl[j]; in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
351 jl[j] = i; /* update jl */ in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
380 jl[k] = jl[i]; in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
381 jl[i] = k; in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
386 PetscCall(PetscFree2(il, jl)); in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
|
| H A D | sbaijfact11.c | 10 PetscInt *a2anew, k, k1, jmin, jmax, *jl, *il, vj, nexti, ili, i, j, mbs = a->mbs; in MatCholeskyFactorNumeric_SeqSBAIJ_7() local
23 PetscCall(PetscMalloc2(mbs, &il, mbs, &jl)); in MatCholeskyFactorNumeric_SeqSBAIJ_7()
25 for (i = 0; i < mbs; i++) jl[i] = mbs; in MatCholeskyFactorNumeric_SeqSBAIJ_7()
87 i = jl[k]; /* first row to be added to k_th row */ in MatCholeskyFactorNumeric_SeqSBAIJ_7()
90 nexti = jl[i]; /* next row to be added to k_th row */ in MatCholeskyFactorNumeric_SeqSBAIJ_7()
387 jl[i] = jl[j]; in MatCholeskyFactorNumeric_SeqSBAIJ_7()
388 jl[j] = i; /* update jl */ in MatCholeskyFactorNumeric_SeqSBAIJ_7()
417 jl[k] = jl[i]; in MatCholeskyFactorNumeric_SeqSBAIJ_7()
418 jl[i] = k; in MatCholeskyFactorNumeric_SeqSBAIJ_7()
423 PetscCall(PetscFree2(il, jl)); in MatCholeskyFactorNumeric_SeqSBAIJ_7()
|
| H A D | sbaijfact2.c | 2477 PetscInt jmin, jmax, nzk, k, j, *jl, prow, *il, nextprow; in MatICCFactorSymbolic_SeqSBAIJ() local
2524 PetscCall(PetscMalloc4(am, &uj_ptr, am, &uj_lvl_ptr, am, &il, am, &jl)); in MatICCFactorSymbolic_SeqSBAIJ()
2526 jl[i] = am; in MatICCFactorSymbolic_SeqSBAIJ()
2553 prow = jl[k]; /* 1st pivot row */ in MatICCFactorSymbolic_SeqSBAIJ()
2556 nextprow = jl[prow]; in MatICCFactorSymbolic_SeqSBAIJ()
2575 jl[prow] = jl[j]; in MatICCFactorSymbolic_SeqSBAIJ()
2576 jl[j] = prow; in MatICCFactorSymbolic_SeqSBAIJ()
2597 jl[k] = jl[i]; in MatICCFactorSymbolic_SeqSBAIJ()
2598 jl[i] = k; in MatICCFactorSymbolic_SeqSBAIJ()
2615 PetscCall(PetscFree4(uj_ptr, uj_lvl_ptr, il, jl)); in MatICCFactorSymbolic_SeqSBAIJ()
[all …]
|
| H A D | sbaij2.c | 1279 … i, j, k, bs = A->rmap->bs, bs2 = a->bs2, k1, mbs = a->mbs, jmin, jmax, nexti, ik, *jl, *il; in MatNorm_SeqSBAIJ() local
1309 PetscCall(PetscMalloc3(bs, &sum, mbs, &il, mbs, &jl)); in MatNorm_SeqSBAIJ()
1310 for (i = 0; i < mbs; i++) jl[i] = mbs; in MatNorm_SeqSBAIJ()
1317 i = jl[k]; /* first |A(i,k)| to be added */ in MatNorm_SeqSBAIJ()
1321 nexti = jl[i]; /* next block row to be added */ in MatNorm_SeqSBAIJ()
1336 jl[i] = jl[j]; in MatNorm_SeqSBAIJ()
1337 jl[j] = i; in MatNorm_SeqSBAIJ()
1359 jl[k] = jl[j]; in MatNorm_SeqSBAIJ()
1360 jl[j] = k; in MatNorm_SeqSBAIJ()
1366 PetscCall(PetscFree3(sum, il, jl)); in MatNorm_SeqSBAIJ()
|
| /petsc/src/mat/graphops/color/impls/minpack/ |
| H A D | numsrt.c | 11 PetscInt jinc, i, j, k, l, jl, ju; in MINPACKnumsrt() local
77 jl = 0; in MINPACKnumsrt()
81 jl = *nmax; in MINPACKnumsrt()
87 for (j = jl; i__2 < 0 ? j >= i__1 : j <= i__1; j += i__2) { in MINPACKnumsrt()
|
| /petsc/src/mat/impls/baij/seq/ |
| H A D | baijfact.c | 755 PetscInt k, jmin, jmax, *jl, *il, col, nexti, ili, nz; in MatCholeskyFactorNumeric_SeqBAIJ_N() local
776 PetscCall(PetscMalloc3(mbs, &rtmp, mbs, &il, mbs, &jl)); in MatCholeskyFactorNumeric_SeqBAIJ_N()
784 jl[i] = mbs; in MatCholeskyFactorNumeric_SeqBAIJ_N()
806 i = jl[k]; /* first row to be added to k_th row */ in MatCholeskyFactorNumeric_SeqBAIJ_N()
809 nexti = jl[i]; /* next row to be added to k_th row */ in MatCholeskyFactorNumeric_SeqBAIJ_N()
825 jl[i] = jl[j]; in MatCholeskyFactorNumeric_SeqBAIJ_N()
826 jl[j] = i; in MatCholeskyFactorNumeric_SeqBAIJ_N()
863 jl[k] = jl[i]; in MatCholeskyFactorNumeric_SeqBAIJ_N()
864 jl[i] = k; in MatCholeskyFactorNumeric_SeqBAIJ_N()
868 PetscCall(PetscFree3(rtmp, il, jl)); in MatCholeskyFactorNumeric_SeqBAIJ_N()
[all …]
|
| /petsc/src/mat/impls/aij/seq/ |
| H A D | aijfact.c | 1879 PetscInt k, jmin, jmax, *jl, *il, col, nexti, ili, nz; in MatCholeskyFactorNumeric_SeqAIJ_inplace() local
1915 PetscCall(PetscMalloc3(mbs, &rtmp, mbs, &il, mbs, &jl)); in MatCholeskyFactorNumeric_SeqAIJ_inplace()
1920 for (i = 0; i < mbs; i++) jl[i] = mbs; in MatCholeskyFactorNumeric_SeqAIJ_inplace()
1945 i = jl[k]; /* first row to be added to k_th row */ in MatCholeskyFactorNumeric_SeqAIJ_inplace()
1948 nexti = jl[i]; /* next row to be added to k_th row */ in MatCholeskyFactorNumeric_SeqAIJ_inplace()
1964 jl[i] = jl[j]; in MatCholeskyFactorNumeric_SeqAIJ_inplace()
1965 jl[j] = i; in MatCholeskyFactorNumeric_SeqAIJ_inplace()
1996 jl[k] = jl[i]; in MatCholeskyFactorNumeric_SeqAIJ_inplace()
1997 jl[i] = k; in MatCholeskyFactorNumeric_SeqAIJ_inplace()
2002 PetscCall(PetscFree3(rtmp, il, jl)); in MatCholeskyFactorNumeric_SeqAIJ_inplace()
[all …]
|
| /petsc/src/dm/impls/da/ |
| H A D | dageometry.c | 28 …const PetscInt il = s.i - da->Xs / da->w, jl = dim > 1 ? s.j - da->Ys : 0, kl = dim > 2 ? s.k - da… in DMDAConvertToCell() local
35 *cell = (kl * my + jl) * mx + il; in DMDAConvertToCell()
|
| /petsc/src/dm/dt/dualspace/impls/lagrange/ |
| H A D | dspacelagrange.c | 1048 PetscInt iT, iF, jT, jF, il, jl; in MatTensorAltV() local
1121 …for (jl = 0; jl < NkFiber; jl++) val += projFstar[il * NkFiber + jl] * PetscRealPart(valsF[jF * Nk… in MatTensorAltV()
1136 …for (jl = 0; jl < NkTrace; jl++) val += projTstar[il * NkTrace + jl] * PetscRealPart(valsT[jT * Nk… in MatTensorAltV()
1146 for (jl = 0; jl < dT; jl++) val += sign * wedgeMat[il * dT + jl] * workT[jl]; in MatTensorAltV()
|
| /petsc/doc/ |
| H A D | index.md | 86 …ub.io/Gridap.jl/stable/) a Julia-based finite element simulation package through the [GridapPETSc]…
|
| /petsc/doc/faq/ |
| H A D | index.md | 748 `PETScBinaryIO.jl` Julia package.
|