22     CeedInt             p = 3, q = 4, dim = 2;  in main()  local
31     CeedVectorCreate(ceed, dim * num_points, &x_points);  in main()
33       CeedScalar x_array[dim * num_points];  in main()
36         for (CeedInt d = 0; d < dim; d++) {  in main()
37           x_array[num_points_per_elem * (e * dim + d) + 0] = 0.25;  in main()
38           x_array[num_points_per_elem * (e * dim + d) + 1] = d == 0 ? -0.25 : 0.25;  in main()
39           x_array[num_points_per_elem * (e * dim + d) + 2] = d == 0 ? 0.25 : -0.25;  in main()
40           x_array[num_points_per_elem * (e * dim + d) + 3] = 0.25;  in main()
50 …CeedElemRestrictionCreateAtPoints(ceed, num_elem, num_points, dim, num_points * dim, CEED_MEM_HOST…  in main()
52 …RestrictionCreateAtPoints(ceed, num_elem, num_points, dim * (dim + 1) / 2, num_points * dim * (dim…  in main()
57     CeedVectorCreate(ceed, dim * num_dofs, &x);  in main()
59       CeedScalar x_array[dim * num_dofs];  in main()
71     CeedVectorCreate(ceed, num_points * dim * (dim + 1) / 2, &q_data);  in main()
83 …CeedElemRestrictionCreate(ceed, num_elem, p * p, dim, num_dofs, dim * num_dofs, CEED_MEM_HOST, CEE…  in main()
88     CeedBasisCreateTensorH1Lagrange(ceed, dim, dim, p, q, CEED_GAUSS, &basis_x);  in main()
89     CeedBasisCreateTensorH1Lagrange(ceed, dim, num_comp, p, q, CEED_GAUSS, &basis_u);  in main()
93     CeedQFunctionAddInput(qf_setup, "dx", dim * dim, CEED_EVAL_GRAD);  in main()
95     CeedQFunctionAddOutput(qf_setup, "q data", dim * (dim + 1) / 2, CEED_EVAL_NONE);  in main()
109     CeedQFunctionAddInput(qf_diff, "du", num_comp * dim, CEED_EVAL_GRAD);  in main()
110     CeedQFunctionAddInput(qf_diff, "q data", dim * (dim + 1) / 2, CEED_EVAL_NONE);  in main()
111     CeedQFunctionAddOutput(qf_diff, "dv", num_comp * dim, CEED_EVAL_GRAD);  in main()