Home
last modified time | relevance | path

Searched refs:ynorm (Results 1 – 25 of 41) sorted by relevance

12

/petsc/src/mat/tests/
H A Dex185.c11 PetscReal xnorm, ynorm, anorm; in main() local
23 PetscCall(VecNorm(Y, NORM_2, &ynorm)); in main()
24 …bsReal(ynorm - 3 * xnorm) <= PETSC_SMALL, PETSC_COMM_WORLD, PETSC_ERR_PLIB, "Expected norm %g actu… in main()
51 PetscCall(VecNorm(Y, NORM_2, &ynorm)); in main()
52 …sReal(ynorm - xnorm / 4) <= PETSC_SMALL, PETSC_COMM_WORLD, PETSC_ERR_PLIB, "Expected norm %g actua… in main()
56 PetscCall(VecNorm(Y, NORM_2, &ynorm)); in main()
57 …sReal(ynorm - xnorm / 4) <= PETSC_SMALL, PETSC_COMM_WORLD, PETSC_ERR_PLIB, "Expected norm %g actua… in main()
61 PetscCall(VecNorm(Y, NORM_2, &ynorm)); in main()
62 …sReal(ynorm - xnorm / 4) <= PETSC_SMALL, PETSC_COMM_WORLD, PETSC_ERR_PLIB, "Expected norm %g actua… in main()
/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()
65 PetscCall(VecNormBegin(Y, NORM_2, &ynorm)); in SNESLineSearchApply_NLEQERR()
67 PetscCall(VecNormEnd(Y, NORM_2, &ynorm)); in SNESLineSearchApply_NLEQERR()
72 if (ynorm == 0.0) { in SNESLineSearchApply_NLEQERR()
80 PetscCall(SNESLineSearchSetNorms(linesearch, xnorm, fnorm, ynorm)); in SNESLineSearchApply_NLEQERR()
93 …PetscViewerASCIIPrintf(monitor, " Line search: norm of Newton step: %14.12e\n", (double)ynorm)); in SNESLineSearchApply_NLEQERR()
103 …err->lambda_prev * (nleqerr->norm_delta_x_prev * nleqerr->norm_bar_delta_x_prev) / (gnorm * ynorm); in SNESLineSearchApply_NLEQERR()
159 ynorm = wnorm = -1.0; in SNESLineSearchApply_NLEQERR()
183 theta = wnorm / ynorm; in SNESLineSearchApply_NLEQERR()
[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
/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()
37 PetscCall(VecNormBegin(Y, NORM_2, &linesearch->ynorm)); in SNESLineSearchApply_Basic()
39 PetscCall(VecNormEnd(Y, NORM_2, &linesearch->ynorm)); in SNESLineSearchApply_Basic()
/petsc/src/snes/linesearch/impls/bt/
H A Dlinesearchbt.c63 PetscReal fnorm, xnorm, ynorm, gnorm; in SNESLineSearchApply_BT() local
91 PetscCall(VecNormBegin(Y, NORM_2, &ynorm)); in SNESLineSearchApply_BT()
93 PetscCall(VecNormEnd(Y, NORM_2, &ynorm)); in SNESLineSearchApply_BT()
96 if (ynorm == 0.0) { 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 … search: Ended due to ynorm < stol*xnorm (%14.12e < %14.12e).\n", (double)ynorm, (double)(stol * x… in SNESLineSearchApply_BT()
245 …=%18.16e, initial slope=%18.16e\n", (double)fnorm, (double)gnorm, (double)ynorm, (double)minlambda… in SNESLineSearchApply_BT()
247 …, lambda=%18.16e, initial slope=%18.16e\n", (double)f, (double)g, (double)ynorm, (double)minlambda… in SNESLineSearchApply_BT()
[all …]
/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()
159 snes->ynorm = ynorm; in SNESSolve_Anderson()
162 PetscCall(SNESConverged(snes, snes->iter, xnorm, ynorm, fnorm)); in SNESSolve_Anderson()
H A Dngmresfunc.c163 …min, PetscReal dminnorm, Vec X, Vec F, Vec Y, PetscReal *xnorm, PetscReal *fnorm, PetscReal *ynorm) in SNESNGMRESSelect_Private() argument
183 *ynorm = yMnorm; in SNESNGMRESSelect_Private()
206 PetscCall(SNESLineSearchGetNorms(ngmres->additive_linesearch, xnorm, fnorm, ynorm)); in SNESNGMRESSelect_Private()
228 *ynorm = yAnorm; in SNESNGMRESSelect_Private()
235 *ynorm = yMnorm; in SNESNGMRESSelect_Private()
243 *ynorm = yAnorm; in SNESNGMRESSelect_Private()
/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()
34 if (PetscAbsScalar(fty_old) < atol * ynorm) { in SNESLineSearchApply_CP()
37 …ause dot(F,Y) = %g < atol*||y|| = %g\n", (double)PetscAbsScalar(fty_old), (double)(atol * ynorm))); in SNESLineSearchApply_CP()
80 if (PetscAbsScalar(fty) < atol * ynorm && i > 0) { in SNESLineSearchApply_CP()
83 …search: abs(fty)/||y|| = %g <= atol = %g\n", (double)(PetscAbsScalar(fty) / ynorm), (double)atol)); in SNESLineSearchApply_CP()
190 PetscCall(SNESLineSearchGetNorms(linesearch, &xnorm, &gnorm, &ynorm)); in SNESLineSearchApply_CP()
/petsc/src/vec/vec/tests/
H A Dex42.c43 PetscReal ynorm; in main() local
45 PetscCall(VecNormBegin(y, NORM_2, &ynorm)); in main()
51 PetscCall(VecNormEnd(y, NORM_2, &ynorm)); in main()
/petsc/src/snes/impls/ncg/
H A Dsnesncg.c26 PetscReal *fnorm, *xnorm, *ynorm; in SNESLineSearchApply_NCGLinear() local
36 ynorm = &linesearch->ynorm; in SNESLineSearchApply_NCGLinear()
57 PetscCall(VecNorm(Y, NORM_2, ynorm)); 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()
271 snes->ynorm = ynorm; 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()
102 snes->ynorm = ynorm; in SNESSolve_NRichardson()
106 PetscCall(SNESConverged(snes, snes->iter, xnorm, ynorm, fnorm)); in SNESSolve_NRichardson()
/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()
237 …8.16e, ynorm=%18.16e, lsreason=%d\n", (double)gnorm, (double)fnorm, (double)ynorm, (int)lsreason)); in SNESSolve_NEWTONLS()
241 if (snes->stol * xnorm > ynorm) { in SNESSolve_NEWTONLS()
269 snes->ynorm = ynorm; in SNESSolve_NEWTONLS()
274 PetscCall(SNESConverged(snes, snes->iter, xnorm, ynorm, fnorm)); in SNESSolve_NEWTONLS()
/petsc/src/snes/linesearch/impls/bisection/
H A Dlinesearchbisection.c9 PetscReal ynorm; in SNESLineSearchApply_Bisection() local
28 PetscCall(VecNorm(Y, NORM_2, &ynorm)); in SNESLineSearchApply_Bisection()
88 if (PetscAbsScalar(fty) <= atol * ynorm) { in SNESLineSearchApply_Bisection()
91 …search: abs(fty)/||y|| = %g <= atol = %g\n", (double)(PetscAbsScalar(fty) / ynorm), (double)atol)); in SNESLineSearchApply_Bisection()
162 …search: fty/||y|| = %g, lambda = %g\n", it, (double)(PetscRealPart(fty) / ynorm), (double)lambda)); in SNESLineSearchApply_Bisection()
/petsc/src/snes/impls/tr/
H A Dtr.c490 PetscReal rho, fnorm, gnorm = 0.0, xnorm = 0.0, delta, ynorm; in SNESSolve_NEWTONTR() local
670 PetscCall(VecNorm(Y, neP->norm, &ynorm)); in SNESSolve_NEWTONTR()
673 if (!use_cauchy && (ynorm > delta || ynorm == 0.0)) { in SNESSolve_NEWTONTR()
674 SNESNewtonTRFallbackType fallback = ynorm > 0.0 ? neP->fallback : SNES_TR_FALLBACK_CAUCHY; in SNESSolve_NEWTONTR()
679 auk = delta / ynorm; in SNESSolve_NEWTONTR()
681 … PetscCall(PetscInfo(snes, "SN evaluated. delta: %g, ynorm: %g\n", (double)delta, (double)ynorm)); in SNESSolve_NEWTONTR()
687 …ed. delta: %g, ynorm: %g, ycnorm: %g, gTBg: %g\n", (double)delta, (double)ynorm, (double)ycnorm, (… in SNESSolve_NEWTONTR()
692 …tside region). delta: %g, ynorm: %g, ycnorm: %g\n", (double)delta, (double)ynorm, (double)ycnorm)); in SNESSolve_NEWTONTR()
726 …, noroots ? "CP" : "DL", (double)tneg, (double)tpos, (double)tau, (double)ynorm, (double)ycnorm, (… in SNESSolve_NEWTONTR()
747 PetscCall(VecNorm(Y, neP->norm, &ynorm)); in SNESSolve_NEWTONTR()
[all …]
/petsc/src/snes/linesearch/interface/
H A Dlinesearch.c188 linesearch->ynorm = 1.0; in SNESLineSearchCreate()
539 PetscReal ynorm, ylastnorm, theta, angle_radians; in SNESLineSearchPreCheckPicard() local
558 PetscCall(VecNorm(Y, NORM_2, &ynorm)); in SNESLineSearchPreCheckPicard()
560 if (ynorm == 0. || ylastnorm == 0.) { in SNESLineSearchPreCheckPicard()
565 …theta = PetscAcosReal((PetscReal)PetscClipInterval(PetscAbsScalar(dot) / (ynorm * ylastnor… in SNESLineSearchPreCheckPicard()
1355 …LineSearchGetNorms(SNESLineSearch linesearch, PetscReal *xnorm, PetscReal *fnorm, PetscReal *ynorm) in SNESLineSearchGetNorms() argument
1361 if (ynorm) *ynorm = linesearch->ynorm; in SNESLineSearchGetNorms()
1383 …NESLineSearchSetNorms(SNESLineSearch linesearch, PetscReal xnorm, PetscReal fnorm, PetscReal ynorm) in SNESLineSearchSetNorms() argument
1389 linesearch->ynorm = ynorm; in SNESLineSearchSetNorms()
1418 PetscCall(VecNorm(linesearch->vec_update, NORM_2, &linesearch->ynorm)); in SNESLineSearchComputeNorms()
[all …]
/petsc/src/snes/tutorials/
H A Dex1f.F90122 PetscReal ynorm, gnorm, xnorm
130 PetscCall(VecNorm(y, NORM_2, ynorm, ierr))
135 PetscCall(VecNorm(y, NORM_2, ynorm, ierr))
136 PetscCall(SNESLineSearchSetNorms(linesearch, xnorm, gnorm, ynorm, ierr))
/petsc/src/snes/impls/vi/ss/
H A Dviss.c200 PetscReal gnorm, xnorm = 0, ynorm; in SNESSolve_VINEWTONSSLS() local
300 ynorm = 1; in SNESSolve_VINEWTONSSLS()
304 PetscCall(SNESLineSearchGetNorms(snes->linesearch, &xnorm, &gnorm, &ynorm)); in SNESSolve_VINEWTONSSLS()
305 …norm=%18.16e, lssucceed=%d\n", (double)vi->phinorm, (double)gnorm, (double)ynorm, (int)lssucceed)); in SNESSolve_VINEWTONSSLS()
330 snes->ynorm = ynorm; in SNESSolve_VINEWTONSSLS()
335 PetscCall(SNESConverged(snes, snes->iter, xnorm, ynorm, vi->phinorm)); in SNESSolve_VINEWTONSSLS()
/petsc/src/snes/impls/ntrdc/
H A Dntrdc.c302 …PetscReal rho, fnorm, gnorm, xnorm = 0, delta, ynorm, temp_xnorm, temp_ynorm; /*… in SNESSolve_NEWTONTRDC() local
451 PetscCall(VecNorm(Y, NORM_2, &ynorm)); /* compute the final ynorm */ in SNESSolve_NEWTONTRDC()
472 temp_ynorm = ynorm; in SNESSolve_NEWTONTRDC()
501 ynorm = temp_ynorm; in SNESSolve_NEWTONTRDC()
511 … PetscCall(SNESTRDC_Converged_Private(snes, snes->iter, xnorm, ynorm, fnorm, &reason, snes->cnvP)); in SNESSolve_NEWTONTRDC()
518 ynorm = temp_ynorm; in SNESSolve_NEWTONTRDC()
549 snes->ynorm = ynorm; in SNESSolve_NEWTONTRDC()
556 PetscUseTypeMethod(snes, converged, snes->iter, xnorm, ynorm, fnorm, &reason, snes->cnvP); in SNESSolve_NEWTONTRDC()
/petsc/src/snes/impls/qn/
H A Dqn.c68 PetscReal fnorm, xnorm, ynorm; in SNESSolve_QN() local
176 ynorm = 1; in SNESSolve_QN()
186 PetscCall(SNESLineSearchGetNorms(snes->linesearch, &xnorm, &fnorm, &ynorm)); in SNESSolve_QN()
188 ….16e, ynorm=%18.16e, lssucceed=%d\n", (double)fnorm, (double)gnorm, (double)ynorm, (int)lsreason)); in SNESSolve_QN()
194 snes->ynorm = ynorm; in SNESSolve_QN()
200 PetscCall(SNESConverged(snes, snes->iter, xnorm, ynorm, fnorm)); in SNESSolve_QN()
/petsc/src/ksp/ksp/tests/
H A Dex14.c90 PetscReal fnorm, ynorm, xnorm; /* various vector norms */ in main() local
220 PetscCall(VecNorm(Y, NORM_2, &ynorm)); /* ynorm = || Y || */ 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 …ntf(comm, "Converged due to small update length: %g < %g * %g\n", (double)ynorm, (double)xtol, (do… in main()
/petsc/src/snes/impls/al/
H A Dal.c287 PetscReal fnorm, xnorm, ynorm, stepSize; in SNESSolve_NEWTONAL() local
483 PetscCall(VecNormBegin(deltaX, NORM_2, &ynorm)); in SNESSolve_NEWTONAL()
486 PetscCall(VecNormEnd(deltaX, NORM_2, &ynorm)); in SNESSolve_NEWTONAL()
494 snes->ynorm = ynorm; in SNESSolve_NEWTONAL()
498 PetscCall(SNESConverged(snes, snes->iter, xnorm, ynorm, fnorm)); in SNESSolve_NEWTONAL()
/petsc/src/snes/impls/patch/
H A Dsnespatch.c260 PetscReal xnorm = 0.0, ynorm = 0.0, fnorm = 0.0; in SNESSolve_Patch() local
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()
311 PetscCall(VecNorm(update, NORM_2, &ynorm)); 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/impls/vi/rs/
H A Dvirs.c301 PetscReal fnorm, gnorm, xnorm = 0, ynorm; in SNESSolve_VINEWTONRSLS() local
535 ynorm = 1; in SNESSolve_VINEWTONRSLS()
542 if (snes->stol * xnorm > ynorm) { in SNESSolve_VINEWTONRSLS()
566 PetscCall(SNESLineSearchGetNorms(snes->linesearch, &xnorm, &gnorm, &ynorm)); in SNESSolve_VINEWTONRSLS()
567 ….16e, ynorm=%18.16e, lssucceed=%d\n", (double)fnorm, (double)gnorm, (double)ynorm, (int)lsreason)); in SNESSolve_VINEWTONRSLS()
575 snes->ynorm = ynorm; in SNESSolve_VINEWTONRSLS()
580 PetscCall(SNESConverged(snes, snes->iter, xnorm, ynorm, fnorm)); in SNESSolve_VINEWTONRSLS()
/petsc/src/tao/leastsquares/impls/pounders/
H A Dgqt.c8 PetscReal e, temp, w, wm, ynorm, znorm, s, sm; in estsv() local
56 PetscCallBLAS("BLASnrm2", ynorm = BLASnrm2_(&blasn, z, &blas1)); in estsv()
64 ynorm *= temp; in estsv()
78 *svmin = ynorm * znorm; in estsv()

12