1570     PetscInt       c, f, fx, fy, fz, e, ex, ey, ez;  in DMPlexCreateCubeMesh_Internal()  local
1580         for (fx = 0; fx < numXEdges; ++fx) {  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()
1601 …if (bdX != DM_BOUNDARY_NONE && fx == numXEdges - 1 && cutLabel) PetscCall(DMLabelSetValue(cutLabel…  in DMPlexCreateCubeMesh_Internal()
1610         for (fx = 0; fx < numXVertices; ++fx) {  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()
1622               if (fx == numXVertices - 1) {  in DMPlexCreateCubeMesh_Internal()
1625               } else if (fx == 0) {  in DMPlexCreateCubeMesh_Internal()
1642       for (fx = 0; fx < numXEdges; ++fx) {  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()
1675       for (fx = 0; fx < numXEdges; ++fx) {  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()
1686             if (bdX == DM_BOUNDARY_TWIST && fx == numXEdges - 1) {  in DMPlexCreateCubeMesh_Internal()
1691               edgeT += numXEdges - 1 - 2 * fx;  in DMPlexCreateCubeMesh_Internal()
1694 …if (bdX != DM_BOUNDARY_NONE && fx == numXEdges - 1 && cutLabel) PetscCall(DMLabelSetValue(cutLabel…  in DMPlexCreateCubeMesh_Internal()