47 int        GetCartesianMeshSize(CeedInt dim, CeedInt degree, CeedInt prob_size, CeedInt num_xyz[dim…
48 int        BuildCartesianRestriction(Ceed ceed, CeedInt dim, CeedInt num_xyz[dim], CeedInt degree, …
50 int        SetCartesianMeshCoords(CeedInt dim, CeedInt num_xyz[dim], CeedInt mesh_degree, CeedVecto…
51 CeedScalar TransformMeshCoords(CeedInt dim, CeedInt mesh_size, CeedVector mesh_coords);
56   CeedInt     dim         = 3;               // dimension of the mesh  in main()  local
73       parse_error = next_arg ? dim = atoi(argv[++ia]), 0 : 1;  in main()
74       num_comp_x                   = dim;  in main()
101     printf("  Mesh dimension         [-d] : %" CeedInt_FMT "\n", dim);  in main()
123 …CeedBasisCreateTensorH1Lagrange(ceed, dim, num_comp_x, mesh_degree + 1, num_qpts, CEED_GAUSS, &mes…  in main()
124   CeedBasisCreateTensorH1Lagrange(ceed, dim, 1, sol_degree + 1, num_qpts, CEED_GAUSS, &sol_basis);  in main()
127   CeedInt num_xyz[dim];  in main()
129   GetCartesianMeshSize(dim, sol_degree, prob_size, num_xyz);  in main()
133     if (dim > 1) printf(", ny = %" CeedInt_FMT, num_xyz[1]);  in main()
134     if (dim > 2) printf(", nz = %" CeedInt_FMT, num_xyz[2]);  in main()
143 …BuildCartesianRestriction(ceed, dim, num_xyz, mesh_degree, num_comp_x, &mesh_size, num_qpts, &mesh…  in main()
144 …BuildCartesianRestriction(ceed, dim, num_xyz, sol_degree, 1 + dim * (dim + 1) / 2, &sol_size, num_…  in main()
145 …BuildCartesianRestriction(ceed, dim, num_xyz, sol_degree, 1, &sol_size, num_qpts, &sol_restriction…  in main()
148     printf("Number of mesh nodes     : %" CeedInt_FMT "\n", mesh_size / dim);  in main()
157   SetCartesianMeshCoords(dim, num_xyz, mesh_degree, mesh_coords);  in main()
160   CeedScalar exact_volume = TransformMeshCoords(dim, mesh_size, mesh_coords);  in main()
166   build_ctx_data.dim = build_ctx_data.space_dim = dim;  in main()
174   CeedQFunctionAddInput(qf_build, "dx", num_comp_x * dim, CEED_EVAL_GRAD);  in main()
176   CeedQFunctionAddOutput(qf_build, "qdata", 1 + dim * (dim + 1) / 2, CEED_EVAL_NONE);  in main()
189   CeedInt    elem_qpts = CeedIntPow(num_qpts, dim);  in main()
192   for (CeedInt d = 0; d < dim; d++) num_elem *= num_xyz[d];  in main()
193   CeedVectorCreate(ceed, num_elem * elem_qpts * (1 + dim * (dim + 1) / 2), &q_data);  in main()
201   CeedQFunctionAddInput(qf_apply, "du", dim, CEED_EVAL_GRAD);  in main()
202   CeedQFunctionAddInput(qf_apply, "qdata", 1 + dim * (dim + 1) / 2, CEED_EVAL_NONE);  in main()
204   CeedQFunctionAddOutput(qf_apply, "dv", dim, CEED_EVAL_GRAD);  in main()
259     CeedScalar tol = (dim == 1 ? 200. * CEED_EPSILON : dim == 2 ? 1E-5 : 1E-5);  in main()
283 int GetCartesianMeshSize(CeedInt dim, CeedInt degree, CeedInt prob_size, CeedInt num_xyz[dim]) {  in GetCartesianMeshSize()  argument
286   CeedInt num_elem = prob_size / CeedIntPow(degree, dim);  in GetCartesianMeshSize()
293   CeedInt r = s % dim;  in GetCartesianMeshSize()
295   for (CeedInt d = 0; d < dim; d++) {  in GetCartesianMeshSize()
296     CeedInt sd = s / dim;  in GetCartesianMeshSize()
307 int BuildCartesianRestriction(Ceed ceed, CeedInt dim, CeedInt num_xyz[dim], CeedInt degree, CeedInt…  in BuildCartesianRestriction()  argument
310   CeedInt num_nodes = CeedIntPow(p, dim);         // number of scalar nodes per element  in BuildCartesianRestriction()
311   CeedInt elem_qpts = CeedIntPow(num_qpts, dim);  // number of qpts per element  in BuildCartesianRestriction()
314   for (CeedInt d = 0; d < dim; d++) {  in BuildCartesianRestriction()
328     for (CeedInt d = 0; d < dim; d++) {  in BuildCartesianRestriction()
337       for (CeedInt d = 0; d < dim; d++) {  in BuildCartesianRestriction()
356 int SetCartesianMeshCoords(CeedInt dim, CeedInt num_xyz[dim], CeedInt mesh_degree, CeedVector mesh_…  in SetCartesianMeshCoords()  argument
360   for (CeedInt d = 0; d < dim; d++) {  in SetCartesianMeshCoords()
375     for (CeedInt d = 0; d < dim; d++) {  in SetCartesianMeshCoords()
391 CeedScalar TransformMeshCoords(CeedInt dim, CeedInt mesh_size, CeedVector mesh_coords) {  in TransformMeshCoords()  argument
396   if (dim == 1) {  in TransformMeshCoords()
403     CeedInt num_nodes = mesh_size / dim;  in TransformMeshCoords()