19   CeedInt             num_elem = 1, p = 4, q = 5, dim = 2;  in main()  local
20 …CeedInt             num_dofs = p * p, num_qpts = num_elem * q * q, q_data_size = dim * (dim + 1) /…  in main()
49   CeedInt strides_u[3] = {1, p * p, p * p};  in main()
50 …CeedElemRestrictionCreateStrided(ceed, num_elem, p * p, 1, num_dofs, strides_u, &elem_restriction_…  in main()
57   CeedBasisCreateTensorH1Lagrange(ceed, dim, 1, p, q, CEED_GAUSS, &basis_u);  in main()
101         u_array[p - 1] = 1.0;  in main()
104         u_array[p * p - p] = 1.0;  in main()
107         u_array[p * p - 1] = 1.0;  in main()
118     S[1 * 4 + i] = -v_array[p - 1];  in main()
119     S[2 * 4 + i] = -v_array[p * p - p];  in main()
120     S[3 * 4 + i] = -v_array[p * p - 1];  in main()
165     CeedScalar  nodes[p];  in main()
168     CeedLobattoQuadrature(p, nodes, NULL);  in main()
170     for (CeedInt i = 0; i < p; i++) {  in main()
171 …for (CeedInt j = 0; j < p; j++) u_array[i * p + j] = -(nodes[i] - 1.0) * (nodes[i] + 1.0) - (nodes…  in main()
186     v_array[p - 1]     = 0.0;  in main()
187     v_array[p * p - p] = 0.0;  in main()
188     v_array[p * p - 1] = 0.0;  in main()
200     w_Pi[1] = w_array[p - 1];  in main()
201     w_Pi[2] = w_array[p * p - p];  in main()
202     w_Pi[3] = w_array[p * p - 1];  in main()
218     v_array[p - 1]     = v_Pi[1];  in main()
219     v_array[p * p - p] = v_Pi[2];  in main()
220     v_array[p * p - 1] = v_Pi[3];  in main()
232     for (CeedInt i = 0; i < p; i++) {  in main()
233       for (CeedInt j = 0; j < p; j++) {  in main()
234         if (fabs(u_array[i * p + j] - w_array[i * p + j]) > 2e-3) {  in main()
236 … ", %" CeedInt_FMT "] Error in inverse: %e != %e\n", i, j, w_array[i * p + j], u_array[i * p + j]);  in main()