30 …CeedInt       num_elem_tet = 6, p_tet = 6, q_tet = 4, num_elem_hex = 6, p_hex = 3, q_hex = 4, dim …  in main()  local
35   CeedScalar    q_ref[dim * q_tet], q_weight[q_tet];  in main()
36   CeedScalar    interp[p_tet * q_tet], grad[dim * p_tet * q_tet];  in main()
41   CeedVectorCreate(ceed, dim * num_dofs, &x);  in main()
43     CeedScalar x_array[dim * num_dofs];  in main()
55   CeedVectorCreate(ceed, num_qpts_tet * dim * (dim + 1) / 2, &q_data_tet);  in main()
56   CeedVectorCreate(ceed, num_qpts_hex * dim * (dim + 1) / 2, &q_data_hex);  in main()
79 …CeedElemRestrictionCreate(ceed, num_elem_tet, p_tet, dim, num_dofs, dim * num_dofs, CEED_MEM_HOST,…  in main()
83   CeedInt strides_q_data_tet[3] = {1, q_tet, q_tet * dim * (dim + 1) / 2};  in main()
84 …CeedElemRestrictionCreateStrided(ceed, num_elem_tet, q_tet, dim * (dim + 1) / 2, dim * (dim + 1) /…  in main()
89 …CeedBasisCreateH1(ceed, CEED_TOPOLOGY_TRIANGLE, dim, p_tet, q_tet, interp, grad, q_ref, q_weight, …  in main()
97   CeedQFunctionAddInput(qf_setup_tet, "dx", dim * dim, CEED_EVAL_GRAD);  in main()
98   CeedQFunctionAddOutput(qf_setup_tet, "rho", dim * (dim + 1) / 2, CEED_EVAL_NONE);  in main()
101   CeedQFunctionAddInput(qf_diff_tet, "rho", dim * (dim + 1) / 2, CEED_EVAL_NONE);  in main()
102   CeedQFunctionAddInput(qf_diff_tet, "u", dim, CEED_EVAL_GRAD);  in main()
103   CeedQFunctionAddOutput(qf_diff_tet, "v", dim, CEED_EVAL_GRAD);  in main()
127 …CeedElemRestrictionCreate(ceed, num_elem_hex, p_hex * p_hex, dim, num_dofs, dim * num_dofs, CEED_M…  in main()
131   CeedInt strides_q_data_hex[3] = {1, q_hex * q_hex, q_hex * q_hex * dim * (dim + 1) / 2};  in main()
132 …mRestrictionCreateStrided(ceed, num_elem_hex, q_hex * q_hex, dim * (dim + 1) / 2, dim * (dim + 1) …  in main()
136   CeedBasisCreateTensorH1Lagrange(ceed, dim, dim, p_hex, q_hex, CEED_GAUSS, &basis_x_hex);  in main()
137   CeedBasisCreateTensorH1Lagrange(ceed, dim, 1, p_hex, q_hex, CEED_GAUSS, &basis_u_hex);  in main()
142   CeedQFunctionAddInput(qf_setup_hex, "dx", dim * dim, CEED_EVAL_GRAD);  in main()
143   CeedQFunctionAddOutput(qf_setup_hex, "rho", dim * (dim + 1) / 2, CEED_EVAL_NONE);  in main()
146   CeedQFunctionAddInput(qf_diff_hex, "rho", dim * (dim + 1) / 2, CEED_EVAL_NONE);  in main()
147   CeedQFunctionAddInput(qf_diff_hex, "u", dim, CEED_EVAL_GRAD);  in main()
148   CeedQFunctionAddOutput(qf_diff_hex, "v", dim, CEED_EVAL_GRAD);  in main()