Lines Matching refs:tau_k
439 PetscReal gamma = 1.0, tau_k, beta; in TaoBNCGStepDirectionUpdate() local
503 PetscCall(TaoBNCGComputeScalarScaling(ynorm2, step * dk_yk, cg->sts, &tau_k, cg->alpha)); in TaoBNCGStepDirectionUpdate()
505 tau_k = 1.0; in TaoBNCGStepDirectionUpdate()
508 PetscCall(VecAXPBY(tao->stepdirection, -tau_k, 0.0, tao->gradient)); in TaoBNCGStepDirectionUpdate()
520 PetscCall(TaoBNCGComputeScalarScaling(ynorm2, step * dk_yk, cg->sts, &tau_k, cg->alpha)); in TaoBNCGStepDirectionUpdate()
521 beta = tau_k * gkp1_yk / dk_yk; in TaoBNCGStepDirectionUpdate()
522 PetscCall(VecAXPBY(tao->stepdirection, -tau_k, beta, tao->gradient)); in TaoBNCGStepDirectionUpdate()
538 …TaoBNCGComputeScalarScaling(ynorm2, step * dk_yk, step * step * dnorm * dnorm, &tau_k, cg->alpha)); in TaoBNCGStepDirectionUpdate()
539 beta = tau_k * gnorm2 / gnorm2_old; in TaoBNCGStepDirectionUpdate()
540 PetscCall(VecAXPBY(tao->stepdirection, -tau_k, beta, tao->gradient)); in TaoBNCGStepDirectionUpdate()
555 … PetscCall(TaoBNCGComputeScalarScaling(ynorm2, step * dk_yk, snorm * snorm, &tau_k, cg->alpha)); in TaoBNCGStepDirectionUpdate()
556 beta = tau_k * gkp1_yk / gnorm2_old; in TaoBNCGStepDirectionUpdate()
557 PetscCall(VecAXPBY(tao->stepdirection, -tau_k, beta, tao->gradient)); in TaoBNCGStepDirectionUpdate()
574 … PetscCall(TaoBNCGComputeScalarScaling(ynorm2, step * dk_yk, snorm * snorm, &tau_k, cg->alpha)); in TaoBNCGStepDirectionUpdate()
575 beta = tau_k * gkp1_yk / gnorm2_old; in TaoBNCGStepDirectionUpdate()
577 PetscCall(VecAXPBY(tao->stepdirection, -tau_k, beta, tao->gradient)); in TaoBNCGStepDirectionUpdate()
594 PetscCall(TaoBNCGComputeScalarScaling(ynorm2, step * dk_yk, cg->yts, &tau_k, cg->alpha)); in TaoBNCGStepDirectionUpdate()
595 beta = tau_k * gnorm2 / (gd - gd_old); in TaoBNCGStepDirectionUpdate()
596 PetscCall(VecAXPBY(tao->stepdirection, -tau_k, beta, tao->gradient)); in TaoBNCGStepDirectionUpdate()
624 PetscCall(TaoBNCGComputeScalarScaling(ynorm2, cg->yts, snorm * snorm, &tau_k, cg->alpha)); 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()
644 PetscCall(VecDot(cg->yk, cg->y_work, &tau_k)); in TaoBNCGStepDirectionUpdate()
645 tau_k = -tau_k * gd / (dk_yk * dk_yk); in TaoBNCGStepDirectionUpdate()
647 beta = gkp1_yk / dk_yk + cg->hz_theta * tau_k; /* HZ; (1.15) from DK 2013 */ in TaoBNCGStepDirectionUpdate()
668 PetscCall(TaoBNCGComputeScalarScaling(ynorm2, cg->yts, snorm * snorm, &tau_k, cg->alpha)); 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()
684 PetscCall(VecDot(cg->yk, cg->y_work, &tau_k)); in TaoBNCGStepDirectionUpdate()
685 tau_k = tau_k * gd / (dk_yk * dk_yk); in TaoBNCGStepDirectionUpdate()
688 beta = gkp1_yk / dk_yk - step * tmp - tau_k; in TaoBNCGStepDirectionUpdate()
714 PetscCall(TaoBNCGComputeScalarScaling(ynorm2, cg->yts, snorm * snorm, &tau_k, cg->alpha)); 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()
739 PetscCall(VecDot(cg->yk, cg->y_work, &tau_k)); in TaoBNCGStepDirectionUpdate()
740 tau_k = tau_k * gd / (dk_yk * dk_yk); in TaoBNCGStepDirectionUpdate()
742 beta = gkp1D_yk / dk_yk - step * gamma - tau_k; in TaoBNCGStepDirectionUpdate()
777 PetscCall(TaoBNCGComputeScalarScaling(cg->yty, cg->yts, cg->sts, &tau_k, cg->alpha)); 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()
806 PetscCall(TaoBNCGComputeScalarScaling(cg->yty, cg->yts, cg->sts, &tau_k, cg->alpha)); in TaoBNCGStepDirectionUpdate()
808 tau_k = cg->dfp_scale * tau_k; in TaoBNCGStepDirectionUpdate()
809 tmp = tau_k * gkp1_yk / cg->yty; in TaoBNCGStepDirectionUpdate()
812 PetscCall(VecAXPBYPCZ(tao->stepdirection, -tau_k, tmp, beta, tao->gradient, cg->yk)); in TaoBNCGStepDirectionUpdate()
840 tau_k = cg->theta * tau_bfgs + (1.0 - cg->theta) * tau_dfp; in TaoBNCGStepDirectionUpdate()
846 PetscCall(VecAXPBYPCZ(tao->stepdirection, -tau_k, tmp, beta, tao->gradient, cg->yk)); in TaoBNCGStepDirectionUpdate()