Home
last modified time | relevance | path

Searched refs:gamma (Results 1 – 25 of 79) sorted by relevance

1234

/petsc/src/ksp/ksp/impls/cg/
H A Dcgls.c32 PetscReal alpha, gamma, oldgamma; in KSPSolve_CGLS() local
54 PetscCall(VecNorm(ss, NORM_2, &gamma)); in KSPSolve_CGLS()
55 KSPCheckNorm(ksp, gamma); in KSPSolve_CGLS()
56 if (ksp->normtype != KSP_NORM_NONE) ksp->rnorm = gamma; in KSPSolve_CGLS()
59 gamma = gamma * gamma; /* gamma = norm2(s)^2 */ in KSPSolve_CGLS()
66 alpha = gamma / alpha; /* alpha = gamma / alpha */ in KSPSolve_CGLS()
70 oldgamma = gamma; /* oldgamma = gamma */ in KSPSolve_CGLS()
71 PetscCall(VecNorm(ss, NORM_2, &gamma)); in KSPSolve_CGLS()
72 KSPCheckNorm(ksp, gamma); in KSPSolve_CGLS()
74 if (ksp->normtype != KSP_NORM_NONE) ksp->rnorm = gamma; in KSPSolve_CGLS()
[all …]
/petsc/src/ksp/ksp/impls/cg/pipecg/
H A Dpipecg.c23 PetscScalar alpha = 0.0, beta = 0.0, gamma = 0.0, gammaold = 0.0, delta = 0.0; in KSPSolve_PIPECG() local
71 PetscCall(VecDotBegin(R, U, &gamma)); /* gamma <- u'*r */ in KSPSolve_PIPECG()
74 PetscCall(VecDotEnd(R, U, &gamma)); in KSPSolve_PIPECG()
75 KSPCheckDot(ksp, gamma); in KSPSolve_PIPECG()
76 dp = PetscSqrtReal(PetscAbsScalar(gamma)); /* dp <- r'*u = r'*B*r = e'*A'*B*A*e */ in KSPSolve_PIPECG()
98 if (!(i == 0 && ksp->normtype == KSP_NORM_NATURAL)) PetscCall(VecDotBegin(R, U, &gamma)); in KSPSolve_PIPECG()
110 if (!(i == 0 && ksp->normtype == KSP_NORM_NATURAL)) PetscCall(VecDotEnd(R, U, &gamma)); in KSPSolve_PIPECG()
114 if (ksp->normtype == KSP_NORM_NATURAL) dp = PetscSqrtReal(PetscAbsScalar(gamma)); in KSPSolve_PIPECG()
125 alpha = gamma / delta; in KSPSolve_PIPECG()
131 beta = gamma / gammaold; in KSPSolve_PIPECG()
[all …]
/petsc/src/ksp/ksp/impls/cr/pipecr/
H A Dpipecr.c23 PetscScalar alpha = 0.0, beta = 0.0, gamma, gammaold = 0.0, delta; in KSPSolve_PIPECR() local
80 PetscCall(VecDotBegin(W, U, &gamma)); in KSPSolve_PIPECR()
87 PetscCall(VecDotEnd(W, U, &gamma)); in KSPSolve_PIPECR()
100 alpha = gamma / delta; in KSPSolve_PIPECR()
105 beta = gamma / gammaold; in KSPSolve_PIPECR()
106 alpha = gamma / (delta - beta / alpha * gamma); in KSPSolve_PIPECR()
114 gammaold = gamma; in KSPSolve_PIPECR()
/petsc/src/snes/impls/ms/
H A Dms.c13 PetscReal *gamma; /* Coefficients of 3S* method */ member
51 PetscReal gamma[3][6] = { in SNESMSRegisterAll() local
58 PetscCall(SNESMSRegister(SNESMSM62, 6, 3, 1.0, &gamma[0][0], delta, betasub)); in SNESMSRegisterAll()
112 PetscCall(PetscFree(t->gamma)); in SNESMSRegisterDestroy()
193 …tscInt nstages, PetscInt nregisters, PetscReal stability, const PetscReal gamma[], const PetscReal… in SNESMSRegister() argument
201 if (gamma || delta) { in SNESMSRegister()
203 PetscAssertPointer(gamma, 5); in SNESMSRegister()
218 if (gamma && delta) { in SNESMSRegister()
219 PetscCall(PetscMalloc1(nstages * nregisters, &t->gamma)); in SNESMSRegister()
221 PetscCall(PetscArraycpy(t->gamma, gamma, nstages * nregisters)); in SNESMSRegister()
[all …]
/petsc/src/ksp/ksp/impls/cg/groppcg/
H A Dgroppcg.c26 PetscScalar alpha, beta = 0.0, gamma, gammaNew, t; in KSPSolve_GROPPCG() local
57 PetscCall(VecDotBegin(r, z, &gamma)); /* gamma <- z'*r */ in KSPSolve_GROPPCG()
60 PetscCall(VecDotEnd(r, z, &gamma)); /* gamma <- z'*r */ in KSPSolve_GROPPCG()
72 KSPCheckDot(ksp, gamma); in KSPSolve_GROPPCG()
73 dp = PetscSqrtReal(PetscAbsScalar(gamma)); /* dp <- r'*z = r'*B*r = e'*A'*B*A*e */ in KSPSolve_GROPPCG()
99 alpha = gamma / t; in KSPSolve_GROPPCG()
133 beta = gammaNew / gamma; in KSPSolve_GROPPCG()
134 gamma = gammaNew; in KSPSolve_GROPPCG()
/petsc/src/ksp/ksp/tests/
H A Dex38.c33 PetscReal h, gamma, beta; /* parameters */ in main() local
46 gamma = 4.0; in main()
49 PetscCall(PetscOptionsGetReal(NULL, NULL, "-gamma", &gamma, NULL)); in main()
51 gamma = gamma / h; in main()
90 co1 = gamma * h * h / 2.0; in main()
H A Dex39.c24 PetscReal h, gamma, beta; /* parameters */ in main() local
38 gamma = 4.0 / h; in main()
41 PetscCall(PetscOptionsGetReal(NULL, NULL, "-gamma", &gamma, NULL)); in main()
62 co1 = gamma * h * h / 2.0; in main()
/petsc/src/ts/impls/implicit/alpha/
H A Dalpha1.c117 PetscReal alpha_m, alpha_f, gamma; in TSAlpha_Restart() local
125 PetscCall(TSAlphaGetParams(ts, &alpha_m, &alpha_f, &gamma)); in TSAlpha_Restart()
172 PetscCall(TSAlphaSetParams(ts, alpha_m, alpha_f, gamma)); in TSAlpha_Restart()
404 PetscReal alpha_m, alpha_f, gamma; in TSAlphaSetRadius_Alpha() local
410 gamma = (PetscReal)0.5 + alpha_m - alpha_f; in TSAlphaSetRadius_Alpha()
411 PetscCall(TSAlphaSetParams(ts, alpha_m, alpha_f, gamma)); in TSAlphaSetRadius_Alpha()
415 …PetscErrorCode TSAlphaSetParams_Alpha(TS ts, PetscReal alpha_m, PetscReal alpha_f, PetscReal gamma) in TSAlphaSetParams_Alpha() argument
419 PetscReal res = ((PetscReal)0.5 + alpha_m - alpha_f) - gamma; in TSAlphaSetParams_Alpha()
424 th->Gamma = gamma; in TSAlphaSetParams_Alpha()
429 …scErrorCode TSAlphaGetParams_Alpha(TS ts, PetscReal *alpha_m, PetscReal *alpha_f, PetscReal *gamma) in TSAlphaGetParams_Alpha() argument
[all …]
H A Dalpha2.c159 PetscReal alpha_m, alpha_f, gamma, beta; in TSAlpha_Restart() local
169 PetscCall(TSAlpha2GetParams(ts, &alpha_m, &alpha_f, &gamma, &beta)); in TSAlpha_Restart()
226 PetscCall(TSAlpha2SetParams(ts, alpha_m, alpha_f, gamma, beta)); in TSAlpha_Restart()
506 PetscReal alpha_m, alpha_f, gamma, beta; in TSAlpha2SetRadius_Alpha() local
512 gamma = (PetscReal)0.5 + alpha_m - alpha_f; in TSAlpha2SetRadius_Alpha()
515 PetscCall(TSAlpha2SetParams(ts, alpha_m, alpha_f, gamma, beta)); in TSAlpha2SetRadius_Alpha()
519 …Alpha2SetParams_Alpha(TS ts, PetscReal alpha_m, PetscReal alpha_f, PetscReal gamma, PetscReal beta) in TSAlpha2SetParams_Alpha() argument
523 PetscReal res = ((PetscReal)0.5 + alpha_m - alpha_f) - gamma; in TSAlpha2SetParams_Alpha()
528 th->Gamma = gamma; in TSAlpha2SetParams_Alpha()
534 …a2GetParams_Alpha(TS ts, PetscReal *alpha_m, PetscReal *alpha_f, PetscReal *gamma, PetscReal *beta) in TSAlpha2GetParams_Alpha() argument
[all …]
/petsc/src/ksp/ksp/impls/cg/pipecgrr/
H A Dpipecgrr.c23 …PetscScalar alpha = 0.0, beta = 0.0, gamma = 0.0, gammaold = 0.0, delta = 0.0, alphap = 0.0, betap… in KSPSolve_PIPECGRR() local
76 PetscCall(VecDotBegin(R, U, &gamma)); /* gamma <- u'*r */ in KSPSolve_PIPECGRR()
80 PetscCall(VecDotEnd(R, U, &gamma)); in KSPSolve_PIPECGRR()
82 KSPCheckDot(ksp, gamma); in KSPSolve_PIPECGRR()
83 dp = PetscSqrtReal(PetscAbsScalar(gamma)); /* dp <- r'*u = r'*B*r = e'*A'*B*A*e */ in KSPSolve_PIPECGRR()
124 if (!(i == 0 && ksp->normtype == KSP_NORM_NATURAL)) PetscCall(VecDotBegin(R, U, &gamma)); in KSPSolve_PIPECGRR()
147 if (!(i == 0 && ksp->normtype == KSP_NORM_NATURAL)) PetscCall(VecDotEnd(R, U, &gamma)); in KSPSolve_PIPECGRR()
162 if (ksp->normtype == KSP_NORM_NATURAL) dp = PetscSqrtReal(PetscAbsScalar(gamma)); in KSPSolve_PIPECGRR()
173 alpha = gamma / delta; in KSPSolve_PIPECGRR()
179 beta = gamma / gammaold; in KSPSolve_PIPECGRR()
[all …]
/petsc/src/ksp/ksp/impls/cg/pipelcg/
H A Dpipelcg.c9 #define gamma(i) (plcg->gamma[i]) macro
23 PetscScalar *gamma, *delta, *alpha; member
224 … gamma(it - l) = (G(it - l, it - l + 1) + sigma(it - l) * G(it - l, it - l)) / G(it - l, it - l); in KSPSolve_InnerLoop_PIPELCG()
226 …gamma(it - l) = (G(it - l, it - l + 1) + sigma(it - l) * G(it - l, it - l) - delta(it - l - 1) * G… in KSPSolve_InnerLoop_PIPELCG()
230 …gamma(it - l) = (G(it - l, it - l) * gamma(it - 2 * l) + G(it - l, it - l + 1) * delta(it - 2 * l)… in KSPSolve_InnerLoop_PIPELCG()
247 PetscCall(VecAXPY(V[0], sigma(0) - gamma(it - l), V[1])); in KSPSolve_InnerLoop_PIPELCG()
249 alpha(0) = sigma(0) - gamma(it - l); in KSPSolve_InnerLoop_PIPELCG()
268 PetscCall(VecAXPY(Q[3 * j], sigma(j + 1) - gamma(it - l), Q[3 * j + 1])); in KSPSolve_InnerLoop_PIPELCG()
270 alpha(0) = sigma(j + 1) - gamma(it - l); in KSPSolve_InnerLoop_PIPELCG()
279 PetscCall(VecAXPY(Z[0], -gamma(it - l), Z[1])); in KSPSolve_InnerLoop_PIPELCG()
[all …]
/petsc/src/ts/tutorials/multirate/
H A Dex7.c174 PetscScalar *x, *f, *r, *min, *alpha, *gamma; in FVRHSFunction() local
190 PetscCall(PetscMalloc4(dof, &r, dof, &min, dof, &alpha, dof, &gamma)); in FVRHSFunction()
207 gamma[0] = 1.0 / 3.0; in FVRHSFunction()
211 …u[j] = x[(i - 1) * dof + j] + PetscMax(0, PetscMin(min[j], alpha[0] + gamma[0] * r[j])) * (x[(i … in FVRHSFunction()
224 gamma[0] = hs * (hs + hs) / (hs + hf) / (hs + hs + hf); in FVRHSFunction()
228 …u[j] = x[(i - 1) * dof + j] + PetscMax(0, PetscMin(min[j], alpha[0] + gamma[0] * r[j])) * (x[(i … in FVRHSFunction()
240 gamma[0] = hf * (hs + hf) / (hf + hf) / (hs + hf + hf); in FVRHSFunction()
244 …u[j] = x[(i - 1) * dof + j] + PetscMax(0, PetscMin(min[j], alpha[0] + gamma[0] * r[j])) * (x[(i … in FVRHSFunction()
256 gamma[0] = 1.0 / 3.0; in FVRHSFunction()
260 …u[j] = x[(i - 1) * dof + j] + PetscMax(0, PetscMin(min[j], alpha[0] + gamma[0] * r[j])) * (x[(i … in FVRHSFunction()
[all …]
/petsc/src/vec/vec/impls/seq/
H A Dbvec1.c120 PetscErrorCode VecAXPBYPCZ_Seq(Vec zin, PetscScalar alpha, PetscScalar beta, PetscScalar gamma, Vec… in VecAXPBYPCZ_Seq() argument
132 for (PetscInt i = 0; i < n; ++i) zz[i] = xx[i] + beta * yy[i] + gamma * zz[i]; in VecAXPBYPCZ_Seq()
133 } else if (gamma == (PetscScalar)1.0) { in VecAXPBYPCZ_Seq()
135 } else if (gamma == (PetscScalar)0.0) { in VecAXPBYPCZ_Seq()
139 for (PetscInt i = 0; i < n; ++i) zz[i] = alpha * xx[i] + beta * yy[i] + gamma * zz[i]; in VecAXPBYPCZ_Seq()
/petsc/src/snes/tutorials/
H A Dex31.c35 PetscReal gamma; member
57 blasius->gamma = 1.4; /* specific heat ratio */ in main()
63 …(PetscOptionsReal("-gamma", "Ratio of specific heats", "", blasius->gamma, &blasius->gamma, NULL)); in main()
190 r[N + 2 + i] = h[2] + Pr * f[0] * h[1] + Pr * (blasius->gamma - 1) * PetscSqr(Ma * f[2]); in FormFunction()
H A Dex74ftube.h1 …ble precision ru1, ru2, ru3, ru4, et1, et2, et3, et4 double precision csubv, gamma, a1, a2, a3, a4,
6 …2, e3, e4, v1, v2, v3, v4, ru1, ru2, ru3, ru4, et1, et2, et3, et4, csubv, gamma, a1, a2, a3, a4, s…
/petsc/src/ksp/ksp/impls/cg/pipecg2/
H A Dpipecg2.c552 PetscScalar alpha[2], beta[2], gamma[2], delta[2], lambda[15]; in KSPSolve_PIPECG2() local
591 PetscCall(PetscMemzero(gamma, 2 * sizeof(PetscScalar))); in KSPSolve_PIPECG2()
609 …PetscCall(VecMergedDot_Private(U, W, R, ksp->normtype, &gamma[0], &delta[0], &dps)); /* gamma <-… in KSPSolve_PIPECG2()
610 lambda[10] = gamma[0]; in KSPSolve_PIPECG2()
632 gamma[0] = lambda[10]; in KSPSolve_PIPECG2()
653 alpha[0] = gamma[0] / delta[0]; in KSPSolve_PIPECG2()
655 …gamma[1] = gamma[0] - alpha[0] * lambda[13] - alpha[0] * delta[0] + alpha[0] * alpha[0] * lambda[1… in KSPSolve_PIPECG2()
658 beta[1] = gamma[1] / gamma[0]; in KSPSolve_PIPECG2()
659 alpha[1] = gamma[1] / (delta[1] - beta[1] / alpha[0] * gamma[1]); in KSPSolve_PIPECG2()
663 beta[0] = gamma[1] / gamma[0]; in KSPSolve_PIPECG2()
[all …]
/petsc/src/ksp/pc/tutorials/
H A Dex4.c81 PetscReal gamma, alpha; member
111 PetscCall(MatScale(ctx->I_plus_gammaUTaDinvU, ctx->gamma)); in SmwSetup()
144 PetscCall(VecScale(y, -ctx->gamma)); in SmwApply()
173 PetscReal gamma = 100, alpha = .01; in main() local
251 PetscCall(MatScale(J, gamma)); in main()
296 ctx.gamma = gamma; in main()
/petsc/src/ts/tutorials/
H A Dex11.h88 PetscReal gamma; member
333 static PetscErrorCode Pressure_PG(const PetscReal gamma, const EulerNode *x, PetscReal *p) in Pressure_PG() argument
339 (*p) = (x->E - 0.5 * ru2 / x->r) * (gamma - 1.0); /* (E - rho V^2/2)(gamma-1) = e rho (gamma-1) */ in Pressure_PG()
343 static PetscErrorCode SpeedOfSound_PG(const PetscReal gamma, const EulerNode *x, PetscReal *c) in SpeedOfSound_PG() argument
348 PetscCall(Pressure_PG(gamma, x, &p)); in SpeedOfSound_PG()
350 (*c) = PetscSqrtReal(gamma * p / x->r); in SpeedOfSound_PG()
368 PetscCall(Pressure_PG(eu->gamma, x, &p)); in EulerFlux()
675 …alar *ul, const PetscScalar *ur, PetscScalar *flux, const PetscReal *nn, int ndim, PetscReal gamma) in godunovflux() argument
748 gaml = gamma; in godunovflux()
749 gamr = gamma; in godunovflux()
[all …]
H A Dex11.c472 int initLinearWave(EulerNode *ux, const PetscReal gamma, const PetscReal coord[],…
479 PetscReal p0, gamma, c = 0.; in PhysicsSolution_Euler() local
486 gamma = eu->gamma; in PhysicsSolution_Euler()
501 p1 = press * (1.0 + 2.0 * gamma / (gamma + 1.0) * (amach * amach - 1.0)); in PhysicsSolution_Euler()
502 gas1 = (gamma - 1.0) / (gamma + 1.0); in PhysicsSolution_Euler()
504 uu->ru[0] = ((uu->r - rho) * PetscSqrtReal(gamma * press / rho) * amach); in PhysicsSolution_Euler()
505 uu->E = p1 / (gamma - 1.0) + .5 / uu->r * uu->ru[0] * uu->ru[0]; in PhysicsSolution_Euler()
508 uu->E = p0 / (gamma - 1.0); in PhysicsSolution_Euler()
512 uu->E = p0 / (gamma - 1.0); in PhysicsSolution_Euler()
518 uu->E = 10. / (gamma - 1.); in PhysicsSolution_Euler()
[all …]
/petsc/src/ts/tutorials/advection-diffusion-reaction/
H A Dreaction_diffusion.c77 f[j][i].u = appctx->D1 * (uxx + uyy) - uc * vc * vc + appctx->gamma * (1.0 - uc); in RHSFunction()
78 f[j][i].v = appctx->D2 * (vxx + vyy) + uc * vc * vc - (appctx->gamma + appctx->kappa) * vc; in RHSFunction()
177 entries[4] = -2.0 * appctx->D1 * (sx + sy) - vc * vc - appctx->gamma; in RHSJacobian()
195 entries[4] = -2.0 * appctx->D2 * (sx + sy) + 2.0 * uc * vc - appctx->gamma - appctx->kappa; in RHSJacobian()
286 … f[j][i].u = udot[j][i].u - (appctx->D1 * (uxx + uyy) - uc * vc * vc + appctx->gamma * (1.0 - uc)); in IFunction()
287 …f[j][i].v = udot[j][i].v - (appctx->D2 * (vxx + vyy) + uc * vc * vc - (appctx->gamma + appctx->kap… in IFunction()
387 entries[4] = 2.0 * appctx->D1 * (sx + sy) + vc * vc + appctx->gamma + a; in IJacobian()
405 … entries[4] = 2.0 * appctx->D2 * (sx + sy) - 2.0 * uc * vc + appctx->gamma + appctx->kappa + a; in IJacobian()
/petsc/src/ts/impls/rosw/
H A Drosw.c961 PetscErrorCode TSRosWRegisterRos4(TSRosWType name, PetscReal gamma, PetscReal a2, PetscReal a3, Pet… in TSRosWRegisterRos4() argument
964 …gamma + gamma * gamma, p42 = one / eight - gamma / three, p43 = one / twelve - gamma / three, p44 … in TSRosWRegisterRos4()
1006 rhs[0] = one / two - gamma; in TSRosWRegisterRos4()
1039 Gamma[0][0] = gamma; in TSRosWRegisterRos4()
1044 Gamma[1][1] = gamma; in TSRosWRegisterRos4()
1049 Gamma[2][2] = gamma; in TSRosWRegisterRos4()
1054 Gamma[3][3] = gamma; in TSRosWRegisterRos4()
1061 …bm[3] = b[3] - e4 * gamma; /* using definition of E4 … in TSRosWRegisterRos4()
1063 bm[1] = (one / two - gamma - beta3p * bm[2] - beta4p * bm[3]) / beta2p; /* second row */ in TSRosWRegisterRos4()
/petsc/src/sys/classes/draw/utils/
H A Dcmap.c254 PetscReal gamma = (beta > 0.0) ? (1 - beta) : (1 / (1 + beta)); in PetscDrawUtilitySetCmap() local
256 PetscReal r = PetscPowReal((PetscReal)R[i] / 255, gamma); in PetscDrawUtilitySetCmap()
257 PetscReal g = PetscPowReal((PetscReal)G[i] / 255, gamma); in PetscDrawUtilitySetCmap()
258 PetscReal b = PetscPowReal((PetscReal)B[i] / 255, gamma); in PetscDrawUtilitySetCmap()
/petsc/src/ksp/ksp/impls/fcg/pipefcg/
H A Dpipefcg.c84 PetscScalar alpha = 0.0, gamma, *betas, *dots; in KSPSolve_PIPEFCG_cycle() local
140 gamma = dots[0]; in KSPSolve_PIPEFCG_cycle()
143 alpha = gamma / delta; in KSPSolve_PIPEFCG_cycle()
164 dp = PetscSqrtReal(PetscAbsScalar(gamma)); /* dp <- sqrt(r'*z) = sqrt(e'*A'*B*A*e) */ in KSPSolve_PIPEFCG_cycle()
226 gamma = betas[j]; in KSPSolve_PIPEFCG_cycle()
242 alpha = gamma / (*eta); /* alpha = gamma/etai */ in KSPSolve_PIPEFCG_cycle()
263 PetscScalar gamma; in KSPSolve_PIPEFCG() local
299 PetscCall(VecXDot(Z, R, &gamma)); in KSPSolve_PIPEFCG()
300 dp = PetscSqrtReal(PetscAbsScalar(gamma)); /* dp <- sqrt(r'*z) = sqrt(e'*A'*B*A*e) */ in KSPSolve_PIPEFCG()
/petsc/src/tao/bound/impls/bncg/
H A Dbncg.c439 PetscReal gamma = 1.0, tau_k, beta; in TaoBNCGStepDirectionUpdate() local
719 gamma = 0.0; in TaoBNCGStepDirectionUpdate()
721 if (gkp1_yk < 0 && cg->neg_xi) gamma = -1.0 * gd / dk_yk; in TaoBNCGStepDirectionUpdate()
724 else gamma = cg->xi * gd / dk_yk; in TaoBNCGStepDirectionUpdate()
727 … PetscCall(VecAXPBYPCZ(tao->stepdirection, -tau_k, gamma * tau_k, beta, tao->gradient, cg->yk)); in TaoBNCGStepDirectionUpdate()
737 gamma = gd / dk_yk; in TaoBNCGStepDirectionUpdate()
742 beta = gkp1D_yk / dk_yk - step * gamma - tau_k; in TaoBNCGStepDirectionUpdate()
747 if (gkp1D_yk / dk_yk < 0) gamma = -1.0 * gd / dk_yk; in TaoBNCGStepDirectionUpdate()
748 else gamma = cg->xi * gd / dk_yk; in TaoBNCGStepDirectionUpdate()
751 gamma = 0.0; in TaoBNCGStepDirectionUpdate()
[all …]
/petsc/src/ksp/ksp/impls/gcr/pipegcr/
H A Dpipegcr.c62 PetscScalar alpha = 0.0, gamma, *betas, *dots; in KSPSolve_PIPEGCR_cycle() local
108 gamma = dots[1]; in KSPSolve_PIPEGCR_cycle()
109 alpha = gamma / delta; in KSPSolve_PIPEGCR_cycle()
181 gamma = dots[j]; in KSPSolve_PIPEGCR_cycle()
195 rnorm = PetscSqrtReal(PetscAbsScalar(gamma)); /* ||r|| <- sqrt(r,w) */ in KSPSolve_PIPEGCR_cycle()
229 alpha = gamma / (*eta); /* alpha = gamma/etai */ in KSPSolve_PIPEGCR_cycle()
254 PetscScalar gamma; in KSPSolve_PIPEGCR() local
279 PetscCall(VecDot(r, w, &gamma)); /* gamma = (r,w) */ in KSPSolve_PIPEGCR()
289 rnorm = PetscSqrtReal(PetscAbsScalar(gamma)); /* ||r|| <- sqrt(r,w) */ in KSPSolve_PIPEGCR()

1234