52 …PetscInt degree = 3, q_extra, *l_size, *xl_size, *g_size, dim = 3, fine_level, mesh_elem[3] = {3, …  in main()  local
156   fine_level = num_levels - 1;  in main()
164       level_degrees[fine_level] = degree;  in main()
193 …PetscCall(SetupDMByDegree(dm[i], level_degrees[fine_level], q_extra, num_comp_u, dim, bp_options[b…  in main()
219   PetscCall(VecDuplicate(X[fine_level], &rhs));  in main()
247 …ype, used_resource, CeedMemTypes[mem_type_backend], P, Q, q_extra, g_size[fine_level] / num_comp_u,  in main()
248 …          l_size[fine_level] / num_comp_u, num_comp_u, CeedElemTopologies[elem_topo], num_levels));  in main()
252   PetscCall(VecDuplicate(X_loc[fine_level], &rhs_loc));  in main()
254   CeedVectorCreate(ceed, xl_size[fine_level], &rhs_ceed);  in main()
261     if (!test_mode && (i == 0 || i == fine_level)) {  in main()
271 …                               ceed_data[i], i == fine_level, i == fine_level, rhs_ceed, &target));  in main()
277   PetscCall(DMLocalToGlobal(dm[fine_level], rhs_loc, ADD_VALUES, rhs));  in main()
284   CeedQFunctionAddInput(qf_error, "qdata", ceed_data[fine_level]->q_data_size, CEED_EVAL_NONE);  in main()
289 …CeedOperatorSetField(op_error, "u", ceed_data[fine_level]->elem_restr_u, ceed_data[fine_level]->ba…  in main()
290 …CeedOperatorSetField(op_error, "true_soln", ceed_data[fine_level]->elem_restr_u_i, CEED_BASIS_NONE…  in main()
291 …ratorSetField(op_error, "qdata", ceed_data[fine_level]->elem_restr_qd_i, CEED_BASIS_NONE, ceed_dat…  in main()
292 …CeedOperatorSetField(op_error, "error", ceed_data[fine_level]->elem_restr_u, ceed_data[fine_level]…  in main()
326   for (PetscInt i = fine_level; i >= 0; i--) {  in main()
387   PetscCall(KSPSetOperators(ksp, mat_O[fine_level], mat_O[fine_level]));  in main()
446     PetscCall(VecZeroEntries(X[fine_level]));  in main()
448     PetscCall(KSPSolve(ksp, rhs, X[fine_level]));  in main()
460   PetscCall(VecZeroEntries(X[fine_level]));  in main()
469   PetscCall(KSPSolve(ksp, rhs, X[fine_level]));  in main()
507 …PetscCall(SetupErrorOperatorCtx(comm, dm[fine_level], ceed, ceed_data[fine_level], X_loc[fine_leve…  in main()
509       PetscCall(ComputeL2Error(X[fine_level], &l2_error, op_error_ctx));  in main()
521 …ec in CG                            : %g (%g) million\n", 1e-6 * g_size[fine_level] * its / rt_max,  in main()
522                             1e-6 * g_size[fine_level] * its / rt_min));  in main()
532     PetscCall(VecView(X[fine_level], vtk_viewer_soln));  in main()