Home
last modified time | relevance | path

Searched refs:lambda (Results 1 – 25 of 291) sorted by relevance

12345678910>>...12

/petsc/src/ksp/ksp/impls/cg/pipecg2/
H A Dpipecg2.c82 … PetscScalar beta0, PetscScalar alpha0, PetscScalar beta1, PetscScalar alpha1, PetscScalar *lambda) in VecMergedOpsShort_Private() argument
112 for (j = 0; j < 15; j++) lambda[j] = 0.0; in VecMergedOpsShort_Private()
152 lambda[0] += ps[j] * PetscConj(pu[j]); in VecMergedOpsShort_Private()
153 lambda[1] += pw[j] * PetscConj(pm[j]); in VecMergedOpsShort_Private()
154 lambda[2] += pw[j] * PetscConj(pq[j]); in VecMergedOpsShort_Private()
155 lambda[4] += ps[j] * PetscConj(pq[j]); in VecMergedOpsShort_Private()
156 lambda[6] += pn[j] * PetscConj(pm[j]); in VecMergedOpsShort_Private()
157 lambda[7] += pn[j] * PetscConj(pq[j]); in VecMergedOpsShort_Private()
158 lambda[9] += pz[j] * PetscConj(pq[j]); in VecMergedOpsShort_Private()
159 lambda[10] += pr[j] * PetscConj(pu[j]); in VecMergedOpsShort_Private()
[all …]
/petsc/src/snes/tutorials/output/
H A Dex58_3.out5 0 Line search: fty/||y|| = -0.033171, lambda = 0.5
6 1 Line search: fty/||y|| = 0.0180532, lambda = 0.25
7 2 Line search: fty/||y|| = 0.0126471, lambda = 0.375
8 3 Line search: fty/||y|| = 0.000436031, lambda = 0.4375
9 4 Line search: fty/||y|| = -0.0166589, lambda = 0.46875
10 5 Line search: fty/||y|| = -0.00698054, lambda = 0.453125
11 6 Line search: fty/||y|| = -0.00295682, lambda = 0.445312
12 7 Line search: fty/||y|| = -0.00118288, lambda = 0.441406
13 8 Line search: fty/||y|| = -0.000354409, lambda = 0.439453
14 9 Line search: fty/||y|| = 4.55115e-05, lambda = 0.438477
[all …]
H A Dex15_nleqerr.out3 Line search: entering iteration with lambda: 1.000000000000e+00
6 Line search: entering iteration with lambda: 2.886019430766e-01
10 Line search: Lipschitz estimate: 4.610151815134e+00; lambda: 1.000000000000e+00
11 Line search: entering iteration with lambda: 1.000000000000e+00
14 Line search: entering iteration with lambda: 2.324329826076e-01
18 Line search: Lipschitz estimate: 1.367676990197e+00; lambda: 1.000000000000e+00
19 Line search: entering iteration with lambda: 1.000000000000e+00
22 Line search: entering iteration with lambda: 2.018780657684e-01
26 Line search: Lipschitz estimate: 1.273457981665e+00; lambda: 1.000000000000e+00
27 Line search: entering iteration with lambda: 1.000000000000e+00
[all …]
H A Dex58_4.out11 0 Line search: fty/||y|| = 0.0106239, lambda = 0.5
12 1 Line search: fty/||y|| = -0.0262441, lambda = 0.75
13 2 Line search: fty/||y|| = 0.000212591, lambda = 0.625
14 3 Line search: fty/||y|| = -0.0104791, lambda = 0.6875
15 4 Line search: fty/||y|| = -0.0044817, lambda = 0.65625
16 5 Line search: fty/||y|| = -0.0019843, lambda = 0.640625
17 6 Line search: fty/||y|| = -0.000850272, lambda = 0.632812
18 7 Line search: fty/||y|| = -0.000310204, lambda = 0.628906
19 8 Line search: fty/||y|| = -4.66799e-05, lambda = 0.626953
20 9 Line search: fty/||y|| = 8.3483e-05, lambda = 0.625977
[all …]
H A Dex58_5.out6 Line search terminated: lambda = 1.33342, fnorms = 0.0764847
14 Line search terminated: lambda = 1.44125, fnorms = 0.0814181
21 Line search terminated: lambda = 1.17396, fnorms = 0.0704441
28 Line search terminated: lambda = 1.34669, fnorms = 0.038434
35 Line search terminated: lambda = 1.50534, fnorms = 0.0328088
42 Line search terminated: lambda = 1.8127, fnorms = 0.0311457
49 Line search terminated: lambda = 1.79072, fnorms = 0.0283698
56 Line search terminated: lambda = 1.30027, fnorms = 0.00935502
63 Line search terminated: lambda = 1.06331, fnorms = 0.00070829
70 Line search terminated: lambda = 1.00323, fnorms = 2.29332e-06
[all …]
H A Dex58_6.out7 Line search terminated: lambda = 0.915104, fnorm = 0.081438
16 Line search terminated: lambda = 0.475511, fnorm = 0.0577764
29 Line search terminated: lambda = 0.669089, fnorm = 0.0360721
40 Line search terminated: lambda = 0.62328, fnorm = 0.0220052
51 Line search terminated: lambda = 0.633595, fnorm = 0.0158549
63 Line search terminated: lambda = 0.508606, fnorm = 0.0126286
74 Line search terminated: lambda = 0.349585, fnorm = 0.011117
85 Line search terminated: lambda = 0.427734, fnorm = 0.00934175
95 Line search terminated: lambda = 0.885085, fnorm = 0.00519364
102 Line search terminated: lambda = 1., fnorm = 0.000645946
[all …]
/petsc/src/snes/linesearch/impls/secant/
H A Dlinesearchsecant.c18 PetscReal lambda, lambda_old, lambda_mid, lambda_update, delLambda; in SNESLineSearchApply_Secant() local
27 PetscCall(SNESLineSearchGetLambda(linesearch, &lambda)); in SNESLineSearchApply_Secant()
44 lambda_mid = 0.5 * (lambda + lambda_old); in SNESLineSearchApply_Secant()
62 PetscCall(VecWAXPY(W, -lambda, Y, X)); in SNESLineSearchApply_Secant()
78 PetscCall(VecWAXPY(W, -lambda, Y, X)); in SNESLineSearchApply_Secant()
86 … search: objective function at lambda = %g is infinity or NaN, cutting lambda\n", (double)lambda)); in SNESLineSearchApply_Secant()
91 if (lambda <= minlambda) { in SNESLineSearchApply_Secant()
94 …ambda_min = %g is infinity or NaN, can not further cut lambda\n", (double)lambda, (double)lambda)); in SNESLineSearchApply_Secant()
102 maxlambda = .95 * lambda; in SNESLineSearchApply_Secant()
105 lambda = .5 * (lambda + lambda_old); in SNESLineSearchApply_Secant()
[all …]
/petsc/src/snes/linesearch/impls/bt/
H A Dlinesearchbt.c64 PetscReal lambda, lambdatemp, lambdaprev, minlambda, initslope, alpha, stol; in SNESLineSearchApply_BT() local
77 PetscCall(SNESLineSearchGetLambda(linesearch, &lambda)); in SNESLineSearchApply_BT()
131 PetscCall(VecWAXPY(W, -lambda, Y, X)); in SNESLineSearchApply_BT()
157 …search: objective function at lambdas = %g is infinity or NaN, cutting lambda\n", (double)lambda)); in SNESLineSearchApply_BT()
160 if (lambda <= minlambda) SNESCheckFunctionDomainError(snes, g); in SNESLineSearchApply_BT()
161 lambda *= .5; in SNESLineSearchApply_BT()
165 …if (g <= f + lambda * alpha * initslope) { /* Sufficient reduction or step tolerance convergence */ in SNESLineSearchApply_BT()
188 lambdaprev = lambda; in SNESLineSearchApply_BT()
192 lambdatemp = -initslope * PetscSqr(lambda) / (2.0 * (g - f - lambda * initslope)); in SNESLineSearchApply_BT()
193 lambda = PetscClipInterval(lambdatemp, .1 * lambda, .5 * lambda); in SNESLineSearchApply_BT()
[all …]
/petsc/src/snes/tutorials/
H A Dex21.c113 PetscScalar *u, *lambda, *w, *fu, *fw, *flambda, d, h; in FormFunction() local
127 PetscCall(DMDAVecGetArray(user->da1, vlambda, &lambda)); in FormFunction()
134 fw[0] = -2. * d * lambda[0]; in FormFunction()
139 if (i == 0) flambda[0] = h * u[0] + 2. * d * lambda[0] - d * lambda[1]; in FormFunction()
140 else if (i == 1) flambda[1] = 2. * h * u[1] + 2. * d * lambda[1] - d * lambda[2]; in FormFunction()
141 else if (i == N - 1) flambda[N - 1] = h * u[N - 1] + 2. * d * lambda[N - 1] - d * lambda[N - 2]; in FormFunction()
142 …else if (i == N - 2) flambda[N - 2] = 2. * h * u[N - 2] + 2. * d * lambda[N - 2] - d * lambda[N - … in FormFunction()
143 else flambda[i] = 2. * h * u[i] - d * (lambda[i + 1] - 2.0 * lambda[i] + lambda[i - 1]); in FormFunction()
157 PetscCall(DMDAVecRestoreArray(user->da1, vlambda, &lambda)); in FormFunction()
169 Vec w, u, lambda, U, F; in Monitor() local
[all …]
H A Dex22.c137 PetscScalar lambda; member
173 fw[0] = -2.0 * d * u_lambda[0].lambda; in ComputeFunction()
178 …if (i == 0) fu_lambda[0].lambda = h * u_lambda[0].u + 2. * d * u_lambda[0].lambda - d * u_lambda[1… in ComputeFunction()
179 …else if (i == 1) fu_lambda[1].lambda = 2. * h * u_lambda[1].u + 2. * d * u_lambda[1].lambda - d * … in ComputeFunction()
180 …(i == N - 1) fu_lambda[N - 1].lambda = h * u_lambda[N - 1].u + 2. * d * u_lambda[N - 1].lambda - d… in ComputeFunction()
181 …== N - 2) fu_lambda[N - 2].lambda = 2. * h * u_lambda[N - 2].u + 2. * d * u_lambda[N - 2].lambda -… in ComputeFunction()
182 …else fu_lambda[i].lambda = 2. * h * u_lambda[i].u - d * (u_lambda[i + 1].lambda - 2.0 * u_lambda[i… in ComputeFunction()
/petsc/src/snes/linesearch/impls/nleqerr/
H A Dlinesearchnleqerr.c39 PetscReal lambda, minlambda, stol; in SNESLineSearchApply_NLEQERR() local
51 PetscCall(SNESLineSearchGetLambda(linesearch, &lambda)); in SNESLineSearchApply_NLEQERR()
104 lambda = PetscMin(1.0, nleqerr->mu_curr); in SNESLineSearchApply_NLEQERR()
108 …earch: Lipschitz estimate: %14.12e; lambda: %14.12e\n", (double)nleqerr->mu_curr, (double)lambda)); in SNESLineSearchApply_NLEQERR()
112 lambda = linesearch->damping; in SNESLineSearchApply_NLEQERR()
122 …SCIIPrintf(monitor, " Line search: entering iteration with lambda: %14.12e\n", (double)lambda)); in SNESLineSearchApply_NLEQERR()
139 if (lambda <= minlambda) { in SNESLineSearchApply_NLEQERR()
146 lambda = 1.0; in SNESLineSearchApply_NLEQERR()
147 PetscCall(VecWAXPY(G, -lambda, Y, X)); in SNESLineSearchApply_NLEQERR()
164 PetscCall(VecWAXPY(W, -lambda, Y, X)); in SNESLineSearchApply_NLEQERR()
[all …]
/petsc/src/tao/unconstrained/impls/bmrm/
H A Dbmrm.c298 PetscReal lambda; in TaoSolve_BMRM() local
313 lambda = bmrm->lambda; in TaoSolve_BMRM()
372 df.Q[i][tao->niter] = df.Q[tao->niter][i] = reg / lambda; in TaoSolve_BMRM()
387 PetscCall(VecAXPY(bmrm->local_w, -df.x[i] / lambda, pgrad->V)); in TaoSolve_BMRM()
392 reg = 0.5 * lambda * reg * reg; in TaoSolve_BMRM()
460 …PetscCall(PetscOptionsReal("-tao_bmrm_lambda", "regulariser weight", "", 100, &bmrm->lambda, NULL)… in TaoSetFromOptions_BMRM()
499 bmrm->lambda = 1.0; in TaoCreate_BMRM()
671 static PetscReal phi(PetscReal *x, PetscInt n, PetscReal lambda, PetscReal *a, PetscReal b, PetscRe… in phi() argument
677 x[i] = -c[i] + lambda * a[i]; in phi()
695 PetscReal lambda, lambdal, lambdau, dlambda, lambda_new; in project() local
[all …]
/petsc/src/snes/linesearch/impls/cp/
H A Dlinesearchcp.c10 PetscReal lambda, lambda_old, lambda_update, delLambda; in SNESLineSearchApply_CP() local
19 PetscCall(SNESLineSearchGetLambda(linesearch, &lambda)); in SNESLineSearchApply_CP()
47 PetscCall(VecWAXPY(W, -lambda, Y, X)); in SNESLineSearchApply_CP()
57 delLambda = lambda - lambda_old; in SNESLineSearchApply_CP()
92 …nitor, " Line search: lambdas = [%g, %g], ftys = [%g, %g]\n", (double)lambda, (double)lambda_ol… in SNESLineSearchApply_CP()
102 PetscCall(VecWAXPY(W, -0.5 * (lambda + lambda_old), Y, X)); in SNESLineSearchApply_CP()
116 PetscCall(VecWAXPY(W, -0.5 * (lambda + lambda_old), Y, X)); in SNESLineSearchApply_CP()
126 PetscCall(VecWAXPY(W, -(lambda + 0.5 * (lambda - lambda_old)), Y, X)); in SNESLineSearchApply_CP()
149 lambda_update = lambda - PetscRealPart(fty / s); in SNESLineSearchApply_CP()
152 if (lambda_update < minlambda) lambda_update = lambda + PetscRealPart(fty / s); in SNESLineSearchApply_CP()
[all …]
/petsc/src/mat/tutorials/
H A Dex20f.F909 PetscReal :: lambda
24 ctxM%lambda = ctxF_pt%lambda
53 ctxF%lambda = 3.14d0
56 print *, 'ctxF%lambda = ', ctxF%lambda
59 print *, 'ctxF_pt%lambda = ', ctxF_pt%lambda
63 print *, 'ctxFcopy_pt%lambda = ', ctxFcopy_pt%lambda
H A Dex6f.F9011 PetscReal :: lambda, kappa
31 ctxF%lambda = 3.14d0
34 print *, 'ctxF%lambda = ', ctxF%lambda
37 print *, 'ctxF_pt%lambda = ', ctxF_pt%lambda
/petsc/doc/manual/
H A Ddt.md20 pdf(x) = \lambda^d P(\lambda x).
27 \int_\Omega P(x) dx &= \int_\Omega P(\lambda s) \lambda^d ds \\
28 &= \int_\Omega P(s) \lambda^{-d} \lambda^d ds \\
44 C(\lambda x) &= \int^{\lambda x}_{x_-} P(s) ds \\
45 &= \int^{x}_{x_-} \lambda^d P(\lambda t) dt \\
/petsc/src/ts/tests/
H A Dex21.c28 PetscReal lambda; member
37 PetscReal lambda, hx, hy; in FormIFunctionLocal() local
41 lambda = app->lambda; in FormIFunctionLocal()
64 f[j][i] = ut - uxx - uyy - lambda * PetscExpScalar(u); in FormIFunctionLocal()
78 PetscScalar v[5], lambda, hx, hy; in FormIJacobianLocal() local
81 lambda = app->lambda; in FormIJacobianLocal()
108 v[k] = shift + 2.0 / (hx * hx) + 2.0 / (hy * hy) - lambda * PetscExpScalar(x[j][i]); in FormIJacobianLocal()
154 app.lambda = 6.8; in main()
155 app.lambda = 6.0; in main()
156 PetscCall(PetscOptionsReal("-lambda", "", "", app.lambda, &app.lambda, NULL)); in main()
/petsc/src/snes/linesearch/impls/bisection/
H A Dlinesearchbisection.c10 PetscReal lambda_left, lambda, lambda_right, lambda_old, fnorm; in SNESLineSearchApply_Bisection() local
18 PetscCall(SNESLineSearchGetLambda(linesearch, &lambda)); in SNESLineSearchApply_Bisection()
32 lambda_right = lambda; in SNESLineSearchApply_Bisection()
43 PetscCall(VecWAXPY(W, -lambda, Y, X)); in SNESLineSearchApply_Bisection()
119 if (PetscAbsReal(lambda - lambda_old) < ltol) { in SNESLineSearchApply_Bisection()
122 … Line search: abs(dlambda) = %g < ltol = %g\n", (double)PetscAbsReal(lambda - lambda_old), (do… in SNESLineSearchApply_Bisection()
131 lambda_right = lambda; in SNESLineSearchApply_Bisection()
133 lambda_left = lambda; in SNESLineSearchApply_Bisection()
140 lambda_old = lambda; in SNESLineSearchApply_Bisection()
141 lambda = 0.5 * (lambda_left + lambda_right); in SNESLineSearchApply_Bisection()
[all …]
/petsc/src/ts/tutorials/hybrid/
H A Dex1adj.c46 Vec *lambda, *mu; in ShiftGradients() local
53 PetscCall(TSGetCostGradients(ts, &numcost, &lambda, &mu)); in ShiftGradients()
78 PetscCall(VecGetArray(lambda[0], &x)); in ShiftGradients()
86 PetscCall(VecRestoreArray(lambda[0], &x)); in ShiftGradients()
88 PetscCall(VecGetArray(lambda[1], &x)); in ShiftGradients()
96 PetscCall(VecRestoreArray(lambda[1], &x)); in ShiftGradients()
204 Vec lambda[2], mu[2]; in main() local
290 PetscCall(MatCreateVecs(A, &lambda[0], NULL)); in main()
291 PetscCall(MatCreateVecs(A, &lambda[1], NULL)); in main()
293 PetscCall(VecZeroEntries(lambda[0])); in main()
[all …]
/petsc/src/snes/tests/
H A Dex1f.F90121 PetscReal lambda
124 common/params/lambda, mx, my, fd_coloring
140 temp1 = lambda/(lambda + one)
185 PetscReal lambda
188 common/params/lambda, mx, my, fd_coloring
243 PetscReal lambda
246 common/params/lambda, mx, my, fd_coloring
263 sc = hx*hy*lambda
316 PetscReal lambda
319 common/params/lambda, mx, my, fd_coloring
[all …]
/petsc/src/binding/petsc4py/demo/legacy/wrap-f2py/
H A DBratu2D.F9056 subroutine InitGuessLocal(grd, x, lambda, ierr) argument
60 PetscReal lambda
69 temp1 = lambda/(lambda + one)
87 subroutine FunctionLocal(grd, x, f, lambda, ierr) argument
92 PetscReal lambda
117 f(i, j) = uxx + uyy - lambda*exp(u)*sc
125 subroutine JacobianLocal(grd, x, Jac, lambda, ierr) argument
130 PetscReal lambda
160 v(3) = two*(hydhx + hxdhy) - lambda*exp(x(i, j))*sc
184 PetscReal lambda
[all …]
/petsc/src/ksp/pc/impls/kaczmarz/
H A Dkaczmarz.c4 PetscReal lambda; /* damping parameter */ member
24 PetscReal lambda = jac->lambda; in PCApply_Kaczmarz() local
43 if (cols[j] >= ys && cols[j] < ye) yarray[cols[j] - ys] += vals[j] * lambda * r / anrm; in PCApply_Kaczmarz()
59 if (cols[j] >= ys && cols[j] < ye) yarray[cols[j] - ys] += vals[j] * lambda * r / anrm; in PCApply_Kaczmarz()
76 …eal("-pc_kaczmarz_lambda", "relaxation factor (0 < lambda)", "", jac->lambda, &jac->lambda, NULL)); in PCSetFromOptions_Kaczmarz()
89 if (isascii) PetscCall(PetscViewerASCIIPrintf(viewer, " lambda = %g\n", (double)jac->lambda)); in PCView_Kaczmarz()
121 jac->lambda = 1.0; in PCCreate_Kaczmarz()
/petsc/src/tao/constrained/tutorials/
H A DtomographyADMM.c19 PetscReal lambda, eps, mumin; member
36 PetscReal lambda, mu; in TaoShellSolve_SoftThreshold() local
48 PetscCall(TaoADMMGetRegularizerCoefficient(admm_tao, &lambda)); in TaoShellSolve_SoftThreshold()
60 PetscCall(TaoSoftThreshold(work, -lambda / mu, lambda / mu, out)); in TaoShellSolve_SoftThreshold()
88 PetscReal lambda; in RegularizerObjectiveAndGradient1() local
102 PetscCall(TaoADMMGetRegularizerCoefficient(admm_tao, &lambda)); in RegularizerObjectiveAndGradient1()
103 *f_reg *= lambda; in RegularizerObjectiveAndGradient1()
106 PetscCall(VecScale(G_reg, lambda)); in RegularizerObjectiveAndGradient1()
114 PetscReal temp, lambda; in RegularizerObjectiveAndGradient2() local
121 PetscCall(TaoADMMGetRegularizerCoefficient(admm_tao, &lambda)); in RegularizerObjectiveAndGradient2()
[all …]
/petsc/src/sys/classes/random/tests/
H A Dex3.c11 static PetscErrorCode PoissonTailProbability(PetscReal lambda, PetscInt Y, PetscReal *prob) in PoissonTailProbability() argument
19 logLambda = PetscLogReal(lambda); in PoissonTailProbability()
21 PetscReal exponent = -lambda + i * logLambda - logFact; in PoissonTailProbability()
37 PetscReal lambda, p; in main() local
58 lambda = PetscPowRealInt(2., 3 * log2n - (2 + log2d * t)); in main()
61 …"Expected spacing collisions per process %g (%g total).\n", (double)lambda, (double)(N * lambda))); in main()
91 PetscCall(PoissonTailProbability(N * lambda, Y, &p)); in main()
/petsc/src/mat/tutorials/output/
H A Dex20f_1.out1 ctxF%lambda = 3.1400000000000001
2 ctxF_pt%lambda = 3.1400000000000001
3 ctxFcopy_pt%lambda = 3.1400000000000001

12345678910>>...12