Home
last modified time | relevance | path

Searched refs:xnorm (Results 1 – 25 of 57) sorted by relevance

123

/petsc/src/mat/tests/
H A Dex185.c11 PetscReal xnorm, ynorm, anorm; in main() local
21 PetscCall(VecNorm(X, NORM_2, &xnorm)); in main()
24 …cAbsReal(ynorm - 3 * xnorm) <= PETSC_SMALL, PETSC_COMM_WORLD, PETSC_ERR_PLIB, "Expected norm %g ac… in main()
52 …scAbsReal(ynorm - xnorm / 4) <= PETSC_SMALL, PETSC_COMM_WORLD, PETSC_ERR_PLIB, "Expected norm %g a… in main()
57 …scAbsReal(ynorm - xnorm / 4) <= PETSC_SMALL, PETSC_COMM_WORLD, PETSC_ERR_PLIB, "Expected norm %g a… in main()
62 …scAbsReal(ynorm - xnorm / 4) <= PETSC_SMALL, PETSC_COMM_WORLD, PETSC_ERR_PLIB, "Expected norm %g a… in main()
/petsc/src/tao/leastsquares/impls/pounders/
H A Dgqt.c220 …anorm, bnorm, parc, parf, parl, pars, par = *retpar, paru, prod, rxnorm, rznorm = 0.0, temp, xnorm; in gqt() local
227 xnorm = 0.0; in gqt()
307 PetscCallBLAS("BLASnrm2", xnorm = BLASnrm2_(&blasn, x, &blas1)); 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()
324 temp = (delta - xnorm) * ((delta + xnorm) / delta); in gqt()
331 if ((rznorm * rznorm + par * xnorm * xnorm) / (delta2) <= par) rednc = 1; in gqt()
341 if (xnorm == 0) { in gqt()
345 temp = 1.0 / xnorm; in gqt()
350 parc = (xnorm - delta) / (delta * temp * temp); in gqt()
[all …]
/petsc/src/ksp/ksp/tests/output/
H A Dex14_1.out2 linear solve iterations = 2, xnorm=1.21872, ynorm=0.228979
4 linear solve iterations = 2, xnorm=1.23797, ynorm=0.0192507
6 linear solve iterations = 2, xnorm=1.23812, ynorm=0.000149569
H A Dex14_1_alt.out2 linear solve iterations = 2, xnorm=1.21872, ynorm=0.228979
4 linear solve iterations = 2, xnorm=1.23797, ynorm=0.0192509
6 linear solve iterations = 2, xnorm=1.23812, ynorm=0.000149565
H A Dex13_qlp_sisc_none.out2 flag rnorm Arnorm Compatible LS Anorm Acond xnorm
14 flag rnorm Arnorm Compatible LS Anorm Acond xnorm
H A Dex13_qlp_sisc.out2 flag rnorm Arnorm Compatible LS Anorm Acond xnorm
14 flag rnorm Arnorm Compatible LS Anorm Acond xnorm
H A Dex42_minres_qlp_nonconsistent.out1 flag rnorm Arnorm Compatible LS Anorm Acond xnorm
14 flag rnorm Arnorm Compatible LS Anorm Acond xnorm
H A Dex42_minres_qlp.out1 flag rnorm Arnorm Compatible LS Anorm Acond xnorm
14 flag rnorm Arnorm Compatible LS Anorm Acond xnorm
/petsc/src/snes/impls/ntrdc/
H A Dntrdc.c63 static PetscErrorCode SNESTRDC_Converged_Private(SNES snes, PetscInt it, PetscReal xnorm, PetscReal… in SNESTRDC_Converged_Private() argument
69 if (neP->delta < xnorm * neP->deltatol) { in SNESTRDC_Converged_Private()
70 … due to too small a trust region %g<%g*%g\n", (double)neP->delta, (double)xnorm, (double)neP->delt… in SNESTRDC_Converged_Private()
302 …PetscReal rho, fnorm, gnorm, xnorm = 0, delta, ynorm, temp_xnorm, temp_ynorm; /*… in SNESSolve_NEWTONTRDC() local
350 PetscCall(VecNorm(X, NORM_2, &xnorm)); /* xnorm <- || X || */ in SNESSolve_NEWTONTRDC()
355 …delta = xnorm ? neP->delta0 * xnorm : neP->delta0; /* initial trust region size scaled by xno… in SNESSolve_NEWTONTRDC()
356 …deltaM = xnorm ? neP->deltaM * xnorm : neP->deltaM; /* maximum trust region size scaled by xno… in SNESSolve_NEWTONTRDC()
393 PetscCall(VecNorm(X, NORM_2, &xnorm)); in SNESSolve_NEWTONTRDC()
395 delta = neP->delta0 * xnorm; in SNESSolve_NEWTONTRDC()
397 delta = neP->delta * xnorm; in SNESSolve_NEWTONTRDC()
[all …]
/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()
38 PetscCall(VecNormBegin(W, NORM_2, &linesearch->xnorm)); in SNESLineSearchApply_Basic()
40 PetscCall(VecNormEnd(W, NORM_2, &linesearch->xnorm)); in SNESLineSearchApply_Basic()
/petsc/src/snes/impls/ngmres/
H A Danderson.c31 PetscReal xnorm, ynorm; in SNESSolve_Anderson() local
122 …e(snes, l, X, F, XM, FM, XA, FA, D, &dnorm, &dminnorm, NULL, NULL, NULL, &xnorm, &fAnorm, &ynorm)); in SNESSolve_Anderson()
128 …_Private(snes, l, X, F, XM, FM, XA, FA, D, NULL, NULL, NULL, NULL, NULL, &xnorm, &fAnorm, &ynorm)); in SNESSolve_Anderson()
134 …_Private(snes, l, X, F, XM, FM, XA, FA, D, NULL, NULL, NULL, NULL, NULL, &xnorm, &fAnorm, &ynorm)); in SNESSolve_Anderson()
158 snes->xnorm = xnorm; in SNESSolve_Anderson()
162 PetscCall(SNESConverged(snes, snes->iter, xnorm, ynorm, fnorm)); in SNESSolve_Anderson()
H A Dngmresfunc.c163 …rm, PetscReal objmin, PetscReal dminnorm, Vec X, Vec F, Vec Y, PetscReal *xnorm, PetscReal *fnorm,… in SNESNGMRESSelect_Private() argument
181 *xnorm = xMnorm; in SNESNGMRESSelect_Private()
206 PetscCall(SNESLineSearchGetNorms(ngmres->additive_linesearch, xnorm, fnorm, ynorm)); in SNESNGMRESSelect_Private()
226 *xnorm = xAnorm; in SNESNGMRESSelect_Private()
233 *xnorm = xMnorm; in SNESNGMRESSelect_Private()
241 *xnorm = xAnorm; in SNESNGMRESSelect_Private()
/petsc/src/snes/linesearch/impls/nleqerr/
H A Dlinesearchnleqerr.c38 PetscReal fnorm, xnorm, ynorm, gnorm, wnorm; in SNESLineSearchApply_NLEQERR() local
50 PetscCall(SNESLineSearchGetNorms(linesearch, &xnorm, &fnorm, &ynorm)); in SNESLineSearchApply_NLEQERR()
66 PetscCall(VecNormBegin(X, NORM_2, &xnorm)); in SNESLineSearchApply_NLEQERR()
68 PetscCall(VecNormEnd(X, NORM_2, &xnorm)); in SNESLineSearchApply_NLEQERR()
80 PetscCall(SNESLineSearchSetNorms(linesearch, xnorm, fnorm, ynorm)); in SNESLineSearchApply_NLEQERR()
250 PetscCall(VecNorm(X, NORM_2, &xnorm)); in SNESLineSearchApply_NLEQERR()
254 PetscCall(SNESLineSearchSetNorms(linesearch, xnorm, fnorm, ynorm < 0 ? PETSC_INFINITY : ynorm)); in SNESLineSearchApply_NLEQERR()
/petsc/src/snes/impls/ncg/
H A Dsnesncg.c26 PetscReal *fnorm, *xnorm, *ynorm; in SNESLineSearchApply_NCGLinear() local
35 xnorm = &linesearch->xnorm; in SNESLineSearchApply_NCGLinear()
56 PetscCall(VecNorm(X, NORM_2, xnorm)); in SNESLineSearchApply_NCGLinear()
181 PetscReal fnorm, ynorm, xnorm, beta = 0.0; in SNESSolve_NCG() local
260 PetscCall(SNESLineSearchGetNorms(linesearch, &xnorm, &fnorm, &ynorm)); in SNESSolve_NCG()
270 snes->xnorm = xnorm; in SNESSolve_NCG()
276 PetscCall(SNESConverged(snes, snes->iter, xnorm, ynorm, fnorm)); in SNESSolve_NCG()
/petsc/src/snes/impls/richardson/
H A Dsnesrichardson.c29 PetscReal xnorm, fnorm, ynorm; in SNESSolve_NRichardson() local
91 PetscCall(SNESLineSearchGetNorms(snes->linesearch, &xnorm, &fnorm, &ynorm)); in SNESSolve_NRichardson()
101 snes->xnorm = xnorm; in SNESSolve_NRichardson()
106 PetscCall(SNESConverged(snes, snes->iter, xnorm, ynorm, fnorm)); in SNESSolve_NRichardson()
/petsc/src/snes/linesearch/impls/bt/
H A Dlinesearchbt.c63 PetscReal fnorm, xnorm, ynorm, gnorm; in SNESLineSearchApply_BT() local
92 PetscCall(VecNormBegin(X, NORM_2, &xnorm)); in SNESLineSearchApply_BT()
94 PetscCall(VecNormEnd(X, NORM_2, &xnorm)); in SNESLineSearchApply_BT()
104 PetscCall(SNESLineSearchSetNorms(linesearch, xnorm, fnorm, ynorm)); in SNESLineSearchApply_BT()
176 if (stol * xnorm > ynorm) { in SNESLineSearchApply_BT()
178 PetscCall(SNESLineSearchSetNorms(linesearch, xnorm, fnorm, ynorm)); in SNESLineSearchApply_BT()
182 …: Ended due to ynorm < stol*xnorm (%14.12e < %14.12e).\n", (double)ynorm, (double)(stol * xnorm))); in SNESLineSearchApply_BT()
344 PetscCall(VecNorm(X, NORM_2, &xnorm)); in SNESLineSearchApply_BT()
345 PetscCall(SNESLineSearchSetNorms(linesearch, xnorm, gnorm, ynorm)); in SNESLineSearchApply_BT()
/petsc/src/snes/impls/ls/
H A Dls.c129 PetscReal fnorm, xnorm, ynorm; in SNESSolve_NEWTONLS() local
236 PetscCall(SNESLineSearchGetNorms(linesearch, &xnorm, &fnorm, &ynorm)); in SNESSolve_NEWTONLS()
241 if (snes->stol * xnorm > ynorm) { in SNESSolve_NEWTONLS()
270 snes->xnorm = xnorm; in SNESSolve_NEWTONLS()
274 PetscCall(SNESConverged(snes, snes->iter, xnorm, ynorm, fnorm)); in SNESSolve_NEWTONLS()
/petsc/src/snes/impls/vi/ss/
H A Dviss.c200 PetscReal gnorm, xnorm = 0, ynorm; in SNESSolve_VINEWTONSSLS() local
238 PetscCall(VecNormBegin(X, NORM_2, &xnorm)); /* xnorm <- ||x|| */ in SNESSolve_VINEWTONSSLS()
239 PetscCall(VecNormEnd(X, NORM_2, &xnorm)); in SNESSolve_VINEWTONSSLS()
304 PetscCall(SNESLineSearchGetNorms(snes->linesearch, &xnorm, &gnorm, &ynorm)); in SNESSolve_VINEWTONSSLS()
329 snes->xnorm = xnorm; in SNESSolve_VINEWTONSSLS()
334 if (snes->ops->converged != SNESConvergedSkip) PetscCall(VecNorm(X, NORM_2, &xnorm)); in SNESSolve_VINEWTONSSLS()
335 PetscCall(SNESConverged(snes, snes->iter, xnorm, ynorm, vi->phinorm)); in SNESSolve_VINEWTONSSLS()
/petsc/src/snes/linesearch/interface/
H A Dlinesearch.c189 linesearch->xnorm = 1.0; in SNESLineSearchCreate()
1355 PetscErrorCode SNESLineSearchGetNorms(SNESLineSearch linesearch, PetscReal *xnorm, PetscReal *fnorm… in SNESLineSearchGetNorms() argument
1359 if (xnorm) *xnorm = linesearch->xnorm; in SNESLineSearchGetNorms()
1383 PetscErrorCode SNESLineSearchSetNorms(SNESLineSearch linesearch, PetscReal xnorm, PetscReal fnorm, … in SNESLineSearchSetNorms() argument
1387 linesearch->xnorm = xnorm; in SNESLineSearchSetNorms()
1417 PetscCall(VecNorm(linesearch->vec_sol, NORM_2, &linesearch->xnorm)); in SNESLineSearchComputeNorms()
1422 PetscCall(VecNormBegin(linesearch->vec_sol, NORM_2, &linesearch->xnorm)); in SNESLineSearchComputeNorms()
1425 PetscCall(VecNormEnd(linesearch->vec_sol, NORM_2, &linesearch->xnorm)); in SNESLineSearchComputeNorms()
/petsc/src/snes/impls/patch/
H A Dsnespatch.c260 PetscReal xnorm = 0.0, ynorm = 0.0, fnorm = 0.0; in SNESSolve_Patch() local
274 PetscCall(VecNorm(state, NORM_2, &xnorm)); in SNESSolve_Patch()
280 PetscUseTypeMethod(snes, converged, its, xnorm, ynorm, fnorm, &snes->reason, snes->cnvP); in SNESSolve_Patch()
282 PetscCall(SNESConvergedSkip(snes, its, xnorm, ynorm, fnorm, &snes->reason, NULL)); in SNESSolve_Patch()
310 PetscCall(VecNorm(state, NORM_2, &xnorm)); in SNESSolve_Patch()
314 PetscUseTypeMethod(snes, converged, its, xnorm, ynorm, fnorm, &snes->reason, snes->cnvP); in SNESSolve_Patch()
316 PetscCall(SNESConvergedSkip(snes, its, xnorm, ynorm, fnorm, &snes->reason, NULL)); in SNESSolve_Patch()
/petsc/src/snes/tutorials/
H A Dex1f.F90122 PetscReal ynorm, gnorm, xnorm
134 PetscCall(VecNorm(x, NORM_2, xnorm, ierr))
136 PetscCall(SNESLineSearchSetNorms(linesearch, xnorm, gnorm, ynorm, ierr))
/petsc/src/ksp/ksp/tests/
H A Dex14.c90 PetscReal fnorm, ynorm, xnorm; /* various vector norms */ in main() local
223 PetscCall(VecNorm(X, NORM_2, &xnorm)); /* xnorm = || X || */ in main()
224 …lve iterations = %" PetscInt_FMT ", xnorm=%g, ynorm=%g\n", lin_its, (double)xnorm, (double)ynorm)); in main()
240 if (ynorm < xtol * (xnorm)) { in main()
241 …onverged due to small update length: %g < %g * %g\n", (double)ynorm, (double)xtol, (double)xnorm)); in main()
/petsc/src/ksp/ksp/impls/minres/
H A Dminres.c104 …PetscReal Axnorm, xnorm, xnorm_tmp, xl2norm = 0.0, pnorm, Anorm = 0.0, gmin = 0.0, gminl = 0.0,… in KSPSolve_MINRES() local
145 PetscCall(VecNorm(X, NORM_2, &xnorm)); in KSPSolve_MINRES()
150 xnorm = 0.0; in KSPSolve_MINRES()
294 xnorm = Norm3(xl2norm, ul, u); in KSPSolve_MINRES()
308 if (xnorm < minres->maxxnorm) { in KSPSolve_MINRES()
396 relres = rnorm / (Anorm * xnorm + beta1); in KSPSolve_MINRES()
401 epsx = Anorm * xnorm * eps; in KSPSolve_MINRES()
408 if (xnorm >= minres->maxxnorm) flag = 6; // xnorm exceeded its limit in KSPSolve_MINRES()
425 …double)Arnorml, (double)relresl, (double)relAresl, (double)Anorml, (double)Acondl, (double)xnorm)); in KSPSolve_MINRES()
460 PetscCall(VecNorm(X, NORM_2, &xnorm)); in KSPSolve_MINRES()
[all …]
/petsc/src/snes/impls/composite/
H A Dsnescomposite.c692 PetscReal fnorm = 0.0, xnorm = 0.0, snorm = 0.0; in SNESSolve_Composite() local
759 PetscCall(VecNormBegin(X, NORM_2, &xnorm)); in SNESSolve_Composite()
762 PetscCall(VecNormEnd(X, NORM_2, &xnorm)); in SNESSolve_Composite()
766 PetscCall(VecNormBegin(X, NORM_2, &xnorm)); in SNESSolve_Composite()
770 PetscCall(VecNormEnd(X, NORM_2, &xnorm)); in SNESSolve_Composite()
775 PetscCall(VecNormBegin(X, NORM_2, &xnorm)); in SNESSolve_Composite()
777 PetscCall(VecNormEnd(X, NORM_2, &xnorm)); in SNESSolve_Composite()
784 snes->xnorm = xnorm; in SNESSolve_Composite()
789 PetscCall(SNESConverged(snes, snes->iter, xnorm, snorm, fnorm)); in SNESSolve_Composite()
/petsc/src/snes/impls/al/
H A Dal.c287 PetscReal fnorm, xnorm, ynorm, stepSize; in SNESSolve_NEWTONAL() local
482 PetscCall(VecNormBegin(X, NORM_2, &xnorm)); in SNESSolve_NEWTONAL()
485 PetscCall(VecNormEnd(X, NORM_2, &xnorm)); in SNESSolve_NEWTONAL()
495 snes->xnorm = xnorm; in SNESSolve_NEWTONAL()
498 PetscCall(SNESConverged(snes, snes->iter, xnorm, ynorm, fnorm)); in SNESSolve_NEWTONAL()

123