59 static void     GlobalNodes(const PetscInt p[3], const PetscInt i_rank[3], PetscInt degree, const P…  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 static PetscInt GlobalStart(const PetscInt p[3], const PetscInt i_rank[3], PetscInt degree, const P…  in GlobalStart()
69         if (i == i_rank[0] && j == i_rank[1] && k == i_rank[2]) return start;  in GlobalStart()
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
428     i_rank[d]            = rank_left / pstride[d];  in main()
429     rank_left -= i_rank[d] * pstride[d];  in main()
432   GlobalNodes(p, i_rank, degree, mesh_elem, m_nodes);  in main()
491           PetscInt ijk_rank[3] = {i_rank[0] + i, i_rank[1] + j, i_rank[2] + k};  in main()
507 …      if ((i_rank[0] == 0 && i == 0) || (i_rank[1] == 0 && j == 0) || (i_rank[2] == 0 && k == 0) ||  in main()
508 …(i_rank[0] + 1 == p[0] && i + 1 == l_nodes[0]) || (i_rank[1] + 1 == p[1] && j + 1 == l_nodes[1]) ||  in main()
509               (i_rank[2] + 1 == p[2] && k + 1 == l_nodes[2])) {  in main()
568 …x_loc[dim * ((i * shape[1] + j) * shape[2] + k) + 0] = 1. * (i_rank[0] * mesh_elem[0] + i) / (p[0]…  in main()
569 …x_loc[dim * ((i * shape[1] + j) * shape[2] + k) + 1] = 1. * (i_rank[1] * mesh_elem[1] + j) / (p[1]…  in main()
570 …x_loc[dim * ((i * shape[1] + j) * shape[2] + k) + 2] = 1. * (i_rank[2] * mesh_elem[2] + k) / (p[2]…  in main()