48 int        GetCartesianMeshSize(CeedInt dim, CeedInt degree, CeedInt prob_size, CeedInt num_xyz[dim…
49 …artesianRestriction(Ceed ceed, CeedInt dim, CeedInt num_xyz[dim], CeedInt degree, CeedInt num_comp…
51 int        SetCartesianMeshCoords(CeedInt dim, CeedInt num_xyz[dim], CeedInt mesh_degree, CeedVecto…
52 CeedScalar TransformMeshCoords(CeedInt dim, CeedInt mesh_size, CeedVector mesh_coords);
57   CeedInt     dim         = 3;               // dimension of the mesh  in main()
58   CeedInt     num_comp_x  = 3;               // number of x components  in main()
59   CeedInt     mesh_degree = 4;               // polynomial degree for the mesh  in main()
60   CeedInt     sol_degree  = 4;               // polynomial degree for the solution  in main()
61   CeedInt     num_qpts    = sol_degree + 2;  // number of 1D quadrature points  in main()
62   CeedInt     prob_size   = -1;              // approximate problem size  in main()
63   CeedInt     help = 0, test = 0, gallery = 0, benchmark = 0;  in main()
130   CeedInt num_xyz[dim];  in main()
143   CeedInt             mesh_size, sol_size;  in main()
199   CeedInt    elem_qpts = CeedIntPow(num_qpts, dim);  in main()
200   CeedInt    num_elem  = 1;  in main()
202   for (CeedInt d = 0; d < dim; d++) num_elem *= num_xyz[d];  in main()
246   for (CeedInt i = 0; i < benchmark; i++) {  in main()
259     for (CeedInt i = 0; i < sol_size; i++) volume += v_array[i];  in main()
294 int GetCartesianMeshSize(CeedInt dim, CeedInt degree, CeedInt prob_size, CeedInt num_xyz[dim]) {  in GetCartesianMeshSize()
297   CeedInt num_elem = prob_size / CeedIntPow(degree, dim);  in GetCartesianMeshSize()
298   CeedInt s        = 0;  // find s: num_elem/2 < 2^s <= num_elem  in GetCartesianMeshSize()
304   CeedInt r = s % dim;  in GetCartesianMeshSize()
306   for (CeedInt d = 0; d < dim; d++) {  in GetCartesianMeshSize()
307     CeedInt sd = s / dim;  in GetCartesianMeshSize()
318 …artesianRestriction(Ceed ceed, CeedInt dim, CeedInt num_xyz[dim], CeedInt degree, CeedInt num_comp…  in BuildCartesianRestriction()
320   CeedInt p         = degree + 1;  in BuildCartesianRestriction()
321   CeedInt num_nodes = CeedIntPow(p, dim);         // number of scalar nodes per element  in BuildCartesianRestriction()
322   CeedInt elem_qpts = CeedIntPow(num_qpts, dim);  // number of qpts per element  in BuildCartesianRestriction()
323   CeedInt nd[3], num_elem = 1, scalar_size = 1;  in BuildCartesianRestriction()
325   for (CeedInt d = 0; d < dim; d++) {  in BuildCartesianRestriction()
334   CeedInt *elem_nodes = malloc(sizeof(CeedInt) * num_elem * num_nodes);  in BuildCartesianRestriction()
336   for (CeedInt e = 0; e < num_elem; e++) {  in BuildCartesianRestriction()
337     CeedInt e_xyz[3] = {1, 1, 1}, re = e;  in BuildCartesianRestriction()
339     for (CeedInt d = 0; d < dim; d++) {  in BuildCartesianRestriction()
343     CeedInt *local_elem_nodes = elem_nodes + e * num_nodes;  in BuildCartesianRestriction()
345     for (CeedInt l_nodes = 0; l_nodes < num_nodes; l_nodes++) {  in BuildCartesianRestriction()
346       CeedInt g_nodes = 0, g_nodes_stride = 1, r_nodes = l_nodes;  in BuildCartesianRestriction()
348       for (CeedInt d = 0; d < dim; d++) {  in BuildCartesianRestriction()
365 int SetCartesianMeshCoords(CeedInt dim, CeedInt num_xyz[dim], CeedInt mesh_degree, CeedVector mesh_…  in SetCartesianMeshCoords()
366   CeedInt p = mesh_degree + 1;  in SetCartesianMeshCoords()
367   CeedInt nd[3], scalar_size = 1;  in SetCartesianMeshCoords()
369   for (CeedInt d = 0; d < dim; d++) {  in SetCartesianMeshCoords()
380   for (CeedInt i = 0; i < p; i++) nodes[i] = 0.5 + 0.5 * nodes[i];  in SetCartesianMeshCoords()
381   for (CeedInt gs_nodes = 0; gs_nodes < scalar_size; gs_nodes++) {  in SetCartesianMeshCoords()
382     CeedInt r_nodes = gs_nodes;  in SetCartesianMeshCoords()
384     for (CeedInt d = 0; d < dim; d++) {  in SetCartesianMeshCoords()
385       CeedInt d_1d = r_nodes % nd[d];  in SetCartesianMeshCoords()
401 CeedScalar TransformMeshCoords(CeedInt dim, CeedInt mesh_size, CeedVector mesh_coords) {  in TransformMeshCoords()
407     for (CeedInt i = 0; i < mesh_size; i++) {  in TransformMeshCoords()
413     CeedInt num_nodes = mesh_size / dim;  in TransformMeshCoords()
415     for (CeedInt i = 0; i < num_nodes; i++) {  in TransformMeshCoords()