4 #include "t591-operator.h"
15 …             num_nodes_coarse = (num_elem_1d * (p_coarse - 1) + 1) * (num_elem_1d * (p_coarse - 1)…  in main()
16 …Int             num_nodes_fine   = (num_elem_1d * (p_fine - 1) + 1) * (num_elem_1d * (p_fine - 1) …  in main()
30     for (CeedInt e = 0; e < num_elem; e++) {  in main()  local
32         x_array[num_points_per_elem * (e * dim + d) + 0] = 0.25;  in main()
33         x_array[num_points_per_elem * (e * dim + d) + 1] = d == 0 ? -0.25 : 0.25;  in main()
34         x_array[num_points_per_elem * (e * dim + d) + 2] = d == 0 ? 0.25 : -0.25;  in main()
35         x_array[num_points_per_elem * (e * dim + d) + 3] = 0.25;  in main()
55     CeedInt p = 2, num_nodes = (num_elem_1d * (p - 1) + 1) * (num_elem_1d * (p - 1) + 1);  in main()
58     for (CeedInt e = 0; e < num_elem; e++) {  in main()  local
61       for (CeedInt d = 0; d < dim; d++) n_d[d] = num_elem_1d * (p - 1) + 1;  in main()
63         CeedInt r_e = e;  in main()
70       CeedInt num_nodes_in_elem = p * p, *elem_nodes = ind_x + e * num_nodes_in_elem;  in main()
76           g_node += (elem_xy[d] * (p - 1) + r_node % p) * g_node_stride;  in main()
104     for (CeedInt e = 0; e < num_elem; e++) {  in main()  local
107       for (CeedInt d = 0; d < dim; d++) n_d[d] = num_elem_1d * (p_coarse - 1) + 1;  in main()
109         CeedInt r_e = e;  in main()
116       CeedInt num_nodes_in_elem = p_coarse * p_coarse, *elem_nodes = ind_u + e * num_nodes_in_elem;  in main()
122           g_node += (elem_xy[d] * (p_coarse - 1) + r_node % p_coarse) * g_node_stride;  in main()
136     for (CeedInt e = 0; e < num_elem; e++) {  in main()  local
139       for (CeedInt d = 0; d < dim; d++) n_d[d] = num_elem_1d * (p_fine - 1) + 1;  in main()
141         CeedInt r_e = e;  in main()
148       CeedInt num_nodes_in_elem = p_fine * p_fine, *elem_nodes = ind_u + e * num_nodes_in_elem;  in main()
154           g_node += (elem_xy[d] * (p_fine - 1) + r_node % p_fine) * g_node_stride;  in main()
216 …if (fabs(sum - num_elem) > 1000. * CEED_EPSILON) printf("Computed Area Coarse Grid: %f != True Are…  in main()
236 …if (fabs(sum - num_elem) > 1000. * CEED_EPSILON) printf("Computed Area Fine Grid: %f != True Area:…  in main()
251 …if (fabs(sum - num_elem) > 1000. * CEED_EPSILON) printf("Computed Area Coarse Grid: %f != True Are…  in main()