127   if (!tao->stepdirection) PetscCall(VecDuplicate(tao->solution, &tao->stepdirection));  in TaoSetUp_BNCG()
407   PetscCall(VecAXPBY(tao->stepdirection, -scaling, 0.0, tao->gradient));  in TaoBNCGResetUpdate()
451     PetscCall(VecDot(cg->yk, tao->stepdirection, &dk_yk));  in TaoBNCGStepDirectionUpdate()
496     PetscCall(VecAXPBY(tao->stepdirection, -1.0, 0.0, cg->g_work));  in TaoBNCGStepDirectionUpdate()
508         PetscCall(VecAXPBY(tao->stepdirection, -tau_k, 0.0, tao->gradient));  in TaoBNCGStepDirectionUpdate()
511         PetscCall(VecAXPBY(tao->stepdirection, -1.0, 0.0, cg->g_work));  in TaoBNCGStepDirectionUpdate()
522         PetscCall(VecAXPBY(tao->stepdirection, -tau_k, beta, tao->gradient));  in TaoBNCGStepDirectionUpdate()
527         PetscCall(VecAXPBY(tao->stepdirection, -1.0, beta, cg->g_work));  in TaoBNCGStepDirectionUpdate()
536       PetscCall(VecDot(cg->yk, tao->stepdirection, &dk_yk));  in TaoBNCGStepDirectionUpdate()
540         PetscCall(VecAXPBY(tao->stepdirection, -tau_k, beta, tao->gradient));  in TaoBNCGStepDirectionUpdate()
546         PetscCall(VecAXPBY(tao->stepdirection, -1.0, beta, cg->g_work));  in TaoBNCGStepDirectionUpdate()
557         PetscCall(VecAXPBY(tao->stepdirection, -tau_k, beta, tao->gradient));  in TaoBNCGStepDirectionUpdate()
563         PetscCall(VecAXPBY(tao->stepdirection, -1.0, beta, cg->g_work));  in TaoBNCGStepDirectionUpdate()
577         PetscCall(VecAXPBY(tao->stepdirection, -tau_k, beta, tao->gradient));  in TaoBNCGStepDirectionUpdate()
584         PetscCall(VecAXPBY(tao->stepdirection, -1.0, beta, cg->g_work));  in TaoBNCGStepDirectionUpdate()
592         PetscCall(VecDot(tao->stepdirection, tao->gradient, &gd));  in TaoBNCGStepDirectionUpdate()
593         PetscCall(VecDot(cg->G_old, tao->stepdirection, &gd_old));  in TaoBNCGStepDirectionUpdate()
596         PetscCall(VecAXPBY(tao->stepdirection, -tau_k, beta, tao->gradient));  in TaoBNCGStepDirectionUpdate()
598         PetscCall(MatMult(cg->B, tao->stepdirection, cg->d_work));  in TaoBNCGStepDirectionUpdate()
601         PetscCall(VecDot(tao->stepdirection, cg->G_old, &gd_old));  in TaoBNCGStepDirectionUpdate()
606         PetscCall(VecWAXPY(tao->stepdirection, -1.0, cg->g_work, cg->d_work));  in TaoBNCGStepDirectionUpdate()
613       PetscCall(VecDot(tao->gradient, tao->stepdirection, &gd));  in TaoBNCGStepDirectionUpdate()
614       PetscCall(VecDot(cg->G_old, tao->stepdirection, &gd_old));  in TaoBNCGStepDirectionUpdate()
631           PetscCall(VecAXPBY(tao->stepdirection, -tau_k, beta, tao->gradient));  in TaoBNCGStepDirectionUpdate()
639           PetscCall(MatSolve(cg->B, tao->stepdirection, cg->d_work));  in TaoBNCGStepDirectionUpdate()
650           PetscCall(VecDot(tao->stepdirection, cg->g_work, &gd));  in TaoBNCGStepDirectionUpdate()
653           PetscCall(VecAXPBY(tao->stepdirection, -1.0, beta, cg->g_work));  in TaoBNCGStepDirectionUpdate()
661       PetscCall(VecDot(tao->gradient, tao->stepdirection, &gd));  in TaoBNCGStepDirectionUpdate()
662       PetscCall(VecDot(cg->G_old, tao->stepdirection, &gd_old));  in TaoBNCGStepDirectionUpdate()
674         PetscCall(VecAXPBYPCZ(tao->stepdirection, -tau_k, 0.0, beta, tao->gradient, cg->yk));  in TaoBNCGStepDirectionUpdate()
681         PetscCall(MatSolve(cg->B, tao->stepdirection, cg->d_work));  in TaoBNCGStepDirectionUpdate()
690         PetscCall(VecDot(cg->y_work, tao->stepdirection, &dk_yk));  in TaoBNCGStepDirectionUpdate()
692         PetscCall(VecDot(tao->stepdirection, cg->g_work, &gd));  in TaoBNCGStepDirectionUpdate()
696         PetscCall(VecAXPBY(tao->stepdirection, -1.0, beta, cg->g_work));  in TaoBNCGStepDirectionUpdate()
703       PetscCall(VecDot(tao->gradient, tao->stepdirection, &gd));  in TaoBNCGStepDirectionUpdate()
704       PetscCall(VecDot(cg->G_old, tao->stepdirection, &gd_old));  in TaoBNCGStepDirectionUpdate()
727 …    PetscCall(VecAXPBYPCZ(tao->stepdirection, -tau_k, gamma * tau_k, beta, tao->gradient, cg->yk));  in TaoBNCGStepDirectionUpdate()
744           PetscCall(VecDot(tao->stepdirection, cg->g_work, &tmp));  in TaoBNCGStepDirectionUpdate()
760           PetscCall(VecAXPBY(tao->stepdirection, -1.0, beta, cg->g_work));  in TaoBNCGStepDirectionUpdate()
761           PetscCall(VecAXPY(tao->stepdirection, gamma, cg->y_work));  in TaoBNCGStepDirectionUpdate()
769       PetscCall(VecDot(tao->gradient, tao->stepdirection, &gd));  in TaoBNCGStepDirectionUpdate()
781 …      PetscCall(VecAXPBYPCZ(tao->stepdirection, -tau_k, tmp * tau_k, beta, tao->gradient, cg->yk));  in TaoBNCGStepDirectionUpdate()
793         PetscCall(VecAXPBYPCZ(tao->stepdirection, -1.0, gamma, beta, cg->g_work, cg->y_work));  in TaoBNCGStepDirectionUpdate()
798       PetscCall(VecDot(tao->gradient, tao->stepdirection, &gd));  in TaoBNCGStepDirectionUpdate()
812         PetscCall(VecAXPBYPCZ(tao->stepdirection, -tau_k, tmp, beta, tao->gradient, cg->yk));  in TaoBNCGStepDirectionUpdate()
824         PetscCall(VecAXPBYPCZ(tao->stepdirection, -1.0, gamma, beta, cg->g_work, cg->y_work));  in TaoBNCGStepDirectionUpdate()
829       PetscCall(VecDot(tao->gradient, tao->stepdirection, &gd));  in TaoBNCGStepDirectionUpdate()
846         PetscCall(VecAXPBYPCZ(tao->stepdirection, -tau_k, tmp, beta, tao->gradient, cg->yk));  in TaoBNCGStepDirectionUpdate()
858         PetscCall(VecAXPBYPCZ(tao->stepdirection, -1.0, gamma, beta, cg->g_work, cg->y_work));  in TaoBNCGStepDirectionUpdate()
892 …ao->linesearch, tao->solution, &cg->f, cg->unprojected_gradient, tao->stepdirection, &step, &ls_st…  in TaoBNCGConductIteration()
917           PetscCall(TaoBNCGBoundStep(tao, cg->as_type, tao->stepdirection));  in TaoBNCGConductIteration()
920 …ao->linesearch, tao->solution, &cg->f, cg->unprojected_gradient, tao->stepdirection, &step, &ls_st…  in TaoBNCGConductIteration()
934           PetscCall(TaoBNCGBoundStep(tao, cg->as_type, tao->stepdirection));  in TaoBNCGConductIteration()
936 …ao->linesearch, tao->solution, &cg->f, cg->unprojected_gradient, tao->stepdirection, &step, &ls_st…  in TaoBNCGConductIteration()
974   PetscCall(VecNorm(tao->stepdirection, NORM_2, &dnorm));  in TaoBNCGConductIteration()
978   PetscCall(TaoBNCGBoundStep(tao, cg->as_type, tao->stepdirection));  in TaoBNCGConductIteration()
986       PetscCall(VecGetSubVector(tao->stepdirection, cg->new_inactives, &cg->inactive_step));  in TaoBNCGConductIteration()
990       PetscCall(VecRestoreSubVector(tao->stepdirection, cg->new_inactives, &cg->inactive_step));  in TaoBNCGConductIteration()
994     PetscCall(VecDot(tao->gradient, tao->stepdirection, &gd));  in TaoBNCGConductIteration()
995     PetscCall(VecNorm(tao->stepdirection, NORM_2, &dnorm));  in TaoBNCGConductIteration()
999       PetscCall(TaoBNCGBoundStep(tao, cg->as_type, tao->stepdirection));  in TaoBNCGConductIteration()