Home
last modified time | relevance | path

Searched refs:gnorm (Results 1 – 25 of 54) sorted by relevance

123

/petsc/src/snes/tests/output/
H A Dex241_1.out6 Line search: Using full step: fnorm 9.999000000000e-01 gnorm 5.624375006250e-01
10 Line search: Using full step: fnorm 5.624375006250e-01 gnorm 1.405843796004e-01
14 Line search: Using full step: fnorm 1.405843796004e-01 gnorm 3.512111267037e-02
18 Line search: Using full step: fnorm 3.512111267037e-02 gnorm 8.755349147747e-03
22 Line search: Using full step: fnorm 8.755349147747e-03 gnorm 2.164119602175e-03
26 Line search: Using full step: fnorm 2.164119602175e-03 gnorm 5.171340824948e-04
30 Line search: Using full step: fnorm 5.171340824948e-04 gnorm 1.083345041472e-04
34 Line search: Using full step: fnorm 1.083345041472e-04 gnorm 1.408355859831e-05
38 Line search: Using full step: fnorm 1.408355859831e-05 gnorm 4.346520770148e-07
42 Line search: Using full step: fnorm 4.346520770148e-07 gnorm 4.702620662805e-10
/petsc/src/snes/linesearch/impls/bt/
H A Dlinesearchbt.c63 PetscReal fnorm, xnorm, ynorm, gnorm; in SNESLineSearchApply_BT() local
145 gnorm = fnorm; in SNESLineSearchApply_BT()
146 PetscCall((*linesearch->ops->vinorm)(snes, G, W, &gnorm)); in SNESLineSearchApply_BT()
148 PetscCall(VecNorm(G, NORM_2, &gnorm)); in SNESLineSearchApply_BT()
150 g = 0.5 * PetscSqr(gnorm); in SNESLineSearchApply_BT()
164 … PetscCall(PetscInfo(snes, "Initial fnorm %14.12e gnorm %14.12e\n", (double)fnorm, (double)gnorm)); in SNESLineSearchApply_BT()
169 … " Line search: Using full step: fnorm %14.12e gnorm %14.12e\n", (double)fnorm, (double)gnorm)); in SNESLineSearchApply_BT()
209 gnorm = fnorm; in SNESLineSearchApply_BT()
210 PetscCall((*linesearch->ops->vinorm)(snes, G, W, &gnorm)); in SNESLineSearchApply_BT()
212 PetscCall(VecNorm(G, NORM_2, &gnorm)); in SNESLineSearchApply_BT()
[all …]
/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
25 PetscCall(VecNorm(tao->gradient, NORM_2, &gnorm)); in TaoSolve_CG()
26 …PetscCheck(!PetscIsInfOrNanReal(f) && !PetscIsInfOrNanReal(gnorm), PetscObjectComm((PetscObject)ta… in TaoSolve_CG()
29 PetscCall(TaoLogConvergenceHistory(tao, f, gnorm, 0.0, tao->ksp_its)); in TaoSolve_CG()
30 PetscCall(TaoMonitor(tao, tao->niter, f, gnorm, 0.0, step)); in TaoSolve_CG()
37 gnorm2 = gnorm * gnorm; in TaoSolve_CG()
136 PetscCall(VecNorm(tao->gradient, NORM_2, &gnorm)); in TaoSolve_CG()
137 …PetscCheck(!PetscIsInfOrNanReal(f) && !PetscIsInfOrNanReal(gnorm), PetscObjectComm((PetscObject)ta… in TaoSolve_CG()
140 gnorm2 = gnorm * gnorm; in TaoSolve_CG()
142 PetscCall(TaoLogConvergenceHistory(tao, f, gnorm, 0.0, tao->ksp_its)); in TaoSolve_CG()
[all …]
/petsc/src/tao/unconstrained/impls/ntr/
H A Dntr.c52 PetscReal f, gnorm; in TaoSolve_NTR() local
92 PetscCall(TaoGradientNorm(tao, tao->gradient, NORM_2, &gnorm)); in TaoSolve_NTR()
93 …PetscCheck(!PetscIsInfOrNanReal(f) && !PetscIsInfOrNanReal(gnorm), PetscObjectComm((PetscObject)ta… in TaoSolve_NTR()
97 PetscCall(TaoLogConvergenceHistory(tao, f, gnorm, 0.0, tao->ksp_its)); in TaoSolve_NTR()
98 PetscCall(TaoMonitor(tao, tao->niter, f, gnorm, 0.0, 1.0)); in TaoSolve_NTR()
123 PetscCall(VecAXPY(tr->W, -tao->trust / gnorm, tao->gradient)); in TaoSolve_NTR()
131 sigma = -tao->trust / gnorm; in TaoSolve_NTR()
137 prered = tao->trust * (gnorm - 0.5 * tao->trust * prered / (gnorm * gnorm)); in TaoSolve_NTR()
145 …tau_1 = tr->theta_i * gnorm * tao->trust / (tr->theta_i * gnorm * tao->trust + (1.0 - tr->theta_… in TaoSolve_NTR()
146 …tau_2 = tr->theta_i * gnorm * tao->trust / (tr->theta_i * gnorm * tao->trust - (1.0 + tr->theta_… in TaoSolve_NTR()
[all …]
/petsc/src/tao/unconstrained/impls/owlqn/
H A Dowlqn.c55 PetscReal f, fold, gdx, gnorm; in TaoSolve_OWLQN() local
68 PetscCall(VecNorm(lmP->GV, NORM_2, &gnorm)); in TaoSolve_OWLQN()
69 …PetscCheck(!PetscIsInfOrNanReal(f) && !PetscIsInfOrNanReal(gnorm), PetscObjectComm((PetscObject)ta… in TaoSolve_OWLQN()
72 PetscCall(TaoLogConvergenceHistory(tao, f, gnorm, 0.0, tao->ksp_its)); in TaoSolve_OWLQN()
73 PetscCall(TaoMonitor(tao, tao->niter, f, gnorm, 0.0, step)); in TaoSolve_OWLQN()
78 delta = 2.0 * PetscMax(1.0, PetscAbsScalar(f)) / (gnorm * gnorm); in TaoSolve_OWLQN()
110 delta = 2.0 * PetscMax(1.0, PetscAbsScalar(f)) / (gnorm * gnorm); in TaoSolve_OWLQN()
156 delta = 2.0 * PetscMax(1.0, PetscAbsScalar(f)) / (gnorm * gnorm); in TaoSolve_OWLQN()
208 PetscCall(VecNorm(lmP->GV, NORM_2, &gnorm)); in TaoSolve_OWLQN()
211 PetscCall(TaoLogConvergenceHistory(tao, f, gnorm, 0.0, tao->ksp_its)); in TaoSolve_OWLQN()
[all …]
/petsc/src/snes/linesearch/impls/secant/
H A Dlinesearchsecant.c13 PetscReal gnorm; in SNESLineSearchApply_Secant() local
26 PetscCall(SNESLineSearchGetNorms(linesearch, &xnorm, &gnorm, &ynorm)); in SNESLineSearchApply_Secant()
39 fnrm_old = gnorm * gnorm; in SNESLineSearchApply_Secant()
55 fnrm_mid = gnorm; in SNESLineSearchApply_Secant()
66 fnrm = gnorm; in SNESLineSearchApply_Secant()
208 PetscCall(SNESLineSearchGetNorms(linesearch, NULL, &gnorm, NULL)); in SNESLineSearchApply_Secant()
209 SNESLineSearchCheckFunctionDomainError(snes, linesearch, gnorm); in SNESLineSearchApply_Secant()
212 …(monitor, " Line search terminated: lambda = %g, fnorm = %g\n", (double)lambda, (double)gnorm)); in SNESLineSearchApply_Secant()
/petsc/src/tao/bound/impls/tron/
H A Dtron.c87 PetscCall(VecNorm(tao->gradient, NORM_2, &tron->gnorm)); in TaoSolve_TRON()
88 …PetscCheck(!PetscIsInfOrNanReal(tron->f) && !PetscIsInfOrNanReal(tron->gnorm), PetscObjectComm((Pe… in TaoSolve_TRON()
92 PetscCall(VecNorm(tao->gradient, NORM_2, &tron->gnorm)); in TaoSolve_TRON()
96 if (tao->trust <= 0) tao->trust = PetscMax(tron->gnorm * tron->gnorm, 1.0); in TaoSolve_TRON()
103 PetscCall(TaoLogConvergenceHistory(tao, tron->f, tron->gnorm, 0.0, tao->ksp_its)); in TaoSolve_TRON()
104 PetscCall(TaoMonitor(tao, tao->niter, tron->f, tron->gnorm, 0.0, tron->stepsize)); in TaoSolve_TRON()
117 PetscCall(VecNorm(tao->gradient, NORM_2, &tron->gnorm)); in TaoSolve_TRON()
132 PetscCall(VecNorm(tao->gradient, NORM_2, &tron->gnorm)); in TaoSolve_TRON()
133 PetscCall(TaoLogConvergenceHistory(tao, tron->f, tron->gnorm, 0.0, tao->ksp_its)); in TaoSolve_TRON()
134 PetscCall(TaoMonitor(tao, tao->niter, tron->f, tron->gnorm, 0.0, delta)); in TaoSolve_TRON()
[all …]
H A Dtron.h25 PetscReal gnorm; member
/petsc/src/tao/unconstrained/impls/nls/
H A Dnls.c49 PetscReal f, fold, gdx, gnorm, pert; in TaoSolve_NLS() local
97 PetscCall(TaoGradientNorm(tao, tao->gradient, NORM_2, &gnorm)); in TaoSolve_NLS()
98 …PetscCheck(!PetscIsInfOrNanReal(f) && !PetscIsInfOrNanReal(gnorm), PetscObjectComm((PetscObject)ta… in TaoSolve_NLS()
101 PetscCall(TaoLogConvergenceHistory(tao, f, gnorm, 0.0, tao->ksp_its)); in TaoSolve_NLS()
102 PetscCall(TaoMonitor(tao, tao->niter, f, gnorm, 0.0, step)); in TaoSolve_NLS()
144 PetscCall(VecAXPY(nlsP->W, -tao->trust / gnorm, tao->gradient)); in TaoSolve_NLS()
151 sigma = -tao->trust / gnorm; in TaoSolve_NLS()
157 prered = tao->trust * (gnorm - 0.5 * tao->trust * prered / (gnorm * gnorm)); in TaoSolve_NLS()
165 …tau_1 = nlsP->theta_i * gnorm * tao->trust / (nlsP->theta_i * gnorm * tao->trust + (1.0 - nlsP->… in TaoSolve_NLS()
166 …tau_2 = nlsP->theta_i * gnorm * tao->trust / (nlsP->theta_i * gnorm * tao->trust - (1.0 + nlsP->… in TaoSolve_NLS()
[all …]
/petsc/src/tao/quadratic/impls/gpcg/
H A Dgpcg.c111 PetscReal actred, f, f_new, gnorm, gdx, stepsize, xtb; in TaoSolve_GPCG() local
135 PetscCall(VecNorm(gpcg->PG, NORM_2, &gpcg->gnorm)); in TaoSolve_GPCG()
141 PetscCall(TaoLogConvergenceHistory(tao, f, gpcg->gnorm, 0.0, tao->ksp_its)); in TaoSolve_GPCG()
142 PetscCall(TaoMonitor(tao, tao->niter, f, gpcg->gnorm, 0.0, tao->step)); in TaoSolve_GPCG()
154 gnorm = gpcg->gnorm; in TaoSolve_GPCG()
196 PetscCall(VecNorm(gpcg->PG, NORM_2, &gnorm)); in TaoSolve_GPCG()
208 gpcg->gnorm = gnorm; in TaoSolve_GPCG()
210 PetscCall(TaoLogConvergenceHistory(tao, f, gpcg->gnorm, 0.0, tao->ksp_its)); in TaoSolve_GPCG()
211 PetscCall(TaoMonitor(tao, tao->niter, f, gpcg->gnorm, 0.0, tao->step)); in TaoSolve_GPCG()
222 PetscReal actred = -1.0, actred_max = 0.0, gAg, gtg = gpcg->gnorm, alpha; in GPCGGradProjections()
[all …]
H A Dgpcg.h12 PetscReal gnorm; member
/petsc/src/tao/bound/impls/blmvm/
H A Dblmvm.c9 PetscReal f, fold, gdx, gnorm, gnorm2; in TaoSolve_BLMVM() local
22 PetscCall(TaoGradientNorm(tao, tao->gradient, NORM_2, &gnorm)); in TaoSolve_BLMVM()
23 …PetscCheck(!PetscIsInfOrNanReal(f) && !PetscIsInfOrNanReal(gnorm), PetscObjectComm((PetscObject)ta… in TaoSolve_BLMVM()
26 PetscCall(TaoLogConvergenceHistory(tao, f, gnorm, 0.0, tao->ksp_its)); in TaoSolve_BLMVM()
27 PetscCall(TaoMonitor(tao, tao->niter, f, gnorm, 0.0, stepsize)); in TaoSolve_BLMVM()
46 gnorm2 = gnorm * gnorm; in TaoSolve_BLMVM()
107 PetscCall(TaoGradientNorm(tao, tao->gradient, NORM_2, &gnorm)); in TaoSolve_BLMVM()
108 …PetscCheck(!PetscIsInfOrNanReal(f) && !PetscIsInfOrNanReal(gnorm), PetscObjectComm((PetscObject)ta… in TaoSolve_BLMVM()
110 PetscCall(TaoLogConvergenceHistory(tao, f, gnorm, 0.0, tao->ksp_its)); in TaoSolve_BLMVM()
111 PetscCall(TaoMonitor(tao, tao->niter, f, gnorm, 0.0, stepsize)); in TaoSolve_BLMVM()
/petsc/src/tao/unconstrained/impls/ntl/
H A Dntl.c39 PetscReal f, fold, gdx, gnorm; in TaoSolve_NTL() local
86 PetscCall(VecNorm(tao->gradient, NORM_2, &gnorm)); in TaoSolve_NTL()
87 …PetscCheck(!PetscIsInfOrNanReal(f) && !PetscIsInfOrNanReal(gnorm), PetscObjectComm((PetscObject)ta… in TaoSolve_NTL()
91 PetscCall(TaoLogConvergenceHistory(tao, f, gnorm, 0.0, tao->ksp_its)); in TaoSolve_NTL()
92 PetscCall(TaoMonitor(tao, tao->niter, f, gnorm, 0.0, step)); in TaoSolve_NTL()
117 PetscCall(VecAXPY(tl->W, -tao->trust / gnorm, tao->gradient)); in TaoSolve_NTL()
125 sigma = -tao->trust / gnorm; in TaoSolve_NTL()
131 prered = tao->trust * (gnorm - 0.5 * tao->trust * prered / (gnorm * gnorm)); in TaoSolve_NTL()
139 …tau_1 = tl->theta_i * gnorm * tao->trust / (tl->theta_i * gnorm * tao->trust + (1.0 - tl->theta_… in TaoSolve_NTL()
140 …tau_2 = tl->theta_i * gnorm * tao->trust / (tl->theta_i * gnorm * tao->trust - (1.0 + tl->theta_… in TaoSolve_NTL()
[all …]
/petsc/src/tao/unconstrained/impls/lmvm/
H A Dlmvm.c10 PetscReal f, fold, gdx, gnorm; in TaoSolve_LMVM() local
20 PetscCall(TaoGradientNorm(tao, tao->gradient, NORM_2, &gnorm)); in TaoSolve_LMVM()
22 …PetscCheck(!PetscIsInfOrNanReal(f) && !PetscIsInfOrNanReal(gnorm), PetscObjectComm((PetscObject)ta… in TaoSolve_LMVM()
25 PetscCall(TaoLogConvergenceHistory(tao, f, gnorm, 0.0, tao->ksp_its)); in TaoSolve_LMVM()
26 PetscCall(TaoMonitor(tao, tao->niter, f, gnorm, 0.0, step)); in TaoSolve_LMVM()
129 PetscCall(TaoGradientNorm(tao, tao->gradient, NORM_2, &gnorm)); in TaoSolve_LMVM()
134 PetscCall(TaoLogConvergenceHistory(tao, f, gnorm, 0.0, tao->ksp_its)); in TaoSolve_LMVM()
135 PetscCall(TaoMonitor(tao, tao->niter, f, gnorm, 0.0, step)); in TaoSolve_LMVM()
/petsc/src/snes/linesearch/impls/basic/
H A Dlinesearchbasic.c9 PetscReal gnorm, xnorm, ynorm, lambda, fnorm = 0.0; in SNESLineSearchApply_Basic() local
13 PetscCall(SNESLineSearchGetNorms(linesearch, &xnorm, &gnorm, &ynorm)); in SNESLineSearchApply_Basic()
43 linesearch->fnorm = gnorm; in SNESLineSearchApply_Basic()
/petsc/src/snes/linesearch/impls/nleqerr/
H A Dlinesearchnleqerr.c38 PetscReal fnorm, xnorm, ynorm, gnorm, wnorm; in SNESLineSearchApply_NLEQERR() local
100 PetscCall(VecNormBegin(G, NORM_2, &gnorm)); in SNESLineSearchApply_NLEQERR()
101 PetscCall(VecNormEnd(G, NORM_2, &gnorm)); in SNESLineSearchApply_NLEQERR()
103 …err->lambda_prev * (nleqerr->norm_delta_x_prev * nleqerr->norm_bar_delta_x_prev) / (gnorm * ynorm); in SNESLineSearchApply_NLEQERR()
186 PetscCall(VecNorm(G, NORM_2, &gnorm)); in SNESLineSearchApply_NLEQERR()
188 mudash = (0.5 * ynorm * lambda * lambda) / gnorm; in SNESLineSearchApply_NLEQERR()
/petsc/src/snes/impls/vi/ss/
H A Dviss.c200 PetscReal gnorm, xnorm = 0, ynorm; in SNESSolve_VINEWTONSSLS() local
301 gnorm = vi->phinorm; in SNESSolve_VINEWTONSSLS()
302 PetscCall(SNESLineSearchApply(snes->linesearch, X, vi->phi, &gnorm, Y)); in SNESSolve_VINEWTONSSLS()
304 PetscCall(SNESLineSearchGetNorms(snes->linesearch, &xnorm, &gnorm, &ynorm)); in SNESSolve_VINEWTONSSLS()
305 …orm=%18.16e, ynorm=%18.16e, lssucceed=%d\n", (double)vi->phinorm, (double)gnorm, (double)ynorm, (i… in SNESSolve_VINEWTONSSLS()
317 PetscCall(SNESVICheckLocalMin_Private(snes, snes->jacobian, vi->phi, X, gnorm, &ismin)); in SNESSolve_VINEWTONSSLS()
323 vi->phinorm = gnorm; in SNESSolve_VINEWTONSSLS()
/petsc/src/tao/interface/
H A Dtaosolver.c1724 PetscReal gnorm; in TaoMonitorDefaultShort() local
1730 gnorm = tao->residual; in TaoMonitorDefaultShort()
1738 if (gnorm >= PETSC_INFINITY) { in TaoMonitorDefaultShort()
1740 } else if (gnorm > 1.e-6) { in TaoMonitorDefaultShort()
1741 PetscCall(PetscViewerASCIIPrintf(viewer, " Residual: %g \n", (double)gnorm)); in TaoMonitorDefaultShort()
1742 } else if (gnorm > 1.e-11) { in TaoMonitorDefaultShort()
2013 PetscReal gnorm = tao->residual, gnorm0 = tao->gnorm0; in TaoDefaultConvergenceTest() local
2030 } else if (gnorm <= gatol && cnorm <= catol) { in TaoDefaultConvergenceTest()
2031 …PetscCall(PetscInfo(tao, "Converged due to residual norm ||g(X)||=%g < %g\n", (double)gnorm, (doub… in TaoDefaultConvergenceTest()
2033 } else if (f != 0 && PetscAbsReal(gnorm / f) <= grtol && cnorm <= crtol) { in TaoDefaultConvergenceTest()
[all …]
/petsc/src/snes/linesearch/impls/cp/
H A Dlinesearchcp.c9 PetscReal xnorm, ynorm, gnorm, minlambda, maxlambda, rtol, atol, ltol; in SNESLineSearchApply_CP() local
17 PetscCall(SNESLineSearchGetNorms(linesearch, &xnorm, &gnorm, &ynorm)); in SNESLineSearchApply_CP()
190 PetscCall(SNESLineSearchGetNorms(linesearch, &xnorm, &gnorm, &ynorm)); in SNESLineSearchApply_CP()
194 …monitor, " Line search terminated: lambda = %g, fnorms = %g\n", (double)lambda, (double)gnorm)); in SNESLineSearchApply_CP()
/petsc/src/snes/tutorials/
H A Dex1f.F90122 PetscReal ynorm, gnorm, xnorm
133 PetscCall(VecNorm(f, NORM_2, gnorm, ierr))
136 PetscCall(SNESLineSearchSetNorms(linesearch, xnorm, gnorm, ynorm, ierr))
/petsc/src/tao/bound/impls/bqnls/
H A Dbqnls.c13 gnorm2 = bnk->gnorm * bnk->gnorm; in TaoBQNLSComputeHessian()
/petsc/src/tao/bound/impls/bnk/
H A Dbnk.c66 PetscCall(TaoGradientNorm(tao, tao->gradient, NORM_2, &bnk->gnorm)); in TaoBNKInitialize()
148 PetscCall(VecAXPY(tao->solution, -tao->trust / bnk->gnorm, tao->gradient)); in TaoBNKInitialize()
162 sigma = -tao->trust / bnk->gnorm; in TaoBNKInitialize()
179 … prered = tao->trust * (bnk->gnorm - 0.5 * tao->trust * prered / (bnk->gnorm * bnk->gnorm)); in TaoBNKInitialize()
187 …tau_1 = bnk->theta_i * bnk->gnorm * tao->trust / (bnk->theta_i * bnk->gnorm * tao->trust + (1.0 … in TaoBNKInitialize()
188 …tau_2 = bnk->theta_i * bnk->gnorm * tao->trust / (bnk->theta_i * bnk->gnorm * tao->trust - (1.0 … in TaoBNKInitialize()
247 PetscCall(TaoGradientNorm(tao, tao->gradient, NORM_2, &bnk->gnorm)); in TaoBNKInitialize()
574 bnk->pert = PetscMin(bnk->imax, PetscMax(bnk->imin, bnk->imfac * bnk->gnorm)); in TaoBNKSafeguardStep()
582 bnk->pert = PetscMin(bnk->pmax, PetscMax(bnk->pgfac * bnk->pert, bnk->pmgfac * bnk->gnorm)); in TaoBNKSafeguardStep()
636 bnk->pert = PetscMin(bnk->imax, PetscMax(bnk->imin, bnk->imfac * bnk->gnorm)); in TaoBNKSafeguardStep()
[all …]
/petsc/src/ts/tutorials/power_grid/
H A Dex3opt_fd.c26 PetscReal f, gnorm, cnorm, xdiff; in monitor() local
32 PetscCall(TaoGetSolutionStatus(tao, &iterate, &f, &gnorm, &cnorm, &xdiff, &reason)); in monitor()
38 …_COMM_WORLD, fp, "%" PetscInt_FMT " %g %.12lf %.12lf\n", iterate, (double)gnorm, (double)PetscReal… in monitor()
/petsc/src/snes/tutorials/output/
H A Dex21_1.out13 Line search: Using full step: fnorm 8.660254037844e-01 gnorm 4.873279917558e-10
/petsc/src/ts/tutorials/
H A Dex24.c102 PetscReal dt, gnorm; in MonitorObjective() local
119 PetscCall(VecNorm(F, NORM_2, &gnorm)); in MonitorObjective()
127 …t_FMT ")\n", step, (double)t, (double)dt, (double)PetscRealPart(f), (double)gnorm, snesit, linit)); in MonitorObjective()

123