2 /// Test 1D mass matrix operator at points with heterogeneous points per element
3 /// \test Test 1D mass matrix operator at points with heterogeneous points per element
13   CeedInt    num_elem = 3, dim = 1, p = 3, q = 5;  in main()
14 …CeedInt    num_nodes_x = num_elem + 1, num_nodes_u = num_elem * (p - 1) + 1, num_points_per_elem =…  in main()
15   CeedInt    ind_x[num_elem * 2], ind_u[num_elem * p], ind_x_points[num_elem + 1 + num_points];  in main()
24   CeedInit(argv[1], &ceed);  in main()
27   for (CeedInt i = 0; i < num_nodes_x; i++) x_array_mesh[i] = (CeedScalar)i / (num_nodes_x - 1);  in main()
30     ind_x[2 * i + 1] = i + 1;  in main()
32 …CeedElemRestrictionCreate(ceed, num_elem, 2, 1, 1, num_nodes_x, CEED_MEM_HOST, CEED_USE_POINTER, i…  in main()
39       ind_u[p * i + j] = i * (p - 1) + j;  in main()
42 …CeedElemRestrictionCreate(ceed, num_elem, p, 1, 1, num_nodes_u, CEED_MEM_HOST, CEED_USE_POINTER, i…  in main()
46     CeedScalar weight_tmp[num_points_per_elem + 1];  in main()
49     // Use num_points_per_elem + 1 to test non-uniform quadrature  in main()
50     CeedGaussQuadrature(num_points_per_elem + 1, x_array_points, weight_tmp);  in main()
51     ind_x_points[0] = num_elem + 1;  in main()
52     for (CeedInt p = 0; p < num_points_per_elem + 1; p++, current_index++) {  in main()
53       ind_x_points[num_elem + 1 + current_index] = current_index;  in main()
56     for (CeedInt e = 1; e < num_elem - 1; e++) {  in main()
58       ind_x_points[e] = num_elem + 1 + current_index;  in main()
60         ind_x_points[num_elem + 1 + current_index] = current_index;  in main()
63     // Use num_points_per_elem - 1 to test non-uniform quadrature  in main()
64     CeedGaussQuadrature(num_points_per_elem - 1, &x_array_points[current_index], weight_tmp);  in main()
65     ind_x_points[num_elem - 1] = num_elem + 1 + current_index;  in main()
66     for (CeedInt p = 0; p < num_points_per_elem - 1; p++, current_index++) {  in main()
67       ind_x_points[num_elem + 1 + current_index] = current_index;  in main()
69     ind_x_points[num_elem] = num_elem + 1 + current_index;  in main()
73 …CeedElemRestrictionCreateAtPoints(ceed, num_elem, num_points, 1, num_points, CEED_MEM_HOST, CEED_C…  in main()
75 …CeedElemRestrictionCreateAtPoints(ceed, num_elem, num_points, 1, num_points, CEED_MEM_HOST, CEED_C…  in main()
84   CeedBasisCreateTensorH1Lagrange(ceed, dim, 1, p, q, CEED_GAUSS, &basis_u);  in main()
87   CeedQFunctionCreateInterior(ceed, 1, setup, setup_loc, &qf_setup);  in main()
88   CeedQFunctionAddInput(qf_setup, "weight", 1, CEED_EVAL_WEIGHT);  in main()
90   CeedQFunctionAddOutput(qf_setup, "rho", 1, CEED_EVAL_NONE);  in main()
101   CeedQFunctionCreateInterior(ceed, 1, mass, mass_loc, &qf_mass);  in main()
102   CeedQFunctionAddInput(qf_mass, "u", 1, CEED_EVAL_INTERP);  in main()
103   CeedQFunctionAddInput(qf_mass, "rho", 1, CEED_EVAL_NONE);  in main()
104   CeedQFunctionAddOutput(qf_mass, "v", 1, CEED_EVAL_INTERP);  in main()
128       if (fabs(v_array[i]) > 1e-14) printf("[%" CeedInt_FMT "] v %g != 0.0\n", i, v_array[i]);  in main()