9 #include "t591-operator.h"
14 …CeedInt num_nodes = (num_elem_1d * (p - 1) + 1) * (num_elem_1d * (p - 1) + 1), num_points_per_elem…  in main()
29     for (CeedInt e = 0; e < num_elem; e++) {  in main()  local
31         x_array[num_points_per_elem * (e * dim + d) + 0] = 0.25;  in main()
32         x_array[num_points_per_elem * (e * dim + d) + 1] = d == 0 ? -0.25 : 0.25;  in main()
33         x_array[num_points_per_elem * (e * dim + d) + 2] = d == 0 ? 0.25 : -0.25;  in main()
34         x_array[num_points_per_elem * (e * dim + d) + 3] = 0.25;  in main()
54     CeedInt p = 2, num_nodes = (num_elem_1d * (p - 1) + 1) * (num_elem_1d * (p - 1) + 1);  in main()
57     for (CeedInt e = 0; e < num_elem; e++) {  in main()  local
60       for (CeedInt d = 0; d < dim; d++) n_d[d] = num_elem_1d * (p - 1) + 1;  in main()
62         CeedInt r_e = e;  in main()
69       CeedInt num_nodes_in_elem = p * p, *elem_nodes = ind_x + e * num_nodes_in_elem;  in main()
75           g_node += (elem_xy[d] * (p - 1) + r_node % p) * g_node_stride;  in main()
103     for (CeedInt e = 0; e < num_elem; e++) {  in main()  local
106       for (CeedInt d = 0; d < dim; d++) n_d[d] = num_elem_1d * (p - 1) + 1;  in main()
108         CeedInt r_e = e;  in main()
115       CeedInt num_nodes_in_elem = p * p, *elem_nodes = ind_u + e * num_nodes_in_elem;  in main()
121           g_node += (elem_xy[d] * (p - 1) + r_node % p) * g_node_stride;  in main()
177       if (fabs(q_data_array[i] - assembled_array[i]) > 1e-9) {  in main()
198 …if (fabs(area - 1.0 * num_elem) > CEED_EPSILON * 5e3) printf("Error: True operator computed area =…  in main()
222 …if (fabs(area - 1.0 * num_elem) > CEED_EPSILON * 5e3) printf("Error: Linearized operator computed …  in main()