10 #include "t502-operator.h"
17 …_elem + 1, num_dofs_u_coarse = num_elem * (p_coarse - 1) + 1, num_dofs_u_fine = num_elem * (p_fine…  in main()
27     for (CeedInt i = 0; i < num_dofs_x; i++) x_array[i] = (CeedScalar)i / (num_dofs_x - 1);  in main()
49     // -- QData  in main()
52     // -- Restrictions  in main()
60     offset = num_elem_sub * i * (p_fine - 1);  in main()
63         ind_u_fine[p_fine * j + k] = offset + j * (p_fine - 1) + k;  in main()
72     // -- Bases  in main()
76     // -- QFunctions  in main()
87     // -- Operators  in main()
98     // -- Create qdata  in main()
101     // -- Composite operators  in main()
104     // -- Cleanup  in main()
127     // -- Fine grid operator  in main()
130     // -- Restrictions  in main()
131     CeedInt offset = num_elem_sub * i * (p_coarse - 1);  in main()
134         ind_u_coarse[p_coarse * j + k] = offset + j * (p_coarse - 1) + k;  in main()
140     // -- Bases  in main()
143     // -- Create multigrid level  in main()
147     // -- Composite operators  in main()
152     // -- Cleanup  in main()
164   // Check output  in main()
174 …if (fabs(sum - 2.) > 1000. * CEED_EPSILON) printf("Computed Area Coarse Grid: %f != True Area: 2.0…  in main()
183   // Check output  in main()
193 …if (fabs(sum - 2.) > 1000. * CEED_EPSILON) printf("Computed Area Fine Grid: %f != True Area: 2.0\n…  in main()
198   // Check output  in main()
208 …if (fabs(sum - 2.) > 1000. * CEED_EPSILON) printf("Computed Area Coarse Grid: %f != True Area: 2.0…  in main()