59 …etscInt p[3], const PetscInt i_rank[3], PetscInt degree, const PetscInt mesh_elem[3], PetscInt m_n…  in GlobalNodes()
60   for (int d = 0; d < 3; d++) m_nodes[d] = degree * mesh_elem[d] + (i_rank[d] == p[d] - 1);  in GlobalNodes()
62 …tart(const PetscInt p[3], const PetscInt i_rank[3], PetscInt degree, const PetscInt mesh_elem[3]) {  in GlobalStart()
70         GlobalNodes(p, ijk_rank, degree, mesh_elem, m_nodes);  in GlobalStart()
77 static PetscErrorCode CreateRestriction(Ceed ceed, const PetscInt mesh_elem[3], CeedInt P, CeedInt …  in CreateRestriction()
78   const PetscInt num_elem = mesh_elem[0] * mesh_elem[1] * mesh_elem[2];  in CreateRestriction()
83   for (int d = 0; d < 3; d++) m_nodes[d] = mesh_elem[d] * (P - 1) + 1;  in CreateRestriction()
85   for (CeedInt i = 0; i < mesh_elem[0]; i++) {  in CreateRestriction()
86     for (CeedInt j = 0; j < mesh_elem[1]; j++) {  in CreateRestriction()
87       for (CeedInt k = 0; k < mesh_elem[2]; k++, idx_p += P * P * P) {  in CreateRestriction()
342 …PetscInt degree, q_extra, local_nodes, local_elem, mesh_elem[3], m_nodes[3], p[3], i_rank[3], l_no…  in main()  local
420     Split3(local_elem, mesh_elem, true);  in main()
421     if (Max3(mesh_elem) / Min3(mesh_elem) <= 2) break;  in main()
432   GlobalNodes(p, i_rank, degree, mesh_elem, m_nodes);  in main()
467 …local_elem, mesh_elem[0], mesh_elem[1], mesh_elem[2], m_nodes[0] * m_nodes[1] * m_nodes[2], m_node…  in main()
474       l_nodes[d] = mesh_elem[d] * degree + 1;  in main()
492           g_start[i][j][k]     = GlobalStart(p, ijk_rank, degree, mesh_elem);  in main()
493           GlobalNodes(p, ijk_rank, degree, mesh_elem, g_m_nodes[i][j][k]);  in main()
553   PetscCall(CreateRestriction(ceed, mesh_elem, P, num_comp_u, &elem_restr_u));  in main()
554   PetscCall(CreateRestriction(ceed, mesh_elem, 2, dim, &elem_restr_x));  in main()
555   CeedInt num_elem = mesh_elem[0] * mesh_elem[1] * mesh_elem[2];  in main()
562 …CeedInt     shape[3] = {mesh_elem[0] + 1, mesh_elem[1] + 1, mesh_elem[2] + 1}, len = shape[0] * sh…  in main()
568 … * shape[1] + j) * shape[2] + k) + 0] = 1. * (i_rank[0] * mesh_elem[0] + i) / (p[0] * mesh_elem[0]…  in main()
569 … * shape[1] + j) * shape[2] + k) + 1] = 1. * (i_rank[1] * mesh_elem[1] + j) / (p[1] * mesh_elem[1]…  in main()
570 … * shape[1] + j) * shape[2] + k) + 2] = 1. * (i_rank[2] * mesh_elem[2] + k) / (p[2] * mesh_elem[2]…  in main()