131   if (!cg->yk) PetscCall(VecDuplicate(tao->gradient, &cg->yk));  in TaoSetUp_BNCG()
141   PetscCall(MatLMVMAllocate(cg->B, cg->sk, cg->yk));  in TaoSetUp_BNCG()
162     PetscCall(VecDestroy(&cg->yk));  in TaoDestroy_BNCG()
448     PetscCall(VecWAXPY(cg->yk, -1.0, cg->G_old, tao->gradient));  in TaoBNCGStepDirectionUpdate()
449     PetscCall(VecNorm(cg->yk, NORM_2, &ynorm));  in TaoBNCGStepDirectionUpdate()
451     PetscCall(VecDot(cg->yk, tao->stepdirection, &dk_yk));  in TaoBNCGStepDirectionUpdate()
519         PetscCall(VecDot(cg->yk, tao->gradient, &gkp1_yk));  in TaoBNCGStepDirectionUpdate()
525         PetscCall(VecDot(cg->yk, cg->g_work, &gkp1_yk));  in TaoBNCGStepDirectionUpdate()
533       PetscCall(VecWAXPY(cg->yk, -1.0, cg->G_old, tao->gradient));  in TaoBNCGStepDirectionUpdate()
534       PetscCall(VecNorm(cg->yk, NORM_2, &ynorm));  in TaoBNCGStepDirectionUpdate()
536       PetscCall(VecDot(cg->yk, tao->stepdirection, &dk_yk));  in TaoBNCGStepDirectionUpdate()
554         PetscCall(VecDot(cg->yk, tao->gradient, &gkp1_yk));  in TaoBNCGStepDirectionUpdate()
561         PetscCall(VecDot(cg->g_work, cg->yk, &gkp1_yk));  in TaoBNCGStepDirectionUpdate()
568       PetscCall(VecWAXPY(cg->yk, -1.0, cg->G_old, tao->gradient));  in TaoBNCGStepDirectionUpdate()
569       PetscCall(VecNorm(cg->yk, NORM_2, &ynorm));  in TaoBNCGStepDirectionUpdate()
573         PetscCall(VecDot(cg->yk, tao->gradient, &gkp1_yk));  in TaoBNCGStepDirectionUpdate()
581         PetscCall(VecDot(cg->g_work, cg->yk, &gkp1_yk));  in TaoBNCGStepDirectionUpdate()
623           PetscCall(VecDot(cg->yk, tao->gradient, &gkp1_yk));  in TaoBNCGStepDirectionUpdate()
638           PetscCall(MatSolve(cg->B, cg->yk, cg->y_work));  in TaoBNCGStepDirectionUpdate()
641           PetscCall(VecDot(cg->yk, cg->g_work, &gkp1_yk));  in TaoBNCGStepDirectionUpdate()
644           PetscCall(VecDot(cg->yk, cg->y_work, &tau_k));  in TaoBNCGStepDirectionUpdate()
667         PetscCall(VecDot(cg->yk, tao->gradient, &gkp1_yk));  in TaoBNCGStepDirectionUpdate()
674         PetscCall(VecAXPBYPCZ(tao->stepdirection, -tau_k, 0.0, beta, tao->gradient, cg->yk));  in TaoBNCGStepDirectionUpdate()
680         PetscCall(MatSolve(cg->B, cg->yk, cg->y_work));  in TaoBNCGStepDirectionUpdate()
683         PetscCall(VecDot(cg->yk, cg->g_work, &gkp1_yk));  in TaoBNCGStepDirectionUpdate()
684         PetscCall(VecDot(cg->yk, cg->y_work, &tau_k));  in TaoBNCGStepDirectionUpdate()
713           PetscCall(VecDot(cg->yk, tao->gradient, &gkp1_yk));  in TaoBNCGStepDirectionUpdate()
727 …    PetscCall(VecAXPBYPCZ(tao->stepdirection, -tau_k, gamma * tau_k, beta, tao->gradient, cg->yk));  in TaoBNCGStepDirectionUpdate()
733           PetscCall(MatSolve(cg->B, cg->yk, cg->y_work));  in TaoBNCGStepDirectionUpdate()
735           PetscCall(VecDot(cg->yk, cg->g_work, &gkp1D_yk));  in TaoBNCGStepDirectionUpdate()
739           PetscCall(VecDot(cg->yk, cg->y_work, &tau_k));  in TaoBNCGStepDirectionUpdate()
776         PetscCall(VecDot(cg->yk, tao->gradient, &gkp1_yk));  in TaoBNCGStepDirectionUpdate()
781 …      PetscCall(VecAXPBYPCZ(tao->stepdirection, -tau_k, tmp * tau_k, beta, tao->gradient, cg->yk));  in TaoBNCGStepDirectionUpdate()
785         PetscCall(MatSolve(cg->B, cg->yk, cg->y_work));  in TaoBNCGStepDirectionUpdate()
787         PetscCall(VecDot(cg->g_work, cg->yk, &gkp1_yk));  in TaoBNCGStepDirectionUpdate()
788         PetscCall(VecDot(cg->y_work, cg->yk, &tmp));  in TaoBNCGStepDirectionUpdate()
807         PetscCall(VecDot(cg->yk, tao->gradient, &gkp1_yk));  in TaoBNCGStepDirectionUpdate()
812         PetscCall(VecAXPBYPCZ(tao->stepdirection, -tau_k, tmp, beta, tao->gradient, cg->yk));  in TaoBNCGStepDirectionUpdate()
816         PetscCall(MatSolve(cg->B, cg->yk, cg->y_work));  in TaoBNCGStepDirectionUpdate()
818         PetscCall(VecDot(cg->g_work, cg->yk, &gkp1_yk));  in TaoBNCGStepDirectionUpdate()
819         PetscCall(VecDot(cg->y_work, cg->yk, &tmp));  in TaoBNCGStepDirectionUpdate()
839         PetscCall(VecDot(cg->yk, tao->gradient, &gkp1_yk));  in TaoBNCGStepDirectionUpdate()
846         PetscCall(VecAXPBYPCZ(tao->stepdirection, -tau_k, tmp, beta, tao->gradient, cg->yk));  in TaoBNCGStepDirectionUpdate()
850         PetscCall(MatSolve(cg->B, cg->yk, cg->y_work));  in TaoBNCGStepDirectionUpdate()
852         PetscCall(VecDot(cg->g_work, cg->yk, &gkp1_yk));  in TaoBNCGStepDirectionUpdate()
853         PetscCall(VecDot(cg->y_work, cg->yk, &tmp));  in TaoBNCGStepDirectionUpdate()