8 static CeedScalar Eval(CeedInt dim, const CeedScalar x[]) {  in Eval()  argument
10   for (CeedInt d = 0; d < dim; d++) {  in Eval()
17 static CeedScalar EvalGrad(CeedInt dim, const CeedScalar x[], CeedInt direction) {  in EvalGrad()  argument
19   for (CeedInt d = 0; d < dim; d++) {  in EvalGrad()
32   for (CeedInt dim = 1; dim <= 3; dim++) {  in main()  local
35 …const CeedInt p = 9, q = 9, num_points = 4, x_dim = CeedIntPow(2, dim), p_dim = CeedIntPow(p, dim);  in main()
37     CeedVectorCreate(ceed, x_dim * dim, &x);  in main()
38     CeedVectorCreate(ceed, p_dim * dim, &x_nodes);  in main()
39     CeedVectorCreate(ceed, num_points * dim, &x_points);  in main()
41     CeedVectorCreate(ceed, num_points * dim, &v);  in main()
44     CeedBasisCreateTensorH1Lagrange(ceed, dim, dim, 2, p, CEED_GAUSS_LOBATTO, &basis_x);  in main()
46       CeedScalar x_array[x_dim * dim];  in main()
48       for (CeedInt d = 0; d < dim; d++) {  in main()
62         CeedScalar coord[dim];  in main()
64         for (CeedInt d = 0; d < dim; d++) coord[d] = x_array[d * p_dim + i];  in main()
65         u_array[i] = Eval(dim, coord);  in main()
72     CeedBasisCreateTensorH1Lagrange(ceed, dim, 1, p, q, CEED_GAUSS, &basis_u);  in main()
86         CeedScalar coord[dim];  in main()
88         for (CeedInt d = 0; d < dim; d++) coord[d] = x_array[d * num_points + i];  in main()
89         for (CeedInt d = 0; d < dim; d++) {  in main()
90           const CeedScalar dfx = EvalGrad(dim, coord, d);  in main()
94 …    printf("[%" CeedInt_FMT "] %f != %f = df(%f", dim, v_array[d * num_points + i], dfx, coord[0]);  in main()
95             for (CeedInt d = 1; d < dim; d++) printf(", %f", coord[d]);  in main()