Home
last modified time | relevance | path

Searched refs:l (Results 1 – 25 of 307) sorted by relevance

12345678910>>...13

/petsc/src/ksp/ksp/impls/cg/pipelcg/
H A Dpipelcg.c4 #define offset(j) PetscMax(((j) - (2 * l)), 0)
6 #define G(i, j) (plcg->G[((j) * (2 * l + 1)) + (shift((i), (j)))])
7 #define G_noshift(i, j) (plcg->G[((j) * (2 * l + 1)) + (i)])
16 PetscInt l; /* pipeline depth */ member
39 PetscInt l = plcg->l, max_it = ksp->max_it; in KSPSetUp_PIPELCG() local
45 …PetscCheck(l >= 1, comm, PETSC_ERR_ARG_OUTOFRANGE, "%s: pipel argument must be positive.", ((Petsc… in KSPSetUp_PIPELCG()
46 …PetscCheck(l <= max_it, comm, PETSC_ERR_ARG_OUTOFRANGE, "%s: pipel argument must be less than max_… in KSPSetUp_PIPELCG()
51 PetscCall(VecDuplicateVecs(plcg->p, PetscMax(3, l + 1), &plcg->Z)); in KSPSetUp_PIPELCG()
54 PetscCall(VecDuplicateVecs(plcg->p, 3 * (l - 1) + 1, &plcg->Q)); in KSPSetUp_PIPELCG()
56 PetscCall(PetscCalloc1(l, &plcg->sigma)); in KSPSetUp_PIPELCG()
[all …]
/petsc/src/ksp/pc/impls/mg/
H A Dmgfunc.c143 PetscErrorCode PCMGSetResidual(PC pc, PetscInt l, PetscErrorCode (*residual)(Mat, Vec, Vec, Vec), M… in PCMGSetResidual() argument
151 if (residual) mglevels[l]->residual = residual; in PCMGSetResidual()
152 if (!mglevels[l]->residual) mglevels[l]->residual = PCMGResidualDefault; in PCMGSetResidual()
153 mglevels[l]->matresidual = PCMGMatResidualDefault; in PCMGSetResidual()
155 PetscCall(MatDestroy(&mglevels[l]->A)); in PCMGSetResidual()
156 mglevels[l]->A = mat; in PCMGSetResidual()
177 PetscErrorCode PCMGSetResidualTranspose(PC pc, PetscInt l, PetscErrorCode (*residualt)(Mat, Vec, Ve… in PCMGSetResidualTranspose() argument
185 if (residualt) mglevels[l]->residualtranspose = residualt; in PCMGSetResidualTranspose()
186 …if (!mglevels[l]->residualtranspose) mglevels[l]->residualtranspose = PCMGResidualTransposeDefault; in PCMGSetResidualTranspose()
187 mglevels[l]->matresidualtranspose = PCMGMatResidualTransposeDefault; in PCMGSetResidualTranspose()
[all …]
H A Dfmg.c8 PetscInt i, l = mglevels[0]->levels; in PCMGFCycle_Private() local
13 for (i = l - 1; i > 0; i--) { in PCMGFCycle_Private()
30 for (i = 0; i < l - 1; i++) { in PCMGFCycle_Private()
37 PetscCall(PCMGMCycle_Private(pc, &mglevels[l - 1], transpose, matapp, NULL)); in PCMGFCycle_Private()
39 PetscCall(PCMGMCycle_Private(pc, &mglevels[l - 1], transpose, matapp, NULL)); in PCMGFCycle_Private()
40 for (i = l - 2; i >= 0; i--) { in PCMGFCycle_Private()
47 for (i = 1; i < l; i++) { in PCMGFCycle_Private()
59 PetscInt i, l = mglevels[0]->levels; in PCMGKCycle_Private() local
63 for (i = l - 1; i > 0; i--) { in PCMGKCycle_Private()
80 for (i = 0; i < l - 1; i++) { in PCMGKCycle_Private()
[all …]
/petsc/src/dm/dt/tests/
H A Dex7.c8 PetscInt Nk, Mk, i, j, l; in CheckPullback() local
40 for (l = 0; l < k; l++) { in CheckPullback()
44 for (j = 0; j < N; j++) sum += L[i * N + j] * x[l * N + j]; in CheckPullback()
45 Lx[l * M + i] = sum; in CheckPullback()
153 PetscInt l; in main() local
156 … for (l = 0; l < k; l++) PetscCall(PetscPrintf(PETSC_COMM_WORLD, " %" PetscInt_FMT, subset[l])); in main()
158 … for (l = k; l < N; l++) PetscCall(PetscPrintf(PETSC_COMM_WORLD, " %" PetscInt_FMT, subset[l])); in main()
190 PetscInt l; in main() local
194 for (l = 0; l < N; l++) wvcheck += w[l] * v[l]; in main()
229 PetscInt l; in main() local
[all …]
H A Dex12.c45 PetscInt i, j, k, l; in TestSetIndex() local
59 for (l = 0; l < 4; ++l) { in TestSetIndex()
60 if ((l == i) || (l == j) || (l == k)) continue; in TestSetIndex()
64 …ll(PetscWeakFormSetIndexResidual(wf, key.label, key.value, key.field, key.part, l, f[l], 0, NULL)); in TestSetIndex()
77 for (l = 0; l < 4; ++l) { in TestSetIndex()
78 if ((l == i) || (l == j) || (l == k)) continue; in TestSetIndex()
82 …ll(PetscWeakFormSetIndexResidual(wf, key.label, key.value, key.field, key.part, 0, NULL, l, f[l])); in TestSetIndex()
95 for (l = 0; l < 4; ++l) { in TestSetIndex()
96 if ((l == i) || (l == j) || (l == k)) continue; in TestSetIndex()
100 …ll(PetscWeakFormSetIndexResidual(wf, key.label, key.value, key.field, key.part, l, f[l], l, f[l])); in TestSetIndex()
[all …]
H A Dex4.c67 PetscInt k, l; in CheckSymmetry() local
75 …for (l = 0; l < dim; l++) valsCopy[kLocal * dim + l] = vals[(offset + k) * dim + l] * (flip ? flip… in CheckSymmetry()
97 …for (l = 0; l < dim; l++) valsCopy2[(offset + k) * dim + l] = valsCopy[kLocal * dim + l] * (flip ?… in CheckSymmetry()
101 for (l = 0; l < dim; l++) { in CheckSymmetry()
102 …l] == vals[dim * k + l], PETSC_COMM_SELF, PETSC_ERR_PLIB, "Symmetry failure: %" PetscInt_FMT "D, %… in CheckSymmetry()
/petsc/src/ts/tutorials/
H A Dex27.c164 PetscInt i, j, l, Mx, My, xs, ys, xm, ym; in FormInitialGuess() local
174 for (l = 0; l < N_SPECIES; l++) { in FormInitialGuess()
176 if (l == 0) x[j][i].sp[l] = (ctx->x_inflow.sp[l] * ((PetscScalar)j) / (My - 1)); in FormInitialGuess()
177 … else if (l == 1) x[j][i].sp[l] = (ctx->x_inflow.sp[l] * (1. - ((PetscScalar)j) / (My - 1))); in FormInitialGuess()
178 else x[j][i].sp[l] = ctx->x_0.sp[l]; in FormInitialGuess()
189 PetscInt i, j, l, m; in FormIFunctionLocal() local
216 for (l = 0; l < N_SPECIES; l++) { in FormIFunctionLocal()
219 sxp = (x[j][i + 1].sp[l] - x[j][i].sp[l]) * dhx; in FormIFunctionLocal()
222 sxp = (x[j][i].sp[l] - x[j][i - 1].sp[l]) * dhx; in FormIFunctionLocal()
225 sxm = (x[j][i + 1].sp[l] - x[j][i].sp[l]) * dhx; in FormIFunctionLocal()
[all …]
/petsc/src/mat/tests/
H A Dex167.c28 …PetscInt n, N, i, j, k, l, nsub, Jlow[2] = {0, 1}, *jlow, Jhigh[2] = {3, 4}, *jhigh, row, c… in main() local
98 for (l = 0; l < nsub; ++l) { in main()
102 for (j = jlow[l]; j < jhigh[l]; ++j) { in main()
107 PetscCall(ISCreateGeneral(PETSC_COMM_SELF, 12, subindices, PETSC_OWN_POINTER, rowis + l)); in main()
109 PetscCall(ISDuplicate(rowis[l], colis + l)); in main()
111 PetscCall(PetscObjectReference((PetscObject)rowis[l])); in main()
112 colis[l] = rowis[l]; in main()
114 if (sort_rows) PetscCall(ISSort(rowis[l])); in main()
115 if (sort_cols) PetscCall(ISSort(colis[l])); in main()
128 for (l = 0; l < nsub; ++l) { in main()
[all …]
/petsc/src/ts/tutorials/power_grid/
H A Dex7.m9 l = size(AA); l = l(2); variable
10 A = cell(1,l-80-1);
11 for i=1:l-80-1; A{1,i} = AA{i+79}; end;
12 l = size(A); l = l(2); variable
14 for i=1:l; A{i} = reshape(A{i},n,n); end
17 for i=1:l
22 for i=1:l
27 a = zeros(l-1,1);
28 for i=1:l-1; a(i) = max(max(A{i+1}))/max(max(A{i})); end
29 b = zeros(l-1,1);
[all …]
/petsc/src/mat/impls/baij/seq/
H A Ddgefa2.c16 PetscInt i__2, i__3, kp1, j, k, l, ll, i, ipvt[2], k3; in PetscKernel_A_gets_inverse_A_2() local
37 l = 1; in PetscKernel_A_gets_inverse_A_2()
42 l = ll + 1; in PetscKernel_A_gets_inverse_A_2()
45 l += k - 1; in PetscKernel_A_gets_inverse_A_2()
46 ipvt[k - 1] = l; in PetscKernel_A_gets_inverse_A_2()
48 if (a[l + k3] == 0.0) { in PetscKernel_A_gets_inverse_A_2()
54 a[l + k3] = shift; in PetscKernel_A_gets_inverse_A_2()
59 if (l != k) { in PetscKernel_A_gets_inverse_A_2()
60 stmp = a[l + k3]; in PetscKernel_A_gets_inverse_A_2()
61 a[l + k3] = a[k4]; in PetscKernel_A_gets_inverse_A_2()
[all …]
H A Ddgefa.c17 PetscInt i__2, i__3, kp1, nm1, j, k, l, ll, kn, knp1, jn1; in PetscLINPACKgefa() local
40 l = 1; in PetscLINPACKgefa()
45 l = ll + 1; in PetscLINPACKgefa()
48 l += k - 1; in PetscLINPACKgefa()
49 ipvt[k] = l; in PetscLINPACKgefa()
51 if (a[l + kn] == 0.0) { in PetscLINPACKgefa()
58 if (l != k) { in PetscLINPACKgefa()
59 t = a[l + kn]; in PetscLINPACKgefa()
60 a[l + kn] = a[knp1]; in PetscLINPACKgefa()
74 t = a[l + jn1]; in PetscLINPACKgefa()
[all …]
H A Ddgefa5.c16 PetscInt i__2, i__3, kp1, j, k, l, ll, i, kb, k3; in PetscKernel_A_gets_inverse_A_5() local
37 l = 1; in PetscKernel_A_gets_inverse_A_5()
42 l = ll + 1; in PetscKernel_A_gets_inverse_A_5()
45 l += k - 1; in PetscKernel_A_gets_inverse_A_5()
46 ipvt[k - 1] = l; in PetscKernel_A_gets_inverse_A_5()
48 if (a[l + k3] == 0.0) { in PetscKernel_A_gets_inverse_A_5()
55 a[l + k3] = shift; in PetscKernel_A_gets_inverse_A_5()
60 if (l != k) { in PetscKernel_A_gets_inverse_A_5()
61 stmp = a[l + k3]; in PetscKernel_A_gets_inverse_A_5()
62 a[l + k3] = a[k4]; in PetscKernel_A_gets_inverse_A_5()
[all …]
H A Ddgefa3.c16 PetscInt i__2, i__3, kp1, j, k, l, ll, i, ipvt[3], kb, k3; in PetscKernel_A_gets_inverse_A_3() local
37 l = 1; in PetscKernel_A_gets_inverse_A_3()
42 l = ll + 1; in PetscKernel_A_gets_inverse_A_3()
45 l += k - 1; in PetscKernel_A_gets_inverse_A_3()
46 ipvt[k - 1] = l; in PetscKernel_A_gets_inverse_A_3()
48 if (a[l + k3] == 0.0) { in PetscKernel_A_gets_inverse_A_3()
55 a[l + k3] = shift; in PetscKernel_A_gets_inverse_A_3()
60 if (l != k) { in PetscKernel_A_gets_inverse_A_3()
61 stmp = a[l + k3]; in PetscKernel_A_gets_inverse_A_3()
62 a[l + k3] = a[k4]; in PetscKernel_A_gets_inverse_A_3()
[all …]
H A Ddgefa4.c16 PetscInt i__2, i__3, kp1, j, k, l, ll, i, ipvt[4], kb, k3; in PetscKernel_A_gets_inverse_A_4() local
37 l = 1; in PetscKernel_A_gets_inverse_A_4()
42 l = ll + 1; in PetscKernel_A_gets_inverse_A_4()
45 l += k - 1; in PetscKernel_A_gets_inverse_A_4()
46 ipvt[k - 1] = l; in PetscKernel_A_gets_inverse_A_4()
48 if (a[l + k3] == 0.0) { in PetscKernel_A_gets_inverse_A_4()
55 a[l + k3] = shift; in PetscKernel_A_gets_inverse_A_4()
60 if (l != k) { in PetscKernel_A_gets_inverse_A_4()
61 stmp = a[l + k3]; in PetscKernel_A_gets_inverse_A_4()
62 a[l + k3] = a[k4]; in PetscKernel_A_gets_inverse_A_4()
[all …]
H A Ddgefa6.c16 PetscInt i__2, i__3, kp1, j, k, l, ll, i, ipvt[6], kb, k3; in PetscKernel_A_gets_inverse_A_6() local
37 l = 1; in PetscKernel_A_gets_inverse_A_6()
42 l = ll + 1; in PetscKernel_A_gets_inverse_A_6()
45 l += k - 1; in PetscKernel_A_gets_inverse_A_6()
46 ipvt[k - 1] = l; in PetscKernel_A_gets_inverse_A_6()
48 if (a[l + k3] == 0.0) { in PetscKernel_A_gets_inverse_A_6()
55 a[l + k3] = shift; in PetscKernel_A_gets_inverse_A_6()
60 if (l != k) { in PetscKernel_A_gets_inverse_A_6()
61 stmp = a[l + k3]; in PetscKernel_A_gets_inverse_A_6()
62 a[l + k3] = a[k4]; in PetscKernel_A_gets_inverse_A_6()
[all …]
H A Ddgefa7.c16 PetscInt i__2, i__3, kp1, j, k, l, ll, i, ipvt[7], kb, k3; in PetscKernel_A_gets_inverse_A_7() local
37 l = 1; in PetscKernel_A_gets_inverse_A_7()
42 l = ll + 1; in PetscKernel_A_gets_inverse_A_7()
45 l += k - 1; in PetscKernel_A_gets_inverse_A_7()
46 ipvt[k - 1] = l; in PetscKernel_A_gets_inverse_A_7()
48 if (a[l + k3] == 0.0) { in PetscKernel_A_gets_inverse_A_7()
55 a[l + k3] = shift; in PetscKernel_A_gets_inverse_A_7()
60 if (l != k) { in PetscKernel_A_gets_inverse_A_7()
61 stmp = a[l + k3]; in PetscKernel_A_gets_inverse_A_7()
62 a[l + k3] = a[k4]; in PetscKernel_A_gets_inverse_A_7()
[all …]
/petsc/config/
H A Dreport_tests.py46 l = line.split()
47 if l[0] == 'failures':
48 if len(l)>1:
49 summary[l[0]] += l[1:]
50 elif l[0] == 'time':
51 if len(l)==1: continue
52 summary[l[0]] += float(l[1])
53 summary['cputime'] += float(l[2])
54 timesummary[cfile]=float(l[1])
55 cputimesummary[cfile]=float(l[2])
[all …]
/petsc/src/vec/vec/impls/mpi/
H A Dcommonmpvec.c17 static PetscErrorCode VecGhostStateSync_Private(Vec g, Vec l) in VecGhostStateSync_Private() argument
23 PetscCall(PetscObjectStateGet((PetscObject)l, &lstate)); in VecGhostStateSync_Private()
25 PetscCall(PetscObjectStateSet((PetscObject)l, PetscMax(gstate, lstate))); in VecGhostStateSync_Private()
66 PetscErrorCode VecGhostGetLocalForm(Vec g, Vec *l) in VecGhostGetLocalForm() argument
72 PetscAssertPointer(l, 2); in VecGhostGetLocalForm()
78 *l = v->localrep; in VecGhostGetLocalForm()
80 *l = g; in VecGhostGetLocalForm()
82 *l = NULL; in VecGhostGetLocalForm()
84 if (*l) PetscCall(VecGhostStateSync_Private(g, *l)); in VecGhostGetLocalForm()
104 PetscErrorCode VecGhostIsLocalForm(Vec g, Vec l, PetscBool *flg) in VecGhostIsLocalForm() argument
[all …]
/petsc/share/petsc/matlab/
H A DPetscBinaryRead.m48 for l=1:nargin-2
49 if ischar(varargin{l}) && strcmpi(varargin{l},'indices')
50 tnargin = min(l,tnargin-1);
51 indices = varargin{l+1};
53 if ischar(varargin{l}) && strcmpi(varargin{l},'precision')
54 tnargin = min(l,tnargin-1);
55 precision = varargin{l+1};
57 if ischar(varargin{l}) && strcmpi(varargin{l},'cell')
58 tnargin = min(l,tnargin-1);
59 arecell = varargin{l+1};
[all …]
H A DPetscBinaryWrite.m33 for l=1:nargin-2
34 if ischar(varargin{l}) && strcmpi(varargin{l},'indices')
35 tnargin = min(l,tnargin-1);
36 indices = varargin{l+1};
38 if ischar(varargin{l}) && strcmpi(varargin{l},'precision')
39 tnargin = min(l,tnargin-1);
40 precision = varargin{l+1};
42 if ischar(varargin{l}) && strcmpi(varargin{l},'complex')
43 tnargin = min(l,tnargin-1);
44 ispetsccomplex = varargin{l+1};
[all …]
/petsc/src/dm/partitioner/impls/multistage/
H A Dmspart.c69 for (PetscInt l = 0; l < p->levels; l++) { in PetscPartitionerView_Multistage() local
70 PetscPartitioner ppart = p->ppart[l]; in PetscPartitionerView_Multistage()
75 if (l) { in PetscPartitionerView_Multistage()
81 if (l) { in PetscPartitionerView_Multistage()
96 …l == p->stage) PetscCall(PetscViewerASCIIPrintf(viewer, "Stage %" PetscInt_FMT " partitioners (%" … in PetscPartitionerView_Multistage()
100 …if (!p->stagedm || l == p->stage) PetscCall(PetscViewerASCIIPrintf(viewer, "Stage %" PetscInt_FMT … in PetscPartitionerView_Multistage()
103 if (format == PETSC_VIEWER_ASCII_INFO_DETAIL && (!p->stagedm || l == p->stage)) { in PetscPartitionerView_Multistage()
124 if (!l) { in PetscPartitionerView_Multistage()
131 if (p->view_tpwgs) PetscCall(PetscSectionView(p->tpwgs[l], pviewer)); in PetscPartitionerView_Multistage()
137 if (l && pgroup != MPI_GROUP_EMPTY) PetscCallMPI(MPI_Group_free(&pgroup)); in PetscPartitionerView_Multistage()
[all …]
/petsc/src/mat/graphops/order/
H A Drcm.c43 PetscInt node, fnbr, lnbr, i, j, k, l, lperm, jstop, jstrt; in SPARSEPACKrcm() local
94 l = k; in SPARSEPACKrcm()
98 if (l < fnbr) goto L500; in SPARSEPACKrcm()
99 lperm = perm[l]; in SPARSEPACKrcm()
101 perm[l + 1] = lperm; in SPARSEPACKrcm()
102 --l; in SPARSEPACKrcm()
105 perm[l + 1] = nbr; in SPARSEPACKrcm()
114 l = *ccsize; in SPARSEPACKrcm()
117 lperm = perm[l]; in SPARSEPACKrcm()
118 perm[l] = perm[i]; in SPARSEPACKrcm()
[all …]
/petsc/src/dm/dt/space/impls/subspace/
H A Dspacesubspace.c96 PetscInt origDim, subDim, origNc, subNc, subNb, origNb, i, j, k, l, m, n, o; in PetscSpaceEvaluate_Subspace() local
139 …for (l = 0; l < origNc * origDim * origDim; l++) phi[l] += inH[(j * origNb + k) * origNc * origDim… in PetscSpaceEvaluate_Subspace()
143 for (l = 0; l < subDim; l++) { in PetscSpaceEvaluate_Subspace()
146 …for (o = 0; o < origDim; o++) psi[(k * subDim + l) * subDim + n] += subsp->Jx[m * subDim + l] * su… in PetscSpaceEvaluate_Subspace()
153 for (l = 0; l < PetscMin(subDim, origDim); l++) { in PetscSpaceEvaluate_Subspace()
154 …; m < PetscMin(subDim, origDim); m++) psi[(k * subDim + l) * subDim + m] += phi[(k * origDim + l) … in PetscSpaceEvaluate_Subspace()
160 for (l = 0; l < origNc; l++) { in PetscSpaceEvaluate_Subspace()
161 …bNb + i) * subNc + k) * subDim * subDim + m] += subsp->Ju[k * origNc + l] * psi[l * subDim * subDi… in PetscSpaceEvaluate_Subspace()
166 …for (l = 0; l < subDim * subDim; l++) H[((j * subNb + i) * subNc + k) * subDim * subDim + l] += ps… in PetscSpaceEvaluate_Subspace()
188 …for (l = 0; l < origNc * origDim; l++) phi[l] += inD[(j * origNb + k) * origNc * origDim + l] * su… in PetscSpaceEvaluate_Subspace()
[all …]
/petsc/src/sys/utils/
H A Dfhost.c75 size_t l, ll; in PetscGetHostName() local
77 PetscCall(PetscStrlen(name, &l)); in PetscGetHostName()
78 if (l == nlen - 1) PetscFunctionReturn(PETSC_SUCCESS); in PetscGetHostName()
79 name[l++] = '.'; in PetscGetHostName()
80 name[l] = 0; in PetscGetHostName()
82 …PetscCheck(!getdomainname(name + l, (int)(nlen - l)), PETSC_COMM_SELF, PETSC_ERR_SYS, "getdomainna… in PetscGetHostName()
93 name[l - 1] = 0; in PetscGetHostName()
/petsc/src/snes/tutorials/
H A Dex48.c712 PetscInt l, ll; in PointwiseNonlinearity() local
719 for (l = 0; l < 8; l++) { in PointwiseNonlinearity()
720 *u += phi[l] * n[l].u; in PointwiseNonlinearity()
721 *v += phi[l] * n[l].v; in PointwiseNonlinearity()
723 du[ll] += dphi[l][ll] * n[l].u; in PointwiseNonlinearity()
724 dv[ll] += dphi[l][ll] * n[l].v; in PointwiseNonlinearity()
733 PetscInt l, ll; in PointwiseNonlinearity2D() local
740 for (l = 0; l < 4; l++) { in PointwiseNonlinearity2D()
741 *u += phi[l] * n[l].u; in PointwiseNonlinearity2D()
742 *v += phi[l] * n[l].v; in PointwiseNonlinearity2D()
[all …]

12345678910>>...13