137     PetscCall(VecDuplicate(tao->solution, &cg->g_work));  in TaoSetUp_BNCG()
158     PetscCall(VecDestroy(&cg->g_work));  in TaoDestroy_BNCG()
495     PetscCall(MatSolve(cg->B, tao->gradient, cg->g_work));  in TaoBNCGStepDirectionUpdate()
496     PetscCall(VecAXPBY(tao->stepdirection, -1.0, 0.0, cg->g_work));  in TaoBNCGStepDirectionUpdate()
510         PetscCall(MatSolve(cg->B, tao->gradient, cg->g_work));  in TaoBNCGStepDirectionUpdate()
511         PetscCall(VecAXPBY(tao->stepdirection, -1.0, 0.0, cg->g_work));  in TaoBNCGStepDirectionUpdate()
524         PetscCall(MatSolve(cg->B, tao->gradient, cg->g_work));  in TaoBNCGStepDirectionUpdate()
525         PetscCall(VecDot(cg->yk, cg->g_work, &gkp1_yk));  in TaoBNCGStepDirectionUpdate()
527         PetscCall(VecAXPBY(tao->stepdirection, -1.0, beta, cg->g_work));  in TaoBNCGStepDirectionUpdate()
542         PetscCall(VecDot(cg->G_old, cg->g_work, &gnorm2_old)); /* Before it's updated */  in TaoBNCGStepDirectionUpdate()
543         PetscCall(MatSolve(cg->B, tao->gradient, cg->g_work));  in TaoBNCGStepDirectionUpdate()
544         PetscCall(VecDot(tao->gradient, cg->g_work, &tmp));  in TaoBNCGStepDirectionUpdate()
546         PetscCall(VecAXPBY(tao->stepdirection, -1.0, beta, cg->g_work));  in TaoBNCGStepDirectionUpdate()
559         PetscCall(VecDot(cg->G_old, cg->g_work, &gnorm2_old));  in TaoBNCGStepDirectionUpdate()
560         PetscCall(MatSolve(cg->B, tao->gradient, cg->g_work));  in TaoBNCGStepDirectionUpdate()
561         PetscCall(VecDot(cg->g_work, cg->yk, &gkp1_yk));  in TaoBNCGStepDirectionUpdate()
563         PetscCall(VecAXPBY(tao->stepdirection, -1.0, beta, cg->g_work));  in TaoBNCGStepDirectionUpdate()
579         PetscCall(VecDot(cg->G_old, cg->g_work, &gnorm2_old)); /* Old gtDg */  in TaoBNCGStepDirectionUpdate()
580         PetscCall(MatSolve(cg->B, tao->gradient, cg->g_work));  in TaoBNCGStepDirectionUpdate()
581         PetscCall(VecDot(cg->g_work, cg->yk, &gkp1_yk));  in TaoBNCGStepDirectionUpdate()
584         PetscCall(VecAXPBY(tao->stepdirection, -1.0, beta, cg->g_work));  in TaoBNCGStepDirectionUpdate()
599         PetscCall(MatSolve(cg->B, tao->gradient, cg->g_work));  in TaoBNCGStepDirectionUpdate()
600         PetscCall(VecDot(cg->g_work, tao->gradient, &gtDg));  in TaoBNCGStepDirectionUpdate()
602         PetscCall(VecDot(cg->d_work, cg->g_work, &dk_yk));  in TaoBNCGStepDirectionUpdate()
606         PetscCall(VecWAXPY(tao->stepdirection, -1.0, cg->g_work, cg->d_work));  in TaoBNCGStepDirectionUpdate()
637           PetscCall(MatSolve(cg->B, tao->gradient, cg->g_work));  in TaoBNCGStepDirectionUpdate()
641           PetscCall(VecDot(cg->yk, cg->g_work, &gkp1_yk));  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()
679         PetscCall(MatSolve(cg->B, tao->gradient, cg->g_work));  in TaoBNCGStepDirectionUpdate()
683         PetscCall(VecDot(cg->yk, cg->g_work, &gkp1_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()
732           PetscCall(MatSolve(cg->B, tao->gradient, cg->g_work));  in TaoBNCGStepDirectionUpdate()
735           PetscCall(VecDot(cg->yk, cg->g_work, &gkp1D_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()
784         PetscCall(MatSolve(cg->B, tao->gradient, cg->g_work));  in TaoBNCGStepDirectionUpdate()
787         PetscCall(VecDot(cg->g_work, cg->yk, &gkp1_yk));  in TaoBNCGStepDirectionUpdate()
793         PetscCall(VecAXPBYPCZ(tao->stepdirection, -1.0, gamma, beta, cg->g_work, cg->y_work));  in TaoBNCGStepDirectionUpdate()
815         PetscCall(MatSolve(cg->B, tao->gradient, cg->g_work));  in TaoBNCGStepDirectionUpdate()
818         PetscCall(VecDot(cg->g_work, cg->yk, &gkp1_yk));  in TaoBNCGStepDirectionUpdate()
824         PetscCall(VecAXPBYPCZ(tao->stepdirection, -1.0, gamma, beta, cg->g_work, cg->y_work));  in TaoBNCGStepDirectionUpdate()
849         PetscCall(MatSolve(cg->B, tao->gradient, cg->g_work));  in TaoBNCGStepDirectionUpdate()
852         PetscCall(VecDot(cg->g_work, cg->yk, &gkp1_yk));  in TaoBNCGStepDirectionUpdate()
858         PetscCall(VecAXPBYPCZ(tao->stepdirection, -1.0, gamma, beta, cg->g_work, cg->y_work));  in TaoBNCGStepDirectionUpdate()