Home
last modified time | relevance | path

Searched refs:kappa (Results 1 – 25 of 36) sorted by relevance

12

/petsc/src/ksp/ksp/impls/cg/stcg/
H A Dstcg.c19 PetscReal alpha, beta, kappa, rz, rzm1; in KSPCGSolve_STCG()
188 PetscCall(VecDot(p, z, &kappa)); /* kappa = p^T Q p */ in KSPCGSolve_STCG()
189 if (PetscIsInfOrNanScalar(kappa)) { in KSPCGSolve_STCG()
197 PetscCall(PetscInfo(ksp, "KSPCGSolve_STCG: bad matrix: kappa=%g\n", (double)kappa)); in KSPCGSolve_STCG()
243 if (kappa <= 0.0) { in KSPCGSolve_STCG()
251 PetscCall(PetscInfo(ksp, "KSPCGSolve_STCG: negative curvature: kappa=%g\n", (double)kappa)); in KSPCGSolve_STCG()
268 cg->o_fcn += step * (0.5 * step * kappa - rz); in KSPCGSolve_STCG()
311 alpha = rz / kappa; in KSPCGSolve_STCG()
343 cg->o_fcn += step * (0.5 * step * kappa - rz); in KSPCGSolve_STCG()
478 PetscCall(VecDot(p, z, &kappa)); /* kappa = p^T Q p */ in KSPCGSolve_STCG()
[all …]
/petsc/src/ksp/ksp/impls/cg/nash/
H A Dnash.c20 PetscReal alpha, beta, kappa, rz, rzm1; in KSPCGSolve_NASH()
191 PetscCall(VecDot(p, z, &kappa)); /* kappa = p^T Q p */ in KSPCGSolve_NASH()
192 if (PetscIsInfOrNanScalar(kappa)) { in KSPCGSolve_NASH()
200 PetscCall(PetscInfo(ksp, "KSPCGSolve_NASH: bad matrix: kappa=%g\n", (double)kappa)); in KSPCGSolve_NASH()
246 if (kappa <= 0.0) { in KSPCGSolve_NASH()
254 PetscCall(PetscInfo(ksp, "KSPCGSolve_NASH: negative curvature: kappa=%g\n", (double)kappa)); in KSPCGSolve_NASH()
271 cg->o_fcn += step * (0.5 * step * kappa - rz); in KSPCGSolve_NASH()
313 alpha = rz / kappa; in KSPCGSolve_NASH()
344 cg->o_fcn += step * (0.5 * step * kappa - rz); in KSPCGSolve_NASH()
479 PetscCall(VecDot(p, z, &kappa)); /* kappa = p^T Q p */ in KSPCGSolve_NASH()
[all …]
/petsc/src/ksp/ksp/tutorials/
H A Dex60.c77 const PetscScalar kappa = 5.0; in diagFunc1() local
78 return 1.0 + (kappa * (PetscScalar)i) / (PetscScalar)(n - 1); in diagFunc1()
83 const PetscScalar kappa = 50.0; in diagFunc2() local
84 return 1.0 + (kappa * (PetscScalar)i) / (PetscScalar)(n - 1); in diagFunc2()
89 const PetscScalar kappa = 10.0; in diagFunc3() local
93 return 1. + (kappa * ((PetscScalar)(i - 1))) / (PetscScalar)(n - 2); in diagFunc3()
/petsc/src/ksp/ksp/impls/cg/gltr/
H A Dgltr.c70 PetscReal alpha, beta, kappa, rz, rzm1; in KSPCGSolve_GLTR()
230 PetscCall(VecDot(p, z, &kappa)); /* kappa = p^T Q p */ in KSPCGSolve_GLTR()
231 if (PetscIsInfOrNanScalar(kappa)) { in KSPCGSolve_GLTR()
238 PetscCall(PetscInfo(ksp, "KSPCGSolve_GLTR: bad matrix: kappa=%g\n", (double)kappa)); in KSPCGSolve_GLTR()
279 cg->diag[t_size] = kappa / rz; in KSPCGSolve_GLTR()
289 if (PetscAbsReal(kappa) <= 0.0) { in KSPCGSolve_GLTR()
293 PetscCall(PetscInfo(ksp, "KSPCGSolve_GLTR: breakdown: kappa=%g\n", (double)kappa)); in KSPCGSolve_GLTR()
304 cg->o_fcn += step * (0.5 * step * kappa - rz); in KSPCGSolve_GLTR()
337 alpha = rz / kappa; in KSPCGSolve_GLTR()
349 PetscCall(PetscInfo(ksp, "KSPCGSolve_GLTR: negative curvature: kappa=%g\n", (double)kappa)); in KSPCGSolve_GLTR()
[all …]
/petsc/src/ts/tutorials/phasefield/
H A Dbiharmonic2.c36 PetscReal kappa; member
62 ctx.kappa = 1.0; in main()
63 PetscCall(PetscOptionsGetReal(NULL, NULL, "-kappa", &ctx.kappa, NULL)); in main()
88 dt = 1.0 / (10. * ctx.kappa * Mx * Mx * Mx * Mx); in main()
138 PetscCall(FormInitialSolution(da, x, ctx.kappa)); in main()
228 f[i].w = x[i].w + ctx->kappa * (x[i - 1].u + x[i + 1].u - 2.0 * x[i].u) * sx; in FormFunction()
259 PetscErrorCode FormInitialSolution(DM da, Vec X, PetscReal kappa) in FormInitialSolution() argument
295 for (i = xs; i < xs + xm; i++) x[i].w = -kappa * (x[i - 1].u + x[i + 1].u - 2.0 * x[i].u) * sx; in FormInitialSolution()
H A Dbiharmonic3.c36 PetscReal kappa; member
62 ctx.kappa = 1.0; in main()
63 PetscCall(PetscOptionsGetReal(NULL, NULL, "-kappa", &ctx.kappa, NULL)); in main()
86 dt = 1.0 / (10. * ctx.kappa * Mx * Mx * Mx * Mx); in main()
137 PetscCall(FormInitialSolution(da, x, ctx.kappa)); in main()
228 f[i].w = x[i].w + ctx->kappa * (x[i - 1].u + x[i + 1].u - 2.0 * x[i].u) * sx; in FormFunction()
269 PetscErrorCode FormInitialSolution(DM da, Vec X, PetscReal kappa) in FormInitialSolution() argument
305 for (i = xs; i < xs + xm; i++) x[i].w = -kappa * (x[i - 1].u + x[i + 1].u - 2.0 * x[i].u) * sx; in FormInitialSolution()
H A Dheat.c37 PetscReal kappa; member
59 ctx.kappa = 1.0; in main()
60 PetscCall(PetscOptionsGetReal(NULL, NULL, "-kappa", &ctx.kappa, NULL)); in main()
73 dt = 1.0 / (ctx.kappa * Mx * Mx); in main()
191 f[i] = ctx->kappa * (x[i - 1] + x[i + 1] - 2.0 * x[i]) * sx; in FormFunction()
327 …yy[0] = PetscRealPart(.25 * ctx->kappa * (u[i - 1] - u[i + 1]) * (u[i - 1] - u[i + 1]) * s… in MyMonitor()
347 yy[0] = PetscRealPart(ctx->kappa * (u[i - 1] + u[i + 1] - 2.0 * u[i]) * sx); in MyMonitor()
387 len = .5 * PetscRealPart(ctx->kappa * (u[i - 1] + u[i + 1] - 2.0 * u[i]) * sx) / max; in MyMonitor()
H A Dbiharmonictest.m7 …rch_minlambda 1.e-30 -ts_type beuler -da_refine ' int2str(i+4) ' -kappa .00001 -ts_time_st…
8 …vi -snes_converged_reason -ts_type cn -da_refine ' int2str(i+4) ' -kappa .00001 -ts_time_st…
H A Dbiharmonic.c56 PetscReal kappa; member
81 ctx.kappa = 1.0; in main()
82 PetscCall(PetscOptionsGetReal(NULL, NULL, "-kappa", &ctx.kappa, NULL)); in main()
109 dt = 1.0 / (10. * ctx.kappa * Mx * Mx * Mx * Mx); in main()
248 f[i] = -ctx->kappa * (l + r - 2.0 * c) * sx; in FormFunction()
352 vals[0] = -ctx->kappa * sx * sx; in FormJacobian()
354 vals[1] = 4.0 * ctx->kappa * sx * sx; in FormJacobian()
356 vals[2] = -6.0 * ctx->kappa * sx * sx; in FormJacobian()
358 vals[3] = 4.0 * ctx->kappa * sx * sx; in FormJacobian()
360 vals[4] = -ctx->kappa * sx * sx; in FormJacobian()
[all …]
H A Dheattest.m6 …system(['./heat -allen-cahn -kappa .001 -ts_monitor -snes_monitor -pc_type lu -snes_converged_…
/petsc/src/snes/tutorials/
H A Dex64.c13 PetscScalar kappa; member
115 const PetscReal kappa = PetscRealPart(constants[3]); in w_1() local
116 for (PetscInt d = 0; d < dim; ++d) f[d] = -kappa * u_x[uOff_x[2] + d]; in w_1()
121 const PetscReal kappa = PetscRealPart(constants[3]); in Jw_1_w1w1() local
123 for (PetscInt d = 0; d < dim; ++d) J[d * dim + d] = -kappa; in Jw_1_w1w1()
128 const PetscReal kappa = PetscRealPart(constants[3]); in pJw_1_w1w1() local
130 for (PetscInt d = 0; d < dim; ++d) J[d * dim + d] = kappa; in pJw_1_w1w1()
137 PetscScalar kappa; member
146 user->kappa = 1.0; in ProcessOptions()
151 PetscCall(PetscOptionsScalar("-kappa", "kappa", NULL, user->kappa, &user->kappa, NULL)); in ProcessOptions()
[all …]
H A Dex15.c74 PetscReal kappa; member
135 user.kappa = 1e-3; in main()
162 …PetscCall(PetscOptionsReal("-kappa", "diffusivity in odd regions", "", user.kappa, &user.kappa, NU… in main()
391 static inline PetscReal kappa(const AppCtx *ctx, PetscReal x, PetscReal y) in kappa() function
393 …return (((PetscInt)(x * ctx->blocks[0])) + ((PetscInt)(y * ctx->blocks[1]))) % 2 ? ctx->kappa : 1.… in kappa()
398 …return kappa(ctx, x, y) * PetscPowScalar(PetscSqr(ctx->epsilon) + 0.5 * (ux * ux + uy * uy), 0.5 *… in eta()
402 …return (ctx->p == 2) ? 0 : kappa(ctx, x, y) * PetscPowScalar(PetscSqr(ctx->epsilon) + 0.5 * (ux * … in deta()
/petsc/src/tao/unconstrained/impls/ntr/
H A Dntr.c50 PetscReal fmin, ftrial, prered, actred, kappa, sigma, beta; in TaoSolve_NTR() local
140 kappa = 1.0; in TaoSolve_NTR()
142 kappa = actred / prered; in TaoSolve_NTR()
150 if (PetscAbsScalar(kappa - (PetscReal)1.0) <= tr->mu1_i) { in TaoSolve_NTR()
161 } else if (PetscAbsScalar(kappa - (PetscReal)1.0) <= tr->mu2_i) { in TaoSolve_NTR()
308 kappa = 1.0; in TaoSolve_NTR()
310 kappa = actred / prered; in TaoSolve_NTR()
314 if (kappa < tr->eta1) { in TaoSolve_NTR()
319 if (kappa < tr->eta2) { in TaoSolve_NTR()
322 } else if (kappa < tr->eta3) { in TaoSolve_NTR()
[all …]
/petsc/src/tao/unconstrained/impls/ntl/
H A Dntl.c37 PetscReal fmin, ftrial, prered, actred, kappa, sigma; in TaoSolve_NTL() local
134 kappa = 1.0; in TaoSolve_NTL()
136 kappa = actred / prered; in TaoSolve_NTL()
144 if (PetscAbsScalar(kappa - (PetscReal)1.0) <= tl->mu1_i) { in TaoSolve_NTL()
159 } else if (PetscAbsScalar(kappa - (PetscReal)1.0) <= tl->mu2_i) { in TaoSolve_NTL()
313 kappa = 1.0; in TaoSolve_NTL()
315 kappa = actred / prered; in TaoSolve_NTL()
319 if (kappa < tl->eta1) { in TaoSolve_NTL()
325 if (kappa < tl->eta2) { in TaoSolve_NTL()
328 } else if (kappa < tl->eta3) { in TaoSolve_NTL()
[all …]
/petsc/src/tao/unconstrained/impls/nls/
H A Dnls.c47 PetscReal fmin, ftrial, f_full, prered, actred, kappa, sigma; in TaoSolve_NLS() local
160 kappa = 1.0; in TaoSolve_NLS()
162 kappa = actred / prered; in TaoSolve_NLS()
170 if (PetscAbsScalar(kappa - (PetscReal)1.0) <= nlsP->mu1_i) { in TaoSolve_NLS()
185 } else if (PetscAbsScalar(kappa - (PetscReal)1.0) <= nlsP->mu2_i) { in TaoSolve_NLS()
577 kappa = 1.0; in TaoSolve_NLS()
579 kappa = actred / prered; in TaoSolve_NLS()
583 if (kappa < nlsP->eta1) { in TaoSolve_NLS()
586 } else if (kappa < nlsP->eta2) { in TaoSolve_NLS()
589 } else if (kappa < nlsP->eta3) { in TaoSolve_NLS()
[all …]
/petsc/src/ts/tutorials/
H A Dex53.c360 const PetscReal kappa = PetscRealPart(constants[4]); in f0_trig_linear_p() local
367 f0[0] -= sum / M - 4 * PetscSqr(PETSC_PI) * kappa * sum * t; in f0_trig_linear_p()
447 PetscScalar kappa = param->k / param->mu_f; /* m^2 / (Pa s) */ in terzaghi_2d_u() local
455 …PetscScalar c = kappa / S; /* m^2 / s, Cheng (B.16) */ in terzaghi_2d_u()
484 PetscScalar kappa = param->k / param->mu_f; /* m^2 / (Pa s) */ in terzaghi_2d_eps() local
492 …PetscScalar c = kappa / S; /* m^2 / s, Cheng (B.16) */ in terzaghi_2d_eps()
521 PetscScalar kappa = param->k / param->mu_f; /* m^2 / (Pa s) */ in terzaghi_2d_p() local
528 … PetscScalar c = kappa / S; /* m^2 / s, Cheng (B.16) */ in terzaghi_2d_p()
559 PetscScalar kappa = param->k / param->mu_f; /* m^2 / (Pa s) */ in terzaghi_2d_u_t() local
567 …PetscScalar c = kappa / S; /* m^2 / s, Cheng (B.16) */ in terzaghi_2d_u_t()
[all …]
/petsc/src/ts/tutorials/advection-diffusion-reaction/
H A Dex4.c33 PetscScalar epsilon, delta, alpha, beta, gamma, kappa, lambda, mu, cstar; member
60 appctx.kappa = .75; in main()
198 kcxrhox = appctx->kappa * (cxx * rho + cx * rhox); in IFunction()
200 …kcxrhox = appctx->kappa * ((u[i + 1].c - u[i].c) * u[i + 1].rho - (u[i].c - u[i - 1].c) * u[i].rho… in IFunction()
H A Dreaction_diffusion.h12 PetscReal D1, D2, gamma, kappa; member
/petsc/src/snes/tutorials/output/
H A Dex64_2d_p2_p1_p2.out3 alpha = 0.5, kappa = 1.
H A Dex64_2d_p2_p1_p2_fetidp.out3 alpha = 0.5, kappa = 1.
/petsc/src/mat/tutorials/
H A Dex6f.F9011 PetscReal :: lambda, kappa
/petsc/src/tao/bound/impls/bnk/
H A Dbnk.c43 PetscReal f_min, ftrial, prered, actred, kappa, sigma, resnorm; in TaoBNKInitialize() local
182 kappa = 1.0; in TaoBNKInitialize()
184 kappa = actred / prered; in TaoBNKInitialize()
192 if (PetscAbsScalar(kappa - (PetscReal)1.0) <= bnk->mu1_i) { in TaoBNKInitialize()
203 } else if (PetscAbsScalar(kappa - (PetscReal)1.0) <= bnk->mu2_i) { in TaoBNKInitialize()
796 PetscReal step, kappa; in TaoBNKUpdateTrustRadius() local
845 kappa = 1.0; in TaoBNKUpdateTrustRadius()
847 kappa = actred / prered; in TaoBNKUpdateTrustRadius()
850 if (kappa < bnk->eta1) { in TaoBNKUpdateTrustRadius()
858 if (kappa < bnk->eta2) { in TaoBNKUpdateTrustRadius()
[all …]
/petsc/src/ts/tests/
H A Dex24.c13 PetscReal D1, D2, gamma, kappa; member
118 entries[4] = -2.0 * appctx->D2 * (sx + sy) + 2.0 * uc * vc - appctx->gamma - appctx->kappa; in RHSJacobian()
198 appctx.kappa = .06; in main()
/petsc/src/ts/adapt/impls/dsp/
H A Dadaptdsp.c38 static PetscReal Limiter(PetscReal value, PetscReal kappa) in Limiter() argument
40 return 1 + kappa * PetscAtanReal((value - 1) / kappa); in Limiter()
/petsc/src/ts/tutorials/autodiff/
H A Dadr_ex5adj_mf.cxx31 PetscReal D1, D2, gamma, kappa; member
63 appctx.kappa = .06; in main()
277 ….v = udot[j][i].v - appctx->D2 * (vxx + vyy) - uc * vc * vc + (appctx->gamma + appctx->kappa) * vc; in IFunctionLocalPassive()
373 ….v = udot[j][i].v - appctx->D2 * (vxx + vyy) - uc * vc * vc + (appctx->gamma + appctx->kappa) * vc; in IFunctionActive()

12