337   Vec               Up, G, Gp;  in FiniteDiff()  local
346   PetscCall(VecDuplicate(U, &Gp));  in FiniteDiff()
354   PetscCall(FormFunctionGradient(tao, Up, &f, Gp, ctx));  in FiniteDiff()
355   PetscCall(VecAXPY(Gp, -1, G));  in FiniteDiff()
356   PetscCall(VecScale(Gp, 1. / eps));  in FiniteDiff()
357   PetscCall(VecGetArray(Gp, &u));  in FiniteDiff()
360   PetscCall(VecRestoreArray(Gp, &u));  in FiniteDiff()
366   PetscCall(FormFunctionGradient(tao, Up, &f, Gp, ctx));  in FiniteDiff()
367   PetscCall(VecAXPY(Gp, -1, G));  in FiniteDiff()
368   PetscCall(VecScale(Gp, 1. / eps));  in FiniteDiff()
369   PetscCall(VecGetArray(Gp, &u));  in FiniteDiff()
372   PetscCall(VecRestoreArray(Gp, &u));  in FiniteDiff()
375   PetscCall(VecDestroy(&Gp));  in FiniteDiff()