59 …[3], const PetscInt i_rank[3], PetscInt degree, const PetscInt mesh_elem[3], PetscInt m_nodes[3]) {  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()
68         PetscInt m_nodes[3], ijk_rank[] = {i, j, k};  in GlobalStart()  local
70         GlobalNodes(p, ijk_rank, degree, mesh_elem, m_nodes);  in GlobalStart()
71         start += m_nodes[0] * m_nodes[1] * m_nodes[2];  in GlobalStart()
79   CeedInt        m_nodes[3], *idx, *idx_p;  in CreateRestriction()  local
83   for (int d = 0; d < 3; d++) m_nodes[d] = mesh_elem[d] * (P - 1) + 1;  in CreateRestriction()
92 … P + jj) * P + kk] = num_comp * (((i * (P - 1) + ii) * m_nodes[1] + (j * (P - 1) + jj)) * m_nodes[…  in CreateRestriction()
94 …P * (jj + P * kk)] = num_comp * (((i * (P - 1) + ii) * m_nodes[1] + (j * (P - 1) + jj)) * m_nodes[…  in CreateRestriction()
104 …emRestrictionCreate(ceed, num_elem, P * P * P, num_comp, 1, m_nodes[0] * m_nodes[1] * m_nodes[2] *…  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
432   GlobalNodes(p, i_rank, degree, mesh_elem, m_nodes);  in main()
437   PetscCall(VecSetSizes(X, m_nodes[0] * m_nodes[1] * m_nodes[2] * num_comp_u, PETSC_DECIDE));  in main()
467 …elem, mesh_elem[0], mesh_elem[1], mesh_elem[2], m_nodes[0] * m_nodes[1] * m_nodes[2], m_nodes[0], …  in main()
468                           m_nodes[2], num_comp_u));  in main()
502 …for (PetscInt i = 0, ir, ii; ir = i >= m_nodes[0], ii = i - ir * m_nodes[0], i < l_nodes[0]; i++) {  in main()
503 …for (PetscInt j = 0, jr, jj; jr = j >= m_nodes[1], jj = j - jr * m_nodes[1], j < l_nodes[1]; j++) {  in main()
504 …for (PetscInt k = 0, kr, kk; kr = k >= m_nodes[2], kk = k - kr * m_nodes[2], k < l_nodes[2]; k++) {  in main()
658 …PetscCall(MatCreateShell(comm, m_nodes[0] * m_nodes[1] * m_nodes[2] * num_comp_u, m_nodes[0] * m_n…  in main()