| /petsc/src/tao/unconstrained/impls/ntr/ |
| H A D | ntr.c | 50 PetscReal fmin, ftrial, prered, actred, kappa, sigma, beta; in TaoSolve_NTR() local
135 PetscCall(VecDot(tao->gradient, tao->stepdirection, &prered)); in TaoSolve_NTR()
137 prered = tao->trust * (gnorm - 0.5 * tao->trust * prered / (gnorm * gnorm)); in TaoSolve_NTR()
139 if ((PetscAbsScalar(actred) <= tr->epsilon) && (PetscAbsScalar(prered) <= tr->epsilon)) { in TaoSolve_NTR()
142 kappa = actred / prered; in TaoSolve_NTR()
145 … * gnorm * tao->trust / (tr->theta_i * gnorm * tao->trust + (1.0 - tr->theta_i) * prered - actred); in TaoSolve_NTR()
146 … * gnorm * tao->trust / (tr->theta_i * gnorm * tao->trust - (1.0 + tr->theta_i) * prered + actred); in TaoSolve_NTR()
289 PetscCall(KSPCGGetObjFcn(tao->ksp, &prered)); in TaoSolve_NTR()
290 if (prered >= 0.0) { in TaoSolve_NTR()
306 prered = -prered; in TaoSolve_NTR()
[all …]
|
| /petsc/src/tao/bound/impls/bnk/ |
| H A D | bntl.c | 112 PetscReal oldTrust, prered, actred, steplen, resnorm; in TaoSolve_BNTL() local
162 PetscCall(TaoBNKRecomputePred(tao, tao->stepdirection, &prered)); in TaoSolve_BNTL()
165 PetscCall(KSPCGGetObjFcn(tao->ksp, &prered)); in TaoSolve_BNTL()
167 prered = -prered; in TaoSolve_BNTL()
173 …PetscCall(TaoBNKUpdateTrustRadius(tao, prered, actred, bnk->update_type, stepType, &stepAccepted)); in TaoSolve_BNTL()
212 … PetscCall(TaoBNKUpdateTrustRadius(tao, prered, actred, BNK_UPDATE_STEP, stepType, &stepAccepted)); in TaoSolve_BNTL()
|
| H A D | bntr.c | 94 PetscReal oldTrust, prered, actred, steplen = 0.0, resnorm; in TaoSolve_BNTR() local
148 PetscCall(TaoBNKRecomputePred(tao, tao->stepdirection, &prered)); in TaoSolve_BNTR()
151 PetscCall(KSPCGGetObjFcn(tao->ksp, &prered)); in TaoSolve_BNTR()
153 prered = -prered; in TaoSolve_BNTR()
160 …PetscCall(TaoBNKUpdateTrustRadius(tao, prered, actred, bnk->update_type, stepType, &stepAccepted)); in TaoSolve_BNTR()
|
| H A D | bnk.c | 43 PetscReal f_min, ftrial, prered, actred, kappa, sigma, resnorm; in TaoBNKInitialize() local
174 PetscCall(VecDot(bnk->X_inactive, bnk->inactive_work, &prered)); in TaoBNKInitialize()
179 … prered = tao->trust * (bnk->gnorm - 0.5 * tao->trust * prered / (bnk->gnorm * bnk->gnorm)); in TaoBNKInitialize()
181 … if ((PetscAbsScalar(actred) <= bnk->epsilon) && (PetscAbsScalar(prered) <= bnk->epsilon)) { in TaoBNKInitialize()
184 kappa = actred / prered; in TaoBNKInitialize()
187 …m * tao->trust / (bnk->theta_i * bnk->gnorm * tao->trust + (1.0 - bnk->theta_i) * prered - actred); in TaoBNKInitialize()
188 …m * tao->trust / (bnk->theta_i * bnk->gnorm * tao->trust - (1.0 + bnk->theta_i) * prered + actred); in TaoBNKInitialize()
523 PetscErrorCode TaoBNKRecomputePred(Tao tao, Vec S, PetscReal *prered) in TaoBNKRecomputePred() argument
541 PetscCall(VecDot(bnk->inactive_work, bnk->X_inactive, prered)); in TaoBNKRecomputePred()
792 PetscErrorCode TaoBNKUpdateTrustRadius(Tao tao, PetscReal prered, PetscReal actred, PetscInt update… in TaoBNKUpdateTrustRadius() argument
[all …]
|
| /petsc/src/tao/unconstrained/impls/ntl/ |
| H A D | ntl.c | 37 PetscReal fmin, ftrial, prered, actred, kappa, sigma; in TaoSolve_NTL() local
129 PetscCall(VecDot(tao->gradient, tao->stepdirection, &prered)); in TaoSolve_NTL()
131 prered = tao->trust * (gnorm - 0.5 * tao->trust * prered / (gnorm * gnorm)); in TaoSolve_NTL()
133 if ((PetscAbsScalar(actred) <= tl->epsilon) && (PetscAbsScalar(prered) <= tl->epsilon)) { in TaoSolve_NTL()
136 kappa = actred / prered; in TaoSolve_NTL()
139 … * gnorm * tao->trust / (tl->theta_i * gnorm * tao->trust + (1.0 - tl->theta_i) * prered - actred); in TaoSolve_NTL()
140 … * gnorm * tao->trust / (tl->theta_i * gnorm * tao->trust - (1.0 + tl->theta_i) * prered + actred); in TaoSolve_NTL()
292 PetscCall(KSPCGGetObjFcn(tao->ksp, &prered)); in TaoSolve_NTL()
293 if (prered >= 0.0) { in TaoSolve_NTL()
311 prered = -prered; in TaoSolve_NTL()
[all …]
|
| /petsc/src/tao/unconstrained/impls/nls/ |
| H A D | nls.c | 47 PetscReal fmin, ftrial, f_full, prered, actred, kappa, sigma; in TaoSolve_NLS() local
155 PetscCall(VecDot(tao->gradient, nlsP->D, &prered)); in TaoSolve_NLS()
157 prered = tao->trust * (gnorm - 0.5 * tao->trust * prered / (gnorm * gnorm)); in TaoSolve_NLS()
159 … if ((PetscAbsScalar(actred) <= nlsP->epsilon) && (PetscAbsScalar(prered) <= nlsP->epsilon)) { in TaoSolve_NLS()
162 kappa = actred / prered; in TaoSolve_NLS()
165 …norm * tao->trust / (nlsP->theta_i * gnorm * tao->trust + (1.0 - nlsP->theta_i) * prered - actred); in TaoSolve_NLS()
166 …norm * tao->trust / (nlsP->theta_i * gnorm * tao->trust - (1.0 + nlsP->theta_i) * prered + actred); in TaoSolve_NLS()
563 PetscCall(KSPCGGetObjFcn(tao->ksp, &prered)); in TaoSolve_NLS()
564 if (prered >= 0.0) { in TaoSolve_NLS()
575 prered = -prered; in TaoSolve_NLS()
[all …]
|
| /petsc/src/tao/linesearch/impls/gpcglinesearch/ |
| H A D | gpcglinesearch.c | 32 PetscReal d1, finit, actred, prered, rho, gdx; in TaoLineSearchApply_GPCG() local
124 PetscCall(VecDot(neP->W1, neP->Gold, &prered)); in TaoLineSearchApply_GPCG()
126 if (PetscAbsReal(prered) < 1.0e-100) prered = 1.0e-12; in TaoLineSearchApply_GPCG()
127 rho = actred / prered; in TaoLineSearchApply_GPCG()
|
| /petsc/src/tao/bound/impls/tron/ |
| H A D | tron.c | 73 PetscReal prered, actred, delta, f, f_new, rhok, gdx, xdiff, stepsize; in TaoSolve_TRON() local
180 PetscCall(VecDot(tao->stepdirection, tron->Work, &prered)); in TaoSolve_TRON()
182 if ((PetscAbsScalar(actred) <= 1e-6) && (PetscAbsScalar(prered) <= 1e-6)) { in TaoSolve_TRON()
185 rhok = PetscAbs(-actred / prered); in TaoSolve_TRON()
|