15     PetscCall(MatMult(gn->Hreg, in, gn->x_work));  in GNHessianProd()
16     PetscCall(VecAXPY(out, gn->lambda, gn->x_work));  in GNHessianProd()
33       PetscCall(MatMultTranspose(gn->D, gn->y_work, gn->x_work)); /* x_work = D'*(diag.*(D*in)) */  in GNHessianProd()
35       PetscCall(VecCopy(gn->y_work, gn->x_work));  in GNHessianProd()
37     PetscCall(VecAXPY(out, gn->lambda, gn->x_work));  in GNHessianProd()
40     PetscCall(VecPointwiseMult(gn->x_work, gn->damping, in));  in GNHessianProd()
41     PetscCall(VecAXPY(out, 1, gn->x_work));  in GNHessianProd()
116     PetscCall((*gn->regularizerobjandgrad)(tao, X, &f_reg, gn->x_work, gn->reg_obj_ctx));  in GNObjectiveGradientEval()
118     PetscCall(VecAXPY(G, gn->lambda, gn->x_work));  in GNObjectiveGradientEval()
129     PetscCall(VecAXPBYPCZ(gn->x_work, 1.0, -1.0, 0.0, X, gn->x_old));  in GNObjectiveGradientEval()
130     PetscCall(VecDot(gn->x_work, gn->x_work, &f_reg));  in GNObjectiveGradientEval()
151       PetscCall(MatMultTranspose(gn->D, gn->y_work, gn->x_work));  in GNObjectiveGradientEval()
153       PetscCall(VecCopy(gn->y_work, gn->x_work));  in GNObjectiveGradientEval()
155     PetscCall(VecAXPY(G, gn->lambda, gn->x_work));  in GNObjectiveGradientEval()
418   if (!gn->x_work) PetscCall(VecDuplicate(tao->solution, &gn->x_work));  in TaoSetUp_BRGN()
486     PetscCall(VecDestroy(&gn->x_work));  in TaoDestroy_BRGN()