| /petsc/src/ksp/ksp/impls/bcgs/fbcgsr/ |
| H A D | fbcgsr.c | 20 PetscScalar tau, sigma, alpha, omega, beta; in KSPSolve_FBCGSR() local
107 PetscCall(VecDot(V, RP, &sigma)); /* sigma <- (v,rp) */ in KSPSolve_FBCGSR()
110 tau = sigma = 0.0; in KSPSolve_FBCGSR()
113 sigma += v[j] * rp[j]; /* sigma <- (v,rp) */ in KSPSolve_FBCGSR()
118 insums[1] = sigma; in KSPSolve_FBCGSR()
123 sigma = outsums[1]; in KSPSolve_FBCGSR()
127 alpha = tau / sigma; in KSPSolve_FBCGSR()
162 if ((xi3 == 0.0) || (sigma == 0.0)) { in KSPSolve_FBCGSR()
171 beta = -xi4 / sigma; in KSPSolve_FBCGSR()
|
| /petsc/src/ksp/ksp/tests/ |
| H A D | ex33.c | 17 PetscScalar v, sigma; in main() local
102 PetscCall(PetscOptionsGetScalar(NULL, NULL, "-sigma", &sigma, &flag)); in main()
104 sigma = -1.0 * sigma; in main()
105 PetscCall(MatAXPY(A, sigma, B, DIFFERENT_NONZERO_PATTERN)); /* A <- A - sigma*B */ in main()
|
| H A D | ex30.c | 38 PetscScalar sigma; in main() local
126 PetscCall(PetscOptionsGetScalar(NULL, NULL, "-mat_sigma", &sigma, &flg)); in main()
136 PetscCall(MatAXPY(A, sigma, B, DIFFERENT_NONZERO_PATTERN)); /* A <- sigma*B + A */ in main()
138 PetscCall(MatShift(A, sigma)); in main()
|
| /petsc/src/tao/linesearch/impls/armijo/ |
| H A D | armijo.c | 42 …scOptionsReal("-tao_ls_armijo_sigma", "acceptance constant", "", armP->sigma, &armP->sigma, NULL)); in TaoLineSearchSetFromOptions_Armijo()
63 PetscCall(PetscViewerASCIIPrintf(pv, "sigma=%g ", (double)armP->sigma)); in TaoLineSearchView_Armijo()
123 } else if ((armP->sigma <= 0) || (armP->sigma >= 0.5)) { in TaoLineSearchApply_Armijo()
124 PetscCall(PetscInfo(ls, "Armijo line search error: sigma (%g) invalid\n", (double)armP->sigma)); in TaoLineSearchApply_Armijo()
187 fact = armP->sigma; in TaoLineSearchApply_Armijo()
189 fact = armP->sigma * gdx; in TaoLineSearchApply_Armijo()
279 armP->sigma = 1e-4; in TaoLineSearchCreate_Armijo()
|
| H A D | armijo.h | 69 PetscReal sigma; /* Acceptance criteria < 1) */ member
|
| /petsc/src/dm/tutorials/ |
| H A D | ex26.c | 14 PetscReal *zeros, *weights, scale, h, sigma = 1 / sqrt(2), g = 0, mu = 0; in main() local
21 PetscCall(PetscOptionsGetReal(NULL, NULL, "-sigma", &sigma, NULL)); in main()
51 …for (s = 0; s < n; s++) g += weights[s] * PetscPowRealInt(sqrt(2) * sigma * zeros[s] + mu, moment); in main()
|
| /petsc/src/tao/linesearch/impls/owarmijo/ |
| H A D | owarmijo.c | 58 …OptionsReal("-tao_ls_OWArmijo_sigma", "acceptance constant", "", armP->sigma, &armP->sigma, NULL)); in TaoLineSearchSetFromOptions_OWArmijo()
78 PetscCall(PetscViewerASCIIPrintf(pv, "sigma=%g ", (double)armP->sigma)); in TaoLineSearchView_OWArmijo()
153 } else if ((armP->sigma <= 0) || (armP->sigma >= 0.5)) { in TaoLineSearchApply_OWArmijo()
154 … PetscCall(PetscInfo(ls, "OWArmijo line search error: sigma (%g) invalid\n", (double)armP->sigma)); in TaoLineSearchApply_OWArmijo()
204 if (armP->nondescending) fact = armP->sigma; in TaoLineSearchApply_OWArmijo()
246 if (!armP->nondescending && *f <= ref + armP->sigma * gdx) break; in TaoLineSearchApply_OWArmijo()
303 armP->sigma = 1e-4; in TaoLineSearchCreate_OWArmijo()
|
| H A D | owarmijo.h | 69 PetscReal sigma; /* Acceptance criteria < 1) */ member
|
| /petsc/src/ksp/ksp/utils/lmvm/rescale/ |
| H A D | symbrdnrescale.c | 12 PetscReal sigma_inv, sigma; in SymBroydenRescaleUpdateScalar() local
19 sigma = 1.0 / sigma_inv; in SymBroydenRescaleUpdateScalar()
60 sigma = ldb->rho * signew + (1.0 - ldb->rho) * sigma; in SymBroydenRescaleUpdateScalar()
61 PetscCall(MatLMVMSetJ0Scale(B, 1.0 / sigma)); in SymBroydenRescaleUpdateScalar()
129 PetscReal sigma; in SymBroydenRescaleUpdateDiagonal() local
259 sigma = ss_sum / ys_sum; in SymBroydenRescaleUpdateDiagonal()
262 sigma = ys_sum / yy_sum; in SymBroydenRescaleUpdateDiagonal()
274 sigma = (2 * c) / denom; in SymBroydenRescaleUpdateDiagonal()
278 sigma = (-b - sqrt_disc) / denom; in SymBroydenRescaleUpdateDiagonal()
284 sigma = (-b + sqrt_disc) / denom; in SymBroydenRescaleUpdateDiagonal()
[all …]
|
| /petsc/src/ksp/ksp/impls/cg/pipelcg/ |
| H A D | pipelcg.c | 11 #define sigma(i) (plcg->sigma[i]) macro
25 PetscReal *sigma; /* base shifts */ member
56 PetscCall(PetscCalloc1(l, &plcg->sigma)); in KSPSetUp_PIPELCG()
66 PetscCall(PetscFree(plcg->sigma)); in KSPReset_PIPELCG()
163 PetscCall(VecAXPY(U[0], -sigma(it), U[1])); in KSPSolve_InnerLoop_PIPELCG()
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()
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()
[all …]
|
| /petsc/src/dm/tests/ |
| H A D | ex30.c | 12 PetscReal alpha = 1, K = 1, rho0 = 1, u0 = 0, sigma = 0.2; in main() local
34 …all(PetscOptionsReal("-sigma", "Width of Gaussian density perturbation", "", sigma, &sigma, NULL)); in main()
130 xx[2 * i] = 0.2 * PetscExpReal(-PetscSqr(xref) / (2 * PetscSqr(sigma))); in main()
|
| /petsc/src/ksp/ksp/tutorials/ |
| H A D | ex21.c | 6 PetscReal sigma; member
16 PetscReal s = rbfctx->sigma; in RBF()
60 PetscCall(PetscRandomGetValueReal(r, &fctx.sigma)); in main()
61 PetscCall(PetscOptionsGetReal(NULL, NULL, "-sigma", &fctx.sigma, NULL)); in main()
|
| /petsc/src/snes/tests/output/ |
| H A D | ex15_particle_hdiv_5.out | 4 sigma: 0.0795775
|
| /petsc/src/ts/tutorials/hamiltonian/output/ |
| H A D | ex3_landau_damping_3d.out | 2 sigma: 1.
|
| H A D | ex3_landau_damping_2d.out | 2 sigma: 1.
|
| H A D | ex2_none_1d.out | 2 sigma: 0.0795773
|
| H A D | ex3_landau_damping_1d.out | 2 sigma: 1.
|
| H A D | ex2_coulomb_1d.out | 2 sigma: 0.0795773
|
| H A D | ex4_landau_damping_1d_bs.out | 2 sigma: 1.
|
| H A D | ex4_landau_damping_1d_dg.out | 2 sigma: 1.
|
| H A D | ex3_sphere_3d.out | 8 sigma: 1.00394
|
| /petsc/src/tao/complementarity/tutorials/ |
| H A D | blackscholes.c | 101 PetscReal sigma, alpha, delta; /* Underlying asset properties */ member
140 user.sigma = 0.40; in main()
157 PetscCall(PetscOptionsGetReal(NULL, NULL, "-sigma", &user.sigma, &flg)); in main()
191 user.d[i] = -0.5 * user.sigma * user.sigma * PetscPowReal(sval, user.alpha); in main()
|
| /petsc/src/snes/impls/ntrdc/ |
| H A D | ntrdcimpl.h | 31 PetscReal sigma; /* used to determine termination */ member
|
| /petsc/src/dm/dt/interface/ |
| H A D | dtprob.c | 159 const PetscReal sigma = scale ? scale[0] : 1.; in PetscPDFGaussian1D() local
160 … = PetscSqrtReal(1. / (2. * PETSC_PI)) * PetscExpReal(-0.5 * PetscSqr(x[0] / sigma)) / sigma; in PetscPDFGaussian1D()
166 const PetscReal sigma = scale ? scale[0] : 1.; in PetscCDFGaussian1D() local
167 p[0] = 0.5 * (1. + PetscErfReal(x[0] / PETSC_SQRT2 / sigma)); in PetscCDFGaussian1D()
|
| /petsc/src/ksp/ksp/impls/cg/gltr/ |
| H A D | gltr.c | 76 PetscInt i, j, max_cg_its, max_lanczos_its, max_newton_its, sigma; in KSPCGSolve_GLTR()
842 sigma = -1; in KSPCGSolve_GLTR()
850 PetscCall(VecScale(d, sigma * t_soln[0] / cg->norm_r[0])); in KSPCGSolve_GLTR()
860 if (alpha >= 0.0) sigma = -sigma; in KSPCGSolve_GLTR()
866 PetscCall(VecAXPY(d, sigma * t_soln[i + 1] / cg->norm_r[i + 1], z)); in KSPCGSolve_GLTR()
877 if (alpha >= 0.0) sigma = -sigma; in KSPCGSolve_GLTR()
883 PetscCall(VecAXPY(d, sigma * t_soln[i + 1] / cg->norm_r[i + 1], z)); in KSPCGSolve_GLTR()
|