1570     PetscInt       c, f, fx, fy, fz, e, ex, ey, ez;  in DMPlexCreateCubeMesh_Internal()  local
1578     for (fz = 0; fz < numZEdges; ++fz) {  in DMPlexCreateCubeMesh_Internal()
1581           PetscInt cell  = (fz * numYEdges + fy) * numXEdges + fx;  in DMPlexCreateCubeMesh_Internal()
1582           PetscInt faceB = firstZFace + (fy * numXEdges + fx) * numZVertices + fz;  in DMPlexCreateCubeMesh_Internal()
1583 …    PetscInt faceT = firstZFace + (fy * numXEdges + fx) * numZVertices + ((fz + 1) % numZVertices);  in DMPlexCreateCubeMesh_Internal()
1584           PetscInt faceF = firstYFace + (fz * numXEdges + fx) * numYVertices + fy;  in DMPlexCreateCubeMesh_Internal()
1585 …    PetscInt faceK = firstYFace + (fz * numXEdges + fx) * numYVertices + ((fy + 1) % numYVertices);  in DMPlexCreateCubeMesh_Internal()
1586           PetscInt faceL = firstXFace + (fz * numYEdges + fy) * numXVertices + fx;  in DMPlexCreateCubeMesh_Internal()
1587 …    PetscInt faceR = firstXFace + (fz * numYEdges + fy) * numXVertices + ((fx + 1) % numXVertices);  in DMPlexCreateCubeMesh_Internal()
1603 …if (bdZ != DM_BOUNDARY_NONE && fz == numZEdges - 1 && cutLabel) PetscCall(DMLabelSetValue(cutLabel…  in DMPlexCreateCubeMesh_Internal()
1608     for (fz = 0; fz < numZEdges; ++fz) {  in DMPlexCreateCubeMesh_Internal()
1611           PetscInt face    = firstXFace + (fz * numYEdges + fy) * numXVertices + fx;  in DMPlexCreateCubeMesh_Internal()
1612           PetscInt edgeL   = firstZEdge + (fy * numXVertices + fx) * numZEdges + fz;  in DMPlexCreateCubeMesh_Internal()
1613 …  PetscInt edgeR   = firstZEdge + (((fy + 1) % numYVertices) * numXVertices + fx) * numZEdges + fz;  in DMPlexCreateCubeMesh_Internal()
1614           PetscInt edgeB   = firstYEdge + (fz * numXVertices + fx) * numYEdges + fy;  in DMPlexCreateCubeMesh_Internal()
1615 …  PetscInt edgeT   = firstYEdge + (((fz + 1) % numZVertices) * numXVertices + fx) * numYEdges + fy;  in DMPlexCreateCubeMesh_Internal()
1641     for (fz = 0; fz < numZEdges; ++fz) {  in DMPlexCreateCubeMesh_Internal()
1644           PetscInt face    = firstYFace + (fz * numXEdges + fx) * numYVertices + fy;  in DMPlexCreateCubeMesh_Internal()
1645           PetscInt edgeL   = firstZEdge + (fy * numXVertices + fx) * numZEdges + fz;  in DMPlexCreateCubeMesh_Internal()
1646 …  PetscInt edgeR   = firstZEdge + (fy * numXVertices + ((fx + 1) % numXVertices)) * numZEdges + fz;  in DMPlexCreateCubeMesh_Internal()
1647           PetscInt edgeB   = firstXEdge + (fz * numYVertices + fy) * numXEdges + fx;  in DMPlexCreateCubeMesh_Internal()
1648 …  PetscInt edgeT   = firstXEdge + (((fz + 1) % numZVertices) * numYVertices + fy) * numXEdges + fx;  in DMPlexCreateCubeMesh_Internal()
1676         for (fz = 0; fz < numZVertices; fz++) {  in DMPlexCreateCubeMesh_Internal()
1677           PetscInt face    = firstZFace + (fy * numXEdges + fx) * numZVertices + fz;  in DMPlexCreateCubeMesh_Internal()
1678           PetscInt edgeL   = firstYEdge + (fz * numXVertices + fx) * numYEdges + fy;  in DMPlexCreateCubeMesh_Internal()
1679 …  PetscInt edgeR   = firstYEdge + (fz * numXVertices + ((fx + 1) % numXVertices)) * numYEdges + fy;  in DMPlexCreateCubeMesh_Internal()
1680           PetscInt edgeB   = firstXEdge + (fz * numYVertices + fy) * numXEdges + fx;  in DMPlexCreateCubeMesh_Internal()
1681 …  PetscInt edgeT   = firstXEdge + (fz * numYVertices + ((fy + 1) % numYVertices)) * numXEdges + fx;  in DMPlexCreateCubeMesh_Internal()
1699               if (fz == numZVertices - 1) {  in DMPlexCreateCubeMesh_Internal()
1702               } else if (fz == 0) {  in DMPlexCreateCubeMesh_Internal()