Home
last modified time | relevance | path

Searched refs:eta (Results 1 – 25 of 50) sorted by relevance

12

/petsc/src/ksp/ksp/tutorials/
H A Dex60.c33 PetscReal eta; member
50 PetscCall(VecScale(xout, ctx->eta * (nrmin / nrmnoise))); in PCApply_Noise()
122 PetscReal eta = 0.1, norm; in main() local
129 PetscCall(PetscOptionsGetReal(NULL, NULL, "-eta", &eta, NULL)); in main()
171 if (eta > 0) { in main()
174 ctx.eta = eta; in main()
H A Dex42.c27 PetscScalar eta[GAUSS_POINTS]; member
105 PetscReal eta = PetscRealPart(_xi[1]); in ShapeFunctionQ13D_Evaluate() local
108 Ni[0] = 0.125 * (1.0 - xi) * (1.0 - eta) * (1.0 - zeta); in ShapeFunctionQ13D_Evaluate()
109 Ni[1] = 0.125 * (1.0 - xi) * (1.0 + eta) * (1.0 - zeta); in ShapeFunctionQ13D_Evaluate()
110 Ni[2] = 0.125 * (1.0 + xi) * (1.0 + eta) * (1.0 - zeta); in ShapeFunctionQ13D_Evaluate()
111 Ni[3] = 0.125 * (1.0 + xi) * (1.0 - eta) * (1.0 - zeta); in ShapeFunctionQ13D_Evaluate()
113 Ni[4] = 0.125 * (1.0 - xi) * (1.0 - eta) * (1.0 + zeta); in ShapeFunctionQ13D_Evaluate()
114 Ni[5] = 0.125 * (1.0 - xi) * (1.0 + eta) * (1.0 + zeta); in ShapeFunctionQ13D_Evaluate()
115 Ni[6] = 0.125 * (1.0 + xi) * (1.0 + eta) * (1.0 + zeta); in ShapeFunctionQ13D_Evaluate()
116 Ni[7] = 0.125 * (1.0 + xi) * (1.0 - eta) * (1.0 + zeta); in ShapeFunctionQ13D_Evaluate()
[all …]
H A Dex43.c72 PetscScalar eta[GAUSS_POINTS]; member
101 for (p = 0; p < GAUSS_POINTS; p++) array[cum++] = props[ej][ei].eta[p]; in glvis_extract_eta()
119 PetscScalar eta = _xi[1]; in ConstructQ12D_Ni() local
121 Ni[0] = 0.25 * (1.0 - xi) * (1.0 - eta); in ConstructQ12D_Ni()
122 Ni[1] = 0.25 * (1.0 - xi) * (1.0 + eta); in ConstructQ12D_Ni()
123 Ni[2] = 0.25 * (1.0 + xi) * (1.0 + eta); in ConstructQ12D_Ni()
124 Ni[3] = 0.25 * (1.0 + xi) * (1.0 - eta); in ConstructQ12D_Ni()
130 PetscScalar eta = _xi[1]; in ConstructQ12D_GNi() local
132 GNi[0][0] = -0.25 * (1.0 - eta); in ConstructQ12D_GNi()
133 GNi[0][1] = -0.25 * (1.0 + eta); in ConstructQ12D_GNi()
[all …]
H A Dex21.c31 PetscReal norm, *coords, eta, scale = 0.5; in main() local
47 eta = 0.6; in main()
48 PetscCall(PetscOptionsGetReal(NULL, NULL, "-eta", &eta, NULL)); in main()
80 …, n, n, PETSC_DECIDE, PETSC_DECIDE, sdim, coords, PETSC_FALSE, RBF, &fctx, eta, leafsize, basisord… in main()
H A Dex70.c47 PetscScalar eta = _xi[1]; in EvaluateBasis_Q1() local
49 N[0] = 0.25 * (1.0 - xi) * (1.0 - eta); in EvaluateBasis_Q1()
50 N[1] = 0.25 * (1.0 + xi) * (1.0 - eta); in EvaluateBasis_Q1()
51 N[2] = 0.25 * (1.0 + xi) * (1.0 + eta); in EvaluateBasis_Q1()
52 N[3] = 0.25 * (1.0 - xi) * (1.0 + eta); in EvaluateBasis_Q1()
58 PetscScalar eta = _xi[1]; in EvaluateBasisDerivatives_Q1() local
60 dN[0][0] = -0.25 * (1.0 - eta); in EvaluateBasisDerivatives_Q1()
61 dN[0][1] = 0.25 * (1.0 - eta); in EvaluateBasisDerivatives_Q1()
62 dN[0][2] = 0.25 * (1.0 + eta); in EvaluateBasisDerivatives_Q1()
63 dN[0][3] = -0.25 * (1.0 + eta); in EvaluateBasisDerivatives_Q1()
[all …]
/petsc/src/dm/tests/
H A Dex36.c88 PetscScalar eta = XX[dim * i + 1]; in DAApplyConformalMapping() local
95 Ni[0] = 0.125 * (1.0 - xi) * (1.0 - eta) * (1.0 - zeta); in DAApplyConformalMapping()
96 Ni[1] = 0.125 * (1.0 + xi) * (1.0 - eta) * (1.0 - zeta); in DAApplyConformalMapping()
97 Ni[2] = 0.125 * (1.0 - xi) * (1.0 + eta) * (1.0 - zeta); in DAApplyConformalMapping()
98 Ni[3] = 0.125 * (1.0 + xi) * (1.0 + eta) * (1.0 - zeta); in DAApplyConformalMapping()
100 Ni[4] = 0.125 * (1.0 - xi) * (1.0 - eta) * (1.0 + zeta); in DAApplyConformalMapping()
101 Ni[5] = 0.125 * (1.0 + xi) * (1.0 - eta) * (1.0 + zeta); in DAApplyConformalMapping()
102 Ni[6] = 0.125 * (1.0 - xi) * (1.0 + eta) * (1.0 + zeta); in DAApplyConformalMapping()
103 Ni[7] = 0.125 * (1.0 + xi) * (1.0 + eta) * (1.0 + zeta); in DAApplyConformalMapping()
172 PetscScalar eta = XX[dim * i + 1]; in DAApplyConformalMapping() local
[all …]
/petsc/src/snes/tutorials/
H A Dex48.c260 PRange eta; member
374 static void THIViscosity(THI thi, PetscReal gam, PetscReal *eta, PetscReal *deta) in THIViscosity() argument
390 *eta = Bd2 * PetscPowReal(eps + gam, exponent); in THIViscosity()
391 *deta = exponent * (*eta) / (eps + gam); in THIViscosity()
547 PRangeClear(&thi->eta); in THICreate()
551 … * units->meter / (3e7 * units->second), gradu = u / (100 * units->meter), eta, deta, rho = 910 * … in THICreate() local
553 THIViscosity(thi, 0.5 * gradu * gradu, &eta, &deta); in THICreate()
558 …, ratio %8.2g\n", (double)u, (double)gradu, (double)eta, (double)(2 * eta * gradu), (double)(2 * e… in THICreate()
559 THIViscosity(thi, 0.5 * PetscSqr(1e-3 * gradu), &eta, &deta); in THICreate()
560 …uble)(1e-3 * u), (double)(1e-3 * gradu), (double)eta, (double)(2 * eta * 1e-3 * gradu), (double)(2… in THICreate()
[all …]
H A Dex15.c396 static inline PetscScalar eta(const AppCtx *ctx, PetscReal x, PetscReal y, PetscScalar ux, PetscSca… in eta() function
430 …[j - 1][i]), e_E = eta(user, xx, yy, ux_E, uy_E), e_W = eta(user, xx, yy, ux_W, uy_W), e_N = eta(u… in FormFunctionLocal()
516 …u = x[j][i], e_E = eta(user, xx, yy, ux_E, uy_E), e_W = eta(user, xx, yy, ux_W, uy_W), e_N = eta(u… in FormJacobianLocal()
784 …* (u - u_S), e_E = eta(user, xx, yy, ux_E, uy_E), e_W = eta(user, xx, yy, ux_W, uy_W), e_N = eta(u… in NonlinearGS()
H A Dex7.c401 Vec eta, etaHat; in TestFreeField() local
424 PetscCall(DMGetGlobalVector(dm, &eta)); in TestFreeField()
425 PetscCall(PetscObjectSetName((PetscObject)eta, "eta")); in TestFreeField()
447 PetscCall(ComputeResidual(dm, psi, eta)); in TestFreeField()
449 PetscCall(VecViewFromOptions(eta, NULL, "-eta_view")); in TestFreeField()
452 PetscCall(ComputeFFT(FT, Nc, eta, etaHat)); in TestFreeField()
517 PetscCall(DMRestoreGlobalVector(dm, &eta)); in TestFreeField()
/petsc/src/binding/petsc4py/docs/source/
H A Dpetsc_options.rst39 eta = OptDB.getReal('eta', 0.014)
48 the options, ``n`` and ``alpha`` will get the values ``50`` and ``8.8``, respectively, while ``eta`…
/petsc/src/ksp/ksp/impls/tfqmr/
H A Dtfqmr.c14 PetscScalar rho, rhoold, a, s, b, eta, etaold, psiold, cf; in KSPSolve_TFQMR() local
83 eta = cm * cm * a; in KSPSolve_TFQMR()
90 PetscCall(VecAXPY(X, eta, D)); in KSPSolve_TFQMR()
102 etaold = eta; in KSPSolve_TFQMR()
/petsc/src/ksp/ksp/impls/bcgs/qmrcgs/
H A Dqmrcgs.c18 PetscScalar eta, rho1, rho2, alpha, eta2, omega, beta, cf, cf1, uu; in KSPSolve_QMRCGS() local
77 eta = 1.0; in KSPSolve_QMRCGS()
117 cf = eta * theta * theta / alpha; in KSPSolve_QMRCGS()
177 eta = c * c * omega; in KSPSolve_QMRCGS()
181 PetscCall(VecWAXPY(X, eta, D, X2)); /* x <- x2 + eta d */ in KSPSolve_QMRCGS()
/petsc/src/ts/tutorials/
H A Dex14.c247 PRange eta; member
424 static void THIViscosity(THI thi, PetscReal gam, PetscReal *eta, PetscReal *deta) in THIViscosity() argument
440 *eta = Bd2 * PetscPowReal(eps + gam, exponent); in THIViscosity()
441 *deta = exponent * (*eta) / (eps + gam); in THIViscosity()
626 PRangeClear(&thi->eta); in THICreate()
630 … * units->meter / (3e7 * units->second), gradu = u / (100 * units->meter), eta, deta, rho = 910 * … in THICreate() local
632 THIViscosity(thi, 0.5 * gradu * gradu, &eta, &deta); in THICreate()
637 …%8.2g, ratio %8.2g\n", (double)u, (double)gradu, (double)eta, (double)(2 * eta * gradu, 2 * eta * … in THICreate()
638 THIViscosity(thi, 0.5 * PetscSqr(1e-3 * gradu), &eta, &deta); in THICreate()
639 …, (double)(1e-3 * u), (double)(1e-3 * gradu), (double)eta, (double)(2 * eta * 1e-3 * gradu, 2 * et… in THICreate()
[all …]
H A Dex48.c50 PetscReal mu, eta; member
67 options->eta = 0.001; in ProcessOptions()
90 PetscCall(PetscOptionsReal("-eta", "Viscosity", "mhd.c", options->eta, &options->eta, NULL)); in ProcessOptions()
96 PetscCall(PetscPrintf(comm, "eta = %g\n", (double)options->eta)); in ProcessOptions()
227 PetscReal eta = PetscRealPart(constants[ETA_CONST]); in f0_psi_4f() local
229 f0[0] = -eta * u[uOff[JZ]]; in f0_psi_4f()
239 PetscReal eta = PetscRealPart(constants[ETA_CONST]); in f0_psi_2f() local
241 f0[0] = -eta * u[uOff[JZ]]; in f0_psi_2f()
262 PetscReal eta = PetscRealPart(constants[ETA_CONST]); in g0_neta() local
264 g0[0] = -eta; in g0_neta()
[all …]
/petsc/src/tao/unconstrained/impls/cg/
H A Dtaocg.h14 PetscReal eta; /* Restart tolerance */ member
H A Dtaocg.c149 if (PetscAbsScalar(ginner) >= cgP->eta * gnorm2) { in TaoSolve_CG()
229 PetscCall(PetscOptionsReal("-tao_cg_eta", "restart tolerance", "", cgP->eta, &cgP->eta, NULL)); in TaoSetFromOptions_CG()
303 cgP->eta = 0.1; in TaoCreate_CG()
/petsc/src/dm/impls/stag/tutorials/
H A Dex4.c657 PetscScalar eta[4], eta_left, eta_right, eta_up, eta_down; in CreateSystem2d() local
693 PetscCall(DMStagVecGetValuesStencil(dm_coefficients, coefficients_local, 4, etaPoint, eta)); in CreateSystem2d()
694 eta_left = eta[0]; in CreateSystem2d()
695 eta_right = eta[1]; in CreateSystem2d()
696 eta_up = eta[2]; in CreateSystem2d()
697 eta_down = eta[3]; in CreateSystem2d()
811 PetscScalar eta[4], eta_left, eta_right, eta_up, eta_down; in CreateSystem2d() local
836 PetscCall(DMStagVecGetValuesStencil(dm_coefficients, coefficients_local, 4, etaPoint, eta)); in CreateSystem2d()
837 eta_left = eta[0]; in CreateSystem2d()
838 eta_right = eta[1]; in CreateSystem2d()
[all …]
/petsc/src/ksp/ksp/impls/fcg/pipefcg/
H A Dpipefcg.c85 PetscReal dp = 0.0, delta, *eta, *etas; in KSPSolve_PIPEFCG_cycle() local
192 eta = pipefcg->etas + idx; in KSPSolve_PIPEFCG_cycle()
229 *eta = 0.; in KSPSolve_PIPEFCG_cycle()
233 *eta -= PetscAbsScalar(betas[j]) * PetscAbsScalar(betas[j]) * etas[kdx]; in KSPSolve_PIPEFCG_cycle()
236 *eta += delta; /* etai = delta -betaik^2 * etak */ in KSPSolve_PIPEFCG_cycle()
237 if (*eta < 0.) { in KSPSolve_PIPEFCG_cycle()
242 alpha = gamma / (*eta); /* alpha = gamma/etai */ in KSPSolve_PIPEFCG_cycle()
/petsc/src/ksp/ksp/impls/cg/pipelcg/
H A Dpipelcg.c142 PetscScalar sum_dummy = 0.0, eta = 0.0, zeta = 0.0, lambda = 0.0; in KSPSolve_InnerLoop_PIPELCG() local
309 eta = gamma(0); in KSPSolve_InnerLoop_PIPELCG()
312 PetscCall(VecScale(p, 1.0 / eta)); in KSPSolve_InnerLoop_PIPELCG()
318 lambda = delta(k - 1) / eta; in KSPSolve_InnerLoop_PIPELCG()
319 eta = gamma(k) - lambda * delta(k - 1); in KSPSolve_InnerLoop_PIPELCG()
321 PetscCall(VecScale(p, -delta(k - 1) / eta)); in KSPSolve_InnerLoop_PIPELCG()
322 PetscCall(VecAXPY(p, 1.0 / eta, V[1])); in KSPSolve_InnerLoop_PIPELCG()
/petsc/src/snes/impls/ntrdc/
H A Dntrdcimpl.h12 PetscReal eta; /* used to compute trust region parameter */ member
/petsc/src/ksp/ksp/tutorials/output/
H A Dex82_1.out34 eta: 10.
83 eta: 10.
128 eta: 10.
/petsc/src/mat/tests/output/
H A Dex66_view.out69 Admissibility parameters: leaf size 17, eta 0.9
H A Dex66_view_cuda.out69 Admissibility parameters: leaf size 17, eta 0.9
/petsc/src/ksp/ksp/impls/gcr/pipegcr/
H A Dpipegcr.c63 PetscReal rnorm = 0.0, delta, *eta, *etas; in KSPSolve_PIPEGCR_cycle() local
141 eta = pipegcr->etas + idx; in KSPSolve_PIPEGCR_cycle()
214 *eta = 0.; in KSPSolve_PIPEGCR_cycle()
218 *eta -= PetscAbsScalar(betas[j]) * PetscAbsScalar(betas[j]) * etas[kdx]; in KSPSolve_PIPEGCR_cycle()
221 *eta += delta; /* etai = delta -betaik^2 * etak */ in KSPSolve_PIPEGCR_cycle()
224 if (*eta < 0.) { in KSPSolve_PIPEGCR_cycle()
229 alpha = gamma / (*eta); /* alpha = gamma/etai */ in KSPSolve_PIPEGCR_cycle()
/petsc/src/dm/impls/da/
H A Ddainterp.c407 PetscScalar *xi, *eta; in DMCreateInterpolation_DA_2D_Q1() local
414 PetscCall(PetscMalloc1(neta, &eta)); in DMCreateInterpolation_DA_2D_Q1()
416 … for (lj = 0; lj < neta; lj++) eta[lj] = -1.0 + (PetscScalar)lj * (2.0 / (PetscScalar)(neta - 1)); in DMCreateInterpolation_DA_2D_Q1()
455 Ni[0] = 0.25 * (1.0 - xi[li]) * (1.0 - eta[lj]); in DMCreateInterpolation_DA_2D_Q1()
456 Ni[1] = 0.25 * (1.0 + xi[li]) * (1.0 - eta[lj]); in DMCreateInterpolation_DA_2D_Q1()
457 Ni[2] = 0.25 * (1.0 - xi[li]) * (1.0 + eta[lj]); in DMCreateInterpolation_DA_2D_Q1()
458 Ni[3] = 0.25 * (1.0 + xi[li]) * (1.0 + eta[lj]); in DMCreateInterpolation_DA_2D_Q1()
470 PetscCall(PetscFree(eta)); in DMCreateInterpolation_DA_2D_Q1()
913 PetscScalar *xi, *eta, *zeta; in DMCreateInterpolation_DA_3D_Q1() local
922 PetscCall(PetscMalloc1(neta, &eta)); in DMCreateInterpolation_DA_3D_Q1()
[all …]

12