Home
last modified time | relevance | path

Searched refs:delta (Results 1 – 25 of 70) sorted by relevance

123

/petsc/src/tao/unconstrained/impls/cg/ !
H A Dtaocg.c17 PetscReal step = 1.0, f, gnorm, gnorm2, delta, gd, ginner, beta; in TaoSolve_CG() local
41 delta = 2.0 * PetscAbsScalar(f) / gnorm2; in TaoSolve_CG()
42 delta = PetscMax(delta, cgP->delta_min); in TaoSolve_CG()
43 delta = PetscMin(delta, cgP->delta_max); in TaoSolve_CG()
45 delta = 2.0 / gnorm2; in TaoSolve_CG()
46 delta = PetscMax(delta, cgP->delta_min); in TaoSolve_CG()
47 delta = PetscMin(delta, cgP->delta_max); in TaoSolve_CG()
68 delta = 2.0 * PetscAbsScalar(f) / gnorm2; in TaoSolve_CG()
69 delta = PetscMax(delta, cgP->delta_min); in TaoSolve_CG()
70 delta = PetscMin(delta, cgP->delta_max); in TaoSolve_CG()
[all …]
/petsc/src/binding/petsc4py/demo/legacy/kspsolve/ !
H A Dpetsc-cg.py17 delta = delta_0
20 delta > delta_0 * eps**2:
22 alpha = delta / d.dot(q)
25 delta_old = delta
26 delta = r.dot(r)
27 beta = delta / delta_old
30 return i, delta**0.5
/petsc/src/ksp/ksp/impls/qcg/ !
H A Dqcg.c20 static PetscErrorCode KSPQCGQuadraticRoots(Vec s, Vec p, PetscReal delta, PetscReal *step1, PetscRe… in KSPQCGQuadraticRoots() argument
28 dsq = delta * delta; in KSPQCGQuadraticRoots()
59 PetscErrorCode KSPQCGSetTrustRegionRadius(KSP ksp, PetscReal delta) in KSPQCGSetTrustRegionRadius() argument
63 …PetscCheck(delta >= 0.0, PetscObjectComm((PetscObject)ksp), PETSC_ERR_ARG_OUTOFRANGE, "Tolerance m… in KSPQCGSetTrustRegionRadius()
64 PetscTryMethod(ksp, "KSPQCGSetTrustRegionRadius_C", (KSP, PetscReal), (ksp, delta)); in KSPQCGSetTrustRegionRadius()
173 …PetscCheck(pcgP->delta > dzero, PetscObjectComm((PetscObject)ksp), PETSC_ERR_ARG_OUTOFRANGE, "Inpu… in KSPSolve_QCG()
222 scal = pcgP->delta / xnorm; in KSPSolve_QCG()
226 PetscCall(KSPQCGQuadraticRoots(W, P, pcgP->delta, &step1, &step2)); in KSPSolve_QCG()
238 pcgP->ltsnrm = pcgP->delta; /* convergence in direction of */ in KSPSolve_QCG()
241 PetscCall(PetscInfo(ksp, "negative curvature: delta=%g\n", (double)pcgP->delta)); in KSPSolve_QCG()
[all …]
/petsc/src/tao/leastsquares/impls/pounders/ !
H A Dgqt.c214 PetscErrorCode gqt(PetscInt n, PetscReal *a, PetscInt lda, PetscReal *b, PetscReal delta, PetscReal… in gqt() argument
216 PetscReal f = 0.0, p001 = 0.001, p5 = 0.5, minusone = -1, delta2 = delta * delta; in gqt()
262 parl = PetscMax(bnorm / delta - parl, pars); in gqt()
264 paru = PetscMax(0.0, bnorm / delta + paru); in gqt()
311 …if (PetscAbs(xnorm - delta) <= rtol * delta || (par == 0 && xnorm <= (1.0 + rtol) * delta)) info =… in gqt()
321 if (xnorm < delta) { in gqt()
323 PetscCallBLAS("BLASdot", prod = BLASdot_(&blasn, z, &blas1, x, &blas1) / delta); in gqt()
324 temp = (delta - xnorm) * ((delta + xnorm) / delta); in gqt()
325 alpha = temp / (PetscAbs(prod) + PetscSqrtScalar(prod * prod + temp / delta)); in gqt()
350 parc = (xnorm - delta) / (delta * temp * temp); in gqt()
[all …]
/petsc/src/snes/impls/ntrdc/ !
H A Dntrdc.c43 if (nrm >= neP->delta) { in SNESTRDC_KSPConverged_Private()
44 …fo(snes, "Ending linear iteration early, delta=%g, length=%g\n", (double)neP->delta, (double)nrm)); in SNESTRDC_KSPConverged_Private()
69 if (neP->delta < xnorm * neP->deltatol) { in SNESTRDC_Converged_Private()
70 …snes, "Diverged due to too small a trust region %g<%g*%g\n", (double)neP->delta, (double)xnorm, (d… in SNESTRDC_Converged_Private()
302 …PetscReal rho, fnorm, gnorm, xnorm = 0, delta, ynorm, temp_xnorm, temp_ynorm; /*… in SNESSolve_NEWTONTRDC() local
355 …delta = xnorm ? neP->delta0 * xnorm : neP->delta0; /* initial trust region size scaled by xno… in SNESSolve_NEWTONTRDC()
357 neP->delta = delta; in SNESSolve_NEWTONTRDC()
395 delta = neP->delta0 * xnorm; in SNESSolve_NEWTONTRDC()
397 delta = neP->delta * xnorm; in SNESSolve_NEWTONTRDC()
410 if (ynnorm <= delta) { /* see if the Newton solution is within the trust region */ in SNESSolve_NEWTONTRDC()
[all …]
/petsc/src/dm/impls/moab/ !
H A Ddmmbfem.cxx727 PetscReal delta[3] = {0.0, 0.0, 0.0}; in DMMoabPToRMapping() local
749 delta[2] = phypts[2] - xphy[2]; in DMMoabPToRMapping()
750 error += (delta[2] * delta[2]); in DMMoabPToRMapping()
752 delta[1] = phypts[1] - xphy[1]; in DMMoabPToRMapping()
753 error += (delta[1] * delta[1]); in DMMoabPToRMapping()
755 delta[0] = phypts[0] - xphy[0]; in DMMoabPToRMapping()
756 error += (delta[0] * delta[0]); in DMMoabPToRMapping()
766 natparam[0] -= ijacobian[0] * delta[0]; in DMMoabPToRMapping()
769 natparam[0] -= ijacobian[0] * delta[0] + ijacobian[1] * delta[1]; in DMMoabPToRMapping()
770 natparam[1] -= ijacobian[2] * delta[0] + ijacobian[3] * delta[1]; in DMMoabPToRMapping()
[all …]
/petsc/src/ksp/ksp/impls/bcgsl/ !
H A Dbcgsl.c65 if (bcgsl->delta > 0.0) { in KSPSolve_BCGSL()
91 if (bcgsl->delta > 0.0) PetscCall(VecAXPY(VX, 1.0, VXR)); in KSPSolve_BCGSL()
137 if (bcgsl->delta > 0.0) { in KSPSolve_BCGSL()
268 if (bcgsl->delta > 0.0) { in KSPSolve_BCGSL()
272 bUpdateX = (PetscBool)(zeta < bcgsl->delta * zeta0 && zeta0 <= rnmax_computed); in KSPSolve_BCGSL()
273 if ((zeta < bcgsl->delta * rnmax_true && zeta0 <= rnmax_true) || bUpdateX) { in KSPSolve_BCGSL()
288 if (bcgsl->delta > 0.0) PetscCall(VecAXPY(VX, 1.0, VXR)); in KSPSolve_BCGSL()
317 PetscErrorCode KSPBCGSLSetXRes(KSP ksp, PetscReal delta) in KSPBCGSLSetXRes() argument
322 PetscValidLogicalCollectiveReal(ksp, delta, 2); in KSPBCGSLSetXRes()
324 if ((delta <= 0 && bcgsl->delta > 0) || (delta > 0 && bcgsl->delta <= 0)) { in KSPBCGSLSetXRes()
[all …]
/petsc/src/ksp/ksp/impls/cg/pipelcg/ !
H A Dpipelcg.c10 #define delta(i) (plcg->delta[i]) macro
23 PetscScalar *gamma, *delta, *alpha; member
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()
228 delta(it - l) = G(it - l + 1, it - l + 1) / G(it - l, it - l); in KSPSolve_InnerLoop_PIPELCG()
230 … - l) * gamma(it - 2 * l) + G(it - l, it - l + 1) * delta(it - 2 * l) - G(it - l - 1, it - l) * de… in KSPSolve_InnerLoop_PIPELCG()
231 delta(it - l) = (G(it - l + 1, it - l + 1) * delta(it - 2 * l)) / G(it - l, it - l); in KSPSolve_InnerLoop_PIPELCG()
250 alpha(1) = -delta(it - l - 1); in KSPSolve_InnerLoop_PIPELCG()
253 PetscCall(VecScale(V[0], 1.0 / delta(it - l))); in KSPSolve_InnerLoop_PIPELCG()
271 alpha(1) = -delta(it - l - 1); in KSPSolve_InnerLoop_PIPELCG()
274 PetscCall(VecScale(Q[3 * j], 1.0 / delta(it - l))); in KSPSolve_InnerLoop_PIPELCG()
[all …]
/petsc/src/snes/impls/ms/ !
H A Dms.c14 PetscReal *delta; /* Coefficients of 3S* method */ member
56 …PetscReal delta[6] = {1.0000000000000000E+00, 5.3275427433201750E-01, 6.0143544663985238E-01, 4.… in SNESMSRegisterAll() local
58 PetscCall(SNESMSRegister(SNESMSM62, 6, 3, 1.0, &gamma[0][0], delta, betasub)); in SNESMSRegisterAll()
113 PetscCall(PetscFree(t->delta)); in SNESMSRegisterDestroy()
193 …nregisters, PetscReal stability, const PetscReal gamma[], const PetscReal delta[], const PetscReal… in SNESMSRegister() argument
201 if (gamma || delta) { in SNESMSRegister()
204 PetscAssertPointer(delta, 6); in SNESMSRegister()
218 if (gamma && delta) { in SNESMSRegister()
220 PetscCall(PetscMalloc1(nstages, &t->delta)); in SNESMSRegister()
222 PetscCall(PetscArraycpy(t->delta, delta, nstages)); in SNESMSRegister()
[all …]
/petsc/src/ts/tutorials/power_grid/stability_9bus/ !
H A Dex9busadj.c97 PetscErrorCode dq2ri(PetscScalar Fd, PetscScalar Fq, PetscScalar delta, PetscScalar *Fr, PetscScala… in dq2ri() argument
100 *Fr = Fd * PetscSinScalar(delta) + Fq * PetscCosScalar(delta); in dq2ri()
101 *Fi = -Fd * PetscCosScalar(delta) + Fq * PetscSinScalar(delta); in dq2ri()
106 PetscErrorCode ri2dq(PetscScalar Fr, PetscScalar Fi, PetscScalar delta, PetscScalar *Fd, PetscScala… in ri2dq() argument
109 *Fd = Fr * PetscSinScalar(delta) - Fi * PetscCosScalar(delta); in ri2dq()
110 *Fq = Fr * PetscCosScalar(delta) + Fi * PetscSinScalar(delta); in ri2dq()
120 PetscScalar Eqp, Edp, delta; in SetInitialGuess() local
150 delta = PetscAtan2Real(Vi + Xq[i] * IGr, Vr - Xq[i] * IGi); /* Machine angle */ in SetInitialGuess()
152 theta = PETSC_PI / 2.0 - delta; in SetInitialGuess()
168 xgen[idx + 2] = delta; in SetInitialGuess()
[all …]
H A Dex9busoptfd.c99 PetscErrorCode dq2ri(PetscScalar Fd, PetscScalar Fq, PetscScalar delta, PetscScalar *Fr, PetscScala… in dq2ri() argument
102 *Fr = Fd * PetscSinScalar(delta) + Fq * PetscCosScalar(delta); in dq2ri()
103 *Fi = -Fd * PetscCosScalar(delta) + Fq * PetscSinScalar(delta); in dq2ri()
108 PetscErrorCode ri2dq(PetscScalar Fr, PetscScalar Fi, PetscScalar delta, PetscScalar *Fd, PetscScala… in ri2dq() argument
111 *Fd = Fr * PetscSinScalar(delta) - Fi * PetscCosScalar(delta); in ri2dq()
112 *Fq = Fr * PetscCosScalar(delta) + Fi * PetscSinScalar(delta); in ri2dq()
122 PetscScalar Eqp, Edp, delta; in SetInitialGuess() local
152 delta = PetscAtan2Real(Vi + Xq[i] * IGr, Vr - Xq[i] * IGi); /* Machine angle */ in SetInitialGuess()
154 theta = PETSC_PI / 2.0 - delta; in SetInitialGuess()
170 xgen[idx + 2] = delta; in SetInitialGuess()
[all …]
H A Dex9bus.c297 PetscErrorCode dq2ri(PetscScalar Fd, PetscScalar Fq, PetscScalar delta, PetscScalar *Fr, PetscScala… in dq2ri() argument
300 *Fr = Fd * PetscSinScalar(delta) + Fq * PetscCosScalar(delta); in dq2ri()
301 *Fi = -Fd * PetscCosScalar(delta) + Fq * PetscSinScalar(delta); in dq2ri()
306 PetscErrorCode ri2dq(PetscScalar Fr, PetscScalar Fi, PetscScalar delta, PetscScalar *Fd, PetscScala… in ri2dq() argument
309 *Fd = Fr * PetscSinScalar(delta) - Fi * PetscCosScalar(delta); in ri2dq()
310 *Fq = Fr * PetscCosScalar(delta) + Fi * PetscSinScalar(delta); in ri2dq()
346 PetscScalar Eqp, Edp, delta; in SetInitialGuess() local
376 delta = PetscAtan2Real(Vi + Xq[i] * IGr, Vr - Xq[i] * IGi); /* Machine angle */ in SetInitialGuess()
378 theta = PETSC_PI / 2.0 - delta; in SetInitialGuess()
394 xgen[idx + 2] = delta; in SetInitialGuess()
[all …]
H A Dex9busopt.c111 PetscErrorCode dq2ri(PetscScalar Fd, PetscScalar Fq, PetscScalar delta, PetscScalar *Fr, PetscScala… in dq2ri() argument
114 *Fr = Fd * PetscSinScalar(delta) + Fq * PetscCosScalar(delta); in dq2ri()
115 *Fi = -Fd * PetscCosScalar(delta) + Fq * PetscSinScalar(delta); in dq2ri()
120 PetscErrorCode ri2dq(PetscScalar Fr, PetscScalar Fi, PetscScalar delta, PetscScalar *Fd, PetscScala… in ri2dq() argument
123 *Fd = Fr * PetscSinScalar(delta) - Fi * PetscCosScalar(delta); in ri2dq()
124 *Fq = Fr * PetscCosScalar(delta) + Fi * PetscSinScalar(delta); in ri2dq()
159 PetscScalar Eqp, Edp, delta; in SetInitialGuess() local
189 delta = PetscAtan2Real(Vi + Xq[i] * IGr, Vr - Xq[i] * IGi); /* Machine angle */ in SetInitialGuess()
191 theta = PETSC_PI / 2.0 - delta; in SetInitialGuess()
207 xgen[idx + 2] = delta; in SetInitialGuess()
[all …]
/petsc/src/ts/tutorials/hybrid/ !
H A Dex1fd.c197 PetscReal delta; in main() local
200 delta = 1e-8; in main()
277 u[0] = 0 + delta; in main()
284 sensi[0] = (u[0] - tmp[0]) / delta; in main()
285 sensi[1] = (u[1] - tmp[1]) / delta; in main()
288 u[1] = 1 + delta; in main()
294 sensi[0] = (u[0] - tmp[0]) / delta; in main()
295 sensi[1] = (u[1] - tmp[1]) / delta; in main()
299 app.lambda1 = app.lambda1 + delta; in main()
305 sensi[0] = (u[0] - tmp[0]) / delta; in main()
[all …]
/petsc/src/vec/vec/tutorials/ !
H A Dex42.m21 delta=zeros(512,1);
24 PetscBinaryWrite(socket,delta);
25 delta = PetscBinaryRead(socket); variable
/petsc/src/binding/petsc4py/test/ !
H A Dtest_ksp_py.py76 delta = delta_0
79 alpha = delta / d.dot(q)
82 delta_old = delta
83 delta = r.dot(r)
84 beta = delta / delta_old
/petsc/src/snes/impls/tr/ !
H A Dtr.c66 if (nrm >= neP->delta) { in SNESTR_KSPConverged_Private()
67 …ation early due to exiting trust region, delta=%g, length=%g\n", (double)neP->delta, (double)nrm)); in SNESTR_KSPConverged_Private()
92 if (neP->delta < neP->deltam) { in SNESTR_Converged_Private()
93 …fo(snes, "Diverged due to too small a trust region %g<%g\n", (double)neP->delta, (double)neP->delt… in SNESTR_Converged_Private()
490 PetscReal rho, fnorm, gnorm = 0.0, xnorm = 0.0, delta, ynorm; in SNESSolve_NEWTONTR() local
550 delta = neP->delta0; in SNESSolve_NEWTONTR()
551 neP->delta = delta; in SNESSolve_NEWTONTR()
634 auk = delta / gfnorm_k; in SNESSolve_NEWTONTR()
636 else tauk = PetscMin(gfnorm * gfnorm * gfnorm_k / (delta * gTBg), 1); in SNESSolve_NEWTONTR()
648 objmin = -neP->kmdc * gnorm * PetscMin(gnorm / beta_k, delta); in SNESSolve_NEWTONTR()
[all …]
/petsc/src/tao/bound/impls/tron/ !
H A Dtron.c73 PetscReal prered, actred, delta, f, f_new, rhok, gdx, xdiff, stepsize; in TaoSolve_TRON() local
121 delta = tao->trust; in TaoSolve_TRON()
134 PetscCall(TaoMonitor(tao, tao->niter, tron->f, tron->gnorm, 0.0, delta)); in TaoSolve_TRON()
156 PetscCall(KSPCGSetRadius(tao->ksp, delta)); in TaoSolve_TRON()
201 delta = PetscMin(xdiff, delta) * tron->sigma1; in TaoSolve_TRON()
203 delta = PetscMin(xdiff, delta) * tron->sigma3; in TaoSolve_TRON()
205 delta = PetscMin(xdiff, delta) * tron->sigma2; in TaoSolve_TRON()
215 } else if (delta <= 1e-30) { in TaoSolve_TRON()
218 delta /= 4.0; in TaoSolve_TRON()
224 tao->trust = delta; in TaoSolve_TRON()
/petsc/src/tao/bound/impls/bqnls/ !
H A Dbqnls.c9 PetscReal gnorm2, delta; in TaoBQNLSComputeHessian() local
15 if (bnk->f == 0.0) delta = 2.0 / gnorm2; in TaoBQNLSComputeHessian()
16 else delta = 2.0 * PetscAbsScalar(bnk->f) / gnorm2; in TaoBQNLSComputeHessian()
17 PetscCall(MatLMVMSymBroydenSetDelta(bqnk->B, delta)); in TaoBQNLSComputeHessian()
/petsc/src/ksp/ksp/tutorials/ !
H A Dexample100.py158 delta = r.dot(z)
162 alpha = delta / z.dot(p)
166 delta_old = delta
167 delta = r.dot(z)
168 beta = delta / delta_old
/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
81 PetscCall(VecDotBegin(M, W, &delta)); in KSPSolve_PIPECR()
88 PetscCall(VecDotEnd(M, W, &delta)); in KSPSolve_PIPECR()
100 alpha = gamma / delta; in KSPSolve_PIPECR()
106 alpha = gamma / (delta - beta / alpha * gamma); in KSPSolve_PIPECR()
/petsc/src/tao/complementarity/impls/ssls/ !
H A Dssils.c39 PetscReal delta, rho; in TaoSolve_SSILS() local
45 delta = ssls->delta; in TaoSolve_SSILS()
81 if (innerd <= delta * PetscPowReal(normd, rho)) { in TaoSolve_SSILS()
121 ssls->delta = 1e-10; in TaoCreate_SSILS()
H A Dssfls.c25 PetscReal delta, rho; in TaoSolve_SSFLS() local
31 delta = ssls->delta; in TaoSolve_SSFLS()
73 if (innerd >= -delta * PetscPowReal(normd, rho)) { in TaoSolve_SSFLS()
131 ssls->delta = 1e-10; in TaoCreate_SSFLS()
/petsc/src/ksp/ksp/utils/lmvm/rescale/ !
H A Dsymbrdnrescale.c352 PETSC_INTERN PetscErrorCode SymBroydenRescaleSetDelta(Mat B, SymBroydenRescale ldb, PetscReal delta) in SymBroydenRescaleSetDelta() argument
358 same = (delta == ldb->delta) ? PETSC_TRUE : PETSC_FALSE; in SymBroydenRescaleSetDelta()
359 ldb->delta = delta; in SymBroydenRescaleSetDelta()
360 ldb->delta = PetscMin(ldb->delta, ldb->delta_max); in SymBroydenRescaleSetDelta()
361 ldb->delta = PetscMax(ldb->delta, ldb->delta_min); in SymBroydenRescaleSetDelta()
378 mctx->delta = bctx->delta; in SymBroydenRescaleCopy()
475 PetscCall(VecSet(invD, ldb->delta)); in SymBroydenRescaleInitializeJ0()
480 PetscCall(MatLMVMSetJ0Scale(B, 1.0 / ldb->delta)); in SymBroydenRescaleInitializeJ0()
545 (*ldb)->delta = 1.0; in SymBroydenRescaleCreate()
/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
99 PetscCall(VecDotBegin(W, U, &delta)); in KSPSolve_PIPECG()
111 PetscCall(VecDotEnd(W, U, &delta)); in KSPSolve_PIPECG()
125 alpha = gamma / delta; in KSPSolve_PIPECG()
132 alpha = gamma / (delta - beta / alpha * gamma); in KSPSolve_PIPECG()

123