31 …PetscCall(VecPointwiseMult(gn->y_work, gn->diag, gn->y)); /* y_work = diag.*(D*in), where diag = e…  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()
142     PetscCall(VecPointwiseMult(gn->y_work, gn->y, gn->y));  in GNObjectiveGradientEval()
143     PetscCall(VecShift(gn->y_work, gn->epsilon * gn->epsilon));  in GNObjectiveGradientEval()
144     PetscCall(VecSqrtAbs(gn->y_work)); /* gn->y_work = sqrt(y.^2+epsilon^2) */  in GNObjectiveGradientEval()
145     PetscCall(VecSum(gn->y_work, &yESum));  in GNObjectiveGradientEval()
149 …PetscCall(VecPointwiseDivide(gn->y_work, gn->y, gn->y_work)); /* reuse y_work = y./sqrt(y.^2+epsil…  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()
193     PetscCall(VecPointwiseMult(gn->y_work, gn->y, gn->y));  in GNComputeHessian()
194     PetscCall(VecShift(gn->y_work, gn->epsilon * gn->epsilon));  in GNComputeHessian()
195     PetscCall(VecCopy(gn->y_work, gn->diag));                    /* gn->diag = y.^2+epsilon^2 */  in GNComputeHessian()
196 …PetscCall(VecSqrtAbs(gn->y_work));                           /* gn->y_work = sqrt(y.^2+epsilon^2) …  in GNComputeHessian()
197 …PetscCall(VecPointwiseMult(gn->diag, gn->y_work, gn->diag)); /* gn->diag = sqrt(y.^2+epsilon^2).^3…  in GNComputeHessian()
435     if (!gn->y_work) PetscCall(VecDuplicate(gn->y, &gn->y_work));  in TaoSetUp_BRGN()
491     PetscCall(VecDestroy(&gn->y_work));  in TaoDestroy_BRGN()