Home
last modified time | relevance | path

Searched refs:mu (Results 1 – 25 of 111) sorted by relevance

12345

/petsc/src/tao/tutorials/output/
H A Dex4_nm_admm_1.out1 Iter 0 : ||x-z||: 1.88745, mu*||z-zold||: 0.
2 Iter 1 : ||x-z||: 1.43502, mu*||z-zold||: 0.
3 Iter 2 : ||x-z||: 1.20971, mu*||z-zold||: 0.0802631
4 Iter 3 : ||x-z||: 0.943163, mu*||z-zold||: 0.522052
5 Iter 4 : ||x-z||: 0.727873, mu*||z-zold||: 0.247286
6 Iter 5 : ||x-z||: 0.525312, mu*||z-zold||: 0.217664
7 Iter 6 : ||x-z||: 0.425896, mu*||z-zold||: 0.128301
8 Iter 7 : ||x-z||: 0.312735, mu*||z-zold||: 0.167137
9 Iter 8 : ||x-z||: 0.239321, mu*||z-zold||: 0.0779858
10 Iter 9 : ||x-z||: 0.182317, mu*||z-zold||: 0.058334
H A Dex4_hessian_admm_1.out1 Iter 0 : ||x-z||: 1.81923, mu*||z-zold||: 0.
2 Iter 1 : ||x-z||: 1.33753, mu*||z-zold||: 0.
3 Iter 2 : ||x-z||: 1.11121, mu*||z-zold||: 0.
4 Iter 3 : ||x-z||: 0.963057, mu*||z-zold||: 0.
5 Iter 4 : ||x-z||: 0.851287, mu*||z-zold||: 0.
6 Iter 5 : ||x-z||: 0.760807, mu*||z-zold||: 0.
7 Iter 6 : ||x-z||: 0.684795, mu*||z-zold||: 0.
8 Iter 7 : ||x-z||: 0.478585, mu*||z-zold||: 0.197965
9 Iter 8 : ||x-z||: 0.354416, mu*||z-zold||: 0.217479
10 Iter 9 : ||x-z||: 0.285869, mu*||z-zold||: 0.220793
H A Dex4_soft_threshold_admm_1.out1 Iter 0 : ||x-z||: 1.81923, mu*||z-zold||: 0.
2 Iter 1 : ||x-z||: 1.33753, mu*||z-zold||: 0.
3 Iter 2 : ||x-z||: 1.02082, mu*||z-zold||: 0.299966
4 Iter 3 : ||x-z||: 0.812309, mu*||z-zold||: 0.356776
5 Iter 4 : ||x-z||: 0.614515, mu*||z-zold||: 0.286523
6 Iter 5 : ||x-z||: 0.460327, mu*||z-zold||: 0.25143
7 Iter 6 : ||x-z||: 0.3454, mu*||z-zold||: 0.180397
8 Iter 7 : ||x-z||: 0.262059, mu*||z-zold||: 0.115619
9 Iter 8 : ||x-z||: 0.200806, mu*||z-zold||: 0.112218
10 Iter 9 : ||x-z||: 0.151443, mu*||z-zold||: 0.0793162
H A Dex4_lmvm_admm_1.out1 Iter 0 : ||x-z||: 1.81923, mu*||z-zold||: 0.
2 Iter 1 : ||x-z||: 1.33753, mu*||z-zold||: 0.
3 Iter 2 : ||x-z||: 1.11121, mu*||z-zold||: 0.
4 Iter 3 : ||x-z||: 0.963057, mu*||z-zold||: 0.
5 Iter 4 : ||x-z||: 0.851287, mu*||z-zold||: 0.
6 Iter 5 : ||x-z||: 0.760807, mu*||z-zold||: 0.
7 Iter 6 : ||x-z||: 0.684795, mu*||z-zold||: 0.
8 Iter 7 : ||x-z||: 2.11936, mu*||z-zold||: 2.50113
9 Iter 8 : ||x-z||: 0.623987, mu*||z-zold||: 1.14322
10 Iter 9 : ||x-z||: 0.15654, mu*||z-zold||: 0.295423
H A Dex4_lmvm_admm_2.out1 Iter 0 : ||x-z||: 0.909616, mu*||z-zold||: 0.909616
2 Iter 1 : ||x-z||: 0.454808, mu*||z-zold||: 0.454808
3 Iter 2 : ||x-z||: 0.227404, mu*||z-zold||: 0.227404
4 Iter 3 : ||x-z||: 0.113702, mu*||z-zold||: 0.113702
5 Iter 4 : ||x-z||: 0.056851, mu*||z-zold||: 0.056851
6 Iter 5 : ||x-z||: 0.0284255, mu*||z-zold||: 0.0284255
7 Iter 6 : ||x-z||: 0.0142128, mu*||z-zold||: 0.0142128
H A Dex4_hessian_admm_2.out1 Iter 0 : ||x-z||: 0.909616, mu*||z-zold||: 0.909616
2 Iter 1 : ||x-z||: 0.454808, mu*||z-zold||: 0.454808
3 Iter 2 : ||x-z||: 0.227404, mu*||z-zold||: 0.227404
4 Iter 3 : ||x-z||: 0.113702, mu*||z-zold||: 0.113702
5 Iter 4 : ||x-z||: 0.056851, mu*||z-zold||: 0.056851
6 Iter 5 : ||x-z||: 0.0284255, mu*||z-zold||: 0.0284255
7 Iter 6 : ||x-z||: 0.0142128, mu*||z-zold||: 0.0142128
H A Dex4_nm_admm_2.out1 Iter 0 : ||x-z||: 1.09045, mu*||z-zold||: 0.876782
2 Iter 1 : ||x-z||: 0.39159, mu*||z-zold||: 0.792021
3 Iter 2 : ||x-z||: 0.248954, mu*||z-zold||: 0.231064
4 Iter 3 : ||x-z||: 0.120882, mu*||z-zold||: 0.132062
5 Iter 4 : ||x-z||: 0.0629571, mu*||z-zold||: 0.061955
6 Iter 5 : ||x-z||: 0.0323657, mu*||z-zold||: 0.0319213
7 Iter 6 : ||x-z||: 0.0158288, mu*||z-zold||: 0.019393
/petsc/src/ts/tutorials/autodiff/
H A Dex16adj.cxx22 PetscReal mu; member
43 f[1] = user->mu * (1. - x[0] * x[0]) * x[1] - x[0]; in RHSFunctionPassive()
71 f_a[1] = user->mu * (1. - x_a[0] * x_a[0]) * x_a[1] - x_a[0]; in RHSFunctionActive()
103 mu_a <<= user->mu; /* Mark independence */ in RHSFunctionActiveP()
141 PetscCall(PetscAdolcComputeRHSJacobianP(3, A, x, &user->mu, user->adctx)); in RHSJacobianP()
179 Vec lambda[2], mu[2], r; in main() local
193 user.mu = 1; in main()
201 PetscCall(PetscOptionsGetReal(NULL, NULL, "-mu", &user.mu, NULL)); in main()
268 …TSC_COMM_WORLD, "mu %g, steps %" PetscInt_FMT ", ftime %g\n", (double)user.mu, steps, (double)ftim… in main()
286 PetscCall(MatCreateVecs(Jacp, &mu[0], NULL)); in main()
[all …]
H A Dex16adj_tl.cxx26 PetscReal mu; member
40 PetscErrorCode EvaluateResidual(const T *x, T mu, T *f) in EvaluateResidual() argument
44 f[1] = mu * (1. - x[0] * x[0]) * x[1] - x[0]; in EvaluateResidual()
60 PetscCall(EvaluateResidual(x, user->mu, f)); in RHSFunctionPassive()
83 mu_a.setValue(user->mu); in RHSJacobian()
136 mu_a.setValue(user->mu); in RHSJacobianP()
199 Vec lambda[2], mu[2]; in main() local
213 user.mu = 1; in main()
221 PetscCall(PetscOptionsGetReal(NULL, NULL, "-mu", &user.mu, NULL)); in main()
278 …TSC_COMM_WORLD, "mu %g, steps %" PetscInt_FMT ", ftime %g\n", (double)user.mu, steps, (double)ftim… in main()
[all …]
/petsc/src/ts/tutorials/eimex/
H A Dct_vdp_imex.c34 PetscReal mu; /*stiffness control coefficient: epsilon*/ member
57 user.mu = 1e0; in main()
58 PetscCall(PetscOptionsReal("-eps", "Stiffness controller", "", user.mu, &user.mu, NULL)); in main()
118 …-2. / 3. + 10. / 81. * (user.mu) - 292. / 2187. * (user.mu) * (user.mu) - 1814. / 19683. * (user.m… in main()
192 f[1] = xdot[1] - ((1. - x[0] * x[0]) * x[1] - x[0]) / user->mu; in IFunction()
202 PetscReal mu = user->mu; in IJacobian() local
211 J[1][0] = (2. * x[0] * x[1] + 1.) / mu; in IJacobian()
212 J[1][1] = a - (1. - x[0] * x[0]) / mu; in IJacobian()
/petsc/src/ts/tutorials/autodiff/output/
H A Dex16opt_ic_1.out1 mu 1., steps 0, ftime 0.5
2 mu 1., steps 7, ftime 0.5
3 mu 1.000000, steps 75, ftime 0.5
5 mu 1.000000, steps 68, ftime 0.5
7 mu 1.000000, steps 68, ftime 0.5
/petsc/src/ts/tutorials/
H A Dex20opt_p.c18 PetscReal mu; member
53 f[1] = user->mu * ((1. - u[0] * u[0]) * u[1] - u[0]); in RHSFunction()
62 PetscReal mu = user->mu; in RHSJacobian() local
71 J[1][0] = -mu * (2.0 * u[1] * u[0] + 1.); in RHSJacobian()
73 J[1][1] = mu * (1.0 - u[0] * u[0]); in RHSJacobian()
100 dJdU[1][0][0] = -2. * user->mu * u[1]; in RHSHessianProductUU()
101 dJdU[1][1][0] = -2. * user->mu * u[0]; in RHSHessianProductUU()
102 dJdU[1][0][1] = -2. * user->mu * u[0]; in RHSHessianProductUU()
192 f[1] = udot[1] - user->mu * ((1.0 - u[0] * u[0]) * u[1] - u[0]); in IFunction()
212 J[1][0] = user->mu * (1.0 + 2.0 * u[0] * u[1]); in IJacobian()
[all …]
H A Dex20adj.c55 PetscReal mu; member
83 f[1] = user->mu * ((1. - u[0] * u[0]) * u[1] - u[0]); in RHSFunction()
93 PetscReal mu = user->mu; in RHSJacobian() local
106 J[1][0] = -mu * (2.0 * u[1] * u[0] + 1.); in RHSJacobian()
107 J[1][1] = mu * (1.0 - u[0] * u[0]); in RHSJacobian()
137 f[1] = udot[1] - user->mu * ((1.0 - u[0] * u[0]) * u[1] - u[0]); in IFunction()
161 J[1][0] = user->mu * (2.0 * u[0] * u[1] + 1.0); in IJacobian()
162 J[1][1] = a - user->mu * (1.0 - u[0] * u[0]); in IJacobian()
257 user.mu = 1.0e3; in main()
262 PetscCall(PetscOptionsGetReal(NULL, NULL, "-mu", &user.mu, NULL)); in main()
[all …]
H A Dex16.c71 PetscReal mu; member
106 f[1] = xdot[1] - user->mu * ((1. - x[0] * x[0]) * x[1] - x[0]); in IFunction()
116 PetscReal mu = user->mu; in IJacobian() local
125 J[1][0] = mu * (2. * x[0] * x[1] + 1.); in IJacobian()
126 J[1][1] = a - mu * (1. - x[0] * x[0]); in IJacobian()
205 user.mu = 1000.0; in main()
209 PetscCall(PetscOptionsGetReal(NULL, NULL, "-mu", &user.mu, NULL)); in main()
239 x_ptr[1] = -2.0 / 3.0 + 10.0 / (81.0 * user.mu) - 292.0 / (2187.0 * user.mu * user.mu); in main()
254 …TSC_COMM_WORLD, "mu %g, steps %" PetscInt_FMT ", ftime %g\n", (double)user.mu, steps, (double)ftim… in main()
H A Dex20.c24 PetscReal mu; member
41 f[1] = user->mu * (1. - x[0] * x[0]) * x[1] - x[0]; in RHSFunction()
58 f[1] = xdot[1] - user->mu * ((1.0 - x[0] * x[0]) * x[1] - x[0]); in IFunction()
76 J[1][0] = user->mu * (2.0 * x[0] * x[1] + 1.0); in IJacobian()
77 J[1][1] = a - user->mu * (1.0 - x[0] * x[0]); in IJacobian()
138 user.mu = 1.0e3; in main()
142 PetscCall(PetscOptionsReal("-mu", "Stiffness parameter", "<1.0e6>", user.mu, &user.mu, NULL)); in main()
177 x_ptr[1] = -2.0 / 3.0 + 10.0 / (81.0 * user.mu) - 292.0 / (2187.0 * user.mu * user.mu); in main()
H A Dex49.c24 PetscReal mu; member
39 f[1] = xdot[1] - user->mu * ((1.0 - x[0] * x[0]) * x[1] - x[0]); in IFunction()
57 J[1][0] = user->mu * (1.0 + 2.0 * x[0] * x[1]); in IJacobian()
58 J[1][1] = a - user->mu * (1.0 - x[0] * x[0]); in IJacobian()
94 user.mu = 1.0e6; in main()
96 PetscCall(PetscOptionsReal("-mu", "Stiffness parameter", "<1.0e6>", user.mu, &user.mu, NULL)); in main()
H A Dex16fwd.c50 PetscReal mu; member
68 f[1] = user->mu * (1. - x[0] * x[0]) * x[1] - x[0]; in RHSFunction()
77 PetscReal mu = user->mu; in RHSJacobian() local
85 J[1][0] = -2. * mu * x[1] * x[0] - 1.; in RHSJacobian()
87 J[1][1] = mu * (1.0 - x[0] * x[0]); in RHSJacobian()
160 user.mu = 1; in main()
163 PetscCall(PetscOptionsGetReal(NULL, NULL, "-mu", &user.mu, NULL)); in main()
219 …TSC_COMM_WORLD, "mu %g, steps %" PetscInt_FMT ", ftime %g\n", (double)user.mu, steps, (double)ftim… in main()
/petsc/src/ts/tutorials/hybrid/
H A Dex1adj.c46 Vec *lambda, *mu; in ShiftGradients() local
53 PetscCall(TSGetCostGradients(ts, &numcost, &lambda, &mu)); in ShiftGradients()
79 PetscCall(VecGetArray(mu[0], &y)); in ShiftGradients()
85 PetscCall(VecRestoreArray(mu[0], &y)); in ShiftGradients()
89 PetscCall(VecGetArray(mu[1], &y)); in ShiftGradients()
95 PetscCall(VecRestoreArray(mu[1], &y)); in ShiftGradients()
204 Vec lambda[2], mu[2]; in main() local
302 PetscCall(MatCreateVecs(Ap, &mu[0], NULL)); in main()
303 PetscCall(MatCreateVecs(Ap, &mu[1], NULL)); in main()
304 PetscCall(VecZeroEntries(mu[0])); in main()
[all …]
/petsc/src/tao/tutorials/
H A Dex4.c21 PetscReal mu; /* the augmented Lagrangian term in ADMM */ member
124 ctx->mu = 1.0; in ConfigureContext()
138 …tionsReal("-mu", "The augmented lagrangian multiplier in ADMM", "ex4.c", ctx->mu, &ctx->mu, NULL)); in ConfigureContext()
222 PetscReal mu, workNorm, misfit; in ObjectiveMisfitADMM() local
226 mu = ctx->mu; in ObjectiveMisfitADMM()
238 *J = misfit + 0.5 * mu * workNorm; in ObjectiveMisfitADMM()
246 PetscReal mu; in GradientMisfitADMM() local
250 mu = ctx->mu; in GradientMisfitADMM()
259 PetscCall(VecAXPY(V, mu, temp)); in GradientMisfitADMM()
270 PetscCall(MatShift(H, ctx->mu)); in HessianMisfitADMM()
[all …]
/petsc/src/snes/tutorials/
H A Dex56.c36 PetscReal trace, mu = s_mu, lambda = s_lambda, rad; in f1_u_3d_alpha() local
44 mu *= s_soft_alpha; in f1_u_3d_alpha()
49 for (j = 0; j < dim; ++j) f1[i * dim + j] = mu * (u_x[i * dim + j] + u_x[j * dim + i]); in f1_u_3d_alpha()
57 PetscReal trace, mu = s_mu, lambda = s_lambda; in f1_u_3d() local
61 for (j = 0; j < dim; ++j) f1[i * dim + j] = mu * (u_x[i * dim + j] + u_x[j * dim + i]); in f1_u_3d()
75 void g3_uu_3d_private(PetscScalar g3[], const PetscReal mu, const PetscReal lambda) in g3_uu_3d_private() argument
79 g3[0] += mu; in g3_uu_3d_private()
80 g3[0] += mu; in g3_uu_3d_private()
83 g3[10] += mu; in g3_uu_3d_private()
84 g3[12] += mu; in g3_uu_3d_private()
[all …]
H A Dex62.c45 PetscScalar mu; /* dynamic shear viscosity */ member
55 const PetscReal mu = PetscRealPart(constants[0]); in f1_u() local
60 for (d = 0; d < dim; ++d) f1[c * dim + d] = mu * (u_x[c * dim + d] + u_x[d * dim + c]); in f1_u()
85 const PetscReal mu = PetscRealPart(constants[0]); in g3_uu() local
91 g3[((c * Nc + c) * dim + d) * dim + d] += mu; /* < \nabla v, \nabla u > */ in g3_uu()
92 g3[((c * Nc + d) * dim + d) * dim + c] += mu; /* < \nabla v, {\nabla u}^T > */ in g3_uu()
99 const PetscReal mu = PetscRealPart(constants[0]); in g0_pp() local
101 g0[0] = 1.0 / mu; in g0_pp()
158 const PetscReal mu = PetscRealPart(constants[0]); in f0_quadratic_u() local
161 f0[0] = (dim - 1) * 4.0 * mu - 1.0; in f0_quadratic_u()
[all …]
H A Dex17.c43 PetscScalar mu; /* shear modulus */ member
128 const PetscReal mu = PetscRealPart(constants[0]); in f0_elas_quadratic_u() local
131 for (PetscInt d = 0; d < dim - 1; ++d) f0[d] += 2.0 * mu; in f0_elas_quadratic_u()
132 f0[dim - 1] += 2.0 * lambda + 4.0 * mu; in f0_elas_quadratic_u()
201 const PetscReal mu = PetscRealPart(constants[0]); in f0_elas_trig_u() local
205 …0; d < dim; ++d) f0[d] += -(2.0 * mu + lambda) * fact * PetscSinReal(2.0 * PETSC_PI * x[d]) - (d <… in f0_elas_trig_u()
215 const PetscReal mu = PetscRealPart(param->mu); in axial_disp_u() local
219 … * lambda * lambda + 8. * lambda * mu + 4 * mu * mu) / (4 * mu * (3 * lambda * lambda + 5. * lambd… in axial_disp_u()
220 u[1] = 0.25 * lambda / mu / (lambda + mu) * N * x[1]; in axial_disp_u()
294 const PetscReal mu = PetscRealPart(constants[0]); in f1_elas_u() local
[all …]
/petsc/src/binding/petsc4py/demo/legacy/ode/
H A Dvanderpol.py39 mu = self.mu_
41 u[1] = -2.0 / 3.0 + 10.0 / (81.0 * mu) - 292.0 / (2187.0 * mu * mu)
45 mu = self.mu_
50 f[1] = mu * ((1.0 - u[0] * u[0]) * u[1] - u[0])
58 mu = self.mu_
65 J[1, 0] = -mu * (2.0 * u[1] * u[0] + 1.0)
66 J[1, 1] = mu * (1.0 - u[0] * u[0])
82 mu = self.mu_
87 f[1] = udot[1] - mu * ((1.0 - u[0] * u[0]) * u[1] - u[0])
95 mu = self.mu_
[all …]
/petsc/src/tao/constrained/impls/admm/
H A Dadmm.c46 am->dualres *= am->mu; in TaoADMMToleranceUpdate()
101 am->muold = am->mu; in AdaptiveADMMPenaltyUpdate()
103 am->mu = PetscSqrtReal(a_k * b_k); in AdaptiveADMMPenaltyUpdate()
105 am->mu = a_k; in AdaptiveADMMPenaltyUpdate()
107 am->mu = b_k; in AdaptiveADMMPenaltyUpdate()
109 if (am->mu > am->muold) am->mu = am->muold; in AdaptiveADMMPenaltyUpdate()
110 if (am->mu < am->mumin) am->mu = am->mumin; in AdaptiveADMMPenaltyUpdate()
189 *f += temp + (am->mu / 2) * temp2; in SubObjGradUpdate()
193 PetscCall(VecAXPY(g, am->mu, tempJR)); in SubObjGradUpdate()
217 *f += temp + (am->mu / 2) * temp2; in RegObjGradUpdate()
[all …]
/petsc/src/ts/tutorials/power_grid/
H A Dex3sa.c37 Vec lambda[1], q, mu[1]; in main() local
147 PetscCall(MatCreateVecs(ctx.Jacp, &mu[0], NULL)); in main()
148 PetscCall(TSSetCostGradients(ts, 1, lambda, mu)); in main()
209 PetscCall(VecGetArray(mu[0], &x_ptr)); in main()
211 PetscCall(VecRestoreArray(mu[0], &x_ptr)); in main()
218 PetscCall(VecView(mu[0], PETSC_VIEWER_STDOUT_WORLD)); in main()
223 PetscCall(ComputeSensiP(lambda[0], mu[0], &ctx)); in main()
224 PetscCall(VecGetArray(mu[0], &x_ptr)); in main()
226 PetscCall(VecRestoreArray(mu[0], &x_ptr)); in main()
228 PetscCall(VecDestroy(&mu[0])); in main()

12345