impls/bncg/bncg.c

441   PetscReal gkp1_yk, gd_old, tau_bfgs, tau_dfp, gkp1D_yk, gtDg;
519         PetscCall(VecDot(cg->yk, tao->gradient, &gkp1_yk));
521         beta = tau_k * gkp1_yk / dk_yk;
525         PetscCall(VecDot(cg->yk, cg->g_work, &gkp1_yk));
526         beta = gkp1_yk / dk_yk;
554         PetscCall(VecDot(cg->yk, tao->gradient, &gkp1_yk));
556         beta = tau_k * gkp1_yk / gnorm2_old;
561         PetscCall(VecDot(cg->g_work, cg->yk, &gkp1_yk));
562         beta = gkp1_yk / gnorm2_old;
573         PetscCall(VecDot(cg->yk, tao->gradient, &gkp1_yk));
575         beta = tau_k * gkp1_yk / gnorm2_old;
581         PetscCall(VecDot(cg->g_work, cg->yk, &gkp1_yk));
582         beta = gkp1_yk / gnorm2_old;
623           PetscCall(VecDot(cg->yk, tao->gradient, &gkp1_yk));
627           beta = tau_k * (gkp1_yk / dk_yk - ynorm2 * gd / (dk_yk * dk_yk));
641           PetscCall(VecDot(cg->yk, cg->g_work, &gkp1_yk));
647           beta = gkp1_yk / dk_yk + cg->hz_theta * tau_k; /* HZ; (1.15) from DK 2013 */
667         PetscCall(VecDot(cg->yk, tao->gradient, &gkp1_yk));
671         beta = tau_k * (gkp1_yk / dk_yk - ynorm2 * gd / (dk_yk * dk_yk) + gd / (dnorm * dnorm)) - step * gd / dk_yk;
683         PetscCall(VecDot(cg->yk, cg->g_work, &gkp1_yk));
688         beta = gkp1_yk / dk_yk - step * tmp - tau_k;
713           PetscCall(VecDot(cg->yk, tao->gradient, &gkp1_yk));
715           beta = tau_k * (gkp1_yk / dk_yk - ynorm2 * gd / (dk_yk * dk_yk)) - step * gd / dk_yk;
721             if (gkp1_yk < 0 && cg->neg_xi) gamma = -1.0 * gd / dk_yk;
776         PetscCall(VecDot(cg->yk, tao->gradient, &gkp1_yk));
779         beta = tau_k * (gkp1_yk / dk_yk - cg->yty * gd / (dk_yk * dk_yk)) - step * tmp;
787         PetscCall(VecDot(cg->g_work, cg->yk, &gkp1_yk));
791         beta = (gkp1_yk / dk_yk - gd * tmp / (dk_yk * dk_yk)) - step * gd / dk_yk;
807         PetscCall(VecDot(cg->yk, tao->gradient, &gkp1_yk));
809         tmp   = tau_k * gkp1_yk / cg->yty;
818         PetscCall(VecDot(cg->g_work, cg->yk, &gkp1_yk));
820         gamma = (gkp1_yk / tmp);
839         PetscCall(VecDot(cg->yk, tao->gradient, &gkp1_yk));
843         tmp  = cg->theta * tau_bfgs * gd / dk_yk + (1 - cg->theta) * tau_dfp * gkp1_yk / cg->yty;
844         beta = cg->theta * tau_bfgs * (gkp1_yk / dk_yk - cg->yty * gd / (dk_yk * dk_yk)) - step * gd / dk_yk;
852         PetscCall(VecDot(cg->g_work, cg->yk, &gkp1_yk));
854         gamma = cg->theta * gd / dk_yk + (1 - cg->theta) * (gkp1_yk / tmp);
856         beta = cg->theta * (gkp1_yk / dk_yk - gd * tmp / (dk_yk * dk_yk)) - step * gd / dk_yk;