436 …l gnorm2, PetscReal step, PetscReal fold, PetscReal gnorm2_old, PetscReal dnorm, PetscBool pcgd_fa…  in TaoBNCGStepDirectionUpdate()  argument
452     if (step * dnorm < PETSC_MACHINE_EPSILON || step * dk_yk < PETSC_MACHINE_EPSILON) {  in TaoBNCGStepDirectionUpdate()
502           cg->sts = step * step * dnorm * dnorm;  in TaoBNCGStepDirectionUpdate()
518         cg->sts = step * step * dnorm * dnorm;  in TaoBNCGStepDirectionUpdate()
538 …PetscCall(TaoBNCGComputeScalarScaling(ynorm2, step * dk_yk, step * step * dnorm * dnorm, &tau_k, c…  in TaoBNCGStepDirectionUpdate()
551       snorm = step * dnorm;  in TaoBNCGStepDirectionUpdate()
616       snorm   = dnorm * step;  in TaoBNCGStepDirectionUpdate()
629 …x(PetscMax(beta, cg->hz_eta * tau_k * gd_old / (dnorm * dnorm)), cg->dk_eta * tau_k * gd / (dnorm …  in TaoBNCGStepDirectionUpdate()
651 …ta = PetscMax(PetscMax(beta, cg->hz_eta * gd_old / (dnorm * dnorm)), cg->dk_eta * gd / (dnorm * dn…  in TaoBNCGStepDirectionUpdate()
664       snorm   = step * dnorm;  in TaoBNCGStepDirectionUpdate()
671 …beta = tau_k * (gkp1_yk / dk_yk - ynorm2 * gd / (dk_yk * dk_yk) + gd / (dnorm * dnorm)) - step * g…  in TaoBNCGStepDirectionUpdate()
672 …x(PetscMax(beta, cg->hz_eta * tau_k * gd_old / (dnorm * dnorm)), cg->dk_eta * tau_k * gd / (dnorm …  in TaoBNCGStepDirectionUpdate()
691         beta += tmp * dk_yk / (dnorm * dnorm); /* projection of y_work onto dk */  in TaoBNCGStepDirectionUpdate()
694 …ta = PetscMax(PetscMax(beta, cg->hz_eta * gd_old / (dnorm * dnorm)), cg->dk_eta * gd / (dnorm * dn…  in TaoBNCGStepDirectionUpdate()
706       snorm   = step * dnorm;  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()
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()
771       snorm   = step * dnorm;  in TaoBNCGStepDirectionUpdate()
800       snorm   = step * dnorm;  in TaoBNCGStepDirectionUpdate()
831       snorm   = step * dnorm;  in TaoBNCGStepDirectionUpdate()
873   PetscReal                    step = 1.0, gnorm2, gd, dnorm = 0.0;  in TaoBNCGConductIteration()  local
914 … PetscCall(TaoBNCGStepDirectionUpdate(tao, gnorm2, step, f_old, gnorm2_old, dnorm, pcgd_fallback));  in TaoBNCGConductIteration()
974   PetscCall(VecNorm(tao->stepdirection, NORM_2, &dnorm));  in TaoBNCGConductIteration()
976   PetscCall(TaoBNCGStepDirectionUpdate(tao, gnorm2, step, f_old, gnorm2_old, dnorm, pcgd_fallback));  in TaoBNCGConductIteration()
995     PetscCall(VecNorm(tao->stepdirection, NORM_2, &dnorm));  in TaoBNCGConductIteration()
996 …if (PetscIsInfOrNanReal(gd) || (gd / (dnorm * dnorm) <= -1e10 || gd / (dnorm * dnorm) >= -1e-10)) {  in TaoBNCGConductIteration()