54   Vec           U, *U_g, *U_loc;     // U: solution, R: residual, F: forcing  in main()  local
203   PetscCall(VecDuplicate(U_g[fine_level], &U));  in main()
234     PetscCall(VecDuplicate(U, &neumann_bcs));  in main()
329     PetscCall(VecGetType(U, &vecType));  in main()
579   PetscCall(PetscObjectSetName((PetscObject)U, ""));  in main()
580   PetscCall(VecSet(U, 0.0));  in main()
584     PetscCall(ViewSolution(comm, app_ctx, U, 0, 0.0));  in main()
618     PetscCall(SNESSolve(snes, F, U));  in main()
622       PetscCall(ViewSolution(comm, app_ctx, U, increment, load_increment));  in main()
639 …  PetscCall(ComputeStrainEnergy(dm_energy, res_ctx, ceed_data[fine_level]->op_energy, U, &energy));  in main()
742     PetscCall(VecDuplicate(U, &error_vec));  in main()
744     PetscCall(VecDuplicate(U, &true_vec));  in main()
755     PetscCall(VecWAXPY(error_vec, -1.0, U, true_vec));  in main()
757     PetscCall(VecNorm(U, NORM_2, &l2_U_norm));  in main()
774   PetscCall(ComputeStrainEnergy(dm_energy, res_ctx, ceed_data[fine_level]->op_energy, U, &energy));  in main()
793 …PetscCall(ViewDiagnosticQuantities(comm, level_dms[fine_level], diagnostic_ctx, app_ctx, U, ceed_d…  in main()
848   PetscCall(VecDestroy(&U));  in main()