Lines Matching refs:dk_yk
440 PetscReal tmp = 1.0, ynorm, ynorm2 = 1.0, snorm = 1.0, dk_yk = 1.0, gd; in TaoBNCGStepDirectionUpdate() local
451 PetscCall(VecDot(cg->yk, tao->stepdirection, &dk_yk)); in TaoBNCGStepDirectionUpdate()
452 if (step * dnorm < PETSC_MACHINE_EPSILON || step * dk_yk < PETSC_MACHINE_EPSILON) { in TaoBNCGStepDirectionUpdate()
503 PetscCall(TaoBNCGComputeScalarScaling(ynorm2, step * dk_yk, cg->sts, &tau_k, cg->alpha)); 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()
526 beta = gkp1_yk / dk_yk; in TaoBNCGStepDirectionUpdate()
536 PetscCall(VecDot(cg->yk, tao->stepdirection, &dk_yk)); in TaoBNCGStepDirectionUpdate()
538 …PetscCall(TaoBNCGComputeScalarScaling(ynorm2, step * dk_yk, step * step * dnorm * dnorm, &tau_k, c… in TaoBNCGStepDirectionUpdate()
555 … PetscCall(TaoBNCGComputeScalarScaling(ynorm2, step * dk_yk, snorm * snorm, &tau_k, cg->alpha)); in TaoBNCGStepDirectionUpdate()
574 … PetscCall(TaoBNCGComputeScalarScaling(ynorm2, step * dk_yk, snorm * snorm, &tau_k, cg->alpha)); in TaoBNCGStepDirectionUpdate()
594 PetscCall(TaoBNCGComputeScalarScaling(ynorm2, step * dk_yk, cg->yts, &tau_k, cg->alpha)); in TaoBNCGStepDirectionUpdate()
602 PetscCall(VecDot(cg->d_work, cg->g_work, &dk_yk)); in TaoBNCGStepDirectionUpdate()
603 dk_yk = dk_yk - gd_old; in TaoBNCGStepDirectionUpdate()
604 beta = gtDg / dk_yk; in TaoBNCGStepDirectionUpdate()
617 cg->yts = step * dk_yk; in TaoBNCGStepDirectionUpdate()
626 tmp = gd / dk_yk; in TaoBNCGStepDirectionUpdate()
627 beta = tau_k * (gkp1_yk / dk_yk - ynorm2 * gd / (dk_yk * dk_yk)); in TaoBNCGStepDirectionUpdate()
643 tmp = gd / dk_yk; 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()
665 cg->yts = dk_yk * step; in TaoBNCGStepDirectionUpdate()
670 tmp = gd / dk_yk; in TaoBNCGStepDirectionUpdate()
671 …beta = tau_k * (gkp1_yk / dk_yk - ynorm2 * gd / (dk_yk * dk_yk) + gd / (dnorm * dnorm)) - step * g… in TaoBNCGStepDirectionUpdate()
685 tau_k = tau_k * gd / (dk_yk * dk_yk); in TaoBNCGStepDirectionUpdate()
686 tmp = gd / dk_yk; in TaoBNCGStepDirectionUpdate()
688 beta = gkp1_yk / dk_yk - step * tmp - tau_k; in TaoBNCGStepDirectionUpdate()
690 PetscCall(VecDot(cg->y_work, tao->stepdirection, &dk_yk)); in TaoBNCGStepDirectionUpdate()
691 beta += tmp * dk_yk / (dnorm * dnorm); /* projection of y_work onto dk */ in TaoBNCGStepDirectionUpdate()
707 cg->yts = dk_yk * step; in TaoBNCGStepDirectionUpdate()
715 beta = tau_k * (gkp1_yk / dk_yk - ynorm2 * gd / (dk_yk * dk_yk)) - step * gd / dk_yk; in TaoBNCGStepDirectionUpdate()
721 if (gkp1_yk < 0 && cg->neg_xi) gamma = -1.0 * gd / dk_yk; in TaoBNCGStepDirectionUpdate()
724 else gamma = cg->xi * gd / dk_yk; in TaoBNCGStepDirectionUpdate()
737 gamma = gd / dk_yk; 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()
747 if (gkp1D_yk / dk_yk < 0) gamma = -1.0 * gd / dk_yk; in TaoBNCGStepDirectionUpdate()
748 else gamma = cg->xi * gd / dk_yk; in TaoBNCGStepDirectionUpdate()
757 } else gamma = cg->xi * gd / dk_yk; in TaoBNCGStepDirectionUpdate()
772 cg->yts = dk_yk * step; in TaoBNCGStepDirectionUpdate()
778 tmp = gd / dk_yk; in TaoBNCGStepDirectionUpdate()
779 beta = tau_k * (gkp1_yk / dk_yk - cg->yty * gd / (dk_yk * dk_yk)) - step * tmp; in TaoBNCGStepDirectionUpdate()
789 gamma = gd / dk_yk; in TaoBNCGStepDirectionUpdate()
791 beta = (gkp1_yk / dk_yk - gd * tmp / (dk_yk * dk_yk)) - step * gd / dk_yk; in TaoBNCGStepDirectionUpdate()
801 cg->yts = dk_yk * step; in TaoBNCGStepDirectionUpdate()
810 beta = -step * gd / dk_yk; in TaoBNCGStepDirectionUpdate()
822 beta = -step * gd / dk_yk; in TaoBNCGStepDirectionUpdate()
832 cg->yts = step * dk_yk; in TaoBNCGStepDirectionUpdate()
837 …PetscCall(TaoBNCGComputeScalarScaling(cg->yty, step * dk_yk, snorm * snorm, &tau_bfgs, cg->bfgs_sc… in TaoBNCGStepDirectionUpdate()
838 …PetscCall(TaoBNCGComputeScalarScaling(cg->yty, step * dk_yk, snorm * snorm, &tau_dfp, cg->dfp_scal… in TaoBNCGStepDirectionUpdate()
843 tmp = cg->theta * tau_bfgs * gd / dk_yk + (1 - cg->theta) * tau_dfp * gkp1_yk / cg->yty; in TaoBNCGStepDirectionUpdate()
844 …beta = cg->theta * tau_bfgs * (gkp1_yk / dk_yk - cg->yty * gd / (dk_yk * dk_yk)) - step * gd / dk_… in TaoBNCGStepDirectionUpdate()
854 gamma = cg->theta * gd / dk_yk + (1 - cg->theta) * (gkp1_yk / tmp); in TaoBNCGStepDirectionUpdate()
856 beta = cg->theta * (gkp1_yk / dk_yk - gd * tmp / (dk_yk * dk_yk)) - step * gd / dk_yk; in TaoBNCGStepDirectionUpdate()