439   PetscReal gamma = 1.0, tau_k, beta;  in TaoBNCGStepDirectionUpdate()  local
719             gamma = 0.0;  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()
727 …    PetscCall(VecAXPBYPCZ(tao->stepdirection, -tau_k, gamma * tau_k, beta, tao->gradient, cg->yk));  in TaoBNCGStepDirectionUpdate()
737           gamma = gd / 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()
751               gamma = 0.0;  in TaoBNCGStepDirectionUpdate()
756               gamma = 0.0;  in TaoBNCGStepDirectionUpdate()
757             } else gamma = cg->xi * gd / dk_yk;  in TaoBNCGStepDirectionUpdate()
761           PetscCall(VecAXPY(tao->stepdirection, gamma, cg->y_work));  in TaoBNCGStepDirectionUpdate()
789         gamma = gd / dk_yk;  in TaoBNCGStepDirectionUpdate()
793         PetscCall(VecAXPBYPCZ(tao->stepdirection, -1.0, gamma, beta, cg->g_work, cg->y_work));  in TaoBNCGStepDirectionUpdate()
820         gamma = (gkp1_yk / tmp);  in TaoBNCGStepDirectionUpdate()
824         PetscCall(VecAXPBYPCZ(tao->stepdirection, -1.0, gamma, beta, cg->g_work, cg->y_work));  in TaoBNCGStepDirectionUpdate()
854         gamma = cg->theta * gd / dk_yk + (1 - cg->theta) * (gkp1_yk / tmp);  in TaoBNCGStepDirectionUpdate()
858         PetscCall(VecAXPBYPCZ(tao->stepdirection, -1.0, gamma, beta, cg->g_work, cg->y_work));  in TaoBNCGStepDirectionUpdate()