151   PetscFE   u, divu, u_sum, divu_sum;  in SetupDiscretization()  local
171 …teDefault(PetscObjectComm((PetscObject)mesh_sum), dim, dim, simplex, "velocity_sum_", -1, &u_sum));  in SetupDiscretization()
172   PetscCall(PetscObjectSetName((PetscObject)u_sum, "velocity_sum"));  in SetupDiscretization()
179   PetscCall(PetscFECopyQuadrature(u_sum, divu_sum));  in SetupDiscretization()
187   PetscCall(DMSetField(mesh_sum, 0, NULL, (PetscObject)u_sum));  in SetupDiscretization()
206   PetscCall(PetscFEDestroy(&u_sum));  in SetupDiscretization()
218   Vec             u, u_sum;  in main()  local
248   PetscCall(DMCreateGlobalVector(dm_sum, &u_sum));  in main()
249   PetscCall(VecSet(u_sum, 0.0));  in main()
250   PetscCall(PetscObjectSetName((PetscObject)u_sum, "solution_sum"));  in main()
253   PetscCall(DMSNESCheckFromOptions(snes_sum, u_sum));  in main()
254   PetscCall(SNESSolve(snes_sum, NULL, u_sum));  in main()
255   PetscCall(SNESGetSolution(snes_sum, &u_sum));  in main()
256   PetscCall(VecViewFromOptions(u_sum, NULL, "-solution_sum_view"));  in main()
259   PetscCall(VecAXPY(u_sum, -1, u));  in main()
260   PetscCall(VecNorm(u_sum, NORM_2, &errNorm));  in main()
264   PetscCall(VecDestroy(&u_sum));  in main()