Lines Matching refs:beta
439 PetscReal gamma = 1.0, tau_k, beta; in TaoBNCGStepDirectionUpdate() local
521 beta = tau_k * gkp1_yk / dk_yk; in TaoBNCGStepDirectionUpdate()
522 PetscCall(VecAXPBY(tao->stepdirection, -tau_k, beta, tao->gradient)); in TaoBNCGStepDirectionUpdate()
526 beta = gkp1_yk / dk_yk; in TaoBNCGStepDirectionUpdate()
527 PetscCall(VecAXPBY(tao->stepdirection, -1.0, beta, cg->g_work)); in TaoBNCGStepDirectionUpdate()
539 beta = tau_k * gnorm2 / gnorm2_old; in TaoBNCGStepDirectionUpdate()
540 PetscCall(VecAXPBY(tao->stepdirection, -tau_k, beta, tao->gradient)); in TaoBNCGStepDirectionUpdate()
545 beta = tmp / gnorm2_old; in TaoBNCGStepDirectionUpdate()
546 PetscCall(VecAXPBY(tao->stepdirection, -1.0, beta, cg->g_work)); in TaoBNCGStepDirectionUpdate()
556 beta = tau_k * gkp1_yk / gnorm2_old; in TaoBNCGStepDirectionUpdate()
557 PetscCall(VecAXPBY(tao->stepdirection, -tau_k, beta, tao->gradient)); in TaoBNCGStepDirectionUpdate()
562 beta = gkp1_yk / gnorm2_old; in TaoBNCGStepDirectionUpdate()
563 PetscCall(VecAXPBY(tao->stepdirection, -1.0, beta, cg->g_work)); in TaoBNCGStepDirectionUpdate()
575 beta = tau_k * gkp1_yk / gnorm2_old; in TaoBNCGStepDirectionUpdate()
576 beta = PetscMax(beta, 0.0); in TaoBNCGStepDirectionUpdate()
577 PetscCall(VecAXPBY(tao->stepdirection, -tau_k, beta, tao->gradient)); in TaoBNCGStepDirectionUpdate()
582 beta = gkp1_yk / gnorm2_old; in TaoBNCGStepDirectionUpdate()
583 beta = PetscMax(beta, 0.0); in TaoBNCGStepDirectionUpdate()
584 PetscCall(VecAXPBY(tao->stepdirection, -1.0, beta, cg->g_work)); in TaoBNCGStepDirectionUpdate()
595 beta = tau_k * gnorm2 / (gd - gd_old); in TaoBNCGStepDirectionUpdate()
596 PetscCall(VecAXPBY(tao->stepdirection, -tau_k, beta, tao->gradient)); in TaoBNCGStepDirectionUpdate()
604 beta = gtDg / dk_yk; in TaoBNCGStepDirectionUpdate()
605 PetscCall(VecScale(cg->d_work, beta)); in TaoBNCGStepDirectionUpdate()
627 beta = tau_k * (gkp1_yk / dk_yk - ynorm2 * gd / (dk_yk * dk_yk)); in TaoBNCGStepDirectionUpdate()
629 …beta = PetscMax(PetscMax(beta, cg->hz_eta * tau_k * gd_old / (dnorm * dnorm)), cg->dk_eta * tau_k … in TaoBNCGStepDirectionUpdate()
631 PetscCall(VecAXPBY(tao->stepdirection, -tau_k, beta, tao->gradient)); in TaoBNCGStepDirectionUpdate()
647 beta = gkp1_yk / dk_yk + cg->hz_theta * tau_k; /* HZ; (1.15) from DK 2013 */ in TaoBNCGStepDirectionUpdate()
651 …beta = PetscMax(PetscMax(beta, cg->hz_eta * gd_old / (dnorm * dnorm)), cg->dk_eta * gd / (dnorm * … in TaoBNCGStepDirectionUpdate()
653 PetscCall(VecAXPBY(tao->stepdirection, -1.0, beta, cg->g_work)); in TaoBNCGStepDirectionUpdate()
671 …beta = tau_k * (gkp1_yk / dk_yk - ynorm2 * gd / (dk_yk * dk_yk) + gd / (dnorm * dnorm)) - step * g… in TaoBNCGStepDirectionUpdate()
672 …beta = PetscMax(PetscMax(beta, cg->hz_eta * tau_k * gd_old / (dnorm * dnorm)), cg->dk_eta * tau_k … in TaoBNCGStepDirectionUpdate()
674 PetscCall(VecAXPBYPCZ(tao->stepdirection, -tau_k, 0.0, beta, tao->gradient, cg->yk)); in TaoBNCGStepDirectionUpdate()
688 beta = gkp1_yk / dk_yk - step * tmp - tau_k; in TaoBNCGStepDirectionUpdate()
691 beta += tmp * dk_yk / (dnorm * dnorm); /* projection of y_work onto dk */ in TaoBNCGStepDirectionUpdate()
694 …beta = PetscMax(PetscMax(beta, cg->hz_eta * gd_old / (dnorm * dnorm)), cg->dk_eta * gd / (dnorm * … in TaoBNCGStepDirectionUpdate()
696 PetscCall(VecAXPBY(tao->stepdirection, -1.0, beta, cg->g_work)); in TaoBNCGStepDirectionUpdate()
715 beta = tau_k * (gkp1_yk / dk_yk - ynorm2 * gd / (dk_yk * dk_yk)) - step * gd / dk_yk; in TaoBNCGStepDirectionUpdate()
716 if (beta < cg->zeta * tau_k * gd / (dnorm * dnorm)) /* 0.1 is KD's zeta parameter */ in TaoBNCGStepDirectionUpdate()
718 beta = cg->zeta * tau_k * gd / (dnorm * dnorm); in TaoBNCGStepDirectionUpdate()
727 … PetscCall(VecAXPBYPCZ(tao->stepdirection, -tau_k, gamma * tau_k, beta, tao->gradient, cg->yk)); in TaoBNCGStepDirectionUpdate()
742 beta = gkp1D_yk / dk_yk - step * gamma - tau_k; in TaoBNCGStepDirectionUpdate()
749 if (beta < cg->zeta * tmp / (dnorm * dnorm)) { in TaoBNCGStepDirectionUpdate()
750 beta = cg->zeta * tmp / (dnorm * dnorm); in TaoBNCGStepDirectionUpdate()
754 if (beta < cg->zeta * tmp / (dnorm * dnorm)) { in TaoBNCGStepDirectionUpdate()
755 beta = cg->zeta * tmp / (dnorm * dnorm); in TaoBNCGStepDirectionUpdate()
760 PetscCall(VecAXPBY(tao->stepdirection, -1.0, beta, cg->g_work)); in TaoBNCGStepDirectionUpdate()
779 beta = tau_k * (gkp1_yk / dk_yk - cg->yty * gd / (dk_yk * dk_yk)) - step * tmp; in TaoBNCGStepDirectionUpdate()
781 … PetscCall(VecAXPBYPCZ(tao->stepdirection, -tau_k, tmp * tau_k, beta, tao->gradient, cg->yk)); in TaoBNCGStepDirectionUpdate()
791 beta = (gkp1_yk / dk_yk - gd * tmp / (dk_yk * dk_yk)) - step * gd / dk_yk; in TaoBNCGStepDirectionUpdate()
793 PetscCall(VecAXPBYPCZ(tao->stepdirection, -1.0, gamma, beta, cg->g_work, cg->y_work)); in TaoBNCGStepDirectionUpdate()
810 beta = -step * gd / dk_yk; in TaoBNCGStepDirectionUpdate()
812 PetscCall(VecAXPBYPCZ(tao->stepdirection, -tau_k, tmp, beta, tao->gradient, cg->yk)); in TaoBNCGStepDirectionUpdate()
822 beta = -step * gd / dk_yk; in TaoBNCGStepDirectionUpdate()
824 PetscCall(VecAXPBYPCZ(tao->stepdirection, -1.0, gamma, beta, cg->g_work, cg->y_work)); in TaoBNCGStepDirectionUpdate()
844 …beta = cg->theta * tau_bfgs * (gkp1_yk / dk_yk - cg->yty * gd / (dk_yk * dk_yk)) - step * gd / dk_… in TaoBNCGStepDirectionUpdate()
846 PetscCall(VecAXPBYPCZ(tao->stepdirection, -tau_k, tmp, beta, tao->gradient, cg->yk)); in TaoBNCGStepDirectionUpdate()
856 beta = cg->theta * (gkp1_yk / dk_yk - gd * tmp / (dk_yk * dk_yk)) - step * gd / dk_yk; in TaoBNCGStepDirectionUpdate()
858 PetscCall(VecAXPBYPCZ(tao->stepdirection, -1.0, gamma, beta, cg->g_work, cg->y_work)); in TaoBNCGStepDirectionUpdate()