Lines Matching refs:neP
33 SNES_NEWTONTRDC *neP = (SNES_NEWTONTRDC *)snes->data; in SNESTRDC_KSPConverged_Private() local
43 if (nrm >= neP->delta) { in SNESTRDC_KSPConverged_Private()
44 …PetscCall(PetscInfo(snes, "Ending linear iteration early, delta=%g, length=%g\n", (double)neP->del… in SNESTRDC_KSPConverged_Private()
65 SNES_NEWTONTRDC *neP = (SNES_NEWTONTRDC *)snes->data; in SNESTRDC_Converged_Private() local
69 if (neP->delta < xnorm * neP->deltatol) { in SNESTRDC_Converged_Private()
70 …ged due to too small a trust region %g<%g*%g\n", (double)neP->delta, (double)xnorm, (double)neP->d… in SNESTRDC_Converged_Private()
297 SNES_NEWTONTRDC *neP = (SNES_NEWTONTRDC *)snes->data; in SNESSolve_NEWTONTRDC() local
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()
357 neP->delta = delta; in SNESSolve_NEWTONTRDC()
361 neP->rho_satisfied = PETSC_FALSE; in SNESSolve_NEWTONTRDC()
383 if (bs > 1 && neP->auto_scale_multiphase) { in SNESSolve_NEWTONTRDC()
386 if (neP->auto_scale_max > 1.0) { in SNESSolve_NEWTONTRDC()
387 if (inorms[j] < 1.0 / neP->auto_scale_max) inorms[j] = 1.0 / neP->auto_scale_max; in SNESSolve_NEWTONTRDC()
395 delta = neP->delta0 * xnorm; in SNESSolve_NEWTONTRDC()
397 delta = neP->delta * xnorm; in SNESSolve_NEWTONTRDC()
399 deltaM = neP->deltaM * xnorm; in SNESSolve_NEWTONTRDC()
408 neP->rho_satisfied = PETSC_FALSE; in SNESSolve_NEWTONTRDC()
412 } else if (neP->use_cauchy) { /* use Cauchy direction if enabled */ in SNESSolve_NEWTONTRDC()
459 if (bs > 1 && neP->auto_scale_multiphase) { in SNESSolve_NEWTONTRDC()
489 if (rho < neP->eta2) { in SNESSolve_NEWTONTRDC()
490 delta *= neP->t1; /* shrink the region */ in SNESSolve_NEWTONTRDC()
491 } else if (rho > neP->eta3) { in SNESSolve_NEWTONTRDC()
492 … delta = PetscMin(neP->t2 * delta, deltaM); /* expand the region, but not greater than deltaM */ in SNESSolve_NEWTONTRDC()
495 neP->delta = delta; in SNESSolve_NEWTONTRDC()
496 if (rho >= neP->eta1) { in SNESSolve_NEWTONTRDC()
498 if (bs > 1 && neP->auto_scale_multiphase) { in SNESSolve_NEWTONTRDC()
499 neP->delta = delta / xnorm; in SNESSolve_NEWTONTRDC()
503 neP->rho_satisfied = PETSC_TRUE; in SNESSolve_NEWTONTRDC()
509 neP->itflag = PETSC_FALSE; in SNESSolve_NEWTONTRDC()
514 if (bs > 1 && neP->auto_scale_multiphase) { in SNESSolve_NEWTONTRDC()
516 neP->delta = delta / xnorm; in SNESSolve_NEWTONTRDC()
554 neP->itflag = PETSC_TRUE; in SNESSolve_NEWTONTRDC()
660 SNES_NEWTONTRDC *neP; in SNESCreate_NEWTONTRDC() local
676 PetscCall(PetscNew(&neP)); in SNESCreate_NEWTONTRDC()
677 snes->data = (void *)neP; in SNESCreate_NEWTONTRDC()
678 neP->eta1 = 0.001; in SNESCreate_NEWTONTRDC()
679 neP->eta2 = 0.25; in SNESCreate_NEWTONTRDC()
680 neP->eta3 = 0.75; in SNESCreate_NEWTONTRDC()
681 neP->t1 = 0.25; in SNESCreate_NEWTONTRDC()
682 neP->t2 = 2.0; in SNESCreate_NEWTONTRDC()
683 neP->sigma = 0.0001; in SNESCreate_NEWTONTRDC()
684 neP->itflag = PETSC_FALSE; in SNESCreate_NEWTONTRDC()
685 neP->rnorm0 = 0.0; in SNESCreate_NEWTONTRDC()
686 neP->ttol = 0.0; in SNESCreate_NEWTONTRDC()
687 neP->use_cauchy = PETSC_TRUE; in SNESCreate_NEWTONTRDC()
688 neP->auto_scale_multiphase = PETSC_FALSE; in SNESCreate_NEWTONTRDC()
689 neP->auto_scale_max = -1.0; in SNESCreate_NEWTONTRDC()
690 neP->rho_satisfied = PETSC_FALSE; in SNESCreate_NEWTONTRDC()
691 neP->delta = 0.0; in SNESCreate_NEWTONTRDC()
692 neP->deltaM = 0.5; in SNESCreate_NEWTONTRDC()
693 neP->delta0 = 0.1; in SNESCreate_NEWTONTRDC()
694 neP->deltatol = 1.e-12; in SNESCreate_NEWTONTRDC()