9 static CeedScalar Eval(CeedInt dim, const CeedScalar x[]) {  in Eval()  argument
11   if (dim > 1) result += (x[1] + 0.2) * (x[1] + 0.2);  in Eval()
12   if (dim > 2) result += -(x[2] + 0.3) * (x[2] + 0.3);  in Eval()
16 static CeedScalar EvalGrad(CeedInt dim, const CeedScalar x[]) {  in EvalGrad()  argument
17   switch (dim) {  in EvalGrad()
27 static CeedScalar GetTolerance(CeedScalarType scalar_type, int dim) {  in GetTolerance()  argument
30     if (dim == 3) tol = 1.e-4;  in GetTolerance()
38 static void VerifyProjectedBasis(CeedBasis basis_project, CeedInt dim, CeedInt p_to_dim, CeedInt p_…  in VerifyProjectedBasis()  argument
46     tol = GetTolerance(scalar_type, dim);  in VerifyProjectedBasis()
56       CeedScalar coord[dim];  in VerifyProjectedBasis()
57       for (CeedInt d = 0; d < dim; d++) coord[d] = x_array[p_from_dim * d + i];  in VerifyProjectedBasis()
58       u_array[i] = Eval(dim, coord);  in VerifyProjectedBasis()
74       CeedScalar coord[dim];  in VerifyProjectedBasis()
75       for (CeedInt d = 0; d < dim; d++) coord[d] = x_array[d * p_to_dim + i];  in VerifyProjectedBasis()
76       const CeedScalar u = Eval(dim, coord);  in VerifyProjectedBasis()
77 …if (fabs(u - u_array[i]) > tol) printf("[%" CeedInt_FMT ", %" CeedInt_FMT "] %f != %f\n", dim, i, …  in VerifyProjectedBasis()
93       CeedScalar coord[dim];  in VerifyProjectedBasis()
95       for (CeedInt d = 0; d < dim; d++) coord[d] = x_array[p_to_dim * d + i];  in VerifyProjectedBasis()
96       for (CeedInt d = 0; d < dim; d++) {  in VerifyProjectedBasis()
101 …nt_FMT ", %" CeedInt_FMT ", %" CeedInt_FMT "] %f != %f\n", dim, i, d, du_array[p_to_dim * (dim - 1…  in VerifyProjectedBasis()
116   for (CeedInt dim = 1; dim <= 3; dim++) {  in main()  local
119 …p_to = 5, q = 6, x_dim = CeedIntPow(2, dim), p_from_dim = CeedIntPow(p_from, dim), p_to_dim = Ceed…  in main()
121     CeedVectorCreate(ceed, x_dim * dim, &x_corners);  in main()
123       CeedScalar x_array[x_dim * dim];  in main()
125       for (CeedInt d = 0; d < dim; d++) {  in main()
130     CeedVectorCreate(ceed, p_from_dim * dim, &x_from);  in main()
131     CeedVectorCreate(ceed, p_to_dim * dim, &x_to);  in main()
136     CeedVectorCreate(ceed, p_to_dim * dim, &du_to);  in main()
140     CeedBasisCreateTensorH1Lagrange(ceed, dim, dim, 2, p_from, CEED_GAUSS_LOBATTO, &basis_x);  in main()
143     CeedBasisCreateTensorH1Lagrange(ceed, dim, dim, 2, p_to, CEED_GAUSS_LOBATTO, &basis_x);  in main()
148     CeedBasisCreateTensorH1Lagrange(ceed, dim, 1, p_from, q, CEED_GAUSS, &basis_from);  in main()
149     CeedBasisCreateTensorH1Lagrange(ceed, dim, 1, p_to, q, CEED_GAUSS, &basis_to);  in main()
152 … VerifyProjectedBasis(basis_project, dim, p_to_dim, p_from_dim, x_to, x_from, u_to, u_from, du_to);  in main()
181 … VerifyProjectedBasis(basis_project, dim, p_to_dim, p_from_dim, x_to, x_from, u_to, u_from, du_to);  in main()
186 … VerifyProjectedBasis(basis_project, dim, p_to_dim, p_from_dim, x_to, x_from, u_to, u_from, du_to);  in main()
191 … VerifyProjectedBasis(basis_project, dim, p_to_dim, p_from_dim, x_to, x_from, u_to, u_from, du_to);  in main()